Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,533,561 T→A 5.9% L129Q (CTG→CAG)  dnaEB → error‑prone DNA polymerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,533,5610TA 5.9% 1111.9 / 37.9 334L129Q (CTG→CAG) dnaEBerror‑prone DNA polymerase
Reads supporting (aligned to +/- strand):  ref base T (180/134);  new base A (7/13);  total (187/147)
Fisher's exact test for biased strand distribution p-value = 6.34e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

Only 100 of 352 total aligned reads displayed.

GTGCGACGGGCCCAAGCTGGTGCTGCTGGTCGAAAACCTGACGGGTTACCAGAACCTGTGCGCGCTGATCACCCGTGCCCGGCGGCGGGCAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGCCTTGTGGGTGGCCGCTGACAGTGGCGATACTGCTACCGGCGCGTGGCTGCGCAGTGTGTTCGCTGAGCGCCTGTGGCTGGCCGTGCACCTGCATCGGGGCAGTGACGATGCCGTGCGCCTGCAACGGTTGCGGGCACTGGCGGCTG  >  NC_002947/3533412‑3533710
                                                                                                                                                     |                                                                                                                                                     
 tgCGACGGGCCCAAGCTGGTGCTGCTGGTCGAAAACCTGACGGGTTACCAGAACCTGTGCGCGCTGATCACCCGTGCCCGGCGGCGGGCAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCctgctg                                                                                                                                                      <  1:608613/150‑1 (MQ=255)
   cGACGGGCCCAAGCTGGTGCTGCTGGTCGAAAACCTGACGGGTTACCAGAACCTGTGCGCGCTGATCACCCGTGCCCGGCGGCGGGCAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGc                                                                                                                                                    >  2:5482130/1‑150 (MQ=255)
     aCGGGCCCAAGCTGGTGCTGCTGGTCGAAAACCTGACGGGTTACCAGAACCTGTGCGCGCTGATCACCCGTGCCCGGCGGCGGGCAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGCCt                                                                                                                                                  <  1:1521080/150‑1 (MQ=255)
     aCGGGCCCAAGCTGGTGCTGCTGGTCGAAAACCTGACGGGTTACCAGAACCTGTGCGCGCTGATCACCCGTGCCCGGCGGCGGGCAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGCCt                                                                                                                                                  >  1:2703043/1‑150 (MQ=255)
       gggCCCAAGCTGGTGCTGCTGGTCGAAAACCTGACGGGTTACCAGAACCTGTGCGCGCTGATCACCCGTGCCCGGCGGCGGGCAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCAGCAGGCCTtg                                                                                                                                                <  2:30967/150‑1 (MQ=255)
        ggCCCAAGCTGGTGCTGCTGGTCGAAAACCTGACGGGTTACCAGAACCTGTGCGCGCTGATCACCCGTGCCCGGCGGCGGGCAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGCCTtgt                                                                                                                                               <  2:4331073/150‑1 (MQ=255)
        ggCCCAAGCTGGTGCTGCTGGTCGAAAACCTGACGGGTTACCAGAACCTGTGCGCGCTGATCACCCGTGCCCGGCGGCGGGCAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGCCTtgt                                                                                                                                               <  2:6948371/150‑1 (MQ=255)
        ggCCCAGCGGGTGCTGCTGGTCGACAACCTGACGGGTTACCTGAACCTGTGCGCGCTGATCACCCGGGCCCGGCGGCGGGCAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACAAGGCCAGCAGGCCTtggtg                                                                                                                                             <  2:6538308/150‑4 (MQ=255)
         gCCCAAGCTGGTGCTGCTGGTCGAAAACCTGACGGGTTACCAGAACCTGTGCGCGCTGATCACCCGTGCCCGGCGGCGGGCAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGCCTtgt                                                                                                                                               <  1:1918881/149‑1 (MQ=255)
          cccAAGCTGGTGCTGCTGGTCGAAAACCTGACGGGTTACCAGAACCTGTGCGCGCTGATCACCCGTGCCCGGCGGCGGGCAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGCCTTgtgg                                                                                                                                             >  1:1044450/1‑150 (MQ=255)
          cccAAGCTGGTGCTGCTGGTCGAAAACCTGACGGGTTACCAGAACCTGTGCGCGCTGATCACCCGTGCCCGGCGGCGGGCAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGCCTTgtgg                                                                                                                                             <  1:5563921/150‑1 (MQ=255)
          cccAAGCTGGTGCTGCTGGTCGAAAACCTGACGGGTTACCAGAACCTGTGCGCGCTGATCACCCGTGCCCGGCGGCGGGCAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGCCTTgtgg                                                                                                                                             >  1:6542258/1‑150 (MQ=255)
          cccAAGCTGGTGCTGCTGGTCGAAAACCTGACGGGTTACCAGAACCTGTGCGCGCTGATCACCCGTGCCCGGCGGCGGGCAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGCCTTgtgg                                                                                                                                             >  1:216068/1‑150 (MQ=255)
          cccAAGCTGGTGCTGCTGGTCGAAAACCTGACGGGTTACCAGAACCTGTGCGCGCTGATCACCCGTGCCCGGCGGCGGGCAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGCCTTgtgg                                                                                                                                             >  2:680416/1‑150 (MQ=255)
               gCTGGTGCTGCTGGTCGAAAACCTGACGGGTTACCAGAACCTGTGCGCGCTGATCACCCGTGCCCGGCGGCGGGCAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGCCTTGTGGGTGGc                                                                                                                                        >  2:5784508/1‑150 (MQ=255)
                    tgctgcTGGTCGAAAACCTGACGGGTTACCAGAACCTGTGCGCGCTGATCACCCGTGCCCGGCGGCGGGCAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGCCTTGTGGGTGGCCGCTg                                                                                                                                   >  1:5189010/1‑150 (MQ=255)
                     gctgctGGTCGAAAACCTGACGGGTTACCAGAACCTGTGCGCGCTGATCACCCGTGCCCGGCGGCGGGCAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGCCTTGTGGGTGGCCGCTGa                                                                                                                                  <  2:1220142/150‑1 (MQ=255)
                       tgctgGTCGAAAACCTGACGGGTTACCAGAACCTGTGCGCGCTGATCACCCGTGCCCGGCGGCGGGCAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGCCTTGTGGGTGGCCGCTGACa                                                                                                                                >  1:6941751/1‑150 (MQ=255)
                        gctgGTCGAAAACCTGACGGGTTACCAGAACCTGTGCGCGCTGATCACCCGTGCCCGGCGGCGGGCAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGCCTTGTGGGTGGCCGCTGACAg                                                                                                                               >  2:749010/1‑150 (MQ=255)
                          tgGTCGAAAACCTGACGGGTTACCAGAACCTGTGCGCGCTGATCACCCGTGCCCGGCGGCGGGCAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGCCTTGTGGGTGGCCGCTGACAGTg                                                                                                                             <  2:3414941/150‑1 (MQ=255)
                               gAAAACCTGACGGGTTACCAGAACCTGTGCGCGCTGATCACCCGTGCCCGGCGGCGGGCAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGCCTTGTGGGTGGCCGCTGACAGTGGCGAt                                                                                                                        >  2:3998718/1‑150 (MQ=255)
                               gAAAACCTGACGGGTTACCAGAACCTGTGCGCGCTGATCACCCGTGCCCGGCGGCGGGCAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGCCTTGTGGGTGGCCGCTGACAGTGGCGAt                                                                                                                        >  2:2209044/1‑150 (MQ=255)
                               gAAAACCTGACGGGTTACCAGAACCTGTGCGCGCTGATCACCCGTGCCCGGCGGCGGGCAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGCCTTGTGGGTGGCCGCTGACAGTGGCGAt                                                                                                                        >  2:5052672/1‑150 (MQ=255)
                                aaaaCCTGACGGGTTACCAGAACCTGTGCGCGCTGATCACCCGTGCCCGGCGGCGGGCAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGCCTTGTGGGTGGCCGCTGACAGTGGCGATa                                                                                                                       <  1:4770777/150‑1 (MQ=255)
                                aaaaCCTGACGGGTTACCAGAACCTGTGCGCGCTGATCACCCGTGCCCGGCGGCGGGCAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACAAGGCCAGCAGGCCTTggtggtg                                                                                                                                          <  2:2363668/130‑1 (MQ=255)
                                aaaaCCTGACGGGTTACCAGAACCTGTGCGCGCTGATCACCCGTGCCCGGCGGCGGGCAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACAAGGCCAGCAGGCCTTggtggtg                                                                                                                                          >  1:2363668/1‑130 (MQ=255)
                                 aaaCCTGACGGGTGACCAGAACCTGTGCGCGCTGATCACCCGTGCCCGGCGGCGGGCAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGCCTTGTGGGTGGCCGCTGACAGTGGCGATAc                                                                                                                      <  1:2498852/150‑1 (MQ=255)
                                  aaCCTGACGGGTTACCAGAACCTGTGCGCGCTGATCACCCGTGCCCGGCGGCGGGCAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGCCTTGTGGGTGGCCGCTGACAGTGGCGATACt                                                                                                                     <  1:4598766/150‑1 (MQ=255)
                                   aCCTGACGGGTTACCAGAACCTGTGCGCGCTGATCACCCGTGCCCGGCGGCGGGCAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGCCTTGTGGGTGGCCGCTGACAGTGGCGATACTg                                                                                                                    >  2:5888132/1‑150 (MQ=255)
                                     cTGACGGGTTACCAGAACCTGTGCGCGCTGATCACCCGTGCCCGGCGGCGGGCAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGCCTGGTGGGTGGCCGCTGACAGTGGCGATACTGCt                                                                                                                  >  1:1790725/1‑150 (MQ=255)
                                      tGACGGGTTACCAGAACCTGTGCGCGCTGATCACCCGTGCCCGGCGGCGGGCAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGCCTTGTGGGTGGCCGCTGACAGTGGCGATACTGCTa                                                                                                                 >  1:3727069/1‑150 (MQ=255)
                                      tGACGGGTTACCAGAACCTGTGCGCGCTGATCACCCGTGCCCGGCGGCGGGCAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGCCTTGTGGGTGGCCGCTGACAGTGGCGATACTGCTa                                                                                                                 >  2:1548983/1‑150 (MQ=255)
                                      tGACGGGTTACCAGAACCTGTGCGCGCTGATCACCCGTGCCCGGCGGCGGGCAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGCCTTGTGGGTGGCCGCTGACAGTGGCGATACTGCTa                                                                                                                 >  2:103572/1‑150 (MQ=255)
                                       gACGGGTTACCAGAACCTGTGCGCGCTGATCACCCGTGCCCGGCGGCGGGCAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGCCTTGTGGGTGGCCGCTGACAGTGGCGATACTGCTAc                                                                                                                >  1:3971559/1‑150 (MQ=255)
                                       gACGGGTTACCAGAACCTGTGCGCGCTGATCACCCGTGCCCGGCGGCGGGCAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGCCTTGTGGGTGGCCGCTGACAGTGGCGATACTGCTAc                                                                                                                >  1:6828569/1‑150 (MQ=255)
                                        aCGGGTTACCAGAACCTGTGCGCGCTGATCACCCGTGCCCGGCGGCGGGCAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGCCTTGTGGGTGGCCGCTGACAGTGGCGATACTGCTAcc                                                                                                               <  2:4537643/150‑1 (MQ=255)
                                             ttACCAGAACCTGTGCGCGCTGATCACCCGTGCCCGGCGGCGGGCAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGCCTTGTGGGTGGCCGCTGACAGTGGCGATACTGCTACCGgcgc                                                                                                          >  2:2013219/1‑150 (MQ=255)
                                               aCCAGAACCTGTGCGCGCTGATCACCCGTGCCCGGCGGCGGGCAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGCCTTGTGGGTGGCCGCTGACAGTGGCGATACTGCTACCGGCGCGt                                                                                                        >  1:4892951/1‑150 (MQ=255)
                                                ccAGAACCTGTGCGCGCTGATCACCCGTGCCCGGCGGCGGGCAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGCCTTGTGGGTGGCCGCTGACAGTGGCGATACTGCTACCGGCGCGTg                                                                                                       >  2:4106800/1‑150 (MQ=255)
                                                ccAGAACCTGTGCGCGCTGATCACCCGTGCCCGGCGGCGGGCAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGCCTTGTGGGTGGCCGCTGACAGTGGCGATACTGCTACCGGCGCGTg                                                                                                       >  2:698452/1‑150 (MQ=255)
                                                 cAGAACCTGTGCGCGCTGATCACCCGTTCCCGGCGGCGGGCAGAGAAGGGGGCGTACCAGCTGTTTAGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGCCTTGTGGGTGGCCGCTGACAGTGGCGATACTGCTACCGGCGCGTgg                                                                                                      >  1:2921399/1‑150 (MQ=255)
                                                 cAGAACCTGTGCGCGCTGATCACCCGTGCCCGGCGGCGGGCAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGCCTTGTGGGTGGCCGCTGACAGTGGCGATACTGCTACCGGCGCGTgg                                                                                                      >  1:1642128/1‑150 (MQ=255)
                                                 cAGAACCTGTGCGCGCTGATCACCCGTGCCCGGCGGCGGGCAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGCCTTGTGGGTGGCCGCTGACAGTGGCGATACTGCTACCGGCGCGTgg                                                                                                      <  2:3861138/150‑1 (MQ=255)
                                                 cAGAACCTGTGCGCGCTGATCACCCGTGCCCGGCGGCGGGCAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGCCTGGTGGGTGGCCGCTGACAGTGGCGATACTGCTACCGGCGCGTgg                                                                                                      >  2:1755020/1‑150 (MQ=255)
                                                 cAGAACCTGTGCGCGCTGATCACCCGTGCCCGGCGGCGGGCAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCCGCACCACCAAGGCCTGCTGGCCTTGTGGGTGGCCGCTGACAGTGGCGATACTGCTACCGGCGCGTgg                                                                                                      >  1:628662/1‑150 (MQ=255)
                                                   gAACCTGTGCGCGCTGATCACCCGTGCCCGGCGGCGGGCAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGCCTTGTGGGTGGCCGCTGACAGTGGCGATACTGCTACCGGCGCGTGGCt                                                                                                    <  1:5780163/150‑1 (MQ=255)
                                                   gAACCTGTGCGCGCTGATCACCCGTGCCCGGCGGCGGGCAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGCCTTGTGGGTGGCCGCTGACAGTGGCGATACTGCTACCGGCGCGTGGCt                                                                                                    >  1:3578014/1‑150 (MQ=255)
                                                        tgtgCGCGCTGATCACCCGTGCCCGGCGGCGGGCAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGCCTTGTGGGTGGCCGCTGACAGTGGCGATACTGCTACCGGCGCGTGGCTGCGCa                                                                                               <  1:539408/150‑1 (MQ=255)
                                                             gcgcTGATCACCCGTGCCCGGCGGCGGGCAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGCCTTGTGGGTGGCCGCTGACAGTGGCGATACTGCTACCGGCGCGTGGCTGCGCAgtgtg                                                                                          >  1:2215958/1‑150 (MQ=255)
                                                               gcTGATCACCCGTGCCCGGCGGCGGGCAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGCCTTGTGGGTGGCCGCTGACAGTGGCGATACTGCTACCGGCGCGTGGCTGCGCAGTGTGtt                                                                                        >  1:3420387/1‑150 (MQ=255)
                                                                cTGATCACCCGTGCCCGGCGGCGGGCAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGCCTTGTGGGTGGCCGCTGACAGTGGCGATACTGCTACCGGCGCGTGGCTGCGCAGTGTGTTc                                                                                       <  1:3144/150‑1 (MQ=255)
                                                                 tGATCACCCGTGCCCGGCGGCGGGCAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGCCTTGTGGGTGGCCGCTGACAGTGGCGATACTGCTACCGGCGCGTGGCTGCGCAGTGTGTTCg                                                                                      >  2:3518701/1‑150 (MQ=255)
                                                                 tGATCACCCGTGCCCGGCGGCGGGCAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGCCTTGTGGGTGGCCGCTGACAGTGGCGATACTGCTACCGGCGCGTGGCTGCGCAGTGTGTTCg                                                                                      >  2:6361555/1‑150 (MQ=255)
                                                                 tGATCACCCGTGCCCGGCGGCGGGCAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGCCTTGTGGGTGGCCGCTGACAGTGGCGATACTGCTACCGGCGCGTGGCTGCGCAGTGTGTTCg                                                                                      >  2:627698/1‑150 (MQ=255)
                                                                     cttctctGCCCGCCGCCGGGCAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGCCTTGTGGGTGGCCGCTGACAGTGGCGATACTGCTACCGGCGCGTGGCTGCGCAGTGTGTTCGCTGa                                                                                  >  2:5669040/7‑150 (MQ=255)
                                                                     cACCCGTGCCCGGCGGCGGGCAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGCCTTGTGGGTGGCCGCTGACAGTGGCGATACTGCTACCGGCGCGTGGCTGCGCAGTGTGTTCGCTGa                                                                                  <  1:5663214/150‑1 (MQ=255)
                                                                       cccGTGCCCGGCGGCGGGCAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGCCTTGTGGGTGGCCGCTGACAGTGGCGATACTGCTACCGGCGCGTGGCTGCGCAGTGTGTTCGCTGAgc                                                                                <  2:663367/150‑1 (MQ=255)
                                                                             cccGGCGGCGGGCAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGCCTTGTGGGTGGCCGCTGACAGTGGCGATACTGCTACCGGCGCGTGGCTGCGCAGTGTGTTCGCTGAGCGCCtgt                                                                          >  2:2264723/1‑150 (MQ=255)
                                                                               cggcggcggGCAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGCCTTGTGGGTGGCCGCTGACAGTGGCGATACTGCTACCGGCGCGTGGCTGCGCAGTGTGTTCGCTGAGCGCCTGtgg                                                                        <  1:3293513/150‑1 (MQ=255)
                                                                                     cggGCAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGCCTTGTGGGTGGCCGCTGACAGTGGCGATACTGCTACCGGCGCGTGGCTGCGCAGTGTGTTCGCTGAGCGCCTGTGGCTGGcc                                                                  <  2:3131952/150‑1 (MQ=255)
                                                                                     cggGCAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGCCTTGTGGGTGGCCGCTGACAGTGGCGATACTGCTACCGGCGCGTGGCTGCGCAGTGTGTTCGCTGAGCGCCTGTGGCTGGcc                                                                  <  1:7219296/150‑1 (MQ=255)
                                                                                     cggGCAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGCCTTGTGGGTGGCCGCTGACAGTGGCGATACTGCTACCGGCGCGTGGCTGCGCAGTGTGTTCGCTGAGCGCCTGTGGCTGGcc                                                                  <  1:680416/150‑1 (MQ=255)
                                                                                       ggCAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGCCTTGTGGGTGGCCGCTGACAGTGGCGATACTGCTACCGGCGCGTGGCTGCGCAGTGTGTTCGCTGAGCGCCTGTGGCTGGCCGt                                                                >  1:6415272/1‑150 (MQ=255)
                                                                                          agagAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGCCTTGTGGGTGGCCGCTGACAGTGGCGATACTGCTACCGGCGCGTGGCTGCGCAGTGTGTTCGCTGAGCGCCTGTGGCTGGCCGTGCa                                                             >  2:2681635/1‑150 (MQ=255)
                                                                                           gagaAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGCCTTGTGGGTGGCCGCTGACAGTGGCGATACTGCTACCGGCGCGTGGCTGCGCAGTGTGTTCGCTGAGCGCCTGTGGCTGGCCGTGCAc                                                            <  2:7364050/150‑1 (MQ=255)
                                                                                             gaAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGCCTTGTGGGTGGCCGCTGACAGTGGCGATACTGCTACCGGCGCGTGGCTGCGCAGTGTGTTCGCTGAGCGCCTGTGGCTGGCCGTGCACCt                                                          <  1:5114404/150‑1 (MQ=255)
                                                                                              aaGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGCCTTGTGGGTGGCCGCTGACAGTGGCGATACTGCTACCGGCGCGTGGCTACGCAGTGTGTTCGCTGAGCGCCTGTGGCTGGCCGTGCACCTg                                                         >  1:4352500/1‑150 (MQ=255)
                                                                                                     cGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGCCTTGTGGGTGGCCGCTGACAGTGGCGATACTGCTACCGGCGCGTGGCTGCGCAGTGTGTTCGCTGAGCGCCTGTGGCTGGCCGTGCACCTGCATCggg                                                  <  2:3267930/150‑1 (MQ=255)
                                                                                                      gTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGCCTTGTGGGTGGCCGCTGACAGTGGCGATACTGCTACCGGCGCGTGGCTGCGCAGTGTGTTCGCTGAGCGCCTGTGGCTGGCCGTGCACCTGCATCgggg                                                 <  1:812703/150‑1 (MQ=255)
                                                                                                              tGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGCCTTGTGGGTGGCCGCTGACAGTGGCGATACTGCTACCGGCGCGTGGCTGCGCAGTGTGTTCGCTGAGCGCCTGTGGCTGGCCGTGCACCTGCATCGGGGCAGTGACg                                         >  1:407224/1‑150 (MQ=255)
                                                                                                                   cGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGCCTTGTGGGTGGCCGCTGACAGTGGCGATACTGCTACCGGCGCGTGGCTGCGCAGTGTGTTCGCTGAGCGCCTGTGGCTGGCCGTGCACCTGCATCGGGGCAGTGACGATGcc                                    <  2:5595531/150‑1 (MQ=255)
                                                                                                                     tgatgaCTTGCTGCTGCACCACCAAGGCCTGCTGGCCTTGTGGGTGGCCGCTGACAGTGGCGATACTGCTACCGGCGCGTGGCTGCGCAGTGTGTTCGCTGAg                                                                                 >  1:2197605/1‑103 (MQ=255)
                                                                                                                     tgatgaCTTGCTGCTGCACCACCAAGGCCTGCTGGCCTTGTGGGTGGCCGCTGACAGTGGCGATACTGCTACCGGCGCGTGGCTGCGCAGTGTGTTCGCTGAGCGCCTGTGGCTGGCCGTGCACCTGCATCGGGGCAGTGACGATGCCGt                                  >  1:7100028/1‑150 (MQ=255)
                                                                                                                     tgatgaCTTGCTGCTGCACCACCAAGGCCTGCTGGCCTTGTGGGTGGCCGCTGACAGTGGCGATACTGCTACCGGCGCGTGGCTGCGCAGTGTGTTCGCTGAGCGCCTGTGGCTGGCCGTGCACCTGCATCGGGGCAGTGACGATGCCGt                                  >  2:3680325/1‑150 (MQ=255)
                                                                                                                        tgaCTTGCTGCTGCACCACCAAGGCCTGCTGGCCTTGTGGGTGGCCGCTGACAGTGGCGATACTGCTACCGGCGCGTGGCTGCGCAGTGTGTTCGCTGAGCGCCTGTGGCTGGCCGTGCACCTGCATCGGGGCAGTGACGATGCCGTgcg                               >  1:5373226/1‑150 (MQ=255)
                                                                                                                         gaCTTGCTGCTGCACCACCAAGGCCTGCTGGCCTTGTGGGTGGCCGCTGACAGTGGCGATACTGCTACCGGCGCGTGGCTGCGCAGTGTGTTCGCTGAGCGCCTGTGGCTGGCCGTGCACCTGCATCGGGGCAGTGACGATGCCGTgcgc                              >  2:2466204/1‑150 (MQ=255)
                                                                                                                         gaCTTGCTGCTGCACCACCAAGGCCTGCTGGCCTTGTGGGTGGCCGCTGACAGTGGCGATACTGCTACCGGCGCGTGGCTGCGCAGTGTGTTCGCTGAGCGCCTGTGGCTGGCCGTGCACCTGCATCGGGGCAGTGACGATGCCGTgcgc                              >  1:552688/1‑150 (MQ=255)
                                                                                                                            ttGCTGCTGCACCACCAAGGCCTGCTGGCCTTGTGGGTGGCCGCTGACAGTGGCGATACTGCTACCGGCGCGTGGCTGCGCAGTGTGTTCGCTGAGCGCCTGTGGCTGGCCGTGCACCTGCATCGGGGCAGTGACGATGCCGTGCGCCTg                           >  1:2638458/1‑150 (MQ=255)
                                                                                                                              gctgctgcACCACCAAGGCCTGCTGGCCTTGTGGGTGGCCGCTGACAGTGGCGATACTGCTACCGGCGCGTGGCTGCGCAGTGTGTTCGCTGAGCGCCTGTGGCTGGCCGTGCACCTGCATCggg                                                  <  2:1319493/125‑1 (MQ=255)
                                                                                                                                 gctgcACCACCAAGGCCTGCTGGCCTTGTGGGTGGCCGCTGACAGTGGCGATACTGCTACCGGCGCGTGGCTGCGCAGTGTGTTCGCTGAGCGCCTGTGGCTGGCCGTGCACCTGCATCGGGGCAGTGACGATGCCGTGCGCCTGCAACg                      <  2:4126228/150‑1 (MQ=255)
                                                                                                                                   tgcACCACCAAGGCCTGCTGGCCTTGTGGGTGGCCGCTGACAGTGGCGATACTGCTACCGGCGCGTGGCTGCGCAGTGTGTTCGCTGAGCGCCTGTGGCTGGCCGTGCACCTGCATCGGGGCAGTGACGATGCCGTGCGCCTGCAACGGt                    >  1:4472891/1‑150 (MQ=255)
                                                                                                                                     caccacCAAGGCCTGCTGGCCTTGTGGGTGGCCGCTGACAGTGGCGATACTGCTACCGGCGCGTGGCTGCGCAGTGTGTTCGCTGAGCGCCTGTGGCTGGCCGTGCACCTGCATCGGGGCAGTGACGATGCCGTGCGCCTGCAACGGTTg                  <  1:4246760/150‑1 (MQ=255)
                                                                                                                                      accaccAAGGCCTGCTGGCCTTGTGGGTGGCCGCTGACAGTGGCGATACTGCTACCGGCGCGt                                                                                                        <  2:5468336/63‑1 (MQ=255)
                                                                                                                                      accaccAAGGCCTGCTGGCCTTGTGGGTGGCCGCTGACAGTGGCGATACTGCTACCGGCGCGTGGCTGCGCAGTGTGTTCGCTGAGCGCCTGTGGCTGGCCGTGCACCTGCATCGGGGCAGTGACGATGCCGTGCGCCTGCAACGGTTGc                 <  1:3126759/150‑1 (MQ=255)
                                                                                                                                       ccaccaAGGCCTGCTGGCCTTGTGGGTGGCCGCTGACAGTGGCGATACTGCTACCGGCGCGTGGCTGCGCAGTGTGTTCGCTGAGCGCCTGTGGCTGGCCGTGCACCTGCATCGGGGCAGTGACGATGCCGTGCGCCTGCAACGGTTGCg                >  1:7403505/1‑150 (MQ=255)
                                                                                                                                        caccaAGGCCTGCTGGCCTTGTGGGTGGCCGCTGACAGTGGCGATACTGCTACCGGCGCGTGGCTGCGCAGTGTGTTCGCTGAGCGCCTGTGGCTGGCCGTGCACCTGCATCGGGGCAGTGACGATGCCGTGCGCCTGCAACGGTTGCgg               <  1:2059591/150‑1 (MQ=255)
                                                                                                                                        caccaAGGCCTGCTGGCCTTGTGGGTGGCCGCTGACAGTGGCGATACTGCTACCGGCGCGTGGCTGCGCAGTGTGTTCGCTGAGCGCCTGTGGCTGGCCGTGCACCTGCATCGGGGCAGTGACGATGCCGTGCGCCTGCAACGGTTGCgg               <  2:2106684/150‑1 (MQ=255)
                                                                                                                                         accaAGGCCTGCTGGCCTTGTGGGTGGCCGCTGACAGTGGCGATACTGCTACCGGCGCGTGGCTGCGCAGTGTGTTCGCTGAGCGCCTGTGGCTGGCCGTGCACCTGCATCGGGGCAGTGACGATGCCGTGCGCCTGCAACGGTTGCggg              <  2:4043013/150‑1 (MQ=255)
                                                                                                                                          acaAGGCCAGCAGGCCTTGTGGGTGGCCGCTGACAGTGGCGATACTGCTACCGGCGCGTGGCTGCGCAGTGTGTTCGCTGAGCGCCt                                                                            >  1:6000712/2‑87 (MQ=255)
                                                                                                                                           caAGGCCTGCTGGCCTTGTGGGTGGCCGCTGACAGTGGCGATACTGCTACCGGCGCGTGGCTGCGCAGTGTGTTCGCTGAGCGCCTGTGGCTGGCCGTGCACCTGCATCGGGGCAGTGACGATGCCGTGCGCCTGCAACGGTTGCGGGCa            >  2:1595622/1‑150 (MQ=255)
                                                                                                                                           caAGGCCAGCAGGCCTTGGTGGTGGCCGCTGACAGTGGCGATACTGCTACCGGCGCGTGGCTGCGCAGTGTGTTCGCTGAGCGCCTGTGGCTGGCCGTGCACCTGCATCGGGGCAGTGACGATGCCGTGCGCCTGCAACGGTTGCGGGCa            >  1:6653009/1‑150 (MQ=255)
                                                                                                                                               gCCTGCTGGCCTTGTGGGTGGCCGCTGACAGTGGCGATACTGCTACCGGCGCGTGGCTGCGCAGTGTGTTCGCTGAGCGCCTGTGGCTGGCCGTGCACCTGCATCGGGGCAGTGACGATGCCGTGCGCCTGCAACGGTTGCGGGCACTgg        >  2:5351466/1‑150 (MQ=255)
                                                                                                                                                 ctgctgGCCTTGTGGGTGGCCGCTGACAGTGGCGATACTGCTACCGGCGCGTGGCTGCGCAGTGTGTTCGCTGAGCGCCTGTGGCTGGCCGTGCACCTGCATCGGGGCAGTGACGATGCCGTGCGCCTGCAACGGTTGCGGGCACTggcg      >  1:6971424/1‑150 (MQ=255)
                                                                                                                                                  tgctgGCCTTGTGGGTGGCCGCTGACAGTGGCGATACTGCTACCGGCGCGTGGCTGCGCAGTGTGTTCGCTGAGCGCCTGTGGCTGGCCGTGCACCTGCATCGGGGCAGTGACGATGCCGTGCGCCTGCAACGGTTGCGGGCACTggcgg     >  2:1680717/1‑150 (MQ=255)
                                                                                                                                                  tgctgGCCTTGTGGGTGGCCGCTGACAGTGGCGATACTGCTACCGGCGCGTGGCTGCGCAGTGTGTTCGCTGAGCGCCTGTGGCTGGCCGTGCACCTGCATCGGGGCAGTGACGATGCCGTGCGCCTGCAACGGTTGCGGGCACTggcgg     >  1:2776701/1‑150 (MQ=255)
                                                                                                                                                  tgctgGCCTTGTGGGTGGCCGCTGACAGTGGCGATACTGCTACCGGCGCGTGGCTGCGCAGTGTGTTCGCTGAGCGCCTGTGGCTGGCCGTGCACCTGCATCGGGGCAGTGACGATGCCGTGCGCCTGCAACGGTTGCGGGCACTggcgg     >  2:6064282/1‑150 (MQ=255)
                                                                                                                                                   gctgGCCTTGTGGGTGGCCGCTGACAGTGGCGATACTGCTACCGGCGCGTGGCTGCGCAGTGTGTTCGCTGAGCGCCTGTGGCTGGCCGTGCACCTGCATCGGGGCAGTGACGATGCCGTGCGCCTGCAACGGTTGCGGGCACTggcggc    <  2:7157184/150‑1 (MQ=255)
                                                                                                                                                    ctgGCCTTGTGGGTGGCCGCTGACAGTGGCGATACTGCTACCGGCGCGTGGCTGCGCAGTGTGTTCGCTGAGCGCCTGTGGCTGGCCGTGCACCTGCATCGGGGCAGTGACGATGCCGTGCGCCTGCAACGGTTGCGGGCACTGGCGGCt   <  2:2908288/150‑1 (MQ=255)
                                                                                                                                                    ctgGCCTTGTGGGTGGCCGCTGACAGTGGCGATACGCTACCGGCGCGTGGCTGCGCAGTGTGTTCGCTGAGCGCCTGTGGCTGGCCGTGCACCTGCATCGGGGCAGTGACGATGCCGTGCGCCTGCAACGGTTGCGGGCa            >  1:1243480/1‑140 (MQ=255)
                                                                                                                                                     tgGCCTTGTGGGTGGCCGCTGACAGTGGCGATACTGCTACCGGCGCGTGGCTGCGCAGTGTGTTCGCTGAGCGCCTGTGGCTGGCCGTGCACCTGCATCGGGGCAGTGACGATGCCGTGCGCCTGCAACGGTTGCGGGCACTGGCGGCTg  <  1:3923442/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GTGCGACGGGCCCAAGCTGGTGCTGCTGGTCGAAAACCTGACGGGTTACCAGAACCTGTGCGCGCTGATCACCCGTGCCCGGCGGCGGGCAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGCCTTGTGGGTGGCCGCTGACAGTGGCGATACTGCTACCGGCGCGTGGCTGCGCAGTGTGTTCGCTGAGCGCCTGTGGCTGGCCGTGCACCTGCATCGGGGCAGTGACGATGCCGTGCGCCTGCAACGGTTGCGGGCACTGGCGGCTG  >  NC_002947/3533412‑3533710

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