Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 208,552 T→A 9.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 9.6% 379.2 / 25.8 115L1143I (TTA→ATA) dnaEDNA polymerase III alpha subunit
Reads supporting (aligned to +/- strand):  ref base T (46/58);  new base A (6/5);  total (52/63)
Fisher's exact test for biased strand distribution p-value = 5.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.44e-01

Only 100 of 120 total aligned reads displayed.

TGATGACCAGCTTTTAAACCGACTCCGTCAGTCTCTGGAACCCCACCGCTCTGGGACAATTCCAGTACATCTCTACTATCAGAGGGCGGATGCACGCGCGCGGTTGCGTTTTGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTATTAAACGATCTCCGTGGCCTCATTGGTTCGGAGCAGGTGGAACTGGAGTTTGACTAATACAGGAATACTATGAGTCTGAATTTCCTTGATTTTGAACAGCCGATTGCAGAGCTGGAAGCGAAAATCGATTCTCTGACTGCGGTTAGCC  >  NC_000913/208404‑208699
                                                                                                                                                    |                                                                                                                                                   
tgatgaCCAGCTTTTAAACCGACTCCGTCAGTCTCTGGAACCCCACCGCTCTGGGACAATTCCAGTACATCTCTACTATCAGAGGGCGGATGCACGCGCGCGGTTGCGTTTTGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTttat                                                                                                                                                     <  4:367691/149‑1 (MQ=255)
    gaCCAGCTTTTAAACCGACTCCGTCAGTCTCTGGAACCCCACCGCTCTGGGACAATTCCAGTACATCTCTACTATCAGAGGGCGGATGCACGCGCGCGGTTGCGTTTTGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTATTaa                                                                                                                                                  <  1:760553/148‑1 (MQ=255)
     aCCAGCTTTTAAACCGACTCCGTCAGTCTCTGGAACCCCACCGCTCTGGGACAATTCCAGTACATCTCTACTATCAGAGGGCGGATGCACGCGCGCGGTTGCGTTTTGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTATTaaa                                                                                                                                                 <  2:737376/148‑1 (MQ=255)
     aCCAGCTTTTAAACCGACTCCGTCAGTCTCTGGAACCCCACCGCTCTGGGACAATTCCAGTACATCTCTACTATCAGAGGGCGGATGCACGCGCGCGGTTGCGTTTTGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTATTAAAc                                                                                                                                                <  3:503102/149‑1 (MQ=255)
      ccAGCTTTTAAACCGACTCCGTCAGTCTCTGGAACCCCACCGCTCTGGGACAATTCCAGTACATCTCTACTATCAGAGGGCGGATGCACGCGCGCGGTTGCGTTTTGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTATTAAACg                                                                                                                                               >  3:432354/1‑149 (MQ=255)
      ccAGCTTTTAAACCGACTCCGTCAGTCTCTGGAACCCCACCGCTCTGGGACAATTCCAGTACATCTCTACTATCAGAGGGCGGATGCACGCGCGCGGTTGCGTTTTGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTATTAAACg                                                                                                                                               >  3:394790/1‑149 (MQ=255)
       cAGCTTTTAAACCGACTCCGTCAGTCTCTGGAACCCCACCGCTCTGGGACAATTCCAGTACATCTCTACTATCAGAGGGCGGATGCACGCGCGCGGTTGCGTTTTGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTATTAAACg                                                                                                                                               <  2:106308/148‑1 (MQ=255)
       cAGCTTTTAAACCGACTCCGTCAGTCTCTGGAACCCCACCGCTCTGGGACAATTCCAGTACATCTCTACTATCAGAGGGCGGATGCACGCGCGCGGTTGCGTTTTGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTATTAAACg                                                                                                                                               <  2:655984/148‑1 (MQ=255)
       cAGCTTTTAAACCGACTCCGTCAGTCTCTGGAACCCCACCGCTCTGGGACAATTCCAGTACATCTCTACTATCAGAGGGCGGATGCACGCGCGCGGTTGCGTTTTGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTATTAAACGa                                                                                                                                              >  4:167547/1‑149 (MQ=255)
        aGCTTTTAAACCGACTCCGTCAGTCTCTGGAACCCCACCGCTCTGGGACAATTCCAGTACATCTCTACTATCAGAGGGCGGATGCACGCGCGCGGTTGCGTTTTGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTATTAAACGa                                                                                                                                              >  2:740206/1‑148 (MQ=255)
        aGCTTTTAAACCGACTCCGTCAGTCTCTGGAACCCCACCGCTCTGGGACAATTCCAGTACATCTCTACTATCAGAGGGCGGATGCACGCGCGCGGTTGCGTTTTGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTAATAAACGa                                                                                                                                              >  1:176048/1‑148 (MQ=255)
         gCTTTTAAACCGACTCCGTCAGTCTCTGGAACCCCACCGCTCTGGGACAATTCCAGTACATCTCTACTATCAGAGGGCGGATGCACGCGCGCGGTTGCGTTTTGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTAATAAACGAt                                                                                                                                             <  2:176048/148‑1 (MQ=255)
            tttAAACCGACTCCGTCAGTCTCTGGAACCCCACCGCTCTGGGACAATTCCAGTACATCTCTACTATCAGAGGGCGGATGCACGCGCGCGGTTGCGTTTTGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTATTaaa                                                                                                                                                 <  4:760118/141‑1 (MQ=255)
            tttAAACCGACTCCGTCAGTCTCTGGAACCCCACCGCTCTGGGACAATTCCAGTACATCTCTACTATCAGAGGGCGGATGCACGCGCGCGGTTGCGTTTTGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTATTAAACGAtctc                                                                                                                                          >  2:851258/1‑148 (MQ=255)
             ttAAACCGACTCCGTCAGTCTCTGGAACCCCACCGCTCTGGGACAATTCCAGTACATCTCTACTATCAGAGGGCGGATGCACGCGCGCGGTTGCGTTTTGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTAATAAACGAtcgct                                                                                                                                         <  1:604738/148‑4 (MQ=255)
                  ccGACTCCGTCAGTCTCTGGAACCCCACCGCTCTGGGACAATTCCAGTACATCTCTACTATCAGAGGGCGGATGCACGCGCGCGGTTGCGTTTTGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTATTAAACGATCTCCGTGGc                                                                                                                                    <  2:145430/148‑1 (MQ=255)
                   cGACTCCGTCAGTCTCTGGAACCCCACCGCTCTGGGACAATTCCAGTACATCTCTACTATCAGAGGGCGGATGCACGCGCGCGGTTGCGTTTTGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTATTAAACGATCTCCGTGGcc                                                                                                                                   >  1:593529/1‑148 (MQ=255)
                      cTCCGTCAGTCTCTGGAACCCCACCGCTCTGGGACAATTCCAGTACATCTCTACTATCAGAGGGCGGATGCACGCGCGCGGTTGCGTTTTGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTATTAAACGATCTCCGTGGCCTCa                                                                                                                                >  2:285507/1‑148 (MQ=255)
                        ccGTCAGTCTCTGGAACCCCACCGCTCTGGGACAATTCCAGTACATCTCTACTATCAGAGGGCGGATGCACGCGCGCGGTTGCGTTTTGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTATTAAACGATCTCCGTGGCCTCAtt                                                                                                                              <  2:758910/148‑1 (MQ=255)
                        ccGTCAGTCTCTGGAACCCCACCGCTCTGGGACAATTCCAGTACATCTCTACTATCAGAGGGCGGATGCACGCGCGCGGTTGCGTTTTGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTATTAAACGATCTCCGTGGCCTCAtt                                                                                                                              >  2:460190/1‑148 (MQ=255)
                          gTCAGTCTCTGGAACCCCACCGCTCTGGGACAATTCCAGTACATCTCTACTATCAGAGGGCGGATGCACGCGCGCGGTTGCGTTTTGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTATTAAACGATCTCCGTGGCCTCATTGGt                                                                                                                           >  3:362488/1‑149 (MQ=255)
                           tCAGTCTCTGGAACCCCACCGCTCTGGGACAATTCCAGTACATCTCTACTATCAGAGGGCGGATGCACGCGCGCGGTTGCGTTTTGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTATTAAACGATCTCCGTGGCCTCATTGGtt                                                                                                                          <  3:522892/149‑1 (MQ=255)
                           tCAGTCTCTGGAACCCCACCGCTCTGGGACAATTCCAGTACATCTCTACTATCAGAGGGCGGATGCACGCGCGCGGTTGCGTTTTGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTATTAAACGATCTCCGTGGCCTCATTGGtt                                                                                                                          <  3:525106/149‑1 (MQ=255)
                            cAGTCTCTGGAACCCCACCGCTCTGGGACAATTCCAGTACATCTCTACTATCAGAGGGCGGATGCACGCGCGCGGTTGCGTTTTGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTATTAAACGATCTCCGTGGCCTCATTGGtt                                                                                                                          <  2:656917/148‑1 (MQ=255)
                            cAGTCTCTGGAACCCCACCGCTCTGGGACAATTCCAGTACATCTCTACTATCAGAGGGCGGATGCACGCGCGCGGTTGCGTTTTGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTATTAAACGATCTCCGTGGCCTCATTGGtt                                                                                                                          >  2:146997/1‑148 (MQ=255)
                            cAGTCTCTGGAACCCCACCGCTCTGGGACAATTCCAGTACATCTCTACTATCAGAGGGCGGATGCACGCGCGCGGTTGCGTTTTGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTATTAAACGATCTCCGTGGCCTCATTGGTTc                                                                                                                         <  4:746478/149‑1 (MQ=255)
                            cAGTCTCTGGAACCCCACCGCTCTGGGACAATTCCAGTACATCTCTACTATCAGAGGGCGGATGCACGCGCGCGGTTGCGTTTTGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTATTAAACGATCTCCGTGGCCTCATTGGTTc                                                                                                                         <  3:642109/149‑1 (MQ=255)
                             aGTCTCTGGAACCCCACCGCTCTGGGACAATTCCAGTACATCTCTACTATCAGAGGGCGGATGCACGCGCGCGGTTGCGTTTTGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTATTAAACGATCTCCGTGGCCTCATTGGTTc                                                                                                                         <  2:35383/148‑1 (MQ=255)
                               tctcTGGAACCCCACCGCTCTGGGACAATTCCAGTACATCTCTACTATCAGAGGGCGGATGCACGCGCGCGGTTGCGTTTTGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTATTAAACGATCTCCGTGGCCTCATTGGTTCGGa                                                                                                                      <  3:607842/149‑1 (MQ=255)
                                  ctGGAACCCCACCGCTCTGGGACAATTCCAGTACATCTCTACTATCAGAGGGCGGATGCACGCGCGCGGTTGCGTTTTGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTATTAAACGATCTCCGTGGCCTCa                                                                                                                                <  2:863831/136‑1 (MQ=255)
                                  ctGGAACCCCACCGCTCTGGGACAATTCCAGTACATCTCTACTATCAGAGGGCGGATGCACGCGCGCGGTTGCGTTTTGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTATTAAACGATCTCCGTGGCCTCa                                                                                                                                >  1:863831/1‑136 (MQ=255)
                                   tGGAACCCCACCGCTCTGGGACAATTCCAGTACATCTCTACTATCAGAGGGCGGATGCACGCGCGCGGTTGCGTTTTGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTATTAAACGATCTCCGTGGCCTCATTGGTTCGGAGCAg                                                                                                                  <  4:684641/149‑1 (MQ=255)
                                    ggAACCCCACCGCTCTGGGACAATTCCAGTACATCTCTACTATCAGAGGGCGGATGCACGCGCGCGGTTGCGTTTTGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTATTAAACGATCTCCGTGGCCTCATTGGTTCGGAGCAgg                                                                                                                 >  4:120675/1‑149 (MQ=255)
                                      aaCCCCACCGCTCTGGGACAATTCCAGTACATCTCTACTATCAGAGGGCGGATGCACGCGCGCGGTTGCGTTTTGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTATTAAACGATCTCCGTGGCCTCATTGGTTCGGAGCAGGt                                                                                                                >  2:548728/1‑148 (MQ=255)
                                      aaCCCCACCGCTCTGGGACAATTCCAGTACATCTCTACTATCAGAGGGCGGATGCACGCGCGCGGTTGCGTTTTGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTATTAAACGATCTCCGTGGCCTCATTGGTTCGGAGCAGGt                                                                                                                >  2:357015/1‑148 (MQ=255)
                                         cccACCGCTCTGGGACAATTCCAGTACATCTCTACTATCAGAGGGCGGATGCACGCGCGCGGTTGCGTTTTGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTATTAAACGATCTCCGTGGCCTCATTGGTTCGGAGCAGGTGGaa                                                                                                            >  3:680025/1‑149 (MQ=255)
                                         cccACCGCTCTGGGACAATTCCAGTACATCTCTACTATCAGAGGGCGGATGCACGCGCGCGGTTGCGTTTTGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTATTAAACGATCTCCGTGGCCTCATTGGTTCGGAGCAGGTGGaa                                                                                                            >  4:46231/1‑149 (MQ=255)
                                            aCCGCTCTGGGACAATTCCAGTACATCTCTACTATCAGAGGGCGGATGCACGCGCGCGGTTGCGTTTTGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTATTAAACGATCTCCGTGGCCTCATTGGTTCGGAGCAGGt                                                                                                                >  1:435600/1‑142 (MQ=255)
                                            aCCGCTCTGGGACAATTCCAGTACATCTCTACTATCAGAGGGCGGATGCACGCGCGCGGTTGCGTTTTGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTATTAAACGATCTCCGTGGCCTCATTGGTTCGGAGCAGGt                                                                                                                <  2:435600/142‑1 (MQ=255)
                                             ccGCTCTGGGACAATTCCAGTACATCTCTACTATCAGAGGGCGGATGCACGCGCGCGGTTGCGTTTTGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTATTAAACGATCTCCGTGGCCTCATTGGTTCGGAGCAGGTGGAACTg                                                                                                         >  2:306972/1‑148 (MQ=255)
                                              cGCTCTGGGACAATTCCAGTACATCTCTACTATCAGAGGGCGGATGCACGCGCGCGGTTGCGTTTTGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTATTAAACGATCTCCGTGGCCTCATTGGTTCGGAGCAGGTGGAACTgg                                                                                                        <  1:740206/148‑1 (MQ=255)
                                                ctctGGGACAATTCCAGTACATCTCTACTATCAGAGGGCGGATGCACGCGCGCGGTTGCGTTTTGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTATTAAACGATCTCCGTGGCCTCATTGGTTCGGAGCAGGTGGAACTGGAg                                                                                                      <  1:695388/148‑1 (MQ=255)
                                                 tctGGGACAATTCCAGTACATCTCTACTATCAGAGGGCGGATGCACGCGCGCGGTTGCGTTTTGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTATTAAACGATCTCCGTGGCCTCATTGGTTCGGAGCAGGTGGAACTGGAGt                                                                                                     <  1:187070/148‑1 (MQ=255)
                                                  ctGGGACAATTCCAGTACATCTCTACTATCAGAGGGCGGATGCACGCGCGCGGTTGCGTTTTGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTATTAAACGATCTCCGTGGCCTCATTGGTTCGGAGCAGGTGGAACTGGAGttt                                                                                                   <  3:113639/149‑1 (MQ=255)
                                                    gggACAATTCCAGTACATCTCTACTATCAGAGGGCGGATGCACGCGCGCGGTTGCGTTTTGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTATTAAACGATCTCCGTGGCCTCATTGGTTCGGAGCAGGTGGAACTGGAGTTTg                                                                                                  <  1:279524/148‑1 (MQ=255)
                                                    gggACAATTCCAGTACATCTCTACTATCAGAGGGCGGATGCACGCGCGCGGTTGCGTTTTGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTATTAAACGATCTCCGTGGCCTCATTGGTTCGGAGCAGGTGGAACTGGAGTTTGa                                                                                                 <  3:173013/149‑1 (MQ=255)
                                                     ggACAATTCCAGTACATCTCTACTATCAGAGGGCGGATGCACGCGCGCGGTTGCGTTTTGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTATTAAACGATCTCCGTGGCCTCATTGGTTCGGAGCAGGTGGAACTGGAGTTTGa                                                                                                 >  1:449521/1‑148 (MQ=255)
                                                      gACAATTCCAGTACATCTCTACTATCAGAGGGCGGATGCACGCGCGCGGTTGCGTTTTGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTATTAAACGATCTCCGTGGCCTCATTGGTTCGGAGCAGGTGGAACTGGAGTTTGACt                                                                                               <  3:803562/149‑1 (MQ=255)
                                                            tCCAGTACATCTCTACTATCAGAGGGCGGATGCACGCGCGCGGTTGCGTTTTGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTATTAAACGATCTCCGTGGCCTCATTGGTTCGGAGCAGGTGGAACTGGAGTTTGACTAATACa                                                                                         >  4:424345/1‑149 (MQ=255)
                                                               aGTACATCTCTACTATCAGAGGGCGGATGCACGCGCGCGGTTGCGTTTTGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTATTAAACGATCTCCGTGGCCTCATTGGTTCGGAGCAGGTGGAACTGGAGTTTGACTAATACAgg                                                                                       >  1:223305/1‑148 (MQ=255)
                                                               aGTACATCTCTACTATCAGAGGGCGGATGCACGCGCGCGGTTGCGTTTTGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTATTAAACGATCTCCGTGGCCTCATTGGTTCGGAGCAGGTGGAACTGGAGTTTGACTAATACAgg                                                                                       >  2:172471/1‑148 (MQ=255)
                                                               aGTACATCTCTACTATCAGAGGGCGGATGCACGCGCGCGGTTGCGTTTTGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTATTAAACGATCTCCGTGGCCTCATTGGTTCGGAGCAGGTGGAACTGGAGTTTGACTAATACAGGa                                                                                      <  4:278441/149‑1 (MQ=255)
                                                                  aCATCTCTACTATCAGAGGGCGGATGCACGCGCGCGGTTGCGTTTTGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTATTAAACGATCTCCGTGGCCTCATTGGTTCGGAGCAGGTGGAACTGGAGTTTGACTAATACa                                                                                         <  2:313771/143‑1 (MQ=255)
                                                                  aCATCTCTACTATCAGAGGGCGGATGCACGCGCGCGGTTGCGTTTTGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTATTAAACGATCTCCGTGGCCTCATTGGTTCGGAGCAGGTGGAACTGGAGTTTGACTAATACa                                                                                         >  1:313771/1‑143 (MQ=255)
                                                                    aTCTCTACTATCAGAGGGCGGATGCACGCGCGCGGTTGCGTTTTGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTATTAAACGATCTCCGTGGCCTCATTGGTTCGGAGCAGGTGGAACTGGAGTTTGACTAATACAGGAATACt                                                                                 >  3:133304/1‑149 (MQ=255)
                                                                              tCAGAGTGCGGATGCACGCGCGCGGTTGCGTTTTGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTATTAAACGATCTCCGTGGCCTCATTGGTTCGGAGCAGGTGGAACTGGAGTTTGACTAATACAGTAATACTATGAGTCTg                                                                        >  1:14383/1‑148 (MQ=255)
                                                                                 gagGGCGGATGCACGCGCGCGGTTGCGTTTTGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTATTAAACGATCTCCGTGGCCTCATTGGTTCGGAGCAGGTGGAACTGGAGTTTGACTAATACAGGAATACTATGAGTCTGAAtt                                                                    >  4:516675/1‑149 (MQ=255)
                                                                                    ggCGGATGCACGCGCGCGGTTGCGTTTTGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTATTAAACGATCTCCGTGGCCTCATTGGTTCGGAGCAGGTGGAACTGGAGTTTGACTAATACAGGAATACTATGAGTCTGAATTTc                                                                  <  1:914938/148‑1 (MQ=255)
                                                                                       ggATGCACGCGCGCGGTTGCGTTTTGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTATTAAACGATCTCCGTGGCCTCATTGGTTCGGAGCAGGTGGAACTGGAGTTTGACTAATACAGGAATACTa                                                                                <  3:291657/131‑1 (MQ=255)
                                                                                       ggATGCACGCGCGCGGTTGCGTTTTGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTATTAAACGATCTCCGTGGCCTCATTGGTTCGGAGCAGGTGGAACTGGAGTTTGACTAATACAGGAATACTa                                                                                >  4:291657/1‑131 (MQ=255)
                                                                                            cACGCGCGCGGTTGCGTTTTGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTATTAAACGATCTCCGTGGCCTCATTGGTTCGGAGCAGGTGGAACTGGAGTTTGACTAATACAGGAATACTATGAGTCTGAATTTCCTTGAttt                                                          <  1:549563/148‑1 (MQ=255)
                                                                                               gcgcgcgGTTGCGTTTTGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTATTAAACGATCTCCGTGGCCTCATTGGTTCGGAGCAGGTGGAACTGGAGTTTGACTAATACAGGAATACTATGAGTCTGAATTTCCTTGATTTTGa                                                       <  1:306972/148‑1 (MQ=255)
                                                                                               gcgcgcgGTTGCGTTTTGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTATAAAACGATCTCCGTGGCCTCATTGGTTCGGAGCAGGTGGAACTGGAGTTTGACTAATACAGGAATACTATGAGTCTGAATTTCCTTGATTTTGa                                                       <  2:414835/148‑1 (MQ=255)
                                                                                                   gcgGTTGCGTTTTGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTATTAAACGATCTCCGTGGCCTCATTGGTTCGGAGCAGGTGGAACTGGAGTTTGACTAATACAGGAATACTATGAGTCTGAATTTCCTTGATTTTGAACAg                                                   <  1:580629/148‑1 (MQ=255)
                                                                                                   gcgGTTGCGTTTTGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTATTAAACGATCTCCGTGGCCTCATTGGTTCGGAGCAGGTGGAACTGGAGTTTGACTAATACAGGAATACTATGAGTCTGAATTTCCTTGATTTTGAACAGc                                                  <  4:394790/149‑1 (MQ=255)
                                                                                                        tGCGTTTTGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTATTAAACGATCTCCGTGGCCTCATTGGTTCGGAGCAGGTGGAACTGGAGTTTGACTAATACAGGAATACTATGAGTCTGAATTTCCTTGATTTTGAACAGCCGAtt                                             <  4:680025/149‑1 (MQ=255)
                                                                                                        tGCGTTTTGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTATTAAACGATCTCCGTGGCCTCATTGGTTCGGAGCAGGTGGAACTGGAGTTTGACTAATACAGGAATACTATGAGTCTGAATTTCCTTGATTTTGAACAGCCGAtt                                             <  3:259542/149‑1 (MQ=255)
                                                                                                        tGCGTTTTGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTATTAAACGATCTCCGTGGCCTCATTGGTTCGGAGCAGGTGGAACTGGAGTTTGACTAATACAGGAATACTATGAGTCTGAATTTCCTTGATTTTGAACAGCCGAtt                                             >  4:325832/1‑149 (MQ=255)
                                                                                                           gTTTTGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTATTAAACGATCTCCGTGGCCTCATTGGTTCGGAGCAGGTGGAACTGGAGTTTGACTAATACAGGAATACTATGAGTCTGAATTTCCTTGATTTTGAACAGCCGATTGc                                           <  1:185358/148‑1 (MQ=255)
                                                                                                           gTTTTGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTATTAAACGATCTCCGTGGCCTCATTGGTTCGGAGCAGGTGGAACTAGAGTTTGACTAATACAGGAATACTATGAGTCTGAATTTCCTTGATTTTGAACAGCCGATTGCa                                          <  3:734952/149‑1 (MQ=255)
                                                                                                            ttttGGCGCGACGTGGCGTGTCTCTCCGGGCGATCGTTTATTAAACGATCTCCGTGGCCTCATTGGTTCGGAGCAGGTGGAACTGGAGTTTGACTAATACAGGAATACTATGAGTTTGAATTTCCTTGATTTTGAACAGCCGATTGCa                                          >  1:548208/1‑148 (MQ=255)
                                                                                                                ggCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTATTAAACGATCTCCGTGGCCTCATTGGTTCGGAGCAGGTGGAACTGGAGTTTGACTAATACAGGAATACTATGAGTCTGAATTTCCTTGATTTTGAACAGCCGATTGCAGAGc                                      <  1:172471/148‑1 (MQ=255)
                                                                                                                ggCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTATTAAACGATCTCCGTGGCCTCATTGGTTCGGAGCAGGTGGAACTGGAGTTTGACTAATACAGGAATACTATGAGTCTGAATTTCCTTGATTTTGAACAGCCGATTGCAGAGCt                                     <  3:424345/149‑1 (MQ=255)
                                                                                                                   gcgACGTGGCGTGTCTCTCCGAGCGATCGTTTATTAAACGATCTCCGTGGCCTCATTGGTTCGGAGCAGGTGGAACTGGAGTTTGACTAATACAGGAATACTATGAGTCTGAATTTCCTTGATTTTGAACAGCCGATTGCAGAGCTgg                                   <  1:357015/148‑1 (MQ=255)
                                                                                                                    cgACGTGGCGTGTCTCTCCGAGCGATCGTTTATTAAACGATCTCCGTGGCCTCATTGGTTCGGAGCAGGTGGAACTGGAGTTTGACTAATACAGGAATACTATGAGTCTGAATTTCCTTGATTTTGAACAGCCGATTGCAGAGCTGGaa                                 <  3:254044/149‑1 (MQ=255)
                                                                                                                    cgACGTGGCGTGTCTCTCCGAGCGATCGTTTATTAAACGATCTCCGTGGCCTCATTGGTTCGGAGCAGGTGGAACTGGAGTTTGACTAATACAGGAATACTATGAGTCTGAATTTCCTTGATTTTGAACAGCCGATTGCAGAGCTGGa                                  >  2:438400/1‑148 (MQ=255)
                                                                                                                    cgACGTGGCGTGTCTCTCCGAGCGATCGTTTATTAAACGATCTCCGTGGCCTCATTGGTTCGGAGCAGGTGGAACTGGAGTTTGACTAATACAGGAATACTATGAGTCTGAATTTCCTTGATTTTGAACAGCCGATTGCAGAGCTGGa                                  <  2:558939/148‑1 (MQ=255)
                                                                                                                                tctctcCGAGCGATCGTTTATTAAACGATCTCCGTGGCCTCATTGGTTCGGAGCAGGTGGAACTGGAGTTTTACTAATACAGGAATACTATGAGTCTGAATTTCCTTGATTTTGAACAGCCGATTGCAGAGCt                                     >  4:217665/1‑133 (MQ=255)
                                                                                                                                tctctcCGAGCGATCGTTTATTAAACGATCTCCGTGGCCTCATTGGTTCGGAGCAGGTGGAACTGGAGTTTGACTAATACAGGAATACTATGAGTCTGAATTTCCTTGATTTTGAACAGCCGATTGCAGAGCt                                     <  3:217665/133‑1 (MQ=255)
                                                                                                                                tctctcCGAGCGATCGTTTATTAAACGATCTCCGTGGCCTCATTGGTTCGGAGCAGGTGGAACTGGAGTTTGACTAATACAGGAATACTATGAGTCTGAATTTCCTTGATTTTGAACAGCCGATTGCAGAGCTGGAAGCGAAAATCGa                      <  1:285507/148‑1 (MQ=255)
                                                                                                                                 ctctcCGAGCGATCGTTTATTAAACGATCTCCGTGGCCTCATTGGTTCGGAGCAGGTGGAACTGGAGTTTGACTAATACAGGAATACTATGAGTCTGAATTTCCTTGATTTTGAACAGCCGATTGCAGAGCTGGAAGCGAAAATCGAt                     <  1:851258/148‑1 (MQ=255)
                                                                                                                                    tcCGAGCGATCGTTTATTAAACGATCTCCGTGGCCTCATTGGTTCGGAGCAGGTGGAACTGGAGTTTGACTAATACAGGAATACTATGAGTCTGAATTTCCTTGATTTTGAACAGCCGATTGCAGAGCTGGaa                                 >  2:796673/1‑133 (MQ=255)
                                                                                                                                    tcCGAGCGATCGTTTATTAAACGATCTCCGTGGCCTCATTGGTTCGGAGCAGGTGGAACTGGAGTTTGACTAATACAGGAATACTATGAGTCTGAATTTCCTTGATTTTGAACAGCCGATTGCAGAGCTGGaa                                 <  1:796673/133‑1 (MQ=255)
                                                                                                                                    tcCGAGCGATCGTTTATTAAACGATCTCCGTGGCCTCATTGGTTCGGAGCAGGTGGAACTGGAGTTTGACTAATACAGGAATACTATGAGTCTGAATTTCCTTGATTTTGAACAGCCGATTGCAGAGCTGGAAGCGAAAATCGATtctc                 >  3:509983/1‑149 (MQ=255)
                                                                                                                                     ccGAGCGATCGTTTATTAAACGATCTCCGTGGCCTCATTGGTTCGGAGCAGGTGGAACTGGAGTTTGACTAATACAGGAATACTATGAGTCTGAATTTCCTTGATTTTGAACAGCCGATTGCAGAGCTGGAAGCGAAAATCGATtctc                 >  1:748471/1‑148 (MQ=255)
                                                                                                                                      cGAGCGATCGTTTATTAAACGATCTCCGTGGCCTCATTGGTTCGGAGCAGGTGGAACTGGAGTTTGACTAATACAGGAATACTATGAGTCTGAATTTCCTTGATTTTGAACAGCCGATTGCAGAGCTgg                                   >  3:797518/1‑129 (MQ=255)
                                                                                                                                           gATCGTTTAATAAACGATCTCCGTGGCCTCATTGGTTCGGAGCAGGTGGAACTGGAGTTTGACTAATACAGGAATACTATGAGTCTGAATTTCCTTGATTTTGAACAGCCGATTGCAGAGCTGGAAGCGAAAATCGATTCTCTGACTGc          >  3:658982/1‑149 (MQ=255)
                                                                                                                                            aTCGTTTATTAAACGATCTCCGTGGCCTCATTGGTTCGGAGCAGGTGGAACTGGAGTTTGACTAATACAGGAATACTATGAGTCTGAATTTCCTTGATTTTGAACAGCCGATTGCAGAGCTGGAAGCGAAAATCGATTCTCTGACTGc          >  2:59950/1‑148 (MQ=255)
                                                                                                                                            aTCGTTTAATAAACGATCTCCGTGGCCTCATTGGTTCGGAGCAGGTGGAACTGGAGTTTGACTAATACAGGAATACTATGAGTCTGAATTTCCTTGATTTTGAACAGCCGATTGCAGAGCTGGAAGCGAAAATCGATTCTCTGACTGc          >  2:509895/1‑148 (MQ=255)
                                                                                                                                             tCGTTTAATAAACGATCTCCGTGGCCTCATTGGTTCGGAGCAGGTGGAACTGGAGTTTGACTAATACAGGAATACTATGAGTCTGAATTTCCTTGATTTTGAACAGCCGATTGCAGAGCTGGAAGCGAAAATCGATTCTCTGACTGCg         >  2:189476/1‑148 (MQ=255)
                                                                                                                                             tCGTTTAATAAACGATCTCCGTGGCCTCATTGGTTCGGAGCAGGTGGAACTGGAGTTTGACTAATACAGGAATACTATGAGTCTGAATTTCCTTGATTTTGAACAGCCGATTGCAGAGCTGGAAGCGAAAATCGATTCTCTGACTGCg         >  2:704727/1‑148 (MQ=255)
                                                                                                                                              cGTTTATTAAACGATCTCCGTGGCCTCATTGGTTCGGAGCAGGTGGAACTGGAGTTTGACTAATACAGGAATACTATGAGTCTGAATTTCCTTGATTTTGAACAGCCGATTGCAGAGCTGGAAGCGAAAATCGATTCTCTGACTGCGGt       <  4:133304/149‑1 (MQ=255)
                                                                                                                                               gTTTATTAAACGATCTCCGTGGCCTCATTGGTTCGGAGCAGGTGGAACTGGAGTTTGACTAATACAGGAATACTATGAGTCTGAATTTCCTTGATTTTGAACAGCCGATTGCAGAGCTGGAAGCGAAAATCGATTCTCTGACTGCGGtt      >  3:767477/1‑149 (MQ=255)
                                                                                                                                                tttATTAAACGATCTCCGTGGCCTCATTGGTTCGGAGCAGGTGGAACTGGAGTTTGACTAATACAGGAATACTATGAGTCTGAATTTCCTTGATTTTGAACAGCCGATTGCAGAGCTGGAAGCGAAAATCGATTCTCTGACTGCGGtt      >  1:775861/1‑148 (MQ=255)
                                                                                                                                                tttATTAAACGATCTCCGTGGCCTCATTGGTTCGGAGCAGGTGGAACTGGAGTTTGACTAATACAGGAATACTATGAGTCTGAATTTCCTTGATTTTGAACAGCCGATTGCAGAGCTGGAAGCGAAAATCGATTCTCTGACTGCGGTTa     >  3:662082/1‑149 (MQ=255)
                                                                                                                                                tttATTAAACGATCTCCGTGGCCTCATTGGTTCGGAGCAGGTGGAACTGGAGTTTGACTAATACAGGAATACTATGAGTCTGAATTTCCTTGATTTTGAACAGCCGATTGCAGAGCTGGAAGCGAAAATCGATTCTCTGACTGCGGTTa     >  3:36431/1‑149 (MQ=255)
                                                                                                                                                 ttattaAAGGATATCCGTGGCCTCATTGGTTCGGAGCAGGTGGAACTGGAGTTTGACTAATACAGGAATACTATGAGTCTGAATTTCCTTGATTTTGAACAGCCGATTGCAGAGCTGGAAGCGAAAATCGATTCTCTGACTGCGGTTa     >  2:689508/1‑148 (MQ=255)
                                                                                                                                                 ttattaAACGATCTCCGTGGCCTCATTGGTTCGGAGCAGGTGGAACTGGAGTTTGACTAATACAGGAATACTATGAGTCTGAATTTCCTTGATTTTGAACAGCCGATTGCAGAGCTGGAAGCGAAAATCGATTCTCTGACTGCGGTTa     >  1:552590/1‑148 (MQ=255)
                                                                                                                                                   attaAACGATCTCCGTGGCCTCATTGGTTCGGAGCAGGTGGAACTGGAGTTTGACTAATACAGGAATACTATGAGTCTGAATTTCCTTGATTTTGAACAGCCGATTGCAGAGCTGGAAGCGAAAATCGATTCTCTGACTGCGGTTAGcc  <  4:662082/149‑1 (MQ=255)
                                                                                                                                                   attaAACGATCTCCGTGGCCTCATTGGTTCGGAGCAGGTGGAACTGGAGTTTGACTAATACAGGAATACTATGAGTCTGAATTTCCTTGATTTTGAACAGCCGATTGCAGAGCTGGAAGCGAAAATCGATTCTCTGACTGCGGTTAGc   <  1:429542/148‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
TGATGACCAGCTTTTAAACCGACTCCGTCAGTCTCTGGAACCCCACCGCTCTGGGACAATTCCAGTACATCTCTACTATCAGAGGGCGGATGCACGCGCGCGGTTGCGTTTTGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTATTAAACGATCTCCGTGGCCTCATTGGTTCGGAGCAGGTGGAACTGGAGTTTGACTAATACAGGAATACTATGAGTCTGAATTTCCTTGATTTTGAACAGCCGATTGCAGAGCTGGAAGCGAAAATCGATTCTCTGACTGCGGTTAGCC  >  NC_000913/208404‑208699

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