Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,313,756 A→T 6.1% S753S (TCT→TCA infB ← translation initiation factor IF‑2

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,313,7560AT 6.1% 292.8 / 5.8 95S753S (TCT→TCAinfBtranslation initiation factor IF‑2
Reads supporting (aligned to +/- strand):  ref base A (43/46);  new base T (4/2);  total (47/48)
Fisher's exact test for biased strand distribution p-value = 4.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.30e-01

CGTCACGAACTTCCGCCAGACCGATAATCTGCTGTTTCAGTTCCGGAGACAGCATACCGCTCATCGCCGCTTTCACTTCGTCAATCAGGTTATAGATGACGGAGTAGTAACGCAGATCCAGGCTTTCCGCTTCAATCACTTTACGTGCAGAGGCATCAGCACGTACGTTAAAGCCAACCAGGATGGCGTTGGACGCCGCAGCCAGGGTGGCGTCGGTTTCGGTGATACCACCTACGCCAGAACCGATGATCTTCACTTTAACTTCGTCAGTAGACAGTTTCAGCAAGGAG  >  NC_000913/3313608‑3313897
                                                                                                                                                    |                                                                                                                                             
cGTCACGAACTTCCGCCAGACCGATAATCTGCTGTTTCAGTTCCGGAGACAGCATACCGCTCATCGCCGCTTTCACTTCGTCAATCAGGTTATAGATGACGGAGTAGTAACGCAGATCCAGGCTTTCCGCTTCAATCACTTTACGTGCa                                                                                                                                               <  2:540567/149‑1 (MQ=255)
 gTCACGAACTTCCGCCAGACCGATAATCTGCTGTTTCAGTTCCGGAGACAGCATACCGCTCATCGCCGCTTTCACTTCGTCAATCAGGTTATAGATGACGGAGTAGTAACGCAGATCCAGGCTTTCCGCTTCAATCACTTTACGTGCag                                                                                                                                              >  2:346277/1‑149 (MQ=255)
  tCACGAACTTCCGCCAGACCGATAATCTGCTGTTTCAGTTCCGGAGACAGCATACCGCTCATCGCCGCTTTCACTTCGTCAATCAGGTTATAGATGACGGAGTAGTAACGCAGATCCAGGCTTTCCGCTTCAATCACTTTACGTGCaga                                                                                                                                             >  2:1286597/1‑149 (MQ=255)
   cACGAACTTCCGCCAGACCGATAATCTGCTGTTTCAGTTCCGGAGACAGCATACCGCTCATCGCCGCTTTCACTTCGTCAATCAGGTTATAGATGACGGAGTAGTAACGCAGATCCAGGCTTTCCGCTTCAATCACTTTACGTGCagag                                                                                                                                            <  2:429039/149‑1 (MQ=255)
    aCGAACTTCCGCCAGACCGATAATCTGCTGTTTCAGTTCCGGAGACAGCATACCGCTCATCGCCGCTTTCACTTCGTCAATCAGGTTATAGATGACGGAGTAGTAACGCAGATCCAGGCTTTCCGCTTCAATCACTTTACGTGCAGAgg                                                                                                                                           >  1:1365647/1‑149 (MQ=255)
     cGAACTTCCTCCAGACCGATAATCTGCTGTTTCAGTTCCGGAGACAGCATACCGCTCATCGCCGCTTTCACTTCGTCAATCAGGTTATAGATGACGGAGTAGTAACGCAGATCCAGGCTTTCCGCTTCAATCACTTTACGTGCAGAGGc                                                                                                                                          <  1:214152/149‑1 (MQ=255)
     cGAACTTCCGCCAGACCGATAATCTGCTGTTTCAGTTCCGGAGACAGCATACCGCTCATCGCCGCTTTCACTTCGTCAATCAGGTTATAGATGACGGAGTAGTAACGCAGATCCAGGCTTTCCGCTTCAATCACTTTACGTGCAGAGGc                                                                                                                                          <  2:1536254/149‑1 (MQ=255)
      gAACTTCCGCCAGACCGATAATCTGCTGTTTCAGTTCCGGAGACAGCATACCGCTCATCGCCGCTTTCACTTCGTCAATCAGGTTATAGATGACGGAGTAGTAACGCAGATCCAGGCTTTCCGCTTCAATCACTTTACGTGCAGAGGCa                                                                                                                                         >  2:135699/1‑149 (MQ=255)
             cGCCAGACCGATAATCTGCTGTTTCAGTTCCGGAGACAGCATACCGCTCATCGCCGCTTTCACTTCGTCAATCAGGTTATAGATGACGGAGTAGTAACGCAGATCCAGGCGTTCCGCTTCAATCACTTTACGTGCAGAGGCATCAGCac                                                                                                                                  <  1:999204/149‑1 (MQ=255)
                   aCCGATAATCTGCTGTTTCAGTTCCGGAGACAGCATACCGCTCATCGCCGCTTTCACTTCGTCAATCAGGTTATAGATGACGGAGTAGTAACGCAGATCCAGGCTTTCCGCTTCAATCACTTTACGTGCAGAGGCATCAGCacgtacgt                                                                                                                            <  1:948346/149‑1 (MQ=255)
                       aTAATCTGCTGTTTCAGTTCCGGAGACAGCATACCGCTCATCGCCGCTTTCACTTCGTCACTCATTTTATAGATGACGTAGTAGTAACGCAGATCCAGGCTTTCCGCTTCAATCACTTTACGTGCAGAGGCATCAGCACGTACGTTaaa                                                                                                                        <  1:416116/149‑1 (MQ=255)
                       aTAATCTGCTGTTTCAGTTCCGGAGACAGCATACCGCTCATCGCCGCTTTCACTTCGTCAATCAGGTTATAGATGACGGAGTAGTAACGCAGATCCAGGCTTTCCGCTTCAATCACTTTACGTGCAGAGGCATCAGCACGTACGTTaaa                                                                                                                        >  2:867509/1‑149 (MQ=255)
                       aTAATCTGCTGTTTCAGTTCCGGAGACAGCATACCGCTCATCGCCGCTTTAACTTCGTCAATCAGGTTATAGATGACGGAGTAGTAACGCAGATCCAGGCTTTCCGCTTCAATCACTTTACGTGCTGATGCCTCTGCacgtaaagtgat                                                                                                                        <  2:256272/149‑8 (MQ=255)
                                  tttCAGTTCCGGAGACAGCATACCGCTCATCGCCGCTTTCACTTCGTCAATCAGGTTATAGATGACGGAGTAGTAACGCAGATCCAGGCTTTCCGCTTCAATCACTTTACGTGCAGAGGCATCAGCACGTACGTTAAAGCCAACCAGGa                                                                                                             >  2:881728/1‑149 (MQ=255)
                                  tttCAGTTCCGGAGACAGCATACCGCTCATCGCCGCTTTCACTTCGTCAATCAGGTTATAGATGACGGAGTAGTAACGCAGATCCAGGCTTTCCGCTTCAATCACTTTACGTGCAGAGGCATCAGCACGTACGTTAAAGCCAACCAGGa                                                                                                             >  2:498718/1‑149 (MQ=255)
                                  tttCAGTTCCGGAGACAGCATACCGCTCATCGCCGCTTTCACTTCGTCAATCAGGTTATAGATGACGGAGTAGTAACGCAGATCCAGGCTTTCCGCTTCAATCACTTTACGTGCAGAGGCATCAGCACGTACGTTAAAGACAACCAGGa                                                                                                             >  1:111480/1‑149 (MQ=255)
                                  tttCAGTTCCGGAGACAGCATACCGCTCATCGCCGCTTTCACTTCGTCAATCAGGTTATAGATGACGGAGTAGTAACGCAGATCCAGGCTTTCCGCTTCAATCACTTTACGTGCAGAGGCATCAGCACGTACGNTAAAGCCAACCAGGa                                                                                                             <  1:505058/149‑1 (MQ=255)
                                   ttCAGTTCCGGAGACAGCATACCGCTCATCGCCGCTTTCACTTCGTCAATCAGGTTATAGATGACGGAGTAGTAACGCAGATCCAGGCTTTCCGCTTCAATCACTTTACGTGCAGAGGCATCAGCACGTACGTTAAAGCCAACCAGGAt                                                                                                            <  2:709358/149‑1 (MQ=255)
                                    tCAGTTCCGGAGACAGCATACCGCTCATCGCCGCTTTCACTTCGTCAATCAGGTTATAGATGACGGAGTAGTAACGCAGATCCAGGCTTTCCGCTGCAATCACTTTACGTGCAGAGGCATCAGCACGTACGTTAAAGCCAACCAGGATg                                                                                                           >  1:246720/1‑149 (MQ=255)
                                    tCAGTTCCGGAGACAGCATACCGCTCATCGCCGCTTTCACTTCGTCAATCAGGTTATAGATGACGGAGTAGTAACGCAGATCCAAGCTTTCCGCTTCAATCACTTTACGTGCAGAGGCATCAGCACGTACGTTAAAGCCAACCAGGATg                                                                                                           <  1:1258248/149‑1 (MQ=255)
                                      aGTTCCGGAGACAGCATACCGCTCATCGCCGCTTTCACTTCGTCAAACAGGTTATAGATGACGGAGTAGTAACGCAGATCCAGGCTTTCCGCTTCAATCACTTTACGTGCAGAGGCATCAGCACGTACGTTAAAGCCAACCAGGATGGc                                                                                                         >  1:77965/1‑149 (MQ=255)
                                        ttCCGGAGACAGCATACCGCTCATCGCCGCTTTCACTTCGTCAATCAGGTTATAGATGACGGAGTAGTAACGCAGATCCAGGCTTTCCGCTTCAATCACTTTACGTGCAGAGGCATCAGCACGTACGTTAAAGCCAACCAGGATGGCGt                                                                                                       <  1:1116373/149‑1 (MQ=255)
                                          ccGGAGACAGCATACCGCTCATCGCCGCTTTCACTTCGTCAATCAGGTTATAGATGACGGAGTAGTAACGCAGATCCAGGCTTTCCGCTTCAATCACTTTACGTGCAGAGGCATCAGCACGTACGTTAAAGCCAACCAGGATGGCGTTg                                                                                                     <  2:993238/149‑1 (MQ=255)
                                           cGGAGACAGCATACCGCTCATCGCCGCTTTCACTTCGTCAATCAGGTTATAGATGACGGAGTAGTAACGCAGATCCAGGCTTTCCGCTTCAATCACTTTACGTGCTGATGCCTCTGCacgt                                                                                                                                >  2:1548729/1‑121 (MQ=255)
                                           cGGAGACAGCATACCGCTCATCGCCGCTTTCACTTCGTCAATCAGGTTATAGATGACGGAGTAGTAACGCAGATCCAGGCTTTCCGCTTCAATCACTTTACGTGCTGATGCCTCTGCacgt                                                                                                                                <  1:1548729/121‑1 (MQ=255)
                                             gagaCAGCATACCGCTCATCGCCGCTTTCACTTCGTCAATCAGGTTATAGATGACGGAGTAGTAACGCAGATCCAGGCTTTCCGCTTCAATCACTTTACGTGCAGAGGCATCAGCACGTACGTTAAAGCCAACCAGGATGGCGTTGGAc                                                                                                  >  1:444450/1‑149 (MQ=255)
                                             gagaCAGCATACCGCTCATCGCCGCTTTCACTTCGTCAATCAGGTTATAGATGACGGAGTAGTAACGCAGATCCAGGCTTTCCGCTTCAATCACTTTACGTGCAGAGGCATCAGCACGTACGTTAAAGCCAACCAGGATGGCGTTGGAc                                                                                                  >  2:557251/1‑149 (MQ=255)
                                               gaCAGCATACCGCTCATCGCCGCTTTCACTTCGTCAATCAGGTTATAGATGACGGAGTAGTAACGCAGATCCAGGCTTTCCGCTTCAATCACTTTACGTGCAGAGGCATCAGCACGTACGTTAAAGCCAACCAGGATGGCGTTGGAcgc                                                                                                <  1:869785/149‑1 (MQ=255)
                                                       aCCGCTCATCGCCGCTTTAACTTAGTAAATCAGGTTTTAGATGAAGGATTAGTAACGCAGATCAAGTCTTTCCGCTTCAATCACTTTACGTGCAGATGCATCAGCACGTACGTTAAAGCCAACCAGGATGGCGTTGGACGCCGCAGCCa                                                                                        <  2:978115/149‑1 (MQ=255)
                                                        ccGCTCATCGCCGCTTTCACTTCGTCAATCAGGTTATAGATGACGGAGTAGTAACGCAGATCCAGGCTTTCCGCTTCAATCACTTTACGTGCAGAGGCATCAGCACGTACGTTAAAGCCAACCAGGATGGCGTTGGACGCCGCAGCCAg                                                                                       >  2:1132878/1‑149 (MQ=255)
                                                        ccGCTCATCGCCGCTTTCACTTCGTCAATCAGGTTATAGATGACGGAGTAGTAACGCAGATCCAGGCTTTCCGCTTCAATCACTTTACGTGCAGAGGCATCAGCACGTACGTTAAAGCCAACCAGGATGGCGTTGGACGCCGCAGCCAg                                                                                       >  2:1582121/1‑149 (MQ=255)
                                                         cGCTCATCGCCGCTTTCACTTCGTCAATCAGGTTATAGATGACGGAGTAGTAACGCAGATCCAGGCTTTCCGCTTCAATCACTTTACGTGCAGAGGCATCAGCACGTACGTTAAAGCCAACCAGGATGGCGTTGGACGCCGCAGCCAgg                                                                                      >  1:1123068/1‑149 (MQ=255)
                                                          gCTCATCGCCGCTTTCACTTCGTCAATCAGGTTATAGATGACGGAGTAGTAACGCAGATCCAGGCTTTCCGCTTCAATCACTTTACGTGCAGAGGCATCAGCACGTACGTTAAAGCCAACCAGGATGGCGTTGGACGCCGCAGCCAggg                                                                                     >  2:416583/1‑149 (MQ=255)
                                                              aTCGCCGCTTTCCCTTCTTCAATCAGGTTATAGATGACGGAGTCGTAACGCAGATCCAGGCTTTCCGCTTCAATCACTTTACGTGCAGAGGCATCAGCACGTACGTTAAAGCCAACCAGGATGGCGTTGGACGCCGCAGCCAGGGTGGc                                                                                 <  1:739571/149‑1 (MQ=255)
                                                              aTCGCCGCTTTCACTTCGTCAATCAGGTTATAGATGACGGAGTAGTAACGCAGATCCAGGCTTTCCGCTTCAATCACTTTACGTGCAGAGGCATCAGCACGTACGTTAAAGCCAACCAGGATGGCGTTGGACGCCGCAGCCAGGGTGGc                                                                                 >  2:1123157/1‑149 (MQ=255)
                                                               tCGCCGCTTTCCCTTCGTCAATGAGGTTATAGATGACGGAGTAGTAACGCAGATCCAGGCTTTCCGCTTCAATCACTTTACGTGCAGAGGCATCAGCACGTACGTTTAAGCCCACCAGGATGGCGTTGGACGCCGCAGCCAGGGTGGCg                                                                                >  2:256918/1‑149 (MQ=255)
                                                                 gccgcTTTCACTTCGTCAATCAGGTTATAGATGACGGAGTAGTAACGCAGATCCAGGCTTTCCGCTTCAATCACTTTACGTGCAGAGGCATCAGCACGTACGTTAAAGCCAACCAGGATGGCGTTGGACGCCGCAGCCAGGGTGGCGTc                                                                              >  1:56109/1‑149 (MQ=255)
                                                                  ccgcTTTCACTTCGTCAATCAGGTTATAGATGACGGAGTAGTAACGCAGATCCAGGCTTTCCGCTTCAATCACTTTACGTGCAGAGGCATCAGCACGTACGTTAAAGCCAACCAGGATGGCGTTGGACGCCGCAGCCAGGGTGGCGTCg                                                                             >  2:708079/1‑149 (MQ=255)
                                                                    gcTTTCACTTCGTCAATCAGGTTATAGATGACGGAGTAGTAACGCAGATCCAGGCTTTCCGCTTCAATCACTTTACGTGCAGAGGCATCAGCACGTACGTTAAAGCCAACCAGGATGGCGTTGGACGCCGCAGCCAGGGTGGCGTCGGt                                                                           <  1:514648/149‑1 (MQ=255)
                                                                      tttCACTTCGTCAATCAGGTTATAGATGACGGAGTAGTAACGCAGATCCAGGCTTTCCGCTTCAATCACTTTACGTGCAGAGGCATCAGCACGTACGTTAAAGCCAACCAGGATGGCGTTGGACGCCGCAGCCAGGGTGGCGTCGGttt                                                                         >  2:514204/1‑149 (MQ=255)
                                                                         cACTTCGTCAATCAGGTTATAGATGACGGAGTAGTAACGCAGATCCAGGCTTTCCGCTTCAATCACTTTACGTGCAGAGGCATCAGCACGTACGTTAAAGCCAACCAGGATGGCGTTGGACGCCGCAGCCAGGGTGGCGTCGGTTTCgg                                                                      <  1:745630/149‑1 (MQ=255)
                                                                         cACTTCGTCAATCAGGTTATAGATGACGGAGTAGTAACGCAGATCCAGGCTTTCCGCTTCAATCACTTTACGTGCAGAGGCATCAGCACGTACGTTAAAGCCAACCAGGATGGCGTTGGACGCCGCAGCCAGGGTGGCGTCGGTTTCgg                                                                      <  1:752105/149‑1 (MQ=255)
                                                                          aCTTCGTCAATCAGGTTATAGATGACGGAGTAGTAACGCAGATCCAGGCTTTCCGCTTCAATCACTTTACGTGCAGAGGCATCAGCACGTACGTTAAAGCCAACCAGGATGGCGTTGGAAGCCGCAGCCAGGGTGGCGTCGGTTTCGGt                                                                     >  1:933167/1‑149 (MQ=255)
                                                                           cTTCGTCAATCAGGTTATAGATGACGGAGTAGTAACGCAGATCCAGGCTTTCCGCTTCAATCACTTTACGTGCAGAGGCATCAGCACGTACGTTAAAGCCAACCAGGATGGCGTTGGACGCCGCAGCCAGGGTGGCGTCGGTTTCGGTg                                                                    >  2:412509/1‑149 (MQ=255)
                                                                            ttCGTCAATCAGGTTATAGATGACGGAGTAGTAACGCAGATCCAGGCTTTCCGCTTCAATCACTTTACGTGCAGAGGCATCAGCACGTACGTTAAAGCCAACCAGGATGGCTTTGGACGCCGCAGCCAGGGTGGCGTCGGTTTCGGTGa                                                                   >  2:163578/1‑149 (MQ=255)
                                                                            ttCGTCAATCAGGTTATAGATGACGGAGTAGTAACGCAGATCCAGGCTTTCCGCTTCAATCACTTTACGTGCAGAGGCATCAGCACGTACGTTAAAGCCAACCAGGATGGCGTTGGACGCCGCAGCCAGGGTGGCGTCGGTTTCGGTGa                                                                   <  1:557251/149‑1 (MQ=255)
                                                                            ttCGTCAATCAGGTTATAGATGACGGAGTAGTAACGCAGATCCAGGCTTTCCGCTTCAATCACTTTACGTGCAGAGGCATCAGCACGTACGTTAAAGCCAACCAGGATGGCGTTGGACGCCGCAGCCAGGGTGGCGTCGGTTTCGGTGa                                                                   >  1:786913/1‑149 (MQ=255)
                                                                              cGTCAATCAGGTTATAGATGACGGCGTAGTAACGCAGATCCCGGCTTTCCGCTTCAATCACTTTACGTGCAGAGGCATCAGCACGTACGTTAAAGCCAACCAGGATGGCGTTGGACGCCGCAGCCAGGGTGGCGTCGGTTTCGGTGATa                                                                 <  2:1010545/149‑1 (MQ=255)
                                                                               gTCAATCAGGTTATAGATGACGGAGTAGTAACGCAGATCCAGGCTTTCCGCTTCAATCACTTTACGTGCAGAGGCATCAGCACGTACGTTAAAGCCAACCAGGATGGCGTTGGACGCCGCCGCCAGGGTGGCGTCGGTTTCGGTGATac                                                                >  2:657287/1‑149 (MQ=255)
                                                                               gTCAATCAGGTTATAGATGACGGAGTAGTAACGCAGATCCAGGCTTTCCGCTTCAATCACTTTACGTGCAGAGGCATCAGCACGTACGTTAAAGCCAACCAGGATGGCGTTGGACGCCGCAGCCAGGGTGGCGTCGGTTTCGGTGATac                                                                >  2:78680/1‑149 (MQ=255)
                                                                                  aaTCAGGTTATAGATGACGGAGTAGTAACGCAGATCCAGGCTTTCCGCTTCAATCACTTTACGTGCAGAGGCATCAGCACGTACGTTAAAGCCAACCAGGATGGCGTTGGACGCCGCAGCCAGGGTGGCGTCGGTTTCGGTGATaccac                                                             >  2:1456164/1‑149 (MQ=255)
                                                                                  aaTCAGGTTATAGATGACGGAGTAGTAACGCAGATCCAGGCTTTCCGCTTCAATCACTTTACGTGCAGAGGCATCAGCACGTACGTTAAAGCCAACCAGGATGGCGTTGGACGCCGCAGCCAGGGTGGCGTCGGTTTCGGTGATaccac                                                             <  2:1192365/149‑1 (MQ=255)
                                                                                  aaTCAGGTTATAGATGACGGAGTAGTAACGCAGATCCAGGCTTTCCGCTTCAATCACTTTACGTGCAGAGGCATCAGCACGTACGTTAAAGCCAACCAGGATGGCGTTGGACGCCGCAGCCAGGGTGGCGTCGGTTTCGGTGATaccac                                                             >  2:1074784/1‑149 (MQ=255)
                                                                                     cAGGTTATAGATGACGGAGTAGTAACGCAGATCCAGGCTTTCCGCTTCAATCACTTTACGTGCAGAGGCATCAGCACGTACGTTAAAGCCATCCAGGATGGCGTTGGACGCCGCAGCCAGGGTGGCGTCGGTTTCGGTGATACCACCTa                                                          <  2:643523/149‑1 (MQ=255)
                                                                                      aGGTTATAGATGACGGAGTAGTAACGCAGATCCAGGCTTTCCGCTTCAATCACTTTACGTGCAGAGGCATCAGCACGTACGTTAAAGCCAACCAGGATGGCGTTGGACGCCGCAGCCAGGGTGGCGTCGGTTTCGGTGATACCACCTAc                                                         <  1:446962/149‑1 (MQ=255)
                                                                                        gTTATAGATGACGGAGTAGTAACGCAGATCCAGGCTTTCCGCTTCAATCACTTTACGTGCAGAGGCATCAGCACGTACGTTAAAGCCAACCAGGATGGCGTTGGACGCCGCAGCCAGGGTGGCGTCGGTTTCGGTGATACCACCTACGc                                                       <  1:346277/149‑1 (MQ=255)
                                                                                          tataGATGACGGAGTAGTAACGCAGATCCAGGCTTTCCGCTTCAATCACTTTACGTGCAGAGGCATCAGCACGTACGTTAAAGCCAACCAGGATGGCGTTGGACGCCGCAGCCAGGGTGGCGTCGGTTTCGGTGATACCACCTACGCCa                                                     >  2:1131816/1‑149 (MQ=255)
                                                                                           ataGATGACGGAGTAGTAACGCAGATCCAGGCTTTCCGCTTCAATCACTTTACGTGCAGAGGCATCAGCACGTACGTTAAAGCCAACCAGGATGGCGTTGGACGCCGCAGCCAGGGTGGCGTCGGTTTCGGTGATACCACCTACGCCAg                                                    >  2:905712/1‑149 (MQ=255)
                                                                                                 ttcGGAGTAGTAAGGCAGATCGAGTCTTTCCTCTTCAATCACTTTACGTGCAGAGGCATCAGCACGTACGTTAAAGCCAACCAGGATGGCGTTGGACGCCGCAGCCAGGGTGGCGTCGGTTTCGGTGATACCACCTACGCCAGAACCga                                              <  1:964928/147‑1 (MQ=255)
                                                                                                   cGGAGTAGTAAGGCTGTTGCAGTCTTTCCGCTTCAATCACTTTCCGTGCAGAGGCATCTGCAGGTAGGTTATAGCCAACCAGGATGGAGTTGGACGCCTCAGCCAGGGTGGGGTCGGTTTCGGTGATACCACCTACGCCAGAACCgatg                                            <  2:256429/149‑1 (MQ=255)
                                                                                                    ggAGTAGTAACGCAGATCCAGGCTTTCCGCGTCAATCACTTTACGTGCAGAGGCATGAGCACGTACGTTAAAGCCAACCAGGATGGCGTTGGACGCCGCAGCCAGGGTGGCGTCGGTTTCGGTGATACCACCTACGCCAGAACCgatga                                           <  2:296201/149‑1 (MQ=255)
                                                                                                      tgtagtaACGCAGATCCAGGCTTTCCGCTTCAATCACTTTCCGTGCAGAGGCATCAGCACGTACGTTAAAGCCAACCAGGAGGGCGTTGGACGCCGCAGCCAGGGTGGCGTCGGTTTCGGTGATACCACCTACGCCAGAACCGATGATc                                         <  1:1272554/148‑1 (MQ=255)
                                                                                                          gtaACGCAGATCCAGGCTTTCCGCTTCAATCACTTTACGTGCAGAGGCATCAGCACGTACGTTAAAGCCAACCAGGATGGCGTTGGACGCCGCAGCCAGGGTGGCGTCGGTTTCGGTGATACCACCTACGCCAGAACCGATGATCTTCa                                     >  1:431171/1‑149 (MQ=255)
                                                                                                          gtaACGCAGATCCAGGCTTTCCGCTTCAATCACTTTACGTGCAGAGGCATCAGCACGTACGTTAAAGCCAACCAGGATGGCGTTGGACGCCGCAGCCAGGGTGGCGTCGGTTTCGGTGATACCACCTACGCCAGAACCGATGATCTTCa                                     <  1:729575/149‑1 (MQ=255)
                                                                                                          gtaACGCAGATCCAGGCTTTCCGCTTCAATCACTTTACGTGCAGAGGCATCAGCACGTACGTTAAAGCCAACCAGGATGGCGTTGGACGCCGCAGCCAGGGTGGCGTCGGTTTCGGTGATACCACCTACGCCAGAACCGATGATCTTCa                                     >  2:957516/1‑149 (MQ=255)
                                                                                                           taACGCAGATCCAGGCTTTCCGCTTCAATCACTTTACGTGCAGAGGCATCAGCACGTACGTTAAAGCCAACCAGGATGGCGTTGGACGCCGCAGCCAGGGTGGCGTCGGTTTCGGTGATACCACCTACGCCAGAACCGATGATCTTCAc                                    >  1:737852/1‑149 (MQ=255)
                                                                                                                cAGATCCAGGCTTTCCGCTTCAATCACTTTACGGGCAGAGGCATCAGCACGTACGTTAAAGCCAACCAGGATGGCGTTGGACGCCGCAGCCAGGGTGGCGTCGGTTTCGGTGATACCACCTACGCCAGAACCGATGATCTTCACTTTaa                               <  1:1131220/149‑1 (MQ=255)
                                                                                                                   aTCCAGGCTTTCCGCTTCAATCACTTTACGTGCAGAGGCATCAGCACGTACGTTAAAGCCAACCAGGATGGCGTTGGACGCCGCAGCCAGGGGGGCGTCGGTTTCGGTGATACCACCTACGCCAGAACCGATGATCTTCACTTTAACtt                            <  1:297017/149‑1 (MQ=255)
                                                                                                                    tCCAGGCTTTCCGCTTCAATCACTTTACGTGCAGAGGCATCAGCACGTACGTTAAAGCCAACCAGGATGGCGTTGGACGCCGCAGCCAGGGTGGCGTCGGTTTCGGTGATACCACCTACGCCAGAACCGATGATCTTCACTTTAACTTc                           <  1:1123157/149‑1 (MQ=255)
                                                                                                                    tCCAGGCTTTCCGCTTCAATCACTTTACGTGCAGAGGCATCAGCACGTACGGTAATGCCAACCAGGATGGCGTTGGACGCCGCAGCCAGGGTGGCGTCGGTTTCGGTGATACCACCTACGCCAGAACCGATGATCTTCACTTTAACTTc                           <  1:1150666/149‑1 (MQ=255)
                                                                                                                      cAGGCTTTCCGCTTCAATCACTTTACGTGCAGAGGCATCAGCACGTACGTTAAAGCCAACCAGGATGGCGTTGGACGCCGCAGCCAGGGTGGCGTCGGTTTCGGTGATACCACCTACGCCAGAACCGATGATCTTCACTTTAACTTCGt                         <  2:157378/149‑1 (MQ=255)
                                                                                                                       aGGCTTTCCGCTTCAATCACTTTACGTGCAGAGGCATCAGCACGTACGTTAAAGCCAACCAGGATGGCGTTGGACGCCGCAGCCAGGGTGGCGTCGGTTTCGGTGATACCACCTACGCCAGAACCGATGATCTTCACTTTAACTTCGTc                        >  2:865637/1‑149 (MQ=255)
                                                                                                                       aGGCTTTCCGCTTCAATCACTTTACGTGCAGAGGCATAAGCACGTACGTTAAAGCCAACCAGGATGGCGTTGGACGCCGCAGCCAGGGTGGCGTCGGTTTCGGTGATACCACCTACGCCAGAACCGATGATCTTCACTTTAACTTCGTc                        >  2:954318/1‑149 (MQ=255)
                                                                                                                        ggCTTTCCGCTTCAATCACTTTACGTGCAGAGGCATCAGCACGTACGTTAAAGCCAACCAGGATGGCGTTGGACGCCGCAGCCAGGGTGGCGTCGGTTTCGGTGATACCACCTACGCCAGAACCGATGATCTTCACTTTAACTTCGTCa                       <  2:558215/149‑1 (MQ=255)
                                                                                                                         gCTTTCCGCTTCAATCACTTTACGTGCAGAGGCATCAGCACGTACGTTAAAGCCAACCAGGATGGCGTTGGACGCCGCAGCCAGGGTGGCGTCGGTTTCGGTGATACCACCTACGCCAGAACCGATGATCTTCACTTTAACTTCGTCAg                      <  2:246720/149‑1 (MQ=255)
                                                                                                                         gCTTTCCGCTTCAATCACTTTACGTGCAGAGGCATCAGCACGTACGTTAAAGCCAACCAGGATGGCGTTGGACGCCGCAGCCAGGGTGGCGTCGGTTTCGGTGATACCACCTACGCCAGAACCGATGATCTTCACTTTAACTTCGTCAg                      >  1:726116/1‑149 (MQ=255)
                                                                                                                           tttCCGCTTCAATCACTTTACGTGCAGAGGCATCAGCACGTACGTTAAAGCCAACCAGGATGTCGTTGGACGCCGCAGCCAGGGTGGCGTCGGTTTCGGTGATACCAGATACGCCAGAACCGATGATCTTCACTTTAACTTCGTCAGTa                    <  2:77965/149‑1 (MQ=255)
                                                                                                                           tttCCGCTTCAATCACTTTACGTGCAGAGGCATCAGCACGTACGTTAAAGCCAACCAGGATGGCGTTGGACGCCGCAGCCAGGGTGGCGTCGGTTTCTGTGATACCACCTACGCCTGAACCGATGATCTTCACTTTAACTTCGTCAGTa                    <  2:1220632/149‑1 (MQ=255)
                                                                                                                           tttCCGCTTCAATCACTTTACGTGCAGAGGCATCAGCACGTACGTTAAAGCCAACCAGGATGGCGTTGGACGCCGCAGCCAGGGTGGCGTCGGTTTCGGTGATACCACCTACGCCAGAACCGATGATCTTCACTTTAACTTCGTCAGTa                    <  2:933167/149‑1 (MQ=255)
                                                                                                                             tCCGCTTCAATCACTTTACGTGCAGAGGCATCAGCACGTACGTTAAAGCCAACCAGGATGGCGTTGGCCGCCGCAGCCAGGGTGGCGTCGGTTTCGGTGATACCACCTACGCCAGAACCGATGATCTTCACTTTAACTTCGTCAGTAGa                  >  1:1007570/1‑149 (MQ=255)
                                                                                                                             tCCGCTTCAATCACTTTACGTGCAGAGGCAGCAGCACGTACGTTAAAGCCAACCAGGATGGCGTTGGACGCCGCAGCCAGGGTGGCGTCGGTTTCGGTGATACCACCTACGCCAGAACCGATGATCTTCACTTTAACTTCGTCAGTAGa                  <  2:1027011/149‑1 (MQ=255)
                                                                                                                               cGCTTCAATCACTTTACGTGCAGAGGCATCAGCACGTACGTTAAAGCCAACCAGGATGGCGTTGGACGCCGCAGCCAGGGTGGCGTCGGTTTCGGTGATACCACCTACGCCAGAACCGATGATCTTCACTTTAACTTCGTCAGTAGACa                <  1:498718/149‑1 (MQ=255)
                                                                                                                               cGCTTCAATCACTTTACGTGCAGAGGCATCAGCACGTACGTTAAAGCCAACCAGGATGGCGTTGGACGCCGCAGCCAGGGTGGCGTCGGTTTCGGTGATACCACCTACGCCAGAACCGATGATCTTCACTTTAACTTCGTCAGTAGACa                <  1:514204/149‑1 (MQ=255)
                                                                                                                                gCTTCAATCACTTTACGTGCAGAGGCATCAGCACGTACGTTAAAGCCAACCAGGATGGCGTTGGACGCCGCAGCCAGGGTGGCGTCGGTTTCGGTGATTCCACCTACGCCAGAACCGAGGATCTTCAGGTTAACTTCGGCAGTAgcgag               >  2:1058139/1‑145 (MQ=255)
                                                                                                                                gCTTCAATCACTTTACGTGCAGAGGCATCAGCACGTACGTTAAAGCCAACCAGGATGGCGTTGGACGCCGCAGCCAGGGTGGCGTCGGTTTCGGTGATACCACCTACGCCAGAACCGATGATCTTCACTTTAACTTCGTCAGTAGACAg               >  1:1530053/1‑149 (MQ=255)
                                                                                                                                  ttCAATCACTTTACGTGCAGAGGCATCAGCACGTACGTTAAAGCCAACCAGGATGGCGTTGGACGCCGCAGCCAGGGTGGCGTCGGTTTCGGTGATACCACCTACGCCAGAACCGATGATCTTCACTTTAACTTCGTCAGTAGACAGtt             >  2:587086/1‑149 (MQ=255)
                                                                                                                                     aaTCACTTTACGTGCAGAGGCATCAGCACGTACGTTAAAGCCAACCAGGATGGCGTTGGACGCCGCAGCCATGGTGGCGTCGGTTTCGGTGATACCACCTACGCCAGAACCGATGATCTTCACTTTAACTTCGTCAGTAGACAGTTTCa          <  2:470275/149‑1 (MQ=255)
                                                                                                                                     aaTCACTTTACGTGCAGAGGCATCAGCACGTACGTTAAAGCCAACCAGGATGGCGTTGGACGCCGCAGCCAGGGTGGCGTCGGTTTCGGTGATACCACCTACGCCAGAACCGATGATCTTCACTTTAACTTCGTCAGTAGACAGTTTCa          >  1:435951/1‑149 (MQ=255)
                                                                                                                                     aaTCACTTTACGTGCAGAGGCATCAGCACGTACGTTAAAGCCAACCAGGATGGCGTTGGACGCCGCAGCCAGGGTGGCGTCGGTTTCGGTGATACCACCTACGCCAGAACCGATGATCTTCACTTTAACTTCGTCAGTAGACAGTTTCa          <  1:867509/149‑1 (MQ=255)
                                                                                                                                      tttaacgtACGTGCTGATGCCTCTGCACGTACGTTAAAGCCAACCAGGATGGCGTTGGACGCCGCAGCCAGGGTGGCGTCGGTTTCGGTGATACCACCTACGCCAGAACCGATGATCTTCACTTTAACTTCGTCAGTAGACAGTTTCAg         >  2:689338/8‑149 (MQ=255)
                                                                                                                                      tttaacgtACGTGCTGATGCCTCTGCACGTACGTTAAAGCCAACCAGGATGGCGTTGGACGCCGCAGCCAGGGTGGCGTCGGTTTCGGTGATACCACCTACGCCAGAACCGATGATCTTCACTTTAACTTCGTCAGTAGACAGTTTCAg         >  1:891162/8‑149 (MQ=255)
                                                                                                                                      aTCACTTTACGTGCAGAGGCATCAGCACGTACGTTAAAGCCAACCAGGATGGCGTTGGACGCCGCAGCCAGGGTGGCGTCGGTTTCGGTGATACCACCTACGCCAGCACCGATGATCTTCACTTTAACTTCGTCAGTAGACAGTTTCAg         >  2:43249/1‑149 (MQ=255)
                                                                                                                                      aTCACTTTACGTGCAGAGGCATCAGCACGTACGTTAAAGCCAACCAGGATGGCGTTGGACGCCGCAGCCAGGGTGGCGTCGGTTTCGGTGATACCACCTACGCCAGAACCGATGATCTTCACTTTAACTTCGTCAGTAGACAGTTTCAg         >  2:417612/1‑149 (MQ=255)
                                                                                                                                         aCTTTACGTGCAGAGGCATCAGCACGTACGTTAAAGCCAACCAGGATGGCGTTGGACGCCGCAGCCAGGGTGGCGTCGGTTTCGGTGATACCACCTACGCCAGAACCGATGATCTTCACTTTAACTTCGTCAGTAGACAGTTTCAGCaa      <  1:503054/149‑1 (MQ=255)
                                                                                                                                         aCTTTACGTGCAGAGGCATCAGCACGTACGTTAAAGCCAACCAGGATGGCGTTGGACGAAGCAGCCAGGGTGGCGTCGGTTTCGGTGATTCCACCTACGCCAGAACCGATGATCTTCACTTTAACTTCGTCAGTAGACAGTTTCAGCaa      <  1:562542/149‑1 (MQ=255)
                                                                                                                                          cTTTACGTGCAGAGGCATCAGCACGTACGTTAAAGCCAACCAGGATGGCGTTGGACGCCGCAGCCAGGGTGGCGTCGGTTTCGGTGATACCACCTACGCCAGATCCGATGATGTTCACTTTAACTTCGTCAGTAGACAGTTTCCGCAAg     >  1:838497/1‑149 (MQ=255)
                                                                                                                                          cTTTACGTGCAGAGGCATCAGCACGTACGTTAAAGCCAACCAGGATGGCGTTGGACGCCGCAGCCAGGGGGGCGTCGGTTTCGGTGATACCACCTACGCCAGAACCGATGATCTTCACTTTAACTTCGTCAGTAGACAGTTTCAGCAAg     <  1:821442/149‑1 (MQ=255)
                                                                                                                                           cgtACGTGCTGATGCCTCTGCACGTACGTTAAAGCCAACCAGGATGGCGTTGGACGCCGCAGCCAGGGTGGCGTCGGTTTCGGTGATACCACCTACGCCAGAACCGATGATCTTCACTTTAACTTCGTCAGTAGACAGTTTCAGCAAgg    >  1:1538092/3‑149 (MQ=255)
                                                                                                                                             tACGTGCAGAGGCATCAGCACGTACGTTAAAGCCAACCAGGATGGCGTTGGACGCCGCAGCCAGGGTGGCGTCGGTTTCGGTGATACCACCTACGCCAGAACCGATGATCTTCACTTTAACTTCGTCAGTAGACAGTTTCAGCAAGGAg  <  1:78680/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                             
CGTCACGAACTTCCGCCAGACCGATAATCTGCTGTTTCAGTTCCGGAGACAGCATACCGCTCATCGCCGCTTTCACTTCGTCAATCAGGTTATAGATGACGGAGTAGTAACGCAGATCCAGGCTTTCCGCTTCAATCACTTTACGTGCAGAGGCATCAGCACGTACGTTAAAGCCAACCAGGATGGCGTTGGACGCCGCAGCCAGGGTGGCGTCGGTTTCGGTGATACCACCTACGCCAGAACCGATGATCTTCACTTTAACTTCGTCAGTAGACAGTTTCAGCAAGGAG  >  NC_000913/3313608‑3313897

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