Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 208,552 T→A 7.6% L1143I (TTA→ATA)  dnaE → DNA polymerase III alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913208,5520TA 7.6% 307.8 / 14.3 92L1143I (TTA→ATA) dnaEDNA polymerase III alpha subunit
Reads supporting (aligned to +/- strand):  ref base T (41/44);  new base A (4/3);  total (45/47)
Fisher's exact test for biased strand distribution p-value = 7.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.97e-03

GACCAGCTTTTAAACCGACTCCGTCAGTCTCTGGAACCCCACCGCTCTGGGACAATTCCAGTACATCTCTACTATCAGAGGGCGGATGCACGCGCGCGGTTGCGTTTTGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTATTAAACGATCTCCGTGGCCTCATTGGTTCGGAGCAGGTGGAACTGGAGTTTGACTAATACAGGAATACTATGAGTCTGAATTTCCTTGATTTTGAACAGCCGATTGCAGAGCTGGAAGCGAAAATCGATTCTCTGACTGCG  >  NC_000913/208408‑208692
                                                                                                                                                |                                                                                                                                            
gaCCAGCTTTTAAACCGGCTCCGTCAGTCTCTGGAACCCCACCGCTCTGGGACAATTCCAGTACATCTCTACTATCAGAGGGCGGATGCACGCGCGCGGTTGCGTTTTGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTATTaaa                                                                                                                                          <  1:452109/149‑1 (MQ=255)
gaCCAGCTTTTAAACCGACTCCGTCAGTCTCTGGAACCCCACCGCTCTGGGACAATTCCAGTACATCTCTACTATCAGAGGGCGGATGCACGCGCGCGGTTGCGTTTTGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTATTaaa                                                                                                                                          <  2:42769/149‑1 (MQ=255)
gaCCAGCTTTTAAACCGACTCCGTCAGTCTCTGGAACCCCACCGCTCTGGGACAATTCCAGTACATCTCTACTATCAGAGGGCGGATGCACGCGCGCGGTTGCGTTTTGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTATTaaa                                                                                                                                          <  4:340691/149‑1 (MQ=255)
gaCCAGCTTTTAAACCGACTCCGTCAGTCTCTGGAACCCCACCGCTCTGGGACAATTCCAGTACATCTCTACTATCAGAGGGCGGATGCACGCGCGCGGTTGCGTTTTGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTATTaaa                                                                                                                                          <  2:570835/149‑1 (MQ=255)
 aCCAGCTTTTAAACCGACTCCGTCAGTCTCTGGAACCCCACCGCTCTGGGACAATTCCAGTACATCTCTACTATCAGAGGGCGGATGCACGCGCGCGGTTGCGTTTTGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTATTAAAc                                                                                                                                         <  3:441209/149‑1 (MQ=255)
  ccAGCTTTTAAACCGACTCCGTCAGTCTCTGGAACCCCACCGCTCTGGGACAATTCCAGTACATCTCTACTATCAGAGGGCGGATGCACGCGCGCGGTTGCGTTTTGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTATTAAACg                                                                                                                                        >  1:476294/1‑149 (MQ=255)
   cAGCTTTTAAACCGACTCCGTCAGTCTCTGGAACCCCACCGCTCTGGGACAATTCCAGTACATCTCTACTATCAGAGGGCGGATGCACGCGCGCGGTTGCGTTTTGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTATTAAACGa                                                                                                                                       >  3:176471/1‑149 (MQ=255)
   cAGCTTTTAAACCGACTCCGTCAGTCTCTGGAACCCCACCGCTCTGGGACAATTCCAGTACATCTCTACTATCAGAGGGCGGATGCACGCGCGCGGTTGCGTTTTGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTATTAAACGa                                                                                                                                       >  1:674512/1‑149 (MQ=255)
      cTTTTAAACCGACTCCGTCAGTCTCTGGAACCCCACCGCTCTGGGACAATTCCAGTACATCTCTACTATCAGAGGGCGGATGCACGCGCGCGGTTGCGTTTTGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTATTAAACGAtct                                                                                                                                    >  1:473066/1‑149 (MQ=255)
      cTTTTAAACCGACTCCGTCAGTCTCTGGAACCCCACCGCTCTGGGACAATTCCAGTACATCTCTACTATCAGAGGGCGGATGCACGCGCGCGGTTGCGTTTTGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTAATAAACGAtcg                                                                                                                                    <  2:640637/149‑2 (MQ=255)
          tAAACCGACTCCGTCAGTCTCTGGAACCCCACCGCTCTGGGACAATTCCAGTACATCTCTACTATCAGAGGGCGGATGCACGCGCGCGGTTGCGTTTTGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTATTAAACGATCTCCGt                                                                                                                                <  4:663921/149‑1 (MQ=255)
           aaaCCGACTCCGTCAGTCTCTGGAACCCCACCGCTCTGGGACAATTCCAGTACATCTCTACTATCAGAGGGCGGATGCACGCGCGCGGTTGCGTTTTGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTATTAAACGATCTCCGTg                                                                                                                               >  2:727021/1‑149 (MQ=255)
             aCCGACTCCGTCAGTCTCTGGAACCCCACCGCTCTGGGACAATTCCAGTACATCTCTACTATCAGAGGGCGGATGCACGCGCGCGGTTGCGTTTTGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTATTAAACGATCTCCGTGGc                                                                                                                             <  2:549485/149‑1 (MQ=255)
              ccGACTCCGTCAGTCTCTGGAACCCCACCGCTCTGGGACAATTCCAGTACATCTCTACTATCAGAGGGCGGATGCACGCGCGCGGTTGCGTTTTGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTATTAAACGATCTCCGTGGnc                                                                                                                            <  3:222586/149‑1 (MQ=255)
              ccGACTCCGTCAGTCTCTGGAACCCCACCGCTCTGGGACAATTCCAGTACATCTCTACTATCAGAGGGCGGATGCACGCGCGCGGTTGCGTTTTGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTATTAAACGATCTCCGTGGcc                                                                                                                            >  3:390772/1‑149 (MQ=255)
                 aCTCCGTCAGTCTCTGGAACCCCACCGCTCTGGGACAATTCCAGTACATCTCTACTATCAGAGGGCGGATGCACGCGCGCGGTTGCGTTTTGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTATTAAACGATCTCCGTGGCCTCa                                                                                                                         >  2:419158/1‑149 (MQ=255)
                   tCCGTCAGTCTCTGGAACCCCACCGCTCTGGGACAATTCCAGTACATCTCTACTATCAGAGGGCGGATGCACGCGCGCGGTTGCGTTTTGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTATTAAACGATCTCCGTGGCCTCAtt                                                                                                                       <  1:487677/149‑1 (MQ=255)
                      gTCAGTCTCTGGAACCCCACCGCTCTGGGACAATTCCAGTACATCTCTACTATCAGAGGGCGGATGCACGCGCGCGGTTGCGTTTTGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTATTAAACGATCTCCGTGGCCTCATTGGt                                                                                                                    >  2:706144/1‑149 (MQ=255)
                      gTCAGTCTCTGGAACCCCACCGCTCTGGGACAATTCCAGTACATCTCTACTATCAGAGGGCGGATGCACGCGCGCGGTTGCGTTTTGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTATTAAACGATCTCCGTGGCCTCATTGGt                                                                                                                    >  1:609336/1‑149 (MQ=255)
                       tCAGTCTCTGGAACCCCACCGCTCTGGGACAATTCCAGTACATCTCTACTATCAGAGGGCGGATGCACGCGCGCGGTTGCGTTTTGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTATTAAACGATCTCCGTGGCCTCATTGGtt                                                                                                                   >  1:419058/1‑149 (MQ=255)
                         aGTCTCTGGAACCCCACCGCTCTGGGACAATTCCAGTACATCTCTACTATCAGAGGGCGGATGCACGCGCGCGGTTGCGTTTTGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTATTAAACGATCTCCGTGGCCTCATTGGTTCg                                                                                                                 <  3:204411/149‑1 (MQ=255)
                           tctcTGGAACCCCACCGCTCTGTGACAATTCCAGTACATCTCTACTATCAGAGGGCGGATGCACGCGCGCGGTTGTGTTTTGGCGCGACGTGGCTCGTCTCTCCGAACGATCGTTTATTAAACGATCTCCGTGGCCTCATTGGTTCGGa                                                                                                               <  2:556939/149‑1 (MQ=255)
                           tctcTGGAACCCCACCGCTCTGGGACAATTCCAGTACATCTCTACTATCAGAGGGCGGATGCACGCGCGCGGTTGCGTTTTGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTATTAAACGATCTCCGTGGCCTCATTGGTTCGGa                                                                                                               <  3:296410/149‑1 (MQ=255)
                           tctcTGGAACCCCACCGCTCTGGGACAATTCCAGTACATCTCTACTATCAGAGGGCGGATGCACGCGCGCGGTTGCGTTTTGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTATTAAACGATCTCCGTGGCCTCATTGGTTCGGa                                                                                                               <  2:26900/149‑1 (MQ=255)
                           tctcTGGAACCCCACCGCTCTGGGACAATTCCAGTACATCTCTACTATCAGAGGGCGGATGCACGCGCGCGGTTGCGTTTTGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTATTAAACGATCTCCGTGGCCTCATTGGTTCGGa                                                                                                               <  4:489333/149‑1 (MQ=255)
                           tctcTGGAACCCCACCGCTCTGGGACAATTCCAGTACATCTCTACTATCAGAGGGCGGATGCACGCGCGCGGTTGCGTTTTGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTATTAAACGATCTCCGTGGCCTCATTGGTTCGGa                                                                                                               <  3:85502/149‑1 (MQ=255)
                           tctcTGGAACCCCACCGCTCTGGGACAATTCCAGTACATCTCTACTATCAGAGGGCGGATGCACGCGCGCGGTTGCGTTTTGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTATTAAACGATCTCCGTGGCCTCATTGGTTCGGa                                                                                                               <  4:709429/149‑1 (MQ=255)
                             tctGGAACCCCACCGCTCTGGGACAATTCCAGTACATCTCTACTATCAGAGGGCGGATGCACGCGCGCGGTTGCGTTTTGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTATTAAACGATCTCCGTGGCCTCATTGGTTCGGAGc                                                                                                             >  4:73016/1‑149 (MQ=255)
                             tctGGAACCCCACCGCTCTGGGACAATTCCAGTACATCTCTACTATCAGAGGGCGGATGCACACGCGCGGTTGCGTTTTGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTATTAAACGATCTCCGTGGCCTCATTGGTTCGGAGc                                                                                                             >  3:603016/1‑149 (MQ=255)
                              ctGGAACCCCACCGCTCTGGGACAATTCCAGTACATCTCTACTATCAGAGGGCGGATGCACGCGCGCGGTTGCGTTTTGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTATTAAACGATCTCCGTGGCCTCATTGGTTCGGAGCa                                                                                                            <  2:473066/149‑1 (MQ=255)
                                ggAACCCCACCGCTCTGGGACAATTCCAGTACATCTCTACTATCAGAGGGCGGATGCACGCGCGCGGTTGCGTTTTGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTATTAAACGATCTCCGTGGCCTCATTGGTTCGGAGCAgg                                                                                                          <  3:356818/149‑1 (MQ=255)
                                ggAACCCCACCGCTCTGGGACAATTCCAGTACATCTCTACTATCAGAGGGCGGATGCACGCGCGCGGTTGCGTTTTGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTATTAAACGATCTCCGTGGCCTCATTGGTTCGGAGCAgg                                                                                                          >  2:572764/1‑149 (MQ=255)
                                 gAACCCCACCGCTCTGGGACAATTCCAGTACATCTCTACTATCAGAGGGCGGATGCACGCGCGCGGTTGCGTTTTGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTATTAAACGATCTCCGTGGCCTCATTGGTTCGGAGCAGGt                                                                                                         <  4:139991/149‑1 (MQ=255)
                                  aaCCCCACCGCTCTGGGACAATTCCAGTACATCTCTACTATCAGAGGGCGGATGCACGCGCGCGGTTGCGTTTTGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTATTAAACGATCTCCGTGGCCTCATTGGTTCGGAGCAGGTg                                                                                                        <  1:63487/149‑1 (MQ=255)
                                     cccACCGCTCTGGGACAATTCCAGTACATCTCTACTATCAGAGGGCGGATGCACGCGCGCGGTTGCGTTTTGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTATTAAACGATCTCCGTGGCCTCATTGGTTCGGAGCAGGTGGaa                                                                                                     >  1:379304/1‑149 (MQ=255)
                                      ccACCGCTCTGGGACAATTCCAGTACATCTCTACTATCAGAGGGCGGATGCACGCGCGCGGTTGCGTTTTGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTATTAAACGATCTCCGTGGCCTCATTGGTTCGGAGCAGGTGGAAc                                                                                                    >  2:632047/1‑149 (MQ=255)
                                        aCCGCTCTGGGACAATTCCAGTACATCTCTACTATCAGAGGGCGGATGCACGCGCGCGGTTGCGTTTTGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTATTAAACGATCTCCGTGGCCTCATTGGTTCGGAGCAGGTGGAACTg                                                                                                  >  2:79590/1‑149 (MQ=255)
                                        aCCGCTCTGGGACAATTCCAGTACATCTCTACTATCAGAGGGCGGATGCACGCGCGCGGTTGCGTTTTGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTATTAAACGATCTCCGTGGCCTCATTGGTTCGGAGCAGGTGGAACTg                                                                                                  >  1:645077/1‑149 (MQ=255)
                                             tctGGGACAATTCCAGTACATCTCTACTATCAGAGGGCGGATGCACGCGCGCGGTTGCGTTTTGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTATTAAACGATCTCCGTGGCCTCATTGGTTCGGAGCAGGTGGAACTGGAGtt                                                                                             <  2:340688/149‑1 (MQ=255)
                                             tctGCGACAATTCAAGTCCATCTCTACTATCAGAAGGCGGATGCACGCGCCCGGCTACGCTTTGGCGCGTCGTGGCGTGTCTCTCCGACCGATCGTCAATTAAACCATCTCCGTGGCCTCATTGGTTCGGAGCAGGTGGGACTGGAGtt                                                                                             >  2:58569/1‑149 (MQ=255)
                                                gggACAATTCCAGTACATCTCTTCTATCAGAGGGCGGATGCACGCGCGCGGTTGCGTTTTGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTATTAAACGATCTCCGTGGCCTCATTGGTTCGGAGCAGGTGGAACTGGAGTTTGa                                                                                          <  3:706564/149‑1 (MQ=255)
                                                gggACAATTCCAGTACATCTCTACTATCAGAGGGCGGATGCACGCGCGCGGTTGCGTTTTGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTATTAAACGATCTCCGTGGCCTCATTGGTTCGGAGCAGGTGGAACTGGAGTTTGa                                                                                          <  1:419158/149‑1 (MQ=255)
                                                 ggACAATTCCAGTACATCTGTACTATCAGAGGGCGGATGCACGCGCGCGGTTGCGTTTTGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTATTAAACGATCTCCGTGGCCTCATTGGTTCGGAGCAGGTGGAACTGGAGTTTGAc                                                                                         >  3:444635/1‑149 (MQ=255)
                                                 ggACAATTCCAGTACATCTCTACTATCAGAGGGCGGATGCACGCGCGCGGTTGCGTTTTGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTATTAAACGATCTCCGTGGCCTCATTGGTTCGGAGCAGGTGGAACTGGAGTTTGAc                                                                                         <  4:176471/149‑1 (MQ=255)
                                                 ggACAATTCCAGTACATCTCTACTATCAGAGGGCGGATGCACACGCGCGGTTGCGTTTTGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTATTAAACGATCTCCGTGGCCTCATTGGTTCGGAGCAGGTGGAACTGGAGTTTGAc                                                                                         <  4:603016/149‑1 (MQ=255)
                                                   aCAATTCCAGTACATCTCTACTATCAGAGGGCGGATGCACGCGCGCGGTTGCGTTTTGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTATTAAACGATCTCCGTGGCCTCATTGGTTCGGAGCAGGTGGAACTGGAGTTTGACTa                                                                                       >  4:622165/1‑149 (MQ=255)
                                                    cAATTCCAGTACATCTCTACTATCAGAGGGCGGATGCACGCGCGCGGTTGCGTTTTGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTATTAAACGATCTCCGTGGCCTCATTGGTTCGGAGCAGGTGGAACTGGAGTTTGACTaa                                                                                      <  2:291316/149‑1 (MQ=255)
                                                     aaTTCCAGTACATCTCTACTATCAGAGGGCGGATGCACGCGCGCGGTTGCGTTTTGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTATTAAACGATCTCCGTGGCCTCATTGGTTCGGAGCAGGTGGAACTGGAGTTTGACTAAt                                                                                     >  1:732549/1‑149 (MQ=255)
                                                       ttCCAGTACATCTCTACTATCAGAGGGCGGATGCACGCGCGCGGTTGCGTTTTGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTATTAAACGATCTCCGTGGCCTCATTGGTTCGGAGCAGGTGGAACTGGAGTTTGACTAATAc                                                                                   >  2:733717/1‑149 (MQ=255)
                                                            gTACATCTCTACTATCAGAGGGCGGATGCACGCGCGCGGTTGCGTTTTGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTATTAAACGATCTCCGTGGCCTCATTGGTTCGGAGCAGGTGGAACTGGAGTTTGACTAATACAGGaa                                                                              >  1:702641/1‑149 (MQ=255)
                                                            gTACATCTCTACTATCAGAGGGCGGATGCACGCGCGCGGTTGCGTTTTGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTATTAAACGATCTCCGTGGCCTCATTGGTTCGGAGCAGGTGGAACTGGAGTTTGACTAATACAGGaa                                                                              >  4:609318/1‑149 (MQ=255)
                                                            gTACATCTCTACTATCAGAGGGCGGATGCACGCGCGCGGTTGCGTTTTGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTATTAAACGATCTCCGTGGCCTCATTGGTTCGGAGCAGGTGGAACTGGAGTTTGACTAATACAGGaa                                                                              >  4:313457/1‑149 (MQ=255)
                                                             tACATCTCTACTATCAGAGGGCGGATGCACGCGCGCGGTTGCGTTTTGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTATTAAACGATCTCCGTGGCCTCATTGGTTCGGAGCAGGTGGAACTGGAGTTTGACTAATACAGGAAt                                                                             >  2:563240/1‑149 (MQ=255)
                                                                 tctcTACTATCAGAGGGCGGATGCACGCGCGCGGTTGCGTTTTGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTATTAAACGATCTCCGTGGCCTCATTGGTTCGGAGCAGGTGGAACTGGAGTTTGACTAATACAGGAATACTa                                                                         <  1:514130/149‑1 (MQ=255)
                                                                  ctctACTATCAGAGGGCGGATGCACGCGCGCGGTTGCGTTTTGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTATTAAACGATCTCCGTGGCCTCATTGGTTCGGAGCAGGTGGAACTGGAGTTTGACTAATACAGGAATACTAt                                                                        <  3:73016/149‑1 (MQ=255)
                                                                  ctctACTATCAGAGGGCAGATGCACGCGCGCGGTTGCGTTTTGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTATTaaa                                                                                                                                          >  2:503618/1‑83 (MQ=255)
                                                                  ctctACTATCAGAGGGCAGATGCACGCGCGCGGTTGCGTTTTGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTATTaaa                                                                                                                                          <  1:503618/83‑1 (MQ=255)
                                                                   tctACTATCAGAGGGCGGATGCACGCGCGCGGTTGCGTTTTGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTATTAAACGATCTCCGTGGCCTCATTGGTTCGGAGCAGGTGGAACTGGAGTTTGACTAATACAGGAATACTATg                                                                       >  1:141828/1‑149 (MQ=255)
                                                                      actaTCAGAGGGCGGATGCACGCGCGCGGTTGCGTTTTGGCGCGACGTGGCGTGTCCCTCCGAGCGATCGTTTATTAAACGATCTCCGTGGCCTCATTGGTTCGGAGCAGGTGGAACTGGAGTTTGACTAATACAGGAATACTATGAGt                                                                    >  1:745084/1‑149 (MQ=255)
                                                                        taTCAGAGGGCGGATGCACGCGCGCGGTTGCGTTTTGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTATTAAACGATCTCCGTGGCCTCATTGGTTCGGAGCAGGTGGAACTGGAGTTTGACTAATACAGGAATACTATGAGTCt                                                                  <  1:601088/149‑1 (MQ=255)
                                                                           cAGAGGGCGGATGCACGCGCGCGGTTGCGTTTTGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTATTAAACGATCTCCGTGGCCTCATTGGTTCGGAGCAGGTGGAACTGGAGTTTGACTAATACAGGAATACTATGAGTCTGaa                                                               <  2:257746/149‑1 (MQ=255)
                                                                            agagGGCGGATGCACGCGCGCGGTTGCGTTTTGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTATTAAACGATCTCCGTGGCCTCATTGGTTCGGAGCAGGTGGAACTGGAGTTTGACTAATACAGGAATACTATGAGTCTGAAt                                                              >  3:150840/1‑149 (MQ=255)
                                                                            agagGGCGGATGCACGCGCGCGGTTGCGTTTTGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTATTAAACGATCTCCGTGGCCTCATTGGTTCGGAGCAGGTGGAACTGGAGTTTGAATAATACAGGAATACTATGAGTTTCAAt                                                              >  4:655078/1‑149 (MQ=255)
                                                                                ggCGGATGCACGCGCGCGGTTGCGTTTTGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTATTAAACGATCTCCGTGGCCTCATTGGTTCGGAGCAGGTGGAACTTGAGTTTGACTAATACAGGAATACTATGAGTCTGAATTTcc                                                          <  4:439419/149‑1 (MQ=255)
                                                                                ggCGGATGCACGCGCGCGGTTGCGTTTTGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTATTAAACGATCTCCGTGGCCTCATTGGTTCGGAGCAGGTGGAACTGGAGTTTGACTAATACAGGAATACTATGAGTCTGAATTTcc                                                          <  3:208671/149‑1 (MQ=255)
                                                                                ggCGGATGCACGCGCGCGGTTGCGTTTTGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTATTAAACGATCTCCGTGGCCTCATTGGTTCGGAGCAGGTGGAACTGGAGTTTGACTAATACAGGAATACTATGAGTCTGAATTTcc                                                          >  4:109422/1‑149 (MQ=255)
                                                                                 gCGGATGCACGCGCGCGGTTGCGTTTTGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTATTAAACGATCTCCGTGGCCTCATTGGTTCGGAGCAGGTGGAACTGGAGTTTGACTAATACAGGAATACTATGAGTCTGAATTTCCt                                                         <  1:706144/149‑1 (MQ=255)
                                                                                     aTGCACGCGCGCGGTTGCGTTTTGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTATTAAACGATCTCCGTGGCCTCATTGGTTCGGAGCAGGTGGAACTGGAGTTTGACTAATACAGGAATACTATGAGTCTGAATTTCCTTGAt                                                     >  2:331558/1‑149 (MQ=255)
                                                                                      tGCACGCGCGCGGTTGCGTTTTGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTATTAAACGATCTCCGTGGCCTCATTGGTTCGGAGCAGGTGGAACTGGAGTTTGACTAATACAGGAATACTATGAGTCTGAATTTCCTTGAtt                                                    <  4:112642/149‑1 (MQ=255)
                                                                                           gcgcgcgGTTGCGTTTTGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTATTAAACGATCTCCGTGGCCTCATTGGTTCGGAGCAGGTGGAACTGGAGTTTGACTAATACAGGAATACTATGAGTCTGAATTTCCTTGATTTTGaa                                               <  3:418358/149‑1 (MQ=255)
                                                                                            cgcgcgGTTGCGTTTTGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTATTAAACGATCTCCGTGGCCTCATTGGTTCGGAGCAGGTGGAACTGGAGTTTGACTAATACAGGAATACTATGAGTCTGAATTTCCTTGATTTTGAAc                                              <  1:54653/149‑1 (MQ=255)
                                                                                                   ttGCGTTTTGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTATTAAACGATCTCCGTGGCCTCATTGGTTCGGAGCAGGTGGAACTGGAGTTTGACTAATACAGGAATACTATGAGTCTGAATTTCCTTGATTTTGAACAGCCGAt                                       >  1:365543/1‑149 (MQ=255)
                                                                                                   ttGCGTTTTGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTATTAAACGATCTCCGTGGCCTCATTGGTTCGGAGCAGGTGGAACTGGAGTTTGACTAATACAGGAATACTATGAGTCTGAATTTCCTTGATTTTGAACAGCCGAt                                       >  2:742972/1‑149 (MQ=255)
                                                                                                       gTTTTGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTATTAAACGATCTCCGTGGCCTCATTGGTTCGGAGCAGGTGGAACTGGAGTTTGACTAATACAGGAATACTATGAGTCTGAATTTCCTTGATTTTGAACAGCCGATTGCa                                   >  2:720051/1‑149 (MQ=255)
                                                                                                        ttttGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTATTAAACGATCTCCGTGGCCTCATTGGTTCGGAGCAGGTGGAACTGGAGTTTGAc                                                                                         >  1:419759/1‑94 (MQ=255)
                                                                                                        ttttGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTATTAAACGATCTCCGTGGCCTCATTGGTTCGGAGCAGGTGGAACTGGAGTTTGAc                                                                                         <  2:419759/94‑1 (MQ=255)
                                                                                                            ggCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTATTAAACGATCTCCGTGGCCTCATTGGTTCGGAGCAGGTGGAACTGGAGTTTGACTAATACAGGAATACTATGAGTCTGAATTTCCTTGATTTTGAACAGCCGATTGCAGAGCt                              <  2:732549/149‑1 (MQ=255)
                                                                                                             gcgcGACGTGGCGTGTCTCTCCGAGCGATCGTTTATTAAACGATCTCCGTGGCCTCATTGGTTCGGAGCAGGTGGAACTGGa                                                                                                <  3:328039/82‑1 (MQ=255)
                                                                                                             gcgcGACGTGGCGTGTCTCTCCGAGCGATCGTTTATTAAACGATCTCCGTGGCCTCATTGGTTCGGAGCAGGTGGAACTGGa                                                                                                >  4:328039/1‑82 (MQ=255)
                                                                                                                 gACGTGGCGTGTCTCTCCGAGCGATCGTTTATTAAACGATCTCCGTGGCCTCATTGGTTCGGAGCAGGTGGAACTGGAGTTTGACTAATACAGGAATACTATGAGTCTGAATTTCCTTGATTTTGAACAGCCGATTGCAGAGCTGGAAg                         >  4:583602/1‑149 (MQ=255)
                                                                                                                                tcCGAGCGATCGTTTATTAAACGATCTCCGTGGCCTCATTGGTTCGGAGCAGGTGGAACTGGAGTTTGACTAATACAGGAATACTATGAGTCTGAATTTCCTTGATTTTGAACAGCCGATTGCAGAGCTGGAAGCGAAAATCGATtctc          >  4:121960/1‑149 (MQ=255)
                                                                                                                                 ccGAGCGATCGTTTATTAAACGATCTCCGTGGCCTCATTGGTTCGGAGCAGGTGGAACTGGAGTTTGACTAATACAGGAATACTATGAGTCTGAATTTCCTTGATTTTGAACAGCCGATTGCAGAGCTGGAAGCGAAAATCGATTctct         <  4:593309/149‑1 (MQ=255)
                                                                                                                                 ccGAGCGATCGTTTATTAAACGATCTCCGTGGCCTCATTGGTACGGAGCAGGTGGAACTGGAGTTTGACTAATACAGGAATACTATGAGTCTGAATTTCCTTGATTTTGAACAGCCGATTGCAGAGCTGGAAGCGAAAATCGATTctct         <  1:632047/149‑1 (MQ=255)
                                                                                                                                   gAGCGATCGTTTATTAAACGATCTCCGTGGCCTCATTGGTTCGGAGCAGGTGGAACTGGAGTTTGACTAATACAGGAATACTATGAGTCTGAATTTCCTTGATTTTGAACAGCCGATTGCAGAGCAGGAAGCGAAAATCGATTCTCTGa       >  4:528159/1‑149 (MQ=255)
                                                                                                                                    ggagATCGTTTAATAAACGATCTCCGTGGCCTCATTGGTTCGGAGCAGGTGGAACTGGAGTTTGACTAATACAGGAATACTATGAGTCTGAATTTCCTTGATTTTGAACAGCCGATTGCAGAGCTGGTAGCGAAAATCGATTCTCTGAc      >  1:625109/4‑149 (MQ=255)
                                                                                                                                      agATCGTTTAATAAACGATCTCCGTGGCCTCATTGGTTCGGAGCAGGTGGAACTGGAGTTTGACTAATACAGGAATACTATGAGTCTGAATTTCCTTGATTTTGAACAGCCGATTGCAGAGCTGGAAGCGAAAATCGATTCTCTGACTg    >  4:223477/2‑149 (MQ=255)
                                                                                                                                       gATCGTTTAATAAACGATCTCCGTGGCCTCATTGGTTCGGAGCAGGTGGAACTGGAGTTTGACTAATACAGGAATACTATGAGTCTGAATTTCCtt                                                        >  4:514790/1‑96 (MQ=255)
                                                                                                                                       gATCGTTTAATAAACGATCTCCGTGGCCTCATTGGTTCGGAGCAGGTGGAACTGGAGTTTGACTAATACAGGAATACTATGAGTCTGAATTTCCtt                                                        <  3:514790/96‑1 (MQ=255)
                                                                                                                                       gATCGTTTAATAAACGATCTCCGTGGCCTCATTGGTTCGGAGCAGGTGGAACTGGAGTTTGACTAATACAGGAATACTATGAGTCTGAATTTCCTTGATTTTGAACAGCCGATTGCAGAGCTGGAAGCGAAAATCGATTCTCTGACTGc   >  1:39779/1‑149 (MQ=255)
                                                                                                                                       gATCGTTTAATAAACGATCTCCGTGGCCTCATTGGTTCGGAGCAGGTGGAACTGGAGTTTGACTAATACAGGAATACTATGAGTCTGAATTTCCTTGATTTTGAACAGCCGATTGCAGAGCTGGAAGCGAAAATCGATTCTCTGACTGc   <  2:39779/149‑1 (MQ=255)
                                                                                                                                        aTCGTTTATTAAACGATCTCCGTGGCCTCATTGGTTCGGAGCAGGTGGAACTGGAGTTTGACTAATACAGGAATACTATGGGTCTGAATTTCCTTGATTTTGAACAGCCGATTGCAGAGCTGGAAGCGAAAATCGATTCTCTGACTGCg  <  3:528159/149‑1 (MQ=255)
                                                                                                                                        aTCGTTTATTAAACGATCTCCGTGGCCTCATTGGTTCGGAGCAGGTGGAACTGGAGTTTGACTAATACAGGAATACTATGAGTCTGAATTTCCTTGATTTTGAACAGCCGATTGCAGAGCTGGAAGCGAAAATCGATTCTCTGACTGCg  <  2:141828/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                            
GACCAGCTTTTAAACCGACTCCGTCAGTCTCTGGAACCCCACCGCTCTGGGACAATTCCAGTACATCTCTACTATCAGAGGGCGGATGCACGCGCGCGGTTGCGTTTTGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTATTAAACGATCTCCGTGGCCTCATTGGTTCGGAGCAGGTGGAACTGGAGTTTGACTAATACAGGAATACTATGAGTCTGAATTTCCTTGATTTTGAACAGCCGATTGCAGAGCTGGAAGCGAAAATCGATTCTCTGACTGCG  >  NC_000913/208408‑208692

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