Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 58,993 (A)7→8 coding (520/651 nt) yabP → putative uncharacterized protein YabP

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091358,9861.A98.4% 428.7 / ‑1.3 121coding (513/651 nt)yabPputative uncharacterized protein YabP
Reads supporting (aligned to +/- strand):  ref base . (1/1);  new base A (55/64);  total (56/65)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.46e-01

Only 100 of 128 total aligned reads displayed.

TTAAGTGTCACTTTTTTCATTGAAAAGAATAAGATATGTACTCAGACGTTGCAGCTGCATAAGGGCCAAGGTAACTTGGGCGTTGATCTTAGAAAGGCTTACCTTCCCGGCGTTGACATGAGGGATTGTTACCTTGGTAAAAAAACAATGAAAGGTAGCAATGATATCCTTTATGAGAGACCTGGGTGGAATGCTAACCTGGGCGTGCTACCCCGGACGGTGCTACCCCGGACGGTGCTAACCCGGACGGTGCTAACCTGGACGGTGCTACCGTGA  >  NC_000913/58849‑59124
                                                                                                                                          |                                                                                                                                          
ttAAGTGTCACTTTTTTCATTGAAAAGAATAAGATATGTACTCAGACGTTGCAGCTGCATAAGGGCCAAGGTAACTTGGGCGTTGATCTTAGAAAGGCTTACCTTCCCGGCGTTGACATGAGGGATTGTTACCTTGGTaa                                                                                                                                          >  2:1020320/1‑140 (MQ=255)
 tAAGTGTCACTTTTTTCATTGAAAAGAATAAGATATGTACTCAGACGTTGCAGCTGCATAAGGGCCAAGGTAACTTGGGCGTTGATCTTAGAAAGGCTTACCTTCCCGGCGTTGACATGAGGGATTGTTACCTTGGTaaa                                                                                                                                         <  2:198802/140‑1 (MQ=255)
    gtgtCACTTTTTTCATTGAAAAGAATAAGATATGTACTCAGACGTTGCAGCTGCATAAGGGCCAAGGTAACTTGGGCGTTGATCTTAGAAAGGCTTACCTTCCCGGCGTTGACATGAGGGATTGTTACCTTGGTaaaaaa                                                                                                                                      <  1:192412/140‑1 (MQ=255)
     tgtCACTTTTTTCATTGAAAAGAATAAGATATGTACTCAGACGTTGCAGCTGCATAAGGGCCAAGGTAACTTGGGCGTTGATCTTAGAAAGGCTTACCTTCCCGGCGTTGACATGAGGGATTGTTACCTTGGTaaaaaaa                                                                                                                                     >  1:36093/1‑140 (MQ=255)
      gtCACTTTTTTCATTGAAAAGAATAAGATATGTACTCAGACGTTGCAGCTGCATAAGGGCCAAGGTAACTTGGGCGTTGATCTTAGAAAGGCTTACCTTCCCGGCGTTGACATGAGGGATTGTTACCTTGGTaaaaaaaa                                                                                                                                    <  1:1097012/140‑2 (MQ=255)
         aCTTTTTTCATTGAAAAGAATAAGATATGTACTCAGACGTTGCAGCTGCATAAGGGCCAAGGTAACTTGGGCGTTGATCTTAGAAAGGCTTACCTTCCCGGCGTTGACATGAGGGATTGTTACCTTGGTaaaaaaaacaa                                                                                                                                 >  2:477340/1‑136 (MQ=255)
              tttCATTGAAAAGAATAAGATATGTACTCAGACGTTGCAGCTGCATAAGGGCCAAGGTAACTTGGGCGTTGATCTTAGAAAGGCTTACCTTCCCGGCGTTGACATGAGGGATTGTTACCTTGGTAAAAAAAACAATGaaa                                                                                                                             <  1:1103093/140‑1 (MQ=255)
              tttCATTGAAAAGAATAAGATATGTACTCAGACGTTGCAGCTGCATAAGGGCCAAGGTAACTTGGGCGTTGATCTTAGAAAGGCTTACCTTCCCGGCGTTGACATGAGGGATTGTTACCTTGGTAAAAAAAACAATGaaa                                                                                                                             <  1:1238908/140‑1 (MQ=255)
                 cATTGAAAAGAATAAGATATGTACTCAGACGTTGCAGCTGCATAAGGGCCAAGGTAACTTGGGCGTTGATCTTAGAAAGGCTTACCTTCCCGGCGTTGACATGAGGGATTGTTACCTTGGTAAAAAAAACAATGAAAGGt                                                                                                                          >  2:310328/1‑140 (MQ=255)
                   ttGAAAAGAATAAGATATGTACTCAGACGTTGCAGCTGCATAAGGGCCAAGGTAACTTGGGCGTTGATCTTAGAAAGGCTTACCTTCCCGGCGTTGACATGAGGGATTGTTACCTTGGTAAAAAAAACAATGAAAGGTAg                                                                                                                        <  2:871118/140‑1 (MQ=255)
                   ttGAAAAGAATAAGATATGTACTCAGACGTTGCAGCTGCATAAGGGCCAAGGTAACTTGGGCGTTGATCTTAGAAAGGCTTACCTTCCCGGCGTTGACATGAGGGATTGTTACCTTGGTAAAAAAAACAATGAAAGGTAg                                                                                                                        <  1:908098/140‑1 (MQ=255)
                   ttGAAAAGAATAAGATATGTACTCAGACGTTGCAGCTGCATAAGGGCCAAGGTAACTTGGGCGTTGATCTTAGAAAGGCTTACCTTCCCGGCGTTGACATGAGGGATTGTTACCTTGGTAAAAAAAACAATGAAAGGTAg                                                                                                                        <  2:1161807/140‑1 (MQ=255)
                   ttGAAAAGAATAAGATAAGTACTCAGACGTTGCAGCTGCATAAGGGCCAAGGTAACTTGGGCGTTGATCTTAGAAAGGCTTACCTTCCCGGCGTTGACATGAGGGATTGTTACCTTGGTAAAAAAAACAATGAAAGGTAg                                                                                                                        <  2:212155/140‑1 (MQ=255)
                     gAAAAGAATAAGATATGTACTCAGACGTTGCAGCTGCATAAGGGCCAAGGTAACTTGGGCGTTGATCTTAGAAAGGCTTACCTTCCCGGCGTTGACATGAGGGATTGTTACCTTGGTAAAAAAAACAATGAAAGGTAGCa                                                                                                                      >  2:381475/1‑140 (MQ=255)
                         aGAATAAGATATGTACTCAGACGTTGCAGCTGCATAAGGGCCAAGGTAACTTGGGCGTTGATCTTAGAAAGGCTTACCTTCCCGGCGTTGACATGAGGGATTGTTACCTTGGTAAAAAAAACAATGAAAGGTAGCAATGa                                                                                                                  <  1:1117862/140‑1 (MQ=255)
                          gAATAAGATATGTACTCAGACGTTGCAGCTGCATAAGGGCCAAGGTAACTTGGGCGTTGATCTTAGAAAGGCTTACCTTCCCGGCGTTGACATGAGGGATTGTTACCTTGGTAAAAAAAACAATGAAAGGTAGCAATGat                                                                                                                 <  1:685080/140‑1 (MQ=255)
                           aaTAAGATATGTACTCAGACGTTGCAGCTGCATAAGGGCCAAGGTAACTTGGGCGTTGATCTTAGAAAGGCTTACCTTCCCGGCGTTGACATGAGGGATTGTTACCTTGGTAAAAAAAACAATGAAAGGTAGCAATGata                                                                                                                <  1:842353/140‑1 (MQ=255)
                                gATATGTACTCAGACGTTGCAGCTGCATAAGGGCCAAGGTAACTTGGGCGTTGATCTTAGAAAGGCTTACCTTCCCGGCGTTGACATGAGGGATTGTTACCTTGGTAAAAAAAACAATGAAAGGTAGCAATGATATCCtt                                                                                                           >  2:793191/1‑140 (MQ=255)
                                     gTACTCAGACGTTGCAGCTGCATAAGGGCCAAGGTAACTTGGGCGTTGATCTTAGAAAGGCTTACCTTCCCGGCGTTGACATGAGGGATTGTTACCTTGGTAAAAAAAACAATGAAAGGTAGCAATGATATCCTTTa                                                                                                         <  2:1344213/137‑1 (MQ=255)
                                     gTACTCAGACGTTGCAGCTGCATAAGGGCCAAGGTAACTTGGGCGTTGATCTTAGAAAGGCTTACCTTCCCGGCGTTGACATGAGGGATTGTTACCTTGGTAAAAAAAACAATGAAAGGTAGCAATGATATCCTTTATga                                                                                                      <  1:931686/140‑1 (MQ=255)
                                     gTACTCAGACGTTGCAGCTGCATAAGGGCCAAGGTAACTTGGGCGTTGATCTTAGAAAGGCTTACCTTCCCGGCGTTGACATGAGGGATTGTTACCTTGGTAAAAAAAACAATGAAAGGTAGCAATGATATCCTTTATga                                                                                                      <  1:793191/140‑1 (MQ=255)
                                     gTACTCAGACGTTGCAGCTGCATAAGGGCCAAGGTAACTTGGGCGTTGATCTTAGAAAGGCTTACCTTCCCGGCGTTGACATGAGGGATTGTTACCTTGGTAAAAAAAACAATGAAAGGTAGCAATGATATCCTTTATga                                                                                                      >  1:844466/1‑140 (MQ=255)
                                     gTACTCAGACGTTGCAGCTGCATAAGGGCCAAGGTAACTTGGGCGTTGATCTTAGAAAGGCTTACCTTCCCGGCGTTGACATGAGGGATTGTTACCTTGGTAAAAAAAACAATGAAAGGTAGCAATGATATCCTTTATga                                                                                                      <  1:280829/140‑1 (MQ=255)
                                     gTACTCAGACGTTGCAGCTGCATAAGGGCCAAGGTAACTTGGGCGTTGATCTTAGAAAGGCTTACCTTCCCGGCGTTGACATGAGGGATTGTTACCTTGGTAAAAAAAACAATGAAAGGTAGCAATGATATCCTTTATga                                                                                                      >  1:958785/1‑140 (MQ=255)
                                     gTACTCAGACGTTGCAGCTGCATAAGGGCCAAGGTAACTTGGGCGTTGATCTTAGAAAGGCTTACCTTCCCGGCGTTGACATGAGGGATTGTTACCTTGGTAAAAAAAACAATGAAAGGTAGCAATGATATCCTTTATga                                                                                                      >  2:679822/1‑140 (MQ=255)
                                     gTACTCAGACGTTGCAGCTGCATAAGGGCCAAGGTAACTTGGGCGTTGATCTTAGAAAGGCTTACCTTCCCGGCGTTGACATGAGGGATTGTTACCTTGGTAAAAAAAACAATGAAAGGTAGCAATGATATCCTTTATga                                                                                                      >  2:280829/1‑140 (MQ=255)
                                     gTACTCAGACGTTGCAGCTGCATAAGGGCCAAGGTAACTTGGGCGTTGATCTTAGAAAGGCTTACCTTCCCGGCGTTGACATGAGGGATTGTTACCTTGGTAAAAAAAACAATGAAAGGTAGCAATGATATCCTTTATga                                                                                                      >  1:123421/1‑140 (MQ=255)
                                       aCTCAGACGTTGCAGCTGCATAAGGGCCAAGGTAACTTGGGCGTTGATCTTAGAAAGGCTTACCTTCCCGGCGTTGACATGAGGGATTGTTACCTTGGTAAAAAAAACAATGAAAGGTAGCAATGATATCCTTTATgaga                                                                                                    <  2:123421/140‑1 (MQ=255)
                                          cAGACGTTGCAGCTGCATAAGGGCCAAGGTAACTTGGGCGTTGATCTTAGAAAGGCTTACCTTCCCGGCGTTGACATGAGGGATTGTTACCTTGGTAAAAAAAACAATGAAAGGTAGCAATGATATCCTTTATGAGAGAc                                                                                                 >  1:783872/1‑140 (MQ=255)
                                           aGACGTTGCAGCTGCATAAGGGCCAAGGTAACTTGGGCGTTGATCTTAGAAAGGCTTACCTTCCCGGCGTTGACATGAGGGATTGTTACCTTGGTAAAAAAAACAATGAAAGGTAGCAATGATATCCTTTATGAGAGAcc                                                                                                >  2:758480/1‑140 (MQ=255)
                                           aGACGTTGCAGCTGCATAAGGGCCAAGGTAACTTGGGCGTTGATCTTAGAAAGGCTTACCTTCCCGGCGTTGACATGAGGGATTGTTACCTTGGTAAAAAAAACAATGAAAGGTAGCAATGATATCCTTTATGAGAGAcc                                                                                                >  1:1658545/1‑140 (MQ=255)
                                            gACGTTGCAGCTGCATAAGGGCCAAGGTAACTTGGGCGTTGATCTTAGAAAGGCTTACCTTCCCGGCGTTGACATGAGGGATTGTTACCTTGGTAAAAAAAACAATGAAAGGTAGCAATGATATCCTTTATGAGAGACCt                                                                                               >  2:14743/1‑140 (MQ=255)
                                            gACGTTGCAGCTGCATAAGGGCCAAGGTAACTTGGGCGTTGATCTTAGAAAGGCTTACCTTCCCGGCGTTGACATGAGGGATTGTTACCTTGGTAAAAAAAACAATGAAAGGTAGCAATGATATCCTTTATGAGAGACCt                                                                                               <  1:1020320/140‑1 (MQ=255)
                                             aCGTTGCAGCTGCATAAGGGCCAAGGTAACTTGGGCGTTGATCTTAGAAAGGCTTACCTTCCCGGCGTTGACATGAGGGATTGTTACCTTGGTAAAAAAAACAATGAAAGGTAGCAATGata                                                                                                                <  2:951655/122‑1 (MQ=255)
                                                ttGCAGCTGCATAAGGGCCAAGGTAACTTGGGCGTTGATCTTAGAAAGGCTTACCTTCCCGGCGTTGACATGAGGGATTGTTACCTTGGTAAAAAAAACAATGAAAGGTAGCAATGTTATCCTTTATGAGAGACCTGGGt                                                                                           >  2:306073/1‑140 (MQ=255)
                                                ttGCAGCTGCATAAGGGCCAAGGTAACTTGGGCGTTGATCTTAGAAAGGCTTACCTTCCCGGCGTTGACATGAGGGATTGTTACCTTGGTAAAAAAAACAATGAAAGGTAGCAATGATATCCTTTATGAGAGACCTGGGt                                                                                           >  1:423573/1‑140 (MQ=255)
                                                ttGCAGCTGCATAAGGGCCAAGGTAACTTGGGCGTTGATCTTAGAAAGGCTTACCTTCCCGGCGTTGACATGAGGGATTGTTACCTTGGTAAAAAAAACAATGAAAGGTAGCAATGATATCCTTTATGAGAGACCTGGGt                                                                                           >  2:551278/1‑140 (MQ=255)
                                                ttGCAGCTGCATAAGGGCCAAGGTAACTTGGGCGTTGATCTTAGAAAGGCTTACCTTCCCGGCGTTGACATGAGGGATTGTTACCTTGGTAAAAAAACAATGAAAGGTAGCAATGATATCCTTTATGAGAGACCTGGGTg                                                                                          <  1:1611526/140‑1 (MQ=255)
                                                 tGCAGCTGCATAAGGGCCAAGGTAACTTGGGCGTTGATCTTAGAAAGGCTTACCTTCCCGGCGTTGACATGAGGGATTGTTACCTTGGTAAAAAAAACAATGAAAGGTAGCAATGATATCCTTTATGAGAGACCTGGGTg                                                                                          <  1:722631/140‑1 (MQ=255)
                                                 tGCAGCTGCATAAGGGCCAAGGTAACTTGGGCGTTGATCTTAGAAAGGCTTACCTTCCCGGCGTTGACATGAGGGATTGTTACCTTGGTAAAAAAAACAATGAAAGGTAGCAATGATATCCTTTATGAGAGACCTGGGTg                                                                                          <  1:310328/140‑1 (MQ=255)
                                                 tGCAGCTGCATAAGGGCCAAGGTAACTTGGGCGTTGATCTTAGAAAGGCTTACCTTCCCGGCGTTGACATGAGGGATTGTTACCTTGGTAAAAAAAACAATGAAAGGTAGCAATGATATCCTTTATGAGAGACCTGGGTg                                                                                          <  2:499627/140‑1 (MQ=255)
                                                 tGCAGCTGCATAAGGGCCAAGGTAACTTGGGCGTTGATCTTAGAAAGGCTTACCTTCCCGGCGTTGACATGAGGGATTGTTACCTTGGTAAAAAAAACAATGAAAGGTAGCAATGATATCCTTTATGAGAGACCTGGGTg                                                                                          <  2:844466/140‑1 (MQ=255)
                                                  gCAGCTGCATAAGGGCCAAGGTAACTTGGGCGTTGATCTTAGAAAGGCTTACCTTCCCGGCGTTGACATGAGGGATTGTTACCTTGGTAAAAAAAACAATGAAAGGTAGCAATGATATCCTTTATGAGAGACCTGGGTgg                                                                                         <  2:726368/140‑1 (MQ=255)
                                                  gCAGCTGCATAAGGGCCAAGGTAACTTGGGCGTTGATCTTAGAAAGGCTTACCTTCCCGGCGTTGACATGAGGGATTGTTACCTTGGTAAAAAAAACAATGAAAGGTAGCAATGATATCCTTTATGAGAGACCTGGGTgg                                                                                         >  2:1187997/1‑140 (MQ=255)
                                                     gCTGCATAAGGGCCAAGGTAACTTGGGCGTTGATCTTAGAAAGGCTTACCTTCCCGGCGTTGACATGAGGGATTGTTACCTTGGTAAAAAAAACAATGAAAGGTAGCAATGATATCCTTTATGAGAGACCTGGGTGGAAt                                                                                      >  1:470866/1‑140 (MQ=255)
                                                       tGCATAAGGGCCAAGGTAACTTGGGCGTTGATCTTAGAAAGGCTTACCTTCCCGGCGTTGACATGAGGGATTGTTACCTTGGTAAAAAAAACAATGAAAGGTAGCAATGATATCCTTTATGAGAGACCTGGGTGGAATGc                                                                                    >  2:417045/1‑140 (MQ=255)
                                                           tAAGGGCCAAGGTAACTTGGGCGTTGATCTTAGAAAGGCTTACCTTCCCGGCGTTGACATGAGGGATTGTTACCTTGGTAAAAAAAACAATGAAAGGTAGCAATGATATCCTTTATGAGAGACCTGGGTg                                                                                          >  1:565545/1‑130 (MQ=255)
                                                           tAAGGGCCAAGGTAACTTGGGCGTTGATCTTAGAAAGGCTTACCTTCCCGGCGTTGACATGAGGGATTGTTACCTTGGTAAAAAAAACAATGAAAGGTAGCAATGATATCCTTTATGAGAGACCTGGGTg                                                                                          <  2:565545/130‑1 (MQ=255)
                                                           tAAGGGCCAAGGTAACTTGGGCGTTGATCTTAGAAAGGCTTACCTTCCCGGCGTTGACATGAGGGATTGTTACCTTGGTAAAAAAAACAATGAAAGGTAGCAATGATATCCTTTATGAGAGACCTGGGTg                                                                                          >  1:204610/1‑130 (MQ=255)
                                                           tAAGGGCCAAGGTAACTTGGGCGTTGATCTTAGAAAGGCTTACCTTCCCGGCGTTGACATGAGGGATTGTTACCTTGGTAAAAAAAACAATGAAAGGTAGCAATGATATCCTTTATGAGAGACCTGGGTGGAATGCTAAc                                                                                <  2:1094268/140‑1 (MQ=255)
                                                           tAAGGGCCAAGGTAACTCGGGCGTTGATCTTAGAAAGGCTTACCTTCCCGGCGTTGACATGAGGGATTGTTACCTTGGTAAAAAAAACAATGAAAGGTAGCAATGATATCCTTTATGAGAGACCTGGGTGGAATGCTAAc                                                                                >  2:1544713/1‑140 (MQ=255)
                                                             aGGGCCAAGGTAACTTGGGCGTTGATCTTAGAAAGGCTTACCTTCCCGGCGTTGACATGAGGGATTGTTACCTTGGTAAAAAAAACAATGAAAGGTAGCAATGATATCCTTTATGAGAGACCTGGGTGGAATGCTAACCt                                                                              <  2:1095210/140‑1 (MQ=255)
                                                              gggCCAAGGTAACTTGGGCGTTGATCTTAGAAAGGCTTACCTTCCCGGCGTTGACATGAGGGATTGTTACCTTGGTAAAAAAAACAATGAAAGGTAGCAATGATATCCTTTATGAGAGACCTGGGTGGAATGCTAACCTg                                                                             <  1:1645869/140‑1 (MQ=255)
                                                                   aaGGTAACTTGGGCGTTGATCTTAGAAAGGCTTACCTTCCCGGCGTTGACATGAGGGATTGTTACCTTGGTAAAAAAAACAATGAAAGGTAGCAATGATATCCTTTATGAGAGACCTGGGTGGAATGCTAACCTGGGCgt                                                                        >  1:1030543/1‑140 (MQ=255)
                                                                          cTTGGGCGTTGATCTTAGAAAGGCTTACCTTCCCGGCGTTGACATGAGGGATTGTTACCTTGGTAAAAAAAACAATGAAAGGTAGCAATGATATCCTTTATGAGAGACCTGGGTGGAATGCTAACCTGGGCgtgctaccc                                                                 >  2:1065878/1‑140 (MQ=255)
                                                                          cTTGGGCGTTGATCTTAGAAAGGCTTACCTTCCCGGCGTTGACATGAGGGATTGTTACCTTGGTAAAAAAAACAATGAAAGGTAGCAATGATATCCTTTATGAGAGACCTGGGTGGAATGCTAACCTGGGCgtgctaccc                                                                 >  1:598537/1‑140 (MQ=255)
                                                                             gggCGTTGATCTTAGAAAGGCTTACCTTCCCGGCGTTGACATGAGGGATTGTTACCTTGGTAAAAAAAACAATGAAAGGTAGCAATGATATCCTTTATGAGAGACCTGGGTGGATGCTAACCTGGGCgtgctaccccgga                                                             <  1:1273222/140‑1 (MQ=255)
                                                                              ggCGTTGATCTTAGAAAGGCTTACCTTCCCGGCGTTGACATGAGGGATTGTTACCTTGGTAAAAAAAACAATGAAAGGTAGCAATGATATCCTTTATGAGAGACCTGGGTGGAATGCTAACCTGGGCgtgctaccccgga                                                             <  1:551278/140‑1 (MQ=255)
                                                                              ggCGTTGATCTTAGAAAGGCTTACCTTCCCGGCGTTGACATGAGGGATTGTTACCTTGGTAAAAAAAACAATGAAAGGTAGCAATGATATCCTTTATGAGAGACCTGGGTGGAATGCTAACCTGGGCgtgctaccccgga                                                             <  2:1206408/140‑1 (MQ=255)
                                                                              ggCGTTGATCTTAGAAAGGCTTACCTTCCCGGCGTTGACATGAGGGATTGTTACCTTGGTAAAAAAAACAATGAAAGGTAGCAATGATATCCTTTATGAGAGACCTGGGTGGAATGCTAACCTGGGCgtgctaccccgga                                                             <  2:423811/140‑1 (MQ=255)
                                                                              ggCGTTGATCTTAGAAAGGCTTACCTTCCCGGCGTTGACATGAGGGATTGTTACCTTGGTAAAAAAAACAATGAAAGGTAGCAATGATATCCTTTATGAGAGACCTGGGTGGAATGCTAACCTGGGCgtgctaccccgga                                                             <  1:1221975/140‑1 (MQ=255)
                                                                               gCGTTGATCTTAGAAAGGCTTACCTTCCCGGCGTTGACATGAGGGATTGTTACCTTGGTAAAAAAAACAATGAAAGGTAGCAATGATATCCTTTATGAGAGACCTGGGTGGAATGCTAACCTGGGCgtgctaccccggac                                                            >  2:1324191/1‑140 (MQ=255)
                                                                               gCGTTGATCTTAGAAAGGCTTACCTTCCCGGCGTTGACATGAGGGATTGTTACCTTGGTAAAAAAAACAATGAAAGGTAGCAATGATATCCTTTATGAGAGACCTGGGTGGAATGCTAACCTGGGCgtgctaccccgga                                                             <  2:502763/139‑1 (MQ=255)
                                                                               gCGTTGATCTTAGAAAGGCTTACCTTCCCGGCGTTGACATGAGGGATTGTTACCTTGGTAAAAAAAACAATGAAAGGTAGCAATGATATCCTTTATGAGAGACCTGGGTGGAATGCTAACCTGGGCgtgctaccccgga                                                             >  1:502763/1‑139 (MQ=255)
                                                                                   tGATCTTAGAAAGGCTTACCTTCCCGGCGTTGACATGAGGGATTGTTACCTTGGTAAAAAAAACAATGAAAGGTAGCAATGATATCCTTTATGAGAGACCTGGGTGGAATGCTAACCTGGGCgtgctaccccggacggtg                                                        <  1:522395/140‑1 (MQ=255)
                                                                                   tGATCTTAGAAAGGCTTACCTTCCCGGCGTTGACATGAGGGATTGTTACCTTGGTAAAAAAAACAATGAAAGGTAGCAATGATATCCTTTATGAGAGACCTGGGTGGAATGCTAACCTGGGCgtgctaccccggacggtg                                                        >  2:243601/1‑140 (MQ=255)
                                                                                             aaGGCTTACCTTCCCGGCGTTGACATGAGGGATTGTTACCTTGGTAAAAAAAACAATGAAAGGTAGCAATGATATCCTTTATGAGAGACCTGGGTGGAATGCTAACGTGGGCgtgctaccccggacggtgctaccccgga                                              <  2:1646362/140‑1 (MQ=255)
                                                                                             aaGGCTTACCTTCCCGGCGTTGACATGAGGGATTGTTACCTTGGTAAAAAAAACAATGAAAGGTAGCAATGATATCCTTTATGAGAGACCTGGGTGGAATGCTAACCTGGGCgtgctaccccggatggtgctaccccgga                                              <  2:184144/140‑1 (MQ=255)
                                                                                             aaGGCTTACCTTCCCGGCGTTGACATGAGGGATTGTTACCTTGGTAAAAAAAACAATGAAAGGTAGCAATGATATCCTTTATGAGAGACCTGGGTGGAATGCTAACCTGGGCgtgctaccccggacggtgctaccccgga                                              <  2:919658/140‑1 (MQ=255)
                                                                                             aaGGCTTACCTTCCCGGCGTTGACATGAGGGATTGTTACCTTGGTAAAAAAAACAATGAAAGGTAGCAATGATATCCTTTATGAGAGACCTGGGTGGAATGCTAACCTGGGCgtgctaccccggacggtgctaccccgga                                              <  2:1030543/140‑1 (MQ=255)
                                                                                             aaGGCTTACCTTCCCGGCGTTGACATGAGGGATTGTTACCTTGGTAAAAAAAACAATGAAAGGTAGCAATGATATCCTTTATGAGAGACCTGGGTGGAATGCTAACCTGGGCgtgctaccccggacggtgctaccccgga                                              <  1:1144682/140‑1 (MQ=255)
                                                                                             aaGGCTTACCTTCCCGGCGTTGACATGAGGGATTGTTACCTTGGTAAAAAAAACAATGAAAGGTAGCAATGATATCCTTTATGAGAGACCTGGGTGGAATGCTAACCTGGGCgtgctaccccggacggtgctaccccgga                                              <  1:1448383/140‑1 (MQ=255)
                                                                                              aGGCTTACCTTCCCGGCGTTGACATGAGGGATTGTTACCTTGGTAAAAAAAACAATGAAAGGTAGCAATGATATCCTTTATGAGAGACCTGGGTGGAATGCTAACCTGGGCgtgctaccccggacggtgctacaccggac                                             <  1:710006/140‑1 (MQ=255)
                                                                                                 cTTACCTTCCCGGCGTTGACATGAGGGATTGTTACCTTGGTAAAAAAAACAATGAAAGGTAGCAATGATATCCTTTATGAGAGACCTGGGTGGAATGCTAACCTGGGCGTgctaccccggacggtgctaccccggacggt                                          >  2:1004450/1‑140 (MQ=255)
                                                                                                 cTTACCTTCCCGGCGTTGACATGAGGGATTGTTACCTTGGTAAAAAAAACAATGAAAGGTAGCAATGATATCCTTTATGAGAGACCTGGGTGGAATGCTAACCTGGGCGTgctaccccggacggtgctaccccggacggt                                          >  2:843037/1‑140 (MQ=255)
                                                                                                 cTTACCTTCCCGGCGTTGACATGAGGGATTGTTACCTTGGTAAAAAAAACAATGAAAGGTAGCAATGATATCCTTTATGAGAGACCTGGGTGGAATGCTAACCTGGGCGTgctaccccggacggtgctaccccggacggt                                          >  2:1332281/1‑140 (MQ=255)
                                                                                                  ttACCTTCCCGGCGTTGACATGAGGGATTGTTACCTTGGTAAAAAAAACAATGAAAGGTAGCAATGATATCCTTTATGAGAGACCTGGGTGGAATGCTAACCTGGGCGTGctaccccggacggtgctaccccggacggtg                                         <  1:1544713/140‑1 (MQ=255)
                                                                                                         cccGGCGTTGACATGAGGGATTGTTACCTTGGTAAAAAAAACAATGAAAGGTAGCAATGATATCCTTTATGAGAGACCTGGGTGGAATGCTAACCTGGGCGTGCTACCCCGGACGGTGCTACcccggacggtgctaaccc                                  >  2:80317/1‑140 (MQ=255)
                                                                                                          ccGGCGTTGACATGAGGGATTGTTACCTTGGTAAAAAAAACAATGAAAGGTAGCAATGATATCCTTTATGAGAGACCTGGGTGGAATGCTAACCTGGGCGTGCTACCCCGGACGGTGCTACcccggacggtgctaacccg                                 >  2:493462/1‑140 (MQ=255)
                                                                                                            ggCGTTGACATGAGGGATTGTTACCTTGGTAAAAAAAACAATGAAAGGTAGCAATGATATCCTTTATGAGAGACCTGGGTGGAATGCTAACCTGGGCGTGCTACCCGGGACGGTGCTACcccggacggtgctaacccgga                               <  2:192382/140‑1 (MQ=255)
                                                                                                            ggCGTTGACATGAGGGATTGTTACCTTGGTAAAAAAAACAATGAAAGGTAGCAATGATATCCTTTATGAGAGACCTGGGTGGAATGCTAACCTGGGCGTGCTACCCCGGACGTTGCTACcccggacggtgctaacccgga                               <  1:306073/140‑1 (MQ=255)
                                                                                                            ggCGTTGACATGAGGGATTGTTACCTTGGTAAAAAAAACAATGAAAGGTAGCAATGATATCCTTTATGAGAGACCTGGGTGGAATGCTAACCTGGGCGTGCTACCCCGGACGGTGCTACcccggacggtgctaacccgga                               <  1:1004450/140‑1 (MQ=255)
                                                                                                                ttGACATGAGGGATTGTTACCTTGGTAAAAAAAACAATGAAAGGTAGCAATGATATCCTTTATGAGAGACCTGGGTGGAATGCTAACCTGGGCGTGCTACCCCGGACGGTGCTACcccggacggtgctaacccggacggt                           >  2:1434696/1‑140 (MQ=255)
                                                                                                                 tGACATGAGGGATTGTTACCTTGGTAAAAAAAACAATGAAAGGTAGCAATGATATCCTTTATGAGAGACCTGGGTGGAATGCTAACCTGGGCGTGCTACCCCGGACGGTGCTACcccggacggtgctaacccggacggtg                          <  1:1472761/140‑1 (MQ=255)
                                                                                                                    cATGAGGGATTGTTACCTTGGTAAAAAAAACAATGAAAGGTAGCAATGATATCCTTTATGAGAGACCTGGGTGGAATGCTAACCTGGGCGTGCTACCCCGGACGGTGCTACcccggacggtgctaacccggacggtgcta                       >  1:523389/1‑140 (MQ=255)
                                                                                                                        aGGGATTGTTACCTTGGTAAAAAAAACAATGAAAGGTAGCAATGATATCCTTTATGAGAGACCTGGGTGGAATGCTAACCTGGGCGTGCTACCCCGGACGGTGCTACcccggacggtgctaacccgga                               >  1:482599/1‑128 (MQ=255)
                                                                                                                           gATTGTTACCTTGGTAAAAAAAACAATGAAAGGTAGCAATGATATCCTTTATGAGAGACCTGGGTGGAATGCTAACCTGGGCGTGCTACCCCGGACGGTGCTACCCCGGACGGTGCTAACCCGGACGGTGCTAACCTGGa                <  2:423573/140‑1 (MQ=255)
                                                                                                                           gATTGTTACCTTGGTAAAAAAAACAATGAAAGGTAGCAATGATATCCTTTATGAGAGACCTGGGTGGAATGCTAACCTGGGCGTGCTACCCCGGACGGTGCTACCCCGGACGGTGCTAACCCGGACGGTGCTAACCTGGa                >  2:1546475/1‑140 (MQ=255)
                                                                                                                           gATTGTTACCTTGGTAAAAAAAACAATGAAAGGTAGCAATGATATCCTTTATGAGAGACCTGGGTGGAATGCTAACCTGGGCGTGCTACCCCGGACGGTGCTACCCCGGACGGTGCTAACCCGGACGGTGCTAACCTGGa                <  1:1065878/140‑1 (MQ=255)
                                                                                                                           gATTGTTACCTTGGTAAAAAAAACAATGAAAGGTAGCAATGATATCCTTTATGAGAGACCTGGGTGGAATGCTAACCTGGGCGTGCTACCCCGGACGGTGCTACCCCGGACGGTGCTAACCCGGACGGTGCTAACCTGGa                <  1:914940/140‑1 (MQ=255)
                                                                                                                           gATTGTTACCTTGGTAAAAAAAACAATGAAAGGTAGCAATGATATCCTTTATGAGAGACCTGGGTGGAATGCTAACCTGGGCGTGCTACCCCGGACGGTGCTACCCCGGACGGTGCTAACCCGGACGGTGCTAACCTGGa                <  1:493462/140‑1 (MQ=255)
                                                                                                                             ttGTTACCTTGGTAAAAAAAACAATGAAAGGTAGCAATGATATCCTTTATGAGAGACCTGGGTGGAATGCTAACCTGGGCGTGCTACCCCGGACGGTGCTACCCCGGACGGTGCTAACCCGGACGGTGCTAACCTGGACg              >  1:210045/1‑140 (MQ=255)
                                                                                                                             ttGTTACCTTGGTAAAAAAAACAATGAAAGGTAGCAATGATATCCTTTATGAGAGACCTGGGTGGAATGCTAACCTGGGCGTGCTACCCCGGACGGTGCTACCCCGGACGGTGCTAACCCGGACGGTGCTAACCTGGACg              >  1:535037/1‑140 (MQ=255)
                                                                                                                             ttGTTACCTTGGTAAAAAAAACAATGAAAGGTAGCAATGATATCCTTTATGAGAGACCTGGGTGGAATGCTAACCTGGGCGTGCTACCCCGGACGGTGCTACCCCGGACGGTGCTAACCCGGACGGTGCTAACCTGGACg              >  1:1426827/1‑140 (MQ=255)
                                                                                                                              tGTTACCTTGGTAAAAAAAACAATGAAAGGTAGCAATGATATCCTTTATGAGAGACCTGGGTGGAATGCTAACCTGGGCGTGCTACCCCGGACGGTGCTACCCCGGACGGTGCTAACCCGGACGGTGCTAACCTGGACgg             <  2:535037/140‑1 (MQ=255)
                                                                                                                              tGTTACCTTGGTAAAAAAAACAATGAAAGGTAGCAATGATATCCTTTATGAGAGACCTGGGTGGAATGCTAACCTGGGCGTGCTACCCCGGACGGTGCTACCCCGGACGGTGCTAACCCGGACGGTGCTAACCTGGACgg             >  2:393257/1‑140 (MQ=255)
                                                                                                                                ttACCTTGGTAAAAAAAACAATGAAAGGTAGCAATGATATCCTTTATGAGAGACCTGGGTGGAATGCTAACCTGGGCGTGCTACCCCGGACGGTGCTACCCCGGACGGTGCTAACCCGGACGGTGCTAACCTGGACGGTg           <  2:598537/140‑1 (MQ=255)
                                                                                                                                ttACCTTGGTAAAAAAAACAATGAAAGGTAGCAATGATATCCTTTATGAGAGACCTGGGTGGAATGCTAACCTGGGCGTGCTACCCCGGACGGTGCTACCCCGGACGGTGCTAACCCGGACGGTGCTAACCTGGACGGTg           <  1:1169885/140‑1 (MQ=255)
                                                                                                                                ttACCTTGGTAAAAAAAACAATGAAAGGTAGCAATGATATCCTTTATGAGAGACCTGGGTGGAATGCTAACCTGGGCGTGCTACCCCGGACGGTGCTACCCCGGACGGTGCTAACCCGGACGGTGCTAACCTGGACGGTg           <  1:1187997/140‑1 (MQ=255)
                                                                                                                                        gTAAAAAAACAATGAAAGGTAGCAATGATATCCTTTATGAGAGACCTGGGTGGAATGCTAACCTGGGCGTGCTACCCCGGACGGTGCTACCCCGGACGGTGCTAACCCGGACGGTGCTAACCTGGACGGTGCTACCGTGa  >  1:1191890/1‑140 (MQ=255)
                                                                                                                                          |                                                                                                                                          
TTAAGTGTCACTTTTTTCATTGAAAAGAATAAGATATGTACTCAGACGTTGCAGCTGCATAAGGGCCAAGGTAACTTGGGCGTTGATCTTAGAAAGGCTTACCTTCCCGGCGTTGACATGAGGGATTGTTACCTTGGTAAAAAAACAATGAAAGGTAGCAATGATATCCTTTATGAGAGACCTGGGTGGAATGCTAACCTGGGCGTGCTACCCCGGACGGTGCTACCCCGGACGGTGCTAACCCGGACGGTGCTAACCTGGACGGTGCTACCGTGA  >  NC_000913/58849‑59124

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: