Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,596,492 C→T 6.4% pseudogene (1595/5422 nt) yneO ← AIDA‑I family autotransporter YneO

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,596,4920CT 6.4% 380.1 / 8.5 123pseudogene (1595/5422 nt)yneOAIDA‑I family autotransporter YneO
Reads supporting (aligned to +/- strand):  ref base C (55/59);  new base T (6/2);  total (61/62)
Fisher's exact test for biased strand distribution p-value = 2.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.52e-01

Only 100 of 123 total aligned reads displayed.

TCTGTGTATCACGCGCTTTAACTATACCTCCGGTGGTAATATTCAGTTCACCCGTGCCTAATACACCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGCATCCATCAGGTTCAGTGTGCCGA  >  NC_000913/1596346‑1596639
                                                                                                                                                  |                                                                                                                                                   
tCTGTGTATCACGCGCTTTAACTATACCTCCGGTGGTAATATTCAGTTCACCCGTGCCTAATACACCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCgg                                                                                                                                                  >  2:732150/1‑150 (MQ=255)
  tgtgTATCACGCGCTTTAACTATACCTCCGGTGGTAATATTCAGTTCACCCGTGCCTAATACACCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAg                                                                                                                                                >  1:2052094/1‑150 (MQ=255)
  tgtgTATCACGCGCTTTAACTATACCTCCGGTGGTAATATTCAGTTCACCCGTGCCTAATACACCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAg                                                                                                                                                >  2:2060353/1‑150 (MQ=255)
  tgtgTATCACGCGCTTTAACTATACCTCCGGTGGTAATATTCAGTTCACCCGTGCCTAATACACCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAg                                                                                                                                                <  1:2302444/150‑1 (MQ=255)
     gtATCACGCGCTTTAACTATACCTCCGGTGGTAATATTCAGTTCACCCGTGCCTAATACACCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCa                                                                                                                                             <  2:2340894/150‑1 (MQ=255)
       aTCACGCGCTTTAACTATACCTCCGGTGGTAATATTCAGTTCACCCGTGCCTAATACACCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACAGGAGGCAtt                                                                                                                                           >  1:718635/1‑150 (MQ=255)
        tCACGCGCTTTAACTATACCTCCGGTGGTAATATTCAGTTCACCCGTGCCTAATACACCGACCTTTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCAttt                                                                                                                                          <  2:2086251/150‑1 (MQ=255)
        tCACGCGCTTTAACTATACCTCCGGTGGTAATATTCAGTTCACCCGTGCCTAATACACCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCAttt                                                                                                                                          <  2:2183857/150‑1 (MQ=255)
        tCACGCGCTTTAACTATACCTCCGGAGGTAATATTCAGTTCACCCGTGCCTAATACACCGACCTGTAGCAATTGTGCGTTAGTAGATGATGGCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCAttt                                                                                                                                          <  2:1028036/150‑1 (MQ=255)
         cACGCGCTTTAACTATACCTCCGGTGGTAATATTCAGTTCACCCGTGCCTAATACACCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTtc                                                                                                                                         >  1:1808923/1‑150 (MQ=255)
          aCGCGCTTTAACTATACCTCCGGTGGTAATATTCAGTTCACCCGTGCCTAATACACCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGCCTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTtcc                                                                                                                                        >  2:1220158/1‑150 (MQ=255)
          aCGCGCTTTAACTATACCTCCGGTGGTAATATTCAGTTCACCCGTGCCTAATACACCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTtcc                                                                                                                                        >  2:2231099/1‑150 (MQ=255)
          aCGCGCTTTAACTATACCTCCGGTGGTAATATTCAGTTCACCCGTGCCTAATACACCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTtcc                                                                                                                                        <  1:1753554/150‑1 (MQ=255)
             cgcTTTAACTATACCTCCGGTGGTAATATTCAGTTCACCCGTGCCTAATACACCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTtcctcc                                                                                                                                     >  2:1620928/1‑150 (MQ=255)
             cgcTTTAACTATACCTCCGGTGGTAATATTCAGTTCACCCGTGCCTAATACACCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTtcctcc                                                                                                                                     <  2:1084434/150‑1 (MQ=255)
              gcTTTAACTATACCTCCGGTGGTAATATTCAGTTCACCCGTGCCTAATACACCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCg                                                                                                                                    <  1:330552/150‑1 (MQ=255)
              gcTTTAACTATACCTCCGGTGGTAATATTCAGTTCACCCGTGCCTAATACACCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCg                                                                                                                                    <  1:2253144/150‑1 (MQ=255)
                 ttAACTATACCTCCGGTTGTAATATTCAGTTCACCCGTGCCTAATACACCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGAAATGCCTCCGGTa                                                                                                                                 <  1:474337/150‑1 (MQ=255)
                 ttAACTATACCTCCGGTGGTAATATTCAGTTCACCCGTGCCTAATACACCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATa                                                                                                                                 <  1:1000630/150‑1 (MQ=255)
                 ttAACTATACCTCCGGTGGTAATATTCAGTTCACCCGTGCCTAATACACCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATa                                                                                                                                 <  1:2285453/150‑1 (MQ=255)
                  tAACTATACCTCCGGTGGTAATATTCAGTTCACCCGTGCCTAATACACCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATaa                                                                                                                                >  1:1441545/1‑150 (MQ=255)
                    aCTATACCTCCGGTGGTAATATTCAGTTCACCCGTGCCTAATACACCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTATCGGAGGAAATGCCTCCGGTaagg                                                                                                                              <  2:645315/150‑3 (MQ=255)
                     cTATACCTCCGGTGGTAATATTCAGTTCACCCGTGCCTAATACACCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATc                                                                                                                             >  1:925719/1‑150 (MQ=255)
                       ataCCTCCGGTGGTAATATTCAGTTCACCCGTGCCTAATACACCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCa                                                                                                                           <  2:1037357/150‑1 (MQ=255)
                           cTCCGGTGGTAATATTCAGTTCACCCGTGCCTAATACACCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGc                                                                                                                       >  2:1967015/1‑150 (MQ=255)
                           cTCCGGTGGTAATATTCAGTTCACCCGTGCCTAATACACCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGc                                                                                                                       >  2:2367617/1‑150 (MQ=255)
                           cTCCGGTGGTAATATTCAGTTCACCCGTGCCTAATACACCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGc                                                                                                                       >  1:623404/1‑150 (MQ=255)
                             ccGGTGGTAATATTCAGTTCACCCGTGCCTAATACACCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCACCTTACCGGAGGCATTTCCTCCGATAAATCCATAGcac                                                                                                                     >  2:417657/1‑148 (MQ=255)
                              cGGTGGTAATATTCAGTTCACCCGTGCCTAATACACCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCg                                                                                                                    <  2:261002/150‑1 (MQ=255)
                               ggtggtAATATTCAGTTCACCCGTGCCTAATACACCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGt                                                                                                                   >  1:518264/1‑150 (MQ=255)
                                 tggtAATATTCAGTTCACCCGTGCCTAATACACCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTg                                                                                                                 <  1:106364/150‑1 (MQ=255)
                                  tgtAATATTCAGTTCACCCGTGCCTAATACACCGAACTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGa                                                                                                                <  2:10163/149‑1 (MQ=255)
                                        aatCAGTCCACCCGCGCCTAATACCCCGACCTGTAGCAATTGCGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTGGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAg                                                                                                          <  2:2073627/148‑1 (MQ=255)
                                         ttCAGTTCACCCGTGCCTAATACACCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAgg                                                                                                         <  2:2343894/150‑1 (MQ=255)
                                          tCAGTTCACCCGTGCCTAATACACCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGATGAATTAAggg                                                                                                        <  1:1967015/150‑1 (MQ=255)
                                           cAGTTCACCCGTGCCTAATACACCGACCTGTAGCAATTGTGCGTTAGTAGATGGTGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGt                                                                                                       <  2:547820/150‑1 (MQ=255)
                                           cAGTTCACCCGTGCCTAATACACCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGt                                                                                                       >  1:1072340/1‑150 (MQ=255)
                                             gTTCACCCGTGCCTAATACACCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTgg                                                                                                     >  1:4921/1‑150 (MQ=255)
                                             gTTCACCCGTGCCTAATACACCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAAATATCCCCTTACCGGAGGCATTTACTCCGATAAATCCATAGACCGATGAATAAAGGGTgg                                                                                                     >  1:2489690/1‑150 (MQ=255)
                                               tCACCCGTGCCTAATACACCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCa                                                                                                   <  1:732150/150‑1 (MQ=255)
                                                  cccGTGCCTAATACACCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCAATAGCCCGTTGAATTAAGGGTGGCACCa                                                                                                >  1:1725235/1‑150 (MQ=255)
                                                   ccGTGCCTAATACACCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCAt                                                                                               >  1:1876293/1‑150 (MQ=255)
                                                       gCCTAATACACCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACt                                                                                           >  1:2222804/1‑150 (MQ=255)
                                                         cTAATACACCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGa                                                                                         >  1:1498714/1‑150 (MQ=255)
                                                         cTAATACACCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGa                                                                                         <  1:233249/150‑1 (MQ=255)
                                                          tAATACACCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGAc                                                                                        >  2:700475/1‑150 (MQ=255)
                                                           aaTACACCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACa                                                                                       >  1:216455/1‑150 (MQ=255)
                                                               cacCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTa                                                                                   <  2:2167450/150‑1 (MQ=255)
                                                                 ccGACCTGTAGCAATTGTGCGTCAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACg                                                                                 <  2:1740451/150‑1 (MQ=255)
                                                                       tGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCg                                                                           >  2:1317488/1‑150 (MQ=255)
                                                                       tGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCg                                                                           >  2:474135/1‑150 (MQ=255)
                                                                          aGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTa                                                                        >  2:2945/1‑150 (MQ=255)
                                                                           gCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTAc                                                                       <  1:289682/150‑1 (MQ=255)
                                                                                 gtgCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAAt                                                                 >  1:2127665/1‑150 (MQ=255)
                                                                                   gCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCCCCATTACTGACATTTACGATACCGCTACCAAAATTa                                                               >  2:955702/1‑150 (MQ=255)
                                                                                    cGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTAc                                                              >  2:2248084/1‑150 (MQ=255)
                                                                                     gTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTAcc                                                             >  1:1054878/1‑150 (MQ=255)
                                                                                        agtagATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACTGCTACCAAAATTACCGAc                                                          >  2:298134/1‑150 (MQ=255)
                                                                                        agtagATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGAc                                                          >  2:1092734/1‑150 (MQ=255)
                                                                                         gtagATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACa                                                         >  1:733479/1‑150 (MQ=255)
                                                                                          tagATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACAt                                                        >  1:2006760/1‑150 (MQ=255)
                                                                                           agATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACAtt                                                       >  2:1567211/1‑150 (MQ=255)
                                                                                           agATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACAtt                                                       >  2:2018302/1‑150 (MQ=255)
                                                                                            gatgatGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTg                                                      >  1:152346/1‑150 (MQ=255)
                                                                                            gatgatGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTg                                                      >  2:1340501/1‑150 (MQ=255)
                                                                                                 tgTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATAcc                                                 >  1:17400/1‑150 (MQ=255)
                                                                                                 tgTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATAcc                                                 <  1:2117251/150‑1 (MQ=255)
                                                                                                   tCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCca                                               <  2:518264/150‑1 (MQ=255)
                                                                                                    cTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCcac                                              <  2:623404/150‑1 (MQ=255)
                                                                                                    cTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCcac                                              <  1:784112/150‑1 (MQ=255)
                                                                                                     tttAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCcaca                                             <  2:1498714/150‑1 (MQ=255)
                                                                                                          aTTCCATAGACTGTCCGTTGAAATATTAACTATCCCCTTATCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTTATACCCACAGTATc                                        <  2:1057313/150‑1 (MQ=255)
                                                                                                          aTTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATc                                        <  2:1951254/150‑1 (MQ=255)
                                                                                                          aTTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATc                                        >  2:2206586/1‑150 (MQ=255)
                                                                                                           ttCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCg                                       <  1:2362756/150‑1 (MQ=255)
                                                                                                            tCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGa                                      >  2:543133/1‑150 (MQ=255)
                                                                                                            tCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGa                                      >  2:1891185/1‑150 (MQ=255)
                                                                                                              cAGAGACTGTCCGTTGAAATATTAACTATACCCTTAACTGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTAATGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAAc                                    <  2:718635/150‑1 (MQ=255)
                                                                                                                gagaCTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGACCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACAAAAATTACCGACATTGATACCCACAGTATCGAACga                                  >  2:841689/1‑150 (MQ=255)
                                                                                                                 tgaCTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACgat                                 <  1:791195/149‑1 (MQ=255)
                                                                                                                  taCTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACgatg                                <  2:1072340/149‑1 (MQ=255)
                                                                                                                   aCTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACgatga                               <  2:1240718/150‑1 (MQ=255)
                                                                                                                      gTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGt                            >  2:1821877/1‑150 (MQ=255)
                                                                                                                        ccGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCg                          <  2:1441545/150‑1 (MQ=255)
                                                                                                                            tGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGcatc                      >  2:407004/1‑150 (MQ=255)
                                                                                                                            tGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGcatc                      >  2:1549316/1‑150 (MQ=255)
                                                                                                                              aaaTATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGcatcca                    <  1:1172485/150‑1 (MQ=255)
                                                                                                                              aaaTATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGcatcca                    <  2:551169/150‑1 (MQ=255)
                                                                                                                                      aCTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGCATCCATCAGGTTc            <  2:1876293/150‑1 (MQ=255)
                                                                                                                                        tATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGCATCCATCAGGTTCAg          <  1:1317488/150‑1 (MQ=255)
                                                                                                                                           ggatttATCGGAGGAAATGCCTCCGGTAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGCATCCATCAGGTTCAgtgt       >  1:201087/5‑150 (MQ=255)
                                                                                                                                            gatttATCGGAGGAAATGCCTCCGGTAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGCATCCATCAGGTTCAGtgtg      >  2:217070/4‑150 (MQ=255)
                                                                                                                                            gatttATCGGAGGAAATGCCTCCGGTAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCAACAGTATCGAACGATGATGTCGCATCCATCAGGTTCAGtgtg      >  2:1236948/4‑150 (MQ=255)
                                                                                                                                            gatttATCGGAGGAAATGCCTCCGGTAAATCCATAGCCCGTTGAATAAACGGTGGCACCATTACTGACATTTACGATACAGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGAAGATGTCGCATCCATCAGGTTCAGtgtg      >  2:619488/4‑150 (MQ=255)
                                                                                                                                            gatttATCGGACGAAATTCCTCCGGTAAATCAATAACCAGTTGAATTAAGGGTGGCACCATTACTTACATTTACGATACAGCTACAAAAAATACCCACATAAATACCCACAGTATCGAACGAAGATGTAGCATCAATCACGTTGATtgtg      >  2:432420/4‑150 (MQ=255)
                                                                                                                                             ccTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGCATCCATCAGGTTCAGTGTGc     <  2:1415465/150‑1 (MQ=255)
                                                                                                                                              cTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGCATCCATCAGGTTCAGTGTGcc    <  2:2489690/150‑1 (MQ=255)
                                                                                                                                              cTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGCATCCATCAGGTTCAGTGTGcc    <  1:2945/150‑1 (MQ=255)
                                                                                                                                              cTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGCATCCATCAGGTTCAGTGTGcc    <  2:1421655/150‑1 (MQ=255)
                                                                                                                                               ttACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGCATCCATCAGGTTCAGTGTGCCg   <  2:546758/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
TCTGTGTATCACGCGCTTTAACTATACCTCCGGTGGTAATATTCAGTTCACCCGTGCCTAATACACCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGCATCCATCAGGTTCAGTGTGCCGA  >  NC_000913/1596346‑1596639

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