Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 1,857,700 A→G 7.3% I225T (ATA→ACA)  USA300HOU_1728 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007301,857,7000AG 7.3% 210.2 / 2.1 91I225T (ATA→ACA) USA300HOU_1728hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (35/46);  new base G (4/3);  total (40/51)
Fisher's exact test for biased strand distribution p-value = 6.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTTGCTAATCTTTCAGCTTGTGACTCTGGTTTTGATACTTCTACCTTTGCTGCACTTGCATCTTTAGTTAATGTGTCACTGTCTTTCTCTGCTTGTGCTTCTTTTTTAGCTTCTTGAACCTCTTGTGCCTCTTGTGATATATCACTTAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCTTTCAATGCCACTGCTTTTGGTTCTTCATTTGATACAGCTGCACTTTTATCTGTTTCTGCTTGTGCTTCTTTTTTAGCTTCTTGAATC  >  CP000730/1857562‑1857840
                                                                                                                                          |                                                                                                                                            
tttGCTAATCTTTCAGCTTGTGACTCTGGTTTTGATACTTCTACCTTTGCTGCACTTGCATCTTTAGTTAATGTGTCACTGTCTTTCTCTGCTTGTGCTTCTTTTTTAGCTTCTTGAACCTCTTGTGCCTCTTGTGatata                                                                                                                                            <  2:457057/141‑1 (MQ=255)
    cTAATCTTTCAGCTTGTGACTCTGGTTTTGATACTTCTACCTTTGCTGCACTTGCATCTTTAGTTAATGTGTCACTGTCTTTCTCTGCTTGTGCTTCTTTTTTAGCTTCTTGAACCTCTTGTGCCTCTTGTGTTTTATccc                                                                                                                                        >  2:41539/1‑139 (MQ=255)
      aaTCTTTCGGCTTGTGACTCTGGTTTTGATACTTCTACCTTTGCTGCACTTGCATCTTTTGTTAATGTGTCACTGTCTTTCTCTGCTTGTGCTTCTTTTTTAGCTTCTTGAACCTCTTGTGCCTCTTGTGATATATCACtt                                                                                                                                      <  2:200860/141‑1 (MQ=255)
      aaTCTTTCAGCTTGTGACTCTGGTTTTGATACTTCTACCTTTGCTGCACTTGCATCTTTAGTTAATGTGTCACTGTCTTTCTCTGCTTGTGCTTCTTTTTTAGCTTCTTGAACCTCTTGTGCCTCTTGTGATATATCACtt                                                                                                                                      <  2:695986/141‑1 (MQ=255)
       aTCTTTCAGCTTGTGACTCTGGTTTTGATACTTCTACCTTTGCTGCACTTGCATCTTTAGTTAATGTGTCACTGTCTTTCTCTGCTTGTGCTTCTTTTTTAGCTTCTTGAACCTCTTGTGCCTCTTGTGATATATCACTTa                                                                                                                                     >  1:326327/1‑141 (MQ=255)
          tttCAGCTTGTGACTCTGGTTTTGATACTTCTACCTTTGCTGCACTTGCATCTTTAGTTAATGTGTCACTGTCTTTCTCTGCTTGTGCTTCTTTTTTAGCTTCTTGAACCTCTTGTGCCTCTTGTGATGTATCCCTTAAat                                                                                                                                  >  2:168946/1‑141 (MQ=255)
          tttCAGCTTGTGACTCCGGTTTTGATACTTCTACCTTTGCTGCACTTGCATCTCTAGTCAATGTGTCACTGTCTTTCTCTGCTTGTGCTTCTTTTTTCGCTTCTTGAACCTCTTGTGCCCCTTGTGATATATCACTTCAat                                                                                                                                  >  2:29306/1‑141 (MQ=38)
             cAGCTTGTGACTCTGGTTTTGATACTTCTACCTTTGCTGCACTTGCATCTTTAGTTAATGTGTCACTGTCTTTCTCTGCTTGTGCTTCTTTTTTAGCTTCTTGAACCTCTTGTGCCTCTTGTGATATATCACTTAAattat                                                                                                                               >  1:166601/1‑141 (MQ=255)
             cAGCTTGTGACTCTGGTTTTGATACTTCTACCTTTGCTGCACTTGCATCTTTAGTTAATGTGTCACTGTCTTTCTCTGCTTGTGCTTCTTTTTTAGCTTCTTGAACCTCTTGTGCCTCTTGTGATATATCACTTAAattat                                                                                                                               >  1:730860/1‑141 (MQ=255)
             cAGCTTGTGACTCTGGTTTTGATACTTCTACCTTTGCTGCACTTGCATCTTTAGTTAATGTGTCACTGTCTTTCTCTGCTTGTGCTTCTTTTTTAGCTTCTTGAACCTCTTGTCCCTCTTGTGATATATCACTTAAattct                                                                                                                               >  1:44398/1‑139 (MQ=255)
             cAGCTTGTGACTCTGGTTTAGATACTTCTACCTTTGCTGCACTTGCATCTTTAGTTAATGTGTCACTGTCTTTCTCTGCTTGTGCTTCTTTTTTAGCTTCTTGAACCTCTTGTGCCTCTTGTGATATATCACTTAAattat                                                                                                                               >  1:527451/1‑141 (MQ=255)
              aGCTTGTGACTCTGGTTTTGATACTTCTACCTTTGCTGCACTTGCATCTTTAGTTAATGTGTCACTGTCTTTCTCTGCTTGTGCTTCTTTTTTAGCTTCTTGAACCTCTTGTGCCTCTTGTGATATATCACTTAAattatt                                                                                                                              >  1:615139/1‑141 (MQ=255)
              aGCTTGTGACTCTGGTTTTGATACTTCTACCTTTGCTGCACTTGCATCTTTAGTTAATGTGTCACTGTCTTTCTCTGCTTGTGCTTCTTTTTTAGCTTCTTGAACCTCTTGTGCCTCTTGTGATATATCACTTAAattatt                                                                                                                              <  1:651011/141‑1 (MQ=255)
               gCTTGTGACTCTGGTTTTGATACTTCTACCTTTGCTGCACTTGCATCTTTAGTTAATGTGTCACTGTCTTTCTCTGCTTGTACTTCTTTTTTAGCTTCTTGAACCTCTTGTGCCTCTTGTGATATATCACTTAAATTAttt                                                                                                                             >  1:77953/1‑141 (MQ=255)
                    tgACTCTGGTTTTGATACTTCTACCTTTGCTGCACTTGCATCTTTAGTTAATGTGTCACTGTCTTTCTCTGCTTGTGCTTCTTTTTTAGCTTCTTGAACCTCTTGTGCCTCTTGTGATATATCACTTAAATTATTTGCACt                                                                                                                        >  1:118609/1‑141 (MQ=255)
                    tgACTCTGGTTTTGATACTTCTACCTTTGCTGCACTTGCATCTTTAGTTAATGTGTCACTGTCTTTCTCTGCTTGTGCTTCTTTTTTAGCTTCTTGAACCTCTTGTGCCTCTTGTGATATATCACTTAAATTATTTGCACt                                                                                                                        >  1:232674/1‑141 (MQ=255)
                     gACTCTGGTTTTGATACTTCTACCTTTGCTGCACTTGCATCTTTAGTTAATGTGTCACTGTCTTTCTCTGCTTGTGCTTCTTTTTTAGCTTCTTGAACCTCTTGTGCCTCTTGTGATATATCACTTAAATTATTTGCACtt                                                                                                                       >  2:660180/1‑141 (MQ=255)
                            gTTTTGAGACTTCTACCTTTGCTGCACTTGCATCTTTAGTTAATATGTACCTGTCTGTCTCTGCGTGTGCTTATTGTTCAGAGAAGGTAACCTCGTGGGCCGCTGGTGATATATCACTTAAATTATTTGCACTTGcttctt                                                                                                                >  2:570459/1‑141 (MQ=255)
                              tttGATACTTCTACCTTTGCTGCACTTGCATCTTTAGTTAATGTGTCACTGTCTTTCTCTGCTTGTGCTTCTTTTTTAGCTTCTTGAACCTCTTGTGCCTCTTGTGATATATCACTTAAATTATTTGCACTTGcttcttct                                                                                                              <  2:172421/141‑1 (MQ=255)
                              tgtGATACTTCTCCCTTTGCTGCACGTGCATCTTTAGTTAATGTGTCACTGTCTTTCTCTGCTTGTGCTTCTTTTTTAGCTTCTTGAACCTCTTGTGCCTCTTGTGATATATCACTTAAATTATTTGCACTTGcttcttct                                                                                                              <  2:279841/139‑1 (MQ=255)
                                 taTACGTCCACCTTTGCTGCACTTGCATCTTTAGTTAATGTGTCACTGTCTTTCTCTGCTTGTGCTTCTTTTTTAGCTTCTTGAACCTCTTGTGCCTCTTGTGATATATCACTTAAATTATTTGCACTTGCTTCTTCtttt                                                                                                           <  2:588289/140‑1 (MQ=255)
                                 gagtcgtCTACCTTGGCTGCACTTGCATCTTTATTTAAGGTGTCACTGTCTTTCTCTGCTTGTGCTTCTTTTTTAGCGTCTTGAACCTCTTGTGCCTCTTGTGATATATCACTTAAATTATTTGCACTTGCTTCTTCtttt                                                                                                           <  2:382296/135‑1 (MQ=21)
                                 gATACTTCTACCTTTGCTGCACTTGCATCTTTAGTTAATGTGTCACTGTCTTTCTCTGCTTGTGCTTCTTTTTTAGCTTCTTGAACCTCTTGTGCCTCTTGTGATATATCACTTAAATTATTTGCACTTGCTTCTTCtttt                                                                                                           >  1:32361/1‑141 (MQ=255)
                                       tCTACCTTTGCTGCACTTGCATCTTTAGTTAATGTGTCACTGTCTTTCTCTGCTTGTGCTTCTTTTTTAGCTTCTTGAACCTCTTGTGCCTCTTGTGATATATCACTTAAATTATTTGCACTTGCTTCTTCTTTTATCgct                                                                                                     <  2:467345/141‑1 (MQ=255)
                                         caCCTTTGTTGCTCTTTCATTTTTTGTTTCTGTGTCTTTGTCTTTCTCTGCTTGTGCTTCTTTTTTAGCTTCTTGAACCTCTTGTGCCTCTTGTGATATATCACTTAAATTATTTGCACTTGCTTCTTCTTTTATCgctgc                                                                                                   <  2:32361/140‑1 (MQ=17)
                                              ttGCTGCACTTGCATCTTTCGTTAATGTGTCACTGTCTTTCTCTGCTTGTGCTTCTTTTTTAGCTTCTTGAACCTCTTGTGCCTCTTGTGATATATCACTTAAATTGTTTGCTCTTGCTTCTTCTTTTATCGCTGCttgtt                                                                                              >  2:101456/1‑141 (MQ=255)
                                              ttGCTGCACTTGCATCTTTAGTTAATGTGTCACTGTCTTTCTCTGCTTGTGCTGCTTTTTTAGCTTCTTGAACCTCTTGTGCCTCTTGTGATATATCACTTAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCttgtt                                                                                              >  1:25992/1‑141 (MQ=255)
                                               tgctgcACTTGCATCTTTAGTTAATGTGTCACTGTCTTTCTCTGCTTGTGCTTCTTTTTTAGCTTCTTGAACCTCTTGTGCCTCTTGTGATATATCACTTAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCttgttg                                                                                             >  1:518735/1‑141 (MQ=255)
                                               tgctgcACTTGCATCTTTAGTTAATGTGTCACTGTCTTTCTCTGCTTGTGCTTCTTTTTTAGCTTCTTGAACCTCTTGTGCCTCTTGTGATATATCACTTAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCttgttg                                                                                             >  1:128763/1‑141 (MQ=255)
                                               tgctgcACTTGCATCTTTAGTTAATGTGTCACTGTCTTTCTCTGCTTGTGCTTCTTTTTTAGCTTCTTGAACCTCTTGTGCCTCTTGTGATATATCACTTAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCtggttg                                                                                             >  1:10196/1‑141 (MQ=255)
                                                ttttgttcttcATTTTTTGTTGATGTGTTACTGTCTTTCTCTGCTTGTGCTTCTTTTTTAGCTTCTTGAATCTCTTGTGCCTCTTGTGATGTATCACTTAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTtgttgt                                                                                            <  2:44398/131‑1 (MQ=14)
                                                tttggtcTTTCATTTTTTGTTACTGTGTCTCTGTCTTTTTCTGCTTGTGCTTCTTTTTTAGCTTCTTGAACCTCTTGTGCCTCTTGTGATGTATCACTCAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTtgttgt                                                                                            <  2:166601/135‑1 (MQ=11)
                                                 ctgcACTTGCATCTTTAGTTAATGTGTCACTGTCTTTCTCTGCTTGTGCTTCTTTTTTAGCTTCTTGAACCTCTTGTGCCTCTTGTGATATATCACTTAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTtgttgtt                                                                                           >  1:79204/1‑140 (MQ=255)
                                                 ctgcACTTGCATCTTTAGTTAATGTGTCACTGTCTTTCTCTGCTTGTGCTTCTTTTTTAGCTTCTTGAACCTCTTGTGCCTCTTGTGATATATCACTTAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTtgtg                                                                                           >  2:706897/1‑141 (MQ=255)
                                                        ttcATCTTTTGTTAATGTGTCACTGTCTTTCTCTGCTTGTGCTTCTTTTTTAGCTTCTTGAAGCTCTTGTGCCTCTTGTGATATATCACTTAAATTATTTGCACTTCCTTCTTCTTTTATCGCTGCTTGTTGTGCTTTCaa                                                                                    <  2:36697/139‑1 (MQ=255)
                                                           ttCTTTAGTTAATGTGTCTCTGTCTTTCTCTGCTTGTGCTTCTTTTTTAGCTTCTTGAACCTCTTGTGCCTCTTGTGATATATCACTTAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCTTTCAATGc                                                                                 <  2:701748/140‑1 (MQ=255)
                                                           aTCTTTAGTTAATTTGTCACTGTCTTTCTCTGCTTGTGCTTCTTTTTTAGCTTCTTGAACCTCTTGTGCCTCTTGTGATATATCACTTAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCTTTCAATGc                                                                                 <  2:639158/141‑1 (MQ=255)
                                                              tttAGTTAATGTGTCACTGTCTTTCTCTGCTTGTGCTTCTTTTTTAGCTTCTTGAACCTCTTGTGCCTCTTGTGATATATCACTTAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCTTTCAATGCCAc                                                                              <  1:794008/141‑1 (MQ=255)
                                                                tAGTTAATGTGTCACTGTCTTTCTCTGCTTGTGCTTCTTTTTTAGCTTCTTGAACCTCTTGTGCCTCTTGTGATATATCACTTAAATTATTTGCACTTGCTTCTTCTTTTATCACTGCTTGTTGTGCTTTCAATGCCACTg                                                                            >  1:268284/1‑141 (MQ=255)
                                                                 tgTTAATGTGTCACTGTCTTTTTCTGCTTGTGCTTCTTTTTTAGCTTCTTGAACCTCTTGTGCCTCTTGTGATATATCACTTAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCTTTCAATGCCACTGc                                                                           <  1:47607/140‑1 (MQ=255)
                                                                  gTTAATGTGTCACTGTCTTTCTCTGCTTGTGCTTCTTTTTTAGCTTCTTGAACCTCTTGTGCCTCTTGTGATATATCACTTAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCTTTCAATGCCACTGCt                                                                          >  1:260472/1‑141 (MQ=255)
                                                                        gtgtCACTGTCTTTCTTTGCTTGTGCTTCTTTTTTAGCTTCTTGACCCTCTTGTGCCTCTTGTGATATATCACTTAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCTTTCCATGCCACTGCTTTTGGt                                                                    <  2:164050/141‑1 (MQ=255)
                                                                        gtgtCACTGTCTTTCTCTGCTTGTGCTTCTTTTTTAGCTTCTTGAACCTCTTGTGCCTCTTGTGATATATCACTTAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCTTTCAATGCCACTGCTTTTGGt                                                                    <  1:41539/141‑1 (MQ=255)
                                                                            cACTGTCTTTCTCTGCTTGTGCTTCTTTTTTAGCTTCTTGAACCTCTTGTGCCTCTTGTGATATATCACTTAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCTTTCAATGCCACTGCTTTTGGttctt                                                                <  1:558992/141‑1 (MQ=255)
                                                                            cACTGTCTTTCTCTGCTTGTGCTTCTTTTTTAGCTTCTTGAACCTCTTGTCCCTCTTGTGATACATCACTTAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCTTTCAATGCCAATGCTTTTGGttctt                                                                <  2:552809/141‑1 (MQ=255)
                                                                             cctgtctgtCTCTTCTTGTGCTTCTTTTTTATCTTCTTGAACCTCTTGTGCCTCTTGTGATATATCACTTAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTCGTTGTGCTTTCAATGCCACTGCTTTTGGttcttc                                                               <  2:529112/140‑1 (MQ=255)
                                                                             aCTGTCTTTCTCTGCTTGTGCTTCTTTTTTAGCTTCTTGAACCTCTTGTGCCTCTTGTGATATATCACTTAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCTTTCAATGCCACTGCTTTTGGttcttc                                                               <  1:722997/141‑1 (MQ=255)
                                                                                gTCTTTCTCTGCTTGTGCTTCTTTTTTAGCTTCTTGAACCTCTTGTGCCTCTTGTGATATATCACTTAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCTTTCAATGCCACTGCTTTTGGTTCTTCAtt                                                            >  1:808576/1‑141 (MQ=255)
                                                                                  cTTTCTCTGCTTGTGCTTCTTTTTTAGCTTCTTGAACCTCTTGTGCCTCTTGTGATGTATCACTTAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCTTTCAATGCCACTGCTTTTGATTCTTTATTTg                                                          >  2:750968/1‑141 (MQ=9)
                                                                                  cTTTCTCTGCTTGTGCTTCTTTTTTAGCTTCTTGAACCTCTTGTGCCTCTTGTGATGTATCACTTAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCCTGTTGTGCTTTCACTGCCACTGCTTTTgcatctaaatcta                                                          >  2:564228/1‑129 (MQ=11)
                                                                                     tctcTGCTTGTGCTTCTTTTTTAGCTTCTTGAACCTCTTGTGCCTCTTGTGATATATCACTTAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCTTTCAATGCCACTGCTTTTGGTTCTTCATTTGATa                                                       <  1:49540/141‑1 (MQ=255)
                                                                                           cTTGTTCTTCTTTTTTAGCTTCTTGAACCTCTTGTGCCTCTTGTGATATATCACTTAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCTTTCAATGCCACTGCTTTTGGTTCTTCATTTGATACAGCTg                                                 <  2:71428/141‑1 (MQ=255)
                                                                                           cTTGTGCTTCTTTTTTTGCTTCTTGAGCCTCTTCTGCCTCTTGTGATATATCACTTAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGTTTTCATTGCCACTGCTTTTGGTTCTTCATTTGATACAGCTg                                                 <  2:78939/141‑1 (MQ=14)
                                                                                           cTTGTGCTTCTTTTTTGGCTTCTTGAACCTTTTGTGCCTCTTGTGATATATCACTTAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCTTTCAATGCCACTGCTTTTGGTTCTTCATTTGATACAGCTg                                                 <  2:411602/141‑1 (MQ=255)
                                                                                            ttGTGCTTCTTTTTTAGCTTCTTGAACCTCTTGTGCCTCTTGTGATGTATCACTTAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCTTTCAATGCCACTGCTTTTGGTTCTTCATTTGATACAGCTGc                                                >  1:434166/1‑141 (MQ=255)
                                                                                              tgtcttcttTTTTAGCTTCTTGAACCTCTTCTTCCTCGTGTGTGTTATCTCTTAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGGTGTGCTTTCTATGCCACTGCTTGTGGTTCTTCATTTGATACAGCTGCAc                                              <  2:366095/138‑1 (MQ=14)
                                                                                              gtgCTTCTTTTTTAGCTTCTTGAACCTCTTGTGCCTCTTGTGATATATCACTTAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCTTTCAATGCCACTGCTTTTGGTTCTTCATTTGATACAGCTGCAc                                              <  2:808576/141‑1 (MQ=255)
                                                                                               tgCTTCTTTTTTTGCTTCTTGAACCTCTTGTGCCTCTTGTGATATATCACTTAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGTTTTCAATGCCACTGCTTTTGGTTCTTCATTTGATACAGCTGCACt                                             <  2:626084/141‑1 (MQ=255)
                                                                                                gCTTCTTTTTTAGCTTCTTGAACCTCTTGTGCCTCTTGTGATGTATCACTTAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCTTTCAAt                                                                                   <  1:205947/102‑1 (MQ=17)
                                                                                                gCTTCTTTTTTAGCTTCTTGAACCTCTTGTGCCTCTTGTGATATATCACTTAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCTTTCAAt                                                                                   >  2:738690/1‑102 (MQ=255)
                                                                                                gCTTCTTTTTTAGCTTCTTGAACCTCTTGTGCCTCTTGTGATATATCACTTAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCTTTCAATGCCACTGCTTTTGGTTCTTCATTTGATACAGCTGCACtt                                            >  1:794326/1‑141 (MQ=255)
                                                                                                gCTTCTTTTTTAGCTTCTTGAACCTCTTGTGCCTCTTGTGATATATCACTTAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCTTTCAATGCCACTGCTTTTGGTTCTTCATTTGATACAGCTGAACtt                                            >  1:500672/1‑141 (MQ=255)
                                                                                                gCTTCTTTTTTAGCTTCTTGAACCTCTTGTGCCTCTTGTGATATATCACTCAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCTTTCAAt                                                                                   >  2:213468/1‑102 (MQ=18)
                                                                                                gCTTCTTTTTTAGCTTCTTGAACCTCTTGTGCCTCTTGTGATATATCACTCAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCTTTCAAt                                                                                   >  2:205947/1‑102 (MQ=18)
                                                                                                    cttTTTTAGCTTCTTGAACCTCTTGTGCCTCTTGTGATATATCACTTAAATTATTTGCTCTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCTTTCAATGCCACTGCTTTTGGTTCTTCATTTGATACAGCTGCACTTTTAt                                        <  2:118609/141‑1 (MQ=255)
                                                                                                    cttTTTTAGCTTCTTGAACCTCTTGTGCCTCTTGTGATATATCACTTAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCTTTCAATGCCACTGCTTTTGGTTCTTCATTTGATACAGCTGCACTTTTAt                                        <  2:232674/141‑1 (MQ=255)
                                                                                                       ttttAGCTTCTTGAACCTCTTGTGCCTCTTGTGATATATCACTTAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCTTTCAATGCCACTGCTTTTGGTTCTTCATTTGATACAGCTGCACTTTTATCTg                                     >  1:289882/1‑141 (MQ=255)
                                                                                                       ttttAGCTTCTTGAACCTCTTGTGCCTCTTGTGATATATCACTTAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCTTTCAATGCCACTGCTTTTGGTTCTTCATTTGATACAGCTGCACTTTTATCTg                                     >  1:509085/1‑141 (MQ=255)
                                                                                                          tAGCTTCTTGAACCTCTTGTGCCTCTTTTGATATATCACTTAAATTATTTGCACTTGCTTCTTCTTTTCTCGCTGCTTGTTGTGCTTTCAATGCCACTGCTTTTGGTTCTTCATTTGATACAGCTGCACTTTTATCTGttt                                  <  2:514709/141‑1 (MQ=255)
                                                                                                           aGCTTCTTGAACCTCTTGTGCCTCTTGTGATATATCACTTAGATAATTTGCGGTTGTTTCTTCTTTTATCGCTGCTTGTTGTGCTTTCAATGCCACTGCTTTTGGTTCTTCATTTGATACAATTGCCATTTTATCTGGTTc                                 >  1:156671/1‑141 (MQ=14)
                                                                                                           aGCTTCTTGAACCTCTTGTGCCTCTTGTGATATATCACTTAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTGGTTGTGCTTTCAATGCCACTGCTTTTGGTTCTTCATTTGATACCGGTCCTCTTTTATCTGTTTc                                 >  1:231994/1‑141 (MQ=17)
                                                                                                                cttGAACCTCTTGTGCCTCTTGTGATATATCACTTAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCTTTCAATGCCACTGCTATTGGTTCTTCATTTGATACAGCTGCACTTTTATCTGTTTCTGCtt                            <  1:556088/141‑1 (MQ=255)
                                                                                                                cctGAACCTCTTGTGCCTCTTGTGATATATCACTTAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCTTTCAATGCCACTGCTTTTGGTTCTTCATTTGATACAGCTGCACTTTTATCTGTTTCTGCtt                            <  1:731948/139‑1 (MQ=255)
                                                                                                                  tGAACCTCTTGTGCCTCTTGTGATATATCACTTAAATTATTTGCACTTGCTTCTTCTTTTAt                                                                                                         >  2:200841/1‑62 (MQ=255)
                                                                                                                  tGAACCTCTTGTGCCTCTTGTGATATATCACTTAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCTTTCAATGCCCCTGCTTTTGGTTCTTCATTTGATACAGCTGCACTTTTATCTGTTTCTGCTtgt                          >  2:561190/1‑141 (MQ=255)
                                                                                                                   ggaCCTCTTGTGGCTCTTGTGATATATCACTTAAATTCTTTGCATTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCTTTCAATGCCACTGCTTTTGGTTCTTCATTTGATACAGCTGCACTTTTATCTGTTTCTGCTtgtg                         <  2:144523/139‑1 (MQ=255)
                                                                                                                   gctcctctttgtcctcttgtgATATATCACTTAAATTTTTTGCACTTGCTTCTTCTTTTTTCGCTGCTTGTTGTGCTTTCAATGCCACTGCTTTTGGTTCTTCATTTGATACAGCTGCACTTTTATCTGTTTCTGCTtgtg                         <  2:120708/129‑1 (MQ=17)
                                                                                                                   gAACCTCTTGTGCCTCTTGTGATATATCACTTAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCTTTCAATGCCACTGCTTTTGGTTCTTCATTTGATACAGCTGCACTTTTATCTGTTTCTGCTtgtg                         <  1:706897/141‑1 (MQ=255)
                                                                                                                    aaCCTCTTGTGCCTCTTGTGATATATCACTTAAATTATTTTCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCTTTCAATGCCACTGCTTTTGGTTCTTCATTTGATACAGCTGCACTTTTATCTGTTTCTGCTTGTGc                        <  2:328480/141‑1 (MQ=255)
                                                                                                                    aaCCTCTTGTGCCTCTTGTGATATATCACTTAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCTTTCAATGCCACTGCTTTTGGTTCTTCATTTGATACAGCTGCACTTTTATCTGTTTCTGCTTGTGc                        <  2:464463/141‑1 (MQ=255)
                                                                                                                          ttgtgccccttgtgATATATCACTTAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCTTTCAATGCCACTGCTTTTGGTTCTTCATTTGATACAGCTGCACTTTTATATGTTTCTGCTTGTGCTTCttt                  <  1:417099/141‑1 (MQ=255)
                                                                                                                            gtgcctcttgttATTGTTCCCTTAAATTTTTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCTTTCAATGCCACTGTTTTTGGTTCTTCATTTGTTACAGCTGCACTTTTATCTGTTTCTGCTTGTGCTTCttttt                <  2:152914/141‑1 (MQ=17)
                                                                                                                              gcctcttgtgATATATCACTTAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCTTTCAATGCCACTGCTTTTGGTTCTTTATTTGATACAGCCACACTTTTATCTGTTTCTGCTTGTGCTTCttttttg              >  2:399993/1‑140 (MQ=255)
                                                                                                                                ctcttgtgATATATCACTTAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCTTTCAATGCCACTGCTTTTGGTTCTTCATTTGATACAGCTGCACTTTTATCTGTTTCTGCTTGTGCTTCTTTTTTAGc            <  2:791636/141‑1 (MQ=255)
                                                                                                                                ctcttgtgATATATCACTTAAATTATTTGCACTTGCTTCTTCTTTTATCACTGCTTGTTGTGCTTTCAATGCCACTGCTTTTGGTTCTTCATTTGATACAGCTGCACTTTTATCTGTTTCTGCTTGTGCTTCTTTTTTAGc            <  2:268284/141‑1 (MQ=255)
                                                                                                                                 ttttgtgATATATCACTTAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCTTTCAATGCCACTGCTTTTGGTTCTTCATTTGATACAGCTGCACTTTTATCTGTTTCTGCTTGTGCTTCTTTTTTAGct           <  2:394582/139‑1 (MQ=255)
                                                                                                                                 tcttgtgATATATCACTTAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCTTTCAATGCCACTGCTTTTGGTTCTTCATTTGATACAGCTGCACTTTTATCTGTTTCTGCTTGTGCTTCTTTTTTAGct           <  1:6249/141‑1 (MQ=255)
                                                                                                                                 tcttgtgATATATCACTTAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCTTTCAATGCCACTGCTTTTGGTTCTTCATTTGATACAGCTGCACTTTTATCTGTTTCTGCTTGTGCTTCTTTTTTAGct           <  1:445964/141‑1 (MQ=255)
                                                                                                                                  cttgtgATATATCACTTAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCTTTCAATGCCACTGCTTTTGGTTCTTCATTTGATACAGCTGGACTTTTATCTGTTTCTGCTTGTGCTTCTTTTTTAGctt          <  1:520827/141‑1 (MQ=255)
                                                                                                                                      tgATATATCACTTAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCTTTCAATGCCACTGCTTTTGGTTCTTCATTTGATACAGCTGCACTTTTATCTGTTTCTGCTTGTGCTTCTTTTTTAGCTTCTTg      <  1:791772/141‑1 (MQ=255)
                                                                                                                                       gATATATCACTTAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCTTTCAATGCCACTGCTTTTGGTTCTTCATTTGATACAGCTGCACTTTTATCTGTTTCTGCTTGTGCTTCTTTTTTAGCTTCTTGa     <  1:711166/141‑1 (MQ=255)
                                                                                                                                       gATATATCACTTAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCTTTCAATGCCACTGCTTTTGGTTCTTCATTTGATACAGCTGCACTTTTATCTGTTTCTGCTTGTGCTTCTTTTTTAGCTTCTTGa     <  1:559803/141‑1 (MQ=255)
                                                                                                                                        atatatCACTTAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCTTTCAATGCCACTGCTTTTGGTTCTTCATTTGATACAGCTGCACTTTTATCTGTTTCTGCTTGTGCTTCTTTTTTAGCTTCTTGaa    <  1:265396/141‑1 (MQ=255)
                                                                                                                                        atatatCACTTAAACTATTTCCACGTGCTTCTTCTTTTATCGCTGCTTGTTGTGCCTTCAATGCCACTGCTTTTGGTTCTTCATTTGATACAGCTGCACTTTTATCTGTTTCTGCTTGTGCTTCTTTTTTAGCTTCTTGaa    <  1:619865/141‑1 (MQ=255)
                                                                                                                                          atatCACTTAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCTTTCAATGCCACTGCTTTTGGTTCTTTATTTGATACAATTGCACTTTTATCTGTTTCTGCTTGTGCTTCTTTTTTAGCTTCTTGTAtc  >  1:626748/1‑141 (MQ=255)
                                                                                                                                          atatCACTTAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCTTTCAATGCCACTGCTTTTGGTTCTTCATTTGATACAGCTGCACTTTTATCTGTTTCTGCTTGTGCTTCTTTTTTAGCTTCTTGAAtc  >  1:424037/1‑141 (MQ=255)
                                                                                                                                          |                                                                                                                                            
TTTGCTAATCTTTCAGCTTGTGACTCTGGTTTTGATACTTCTACCTTTGCTGCACTTGCATCTTTAGTTAATGTGTCACTGTCTTTCTCTGCTTGTGCTTCTTTTTTAGCTTCTTGAACCTCTTGTGCCTCTTGTGATATATCACTTAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCTTTCAATGCCACTGCTTTTGGTTCTTCATTTGATACAGCTGCACTTTTATCTGTTTCTGCTTGTGCTTCTTTTTTAGCTTCTTGAATC  >  CP000730/1857562‑1857840

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: