Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,517,056 A→T 5.8% N125I (AAT→ATT)  PP_2210 → LysR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,517,0560AT 5.8% 220.4 / 2.6 69N125I (AAT→ATT) PP_2210LysR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (31/34);  new base T (2/2);  total (33/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

CGAAGAAGCGCAGGCCATGCTTGAGGCCAACCGCCAGTTGCTGACCGGGCGCCTGGAGGTAGCCGCCGACGGCCCGCATCGGGTCATGCCGATGCTGGCGCTGCTGCGTGAGCGGTACCCTGGTGTCACCGTCAACTTGCGCCTGGGCAATGCTCAGGAGACCTTGGCTGCGCTGCTCGGCGAGCATGCCGACGTGGCGGTGCTGACCGAAATCGAGCCGCGCAAGGGGCTGCATCTGCAAAATCTCTGTGAGTCGCGCCTGTGTGCGCTGTTACCGGTGGGGCACGCCTGGGCCA  >  NC_002947/2516907‑2517202
                                                                                                                                                     |                                                                                                                                                  
cGAAGAAGCGCAGGCCATGCTTGAGGCCAACCGCCAGTTGCTGACCGGGCGCCTGGAGGTAGCCGCCGACGGCCCGCATCGGGTCATGCCGATGCTGGCGCTGCTGCGTGAGCGGTACCCTGGTGTCACCGTCAACTTGCGCCTGGGCaa                                                                                                                                                    >  2:497218/1‑150 (MQ=255)
cGAAGAAGCGCAGGCCATGCTTGAGGCCAACCGCCAGTTGCTGACCGGGCGCCTGGAGGTAGCCGCCGACGGCCCGCATCGGGTCATGCCGATGCTGGCGCTGCTGCGTGAGCGGTACCCTGGTGTCACCGTCAACTTGCGCCTGGGCaa                                                                                                                                                    <  2:1095872/150‑1 (MQ=255)
cGAAGAAGCGCAGGCCATGCTTGAGGCCAACCGCCAGTTGCTGACCGGGCGCCTGGAGGTAGCCGACGACGGCCCGCATCGGGTCATGCCGATGCTGGCGCTGCTGCGTGAGCGGTACCCTGGTGTCACCGTCAACTTGCGCCTGGGCaa                                                                                                                                                    >  1:1084882/1‑150 (MQ=255)
 gaagaaGCGCAGGCCATGCTTGAGGCCAACCGCCAGTTGCTGACCGGGCGCCTGGAGGTAGCCGCCGACGGCCCGCATCGGGTCATGCCGATGCTGGCGCTGCTGCGTGAGCGGTACCCTGGTGTCACCGTCAACTTGCGCCTGGGCAAt                                                                                                                                                   <  1:1334059/150‑1 (MQ=255)
      agtgcAGACCATGCTTGAGGCCAACCGCCAGTTGCTGACCGGGCGCCTGGAGGTAGCCGCCGACGGCCCGCATCGGGTCATGCCGATGCTGGCGCTGCTGCGTGAGCGGTACCCTGGTGTCACCGTCAACTTGCGCCTGGGCAATGCTCa                                                                                                                                              <  2:10711/147‑1 (MQ=255)
       gcgcAGGCCATGCTTGAGGCCAACCGCCAGTTGCTGACCGGGCGCCTGGAGGTAGCCGCCGACGGCCCGCATCGGGTCATGCCGATGCTGGCGCTGCTGCGTGAGCGGTACCCTGGTGTCACCGTCAACTTGCGCCTGGGCaa                                                                                                                                                    <  2:935712/143‑1 (MQ=255)
       gcgcAGGCCATGCTTGAGGCCAACCGCCAGTTGCTGACCGGGCGCCTGGAGGTAGCCGCCGACGGCCCGCATCGGGTCATGCCGATGCTGGCGCTGCTGCGTGAGCGGTACCCTGGTGTCACCGTCAACTTGCGCCTGGGCaa                                                                                                                                                    >  1:935712/1‑143 (MQ=255)
         gcAGGCCATGCTTGAGGCCAACCGCCAGTTGCTGACCGGGCGCCTGGAGGTAGCCGCCGACGGCCCGCATCGGGTCATGCCGATGCTGGCGCTGCTGCGTGAGCGGTACCCTGGTGTCACCGTCAACTTGCGCCTGAGCATTGCCCAggc                                                                                                                                           <  1:1766876/150‑2 (MQ=255)
          cAGGCCATGCTTGAGGCCAACCGCCAGTTGCTGACCGGGCGCCTGGAGGTAGCCGCCGACGGCCCGCATCGGGTCATGCCGATGCTGGCGCTGCTGCGTGAGCGGTACCCTGGTGTCACCGTCAACTTGCGCCTGAGCATTGCCCAggcg                                                                                                                                          >  2:588812/1‑148 (MQ=255)
             gCCATGCTTGAGGCCAACCGCCAGTTGCTGACCGGGCGCCTGGAGGTAGCCGCCGACGGCCCGCATCGGGTCATGCCGATGCTGGCGCTGCTGCGTGAGCGGTACCCTGGTGTCACCGTCAACTTGCGCCTGGGCAATGCTCAGGAGAcc                                                                                                                                       >  2:1684175/1‑150 (MQ=255)
             gCCATGCTTGAGGCCAACCGCCAGTTGCTGACCGGGCGCCTGGAGGTAGCCGCCGACGGCCCGCATCGGGTCATGCCGATGCTGGCGCTGCTGCGTGAGCGGTACCCTGGTGTCACCGTCAACTTGCGCCTGAGCATTGCCCAggcgcaa                                                                                                                                       <  1:588812/150‑6 (MQ=255)
                   cTTGAGGCCAACCGCCAGTTGCTGACCGGGCGCCTGGAGGTAGCCGCCGACGGCCCGCATCGGGTCATGCCGATGCTGGCGCTGCTGCGTGAGCGGTACCCTGGTGTCACCGTCAACTTGCGCCTGGGCAATGCTCAGGAGACCTTGgct                                                                                                                                 <  1:1366801/150‑1 (MQ=255)
                      gAGGCCAACCGCCAGTTGCTGACCGGGCGCCTGGAGGTAGCCGCCGACGGCCCGCATCGGGTCATGCCGATGCTGGCGCTGCTGCGTGAGCGGTACCCTGGTGTCACCGTCAACTTGCGCCTGGGCAATGCTCAGGAGACCTTGgctgcg                                                                                                                              >  2:386795/1‑150 (MQ=255)
                      gAGGCCAACCGCCAGTTGCTGACCGGGCGCCTGGAGGTAGCCGCCGACGGCCCGCATCGGGTCATGCCGATGCTGGCGCTGCTGCGTGAGCGGTACCCTGGTGTCACCGTCAACTTGCGCCTGGGCAATGCTCAGGAGACCTTGgctgcg                                                                                                                              >  2:386419/1‑150 (MQ=255)
                         gCCAACCGCCAGTTGCTGACCGGGCGCCTGGAGGTAGCCGCCGACGGCCCGCATCGGGTCATGCCGATGCTGGCGCTGCTGCGTGAGCGGTACCCTGGTGTCACCGTCAACTTGCGCCTGGGCAATGCTCAGGAGACCTTGgctgcgctg                                                                                                                           <  1:755755/150‑1 (MQ=255)
                           cAACCGCCAGTTGCTGACCGGGCGCCTGGAGGTAGCCGCCGACGGCCCGCATCGGGTCATGCCGATGCTGGCGCTGCTGCGTGAGCGGTACCCTGGTGTCACCGTCAACTTGCGCCTGGGCAATGCTCAg                                                                                                                                             >  1:886479/1‑130 (MQ=255)
                           cAACCGCCAGTTGCTGACCGGGCGCCTGGAGGTAGCCGCCGACGGCCCGCATCGGGTCATGCCGATGCTGGCGCTGCTGCGTGAGCGGTACCCTGGTGTCACCGTCAACTTGCGCCTGGGCAATGCTCAg                                                                                                                                             <  2:886479/130‑1 (MQ=255)
                                 ccAGTTGCTGACCGGGCGCCTGGAGGTAGCCGCCGACGGCCCGCATCGGGTCATGCCGATGCTGGCGCTGCTGCGTGAGCGGTACCCTGGTGTCACCGTCAACTTGCGCCTGGGCAATGCTCAGGAGACCTTGGCTGCGCTGCTCGGCGa                                                                                                                   >  1:767667/1‑150 (MQ=255)
                                   aGTTGCTGACCGGGCGCCTGGAGGTAGCCGCCGACGGCCCGCATCGGGTCATGCCGATGCTGGCGCTGCTGCGTGAGCGGTACCCTGGTGTCACCGTCAACTTGCGCCTGGGCAATGCTCAGGAGCCCTTGGCTGCGCTGCTCGGCGAGc                                                                                                                 >  2:1900049/1‑150 (MQ=255)
                                       gCTGACCGGGCGCCTGGAGGTAGCCGCCGACGGCCCGCATCGGGTCATGCCGATGCTGGCGCTGCTGCGTGAGCGGTACCCTGGTGTCACCGTCAGCTTGCGCCTGGGCAATGCTCAGGAGACCTTGGCTGCGCTGCTCGGCGAGCATGc                                                                                                             <  2:1225816/150‑1 (MQ=255)
                                        cTGACCGGGCGCCTGGAGGTAGCCGCCGACGGCCCGCATCGGGTCATGCCGATGCTGGCGCTGCTGCGTGAGCGGTACCCTGGTGTCACCGTCAACTTGCGCCTGGGCAATGCTCAGGAGACCTTGGCTGCGCTGCTCGGCGAGCATGcc                                                                                                            >  2:1913392/1‑150 (MQ=255)
                                         tGACCGGGCGCCTGGAGGTAGCCGCCGACGGCCCGCATCGGGTCATGCCGATGCTGGCGCTGCTGCGTGAGCGGTACCCTGGTGTCACCGTCAACTTGCGCCTGGGCAATGCTCAGGAGACCTTGGCTGCGCTGCTCGGCGAGCATGCCg                                                                                                           >  1:1707774/1‑150 (MQ=255)
                                                 cgcCTGGAGGTAGCCGCCGACGGCCCGCATCGGGTCATGCCGATGCTGGCGCTGCTGCGTGAGCGGTACCCTGGTGTCACCGTCAACTTGCGCCTGGGCAATGCTCAGGAGACCTTGGCTGCGCTGCTCGGCGAGCATGCCGACGTGGCg                                                                                                   <  2:658688/150‑1 (MQ=255)
                                                  gcCTGGAGGTAGCCGCCGACGGCCCGCATCGGGTCATGCCGATGCTGGCGCTGCTGCGTGAGCGGTACCCTGGTGTCACCGTCAACTTGCGCCTGGGCAATGCTCAGGAGACCTTGGCTGCGCTGCTCGGCGAGCATGCCGACGTGGCgg                                                                                                  <  2:767667/150‑1 (MQ=255)
                                                        aGGTAGCCGCCGACGGCCCGCATCGGGTCATGCCGATGCTGGCGCTGCTGCGTGAGCGGTACCCTGGTGTCACCGTCAACTTGCGCCTGGGCAATGCTCAGGAGACCTTGGCTGCGCTGCTCGGCGAGCATGCCGACGTGGCGGTGCTGa                                                                                            <  1:321007/150‑1 (MQ=255)
                                                        aGGTAGCCGCCGACGGCCCGCATCGGGTCATGCCGATGCTGGCGCTGCTGCGTGAGCGGTACCCTGGTGTCACCGTCAACTTGCGCCTGGGCAATGCTCAGGAGACCTTGGCTGCGCTGCTCGGCGAGCATGCCGACGTGGCGGTGCTGa                                                                                            <  1:386419/150‑1 (MQ=255)
                                                        aGGTAGCCGCCGACGGCCCGCATCGGGTCATGCCGATGCTGGCGCTGCTGCGTGAGCGGTACCCTGGTGTCACCGTCAACTTGCGCCTGGGCAATGCTCAGGAGACCTTGGCTGCGCTGCTCGGCGAGCATGCCGACGTGGCGGTGCTGa                                                                                            <  1:386795/150‑1 (MQ=255)
                                                              ccgccgACGGCCCGCATCGGGTCATGCCGATGCTGGCGCTGCTGCGTGAGCGGTACCCTGGTGTCACCGTCAACTTGCGCCTGGGCAATGCTCAGGAGACCTTGGCTGCGCTGCTCGGCGAGCATGCCGACGTGGCGGTGCTGACCGaaa                                                                                      >  2:1016259/1‑150 (MQ=255)
                                                              ccgccgACGGCCCGCATCGGGTCATGCCGATGCTGGCGCTGCTGCGTGAGCGGTACCCTGGTGTCACCGTCAACTTGCGCCTGGGCAATGCTCAGGAGACCTTGGCTGCGCTGCTCGGCGAGCATGCCGACGTGGCGGTGCTGACCGaaa                                                                                      >  1:1116701/1‑150 (MQ=255)
                                                                gccgACGGCCCGCATCGGGTCATGCCGATGCTGGCGCTGCTGCGTGAGCGGTACCCTGGTGTCACCGTCAACTTGCGCCTGGGCAATGCTCAGGAGACCTTGGCTGCGCTGCTCGGCGAGCATGCCGACGTGGCGGTGCTGACCGAAATc                                                                                    >  1:1942156/1‑150 (MQ=255)
                                                                  cgACGGCCCGCATCGGGTCATGCCGATGCTGGCGCTGCTGCGTGAGCGGTACCCTGGTGTCACCGTCAACTTGCGCCTGGGCAATGCTCAGGAGACCTTGGCTGCGCTGCTCGGCGAGCATGCCGACGTGGCGGTGCTGACCGAAATCGa                                                                                  <  2:1857904/150‑1 (MQ=255)
                                                                   gACGGCCCGCATCGGGTCATGCCGATGCTGGTGCTGCTGCGTGAGCGGTACCCTGGTGTCACCGTCAACTTGCGCCTGGGCAATGCTCAGGAGACCTTGGCTGCGCTGCTCGGCGAGc                                                                                                                 <  1:1796187/118‑1 (MQ=255)
                                                                   gACGGCCCGCATCGGGTCATGCCGATGCTGGTGCTGCTGCGTGAGCGGTACCCTGGTGTCACCGTCAACTTGCGCCTGGGCAATGCTCAGGAGACCTTGGCTGCGCTGCTCGGCGAGc                                                                                                                 >  2:1796187/1‑118 (MQ=255)
                                                                    aCGGCCCGCATCGGGTCATGCCGATGCTGGCGCTGCTGCGTGAGCGGTACCCTGGTGTCACCGTCAACTTGCGCCTGGGCAATGCTCAGGAGACCTTGGCTGCGCTGCTCGGCGAGCATGCCGACGTGGCGGTGCTGACCGAAATCGAGc                                                                                <  2:1707774/150‑1 (MQ=255)
                                                                        cccGCATCGGGTCATGCCGATGCTGGCGCTGCTGCGTGAGCGGTACCCTGGTGTCACCGTCAACTTGCGCCTGGGCAATGCTCAGGAGACCTTGGCTGCGCTGCTCGGCGAGCATGCCGACGTGGCGGTGCTGACCGAAATCGAGCcgcg                                                                            <  1:1900049/150‑1 (MQ=255)
                                                                         ccGCATCGGGTCATGCCGATGCTGGCGCTGCTGCGTGAGCGGTACCCTGGTGTCACCGTCAACTTGCGCCTGGGCAATGCTCAGGAGACCTTGGCTGCGCTGCTCGGCGAGCATGCCGACGTGGCGGTGCTGACCGAAATCGAGCCgcgc                                                                           >  2:1757179/1‑150 (MQ=255)
                                                                         ccGAATCGGGTCATGCCGATGCTGGCGCTGCTGCGTGAGCGGTACCCTGGTGTCACCGTCAACTTGCGCCTGGGCAATGCTCAGGAGACCTTGGCTGCGCTGCTCGGCGAGAATGCCGACGTGGCGGTGCTGACCGAAATCGAGCCgcgc                                                                           >  1:1136182/1‑150 (MQ=255)
                                                                          cGCATCGGGTCATGCCGATGCTGGCGCTGCTGCGTGAGCGGTACCCTGGTGTCACCGTCAACTTGCGCCTGGGCAATGCTCAGAAAACCTTGGCTGCGCTGCTCGGCGAGCATGCCGACGTGGCGGTGCTGACCGAAATCGAGCCGCGCa                                                                          >  2:551220/1‑150 (MQ=255)
                                                                              tCGGGTCATGCCGATGCTGGCGCTGCTGCGTGAGCGGTACCCTGGTGTCACCGTCAACTTGCGCCTGGGCAATGCTCAGGAGACCTTGGCTGCGCTGCTCGGCGAGCATGCCGACGTGGCGGTGCTGACCGAAATCGAGCCGCGCAAggg                                                                      <  2:911166/150‑1 (MQ=255)
                                                                               cGGGTAAAGCAGATGCTGGTGCTGCTGCGTGAGCGGTACCCTGGTGTCACCGTCAACTTGCGCCTCGGCAATGCTCAGGAGACCTTGGCTGCGCTGCTCGGCGAGCATGCCGACGTGGCGGTGCTGACCGAAATCGAGCCGCGCAAgggg                                                                     >  2:1392285/1‑150 (MQ=255)
                                                                                 ggTCATGCCGATGCTGGCGCTGCTGCGTGAGCGGTACCCTGGTGTAACCGTCACCTTGCGCTTTGGCAATGCTCAGGAGACCTTGGCTGCGCTGCTCGGCGAGCATGCCGACGTGGCGGTGCTGACCGAAATCGAGCCGCGCAAGGGGCt                                                                   <  1:1684175/150‑1 (MQ=255)
                                                                                   tCATGCCGATGCTGGCGCTGCTGCGTGAGCGGTACCCTGGTGTCACCGTCAACTTGCGCCTGGGCAATGCTCAGGAGACCTTGGCTGCGCTGCTCGGCGAGCATGCCGACGTGGCGGTGCTGACCGAAATCGAGcc                                                                               >  2:326184/1‑136 (MQ=255)
                                                                                   tCATGCCGATGCTGGCGCTGCTGCGTGAGCGGTACCCTGGTGTCACCGTCAACTTGCGCCTGGGCAATGCTCAGGAGACCTTGGCTGCGCTGCTCGGCGAGCATGCCGACGTGGCGGTGCTGACCGAAATCGAGcc                                                                               <  1:326184/136‑1 (MQ=255)
                                                                                   tCATGCCGATGCTGGCGCTGCTGCGTGAGCGGTACCCTGGTGTCACCGTCAACTTGCGCCTGGGCAATGCTCAGGAGACCTTGGCTGCGCTGCTCGGCGAGCATGCCGACGTGGCGGTGCTGACCGAAATCGAGCCGCGCAAGGGGCTGc                                                                 >  1:1647808/1‑150 (MQ=255)
                                                                                              cTGGCGCTGCTGCGTGAGCGGTACCCTGGTGTCACCGTCAACTTGCGCCTGGGCAATGCTCAGGAGACCTTGGCTGCGCTGCTCGGCGAGCATGCCGACGTGGCGGTGCTGACCGAAATCGAGCCGCGCAAGGGGCTGCATCTGCAAAAt                                                      >  2:350270/1‑150 (MQ=255)
                                                                                                  cgcTGCTGCGTGAGCGGTACCCTGGTGTCACCGTCAACTTGCGCCTGGGCAATGCTCAGGAGACCTTGGCTGCGCTGCTCGGCGAGCATGCCGACGTGGCGGTGCTGACCGAAATCGAGCCGCGCAAGGGGCTGCATCTGCAAAATctct                                                  <  1:497218/150‑1 (MQ=255)
                                                                                                    ctgctgCGTGAGCGGTACCCTGGTGTCACCGTCAACTTGCGCCTGGGCAATGCTCAGGAGACCTTGGCTGCGCTGCTCGGCGAGCATGCCGACGTGGCGGTGCTGACCGAAATCGAGCCGCGCAAGGGGCTGCATCTGCAAAATCTCtgt                                                <  2:1647808/150‑1 (MQ=255)
                                                                                                    ctgctgCGTGAGCGGTACCCTGGTGTCACCGTCAACTTGCGCCTGGGCAATGCTCAGGAGACCTTGGCTGCGCTGCTCGGCGAGCATGCCGACGTGGCGGTGCTGACCGAAATCGAGCCGCGCAAGGGGCTGCATCTGCAAAATCTCtgt                                                <  1:350270/150‑1 (MQ=255)
                                                                                                    ctgctgCGTGAGCGGTACCCTGGTGTCACCGTCAACTTGCGCCTGGGCAATGCTCAGGAGACCTTGGCTGCGCTGCTCGGCGAGCATGCCGACGTGGCGGTGCTGACCGAAATCGAGCCGCGCAAGGGGCTGCATCTGCAAAATCTCtgt                                                <  1:1786273/150‑1 (MQ=255)
                                                                                                         gcGTGAGCGGTACCCTGGTGTCACCGTCAACTTGCGCCTGGGCAATGCTCAGGAGACCTTGGCTGCGCTGCTCGGCGAGCATGCCGACGTGGCGGTGCTGACCGAAATCGAGCCGCGCAAGGGGCTGCATCTGCAAAATCTCTGTGAGTc                                           <  2:1116701/150‑1 (MQ=255)
                                                                                                               gCGGTACCCTGGTGTCACCGTCAACTTGCGCCTGGGCAATGCTCCGGAGACCTTGGCTGCGCTGCTCGGCGAGCATGCCGACGTGGCGGTGCTGACCGAAATCGAGCCGCGCAAGGGGCTGCATCTGCAAAATCTCTGTGAGTCGCGCCt                                     <  2:1399421/150‑1 (MQ=255)
                                                                                                                   tACCCTGGTGTCACCGTCAACTTGCGCCTGGGCAATGCTCAGGAGACCTTGGCTGCGCTGCTCGGCGAGCATGCCGACGTGGCGGTGCTGACCGAAATCGAGCCGCGCAAGGGGCTGCATCTGCAAAATCTCTGTGAGTCGCGCCtgtgt                                 >  2:894254/1‑150 (MQ=255)
                                                                                                                          gtgtCACCGTCAACTTGCGCCTGGGCAATGCTCAGGACACCTTGGCTGCGCTGCTCGGCGAGCATGCCGACGTGGCGGTGCTGACCGAAATCGAGCCGCGCAAGGGGCTGCATCTGCAAAATCTCTGTGAGTCGCGCCTGTGTGCGCTGt                          <  2:940869/150‑1 (MQ=255)
                                                                                                                             tCACCGTCAACTTGCGCCTGGGCAATGCTCAGGAGACCTTGGCTGCGCTGCTCGGCGAGCATGCCGACGTGGCGGTGCTGACCGAAATCGAGCCGCGCAAGGGGCTGCATCTGCAAAATCTCTGTGAGTCGCGCCTGTGTGCGCTGTTAc                       >  2:936870/1‑150 (MQ=255)
                                                                                                                             tCACCGTCAACTTGCGCCTGGGCAATGCTCAGGAGACCTTGGCTGCGCTGCTCGGCGAGCATGCCGACGTGGCGGTGCTGACCGAAATCGAGCCGCGCAAGGGGCTGCATCTGCAAAATCTCTGTGAGTCGCGCCTGTGTGCGCTGTTAc                       <  1:894254/150‑1 (MQ=255)
                                                                                                                               aCCGTCAACTTGCGCCTGGGCAATGCTCAGGAGACCTTGGCTGCGCTGCTCGGCGAGCATGCCGACGTGGCGGTGCTGACCGAAATCGAGCCGCGCAAGGGGCTGCATc                                                              >  1:777813/1‑109 (MQ=255)
                                                                                                                               aCCGTCAACTTGCGCCTGGGCAATGCTCAGGAGACCTTGGCTGCGCTGCTCGGCGAGCATGCCGACGTGGCGGTGCTGACCGAAATCGAGCCGCGCAAGGGGCTGCATc                                                              <  2:777813/109‑1 (MQ=255)
                                                                                                                                  gTCAACTTGCGCCTGGGCAATGCTCAGGAGACCTTGGCTGCGCTGCTCGGCGAGCATGCCGACGTGGCGGTGCTGACCGAAATCGAGCCGCGCAAGGGGCTGCATCTGCAAAATCTCTGTGAGTCGCGCCtgtgt                                 <  1:1774037/135‑1 (MQ=255)
                                                                                                                                  gTCAACTTGCGCCTGGGCAATGCTCAGGAGACCTTGGCTGCGCTGCTCGGCGAGCATGCCGACGTGGCGGTGCTGACCGAAATCGAGCCGCGCAAGGGGCTGCATCTGCAAAATCTCTGTGAGTCGCGCCtgtgt                                 >  2:1774037/1‑135 (MQ=255)
                                                                                                                                  gTCAACTTGCGCCTGGGCAATGCTCAGGAGACCTTGGCTGCGCTGCTCGGCGAGCATGCCGACGTGGCGGTGCTGACCGAAATCGAGCCGCGCAAGGGGCTGCATCTGCAAAATCTCTGTGAGTCGCGCCTGTGTGCGCTGTTACCGGTg                  <  2:1942156/150‑1 (MQ=255)
                                                                                                                                   tCAACTTGCGCCTGGGCAATGCTCAGGAGACCTTGGCTGCGCTGCTCGGCGAGCATGCCGAc                                                                                                         <  1:743460/62‑1 (MQ=255)
                                                                                                                                   tCAACTTGCGCCTGGGCAATGCTCAGGAGACCTTGGCTGCGCTGCTCGGCGAGCATGCCGAc                                                                                                         >  2:743460/1‑62 (MQ=255)
                                                                                                                                     aaCTTGCGCCTGGGCAATGCTCAGGAGACCTTGGCTGCGCTGCTCGGCGAGCATGCCGACGTGGCGGTGCTGACCGAAATCGAGCCGCGCAAGGGGCTGCATCTGCAAAATCTCtgt                                                <  1:1302890/117‑1 (MQ=255)
                                                                                                                                     aaCTTGCGCCTGGGCAATGCTCAGGAGACCTTGGCTGCGCTGCTCGGCGAGCATGCCGACGTGGCGGTGCTGACCGAAATCGAGCCGCGCAAGGGGCTGCATCTGCAAAATCTCtgt                                                >  2:1302890/1‑117 (MQ=255)
                                                                                                                                      aCTTGCGCCTGGGCAATGCTCAGGAGACCTTGGCTGCGCTGCTCGGCGAGCATGCCGACGTGGCGGTGCTGACCGAAATCGAGCCGCGCAAGGGGCTGCATCTGCAAAATCTCTGTGAGTCGCGCCTGTGTGCGCTGTTACCGGTGgcgc              <  1:1127698/150‑1 (MQ=255)
                                                                                                                                        ttGCGCCTGGGCAATGCTCAGGAGACCTTGGCTGCGCTGCTCGGCGAGCATGCCGACGTGGCGGTGCTGACCGAAATCGAGCCGCGCAAGGGGCTGCATCTGCAAAATCTCTGTGAGTCGCGCCTGTGTGCGCTGTTACCGGTGGGGCAc            >  1:1203390/1‑150 (MQ=255)
                                                                                                                                        ggtctccTGAGCATTGCCCAGGAGACCTTGGCTGCGCTGCTCGGCGAGCATGCCGACGTGGCGGTGCTGACCGAAATCGAGCCGCGCAAGGGGCTGCATCTGCAAAATCTCTGTGAGTCGCGCCTGTGTGCGCTGTTACCGGTGGGGCAc            >  2:1035117/6‑150 (MQ=255)
                                                                                                                                              cTGGGCAATGCTCAGGAGACCTTGGCTGCGCTGCTCGGCGAGCATGCCGACGTGGCGGTGCTGACCGAAATCGAGCCGCGCAAGGGGCTGCATCTGCAAAATCTCTGTGAGTCGCGCCTGTGTGCGCTGTTACCGGTGGGGCACGCCTgg      <  1:395158/150‑1 (MQ=255)
                                                                                                                                              cTGGGCAATGCTCAGGAGACCTTGGCTGCGCTGCTCGGCGAGCATGCCGACGTGGCGGTGCTGACCGAAATCGAGCCGCGCAAGGGGCTGCATCTGCAAAATCTCTGTGAGTCGCGCCTGTGTGCGCTGTTACCGGTGGGGCACGCCTgg      >  1:546807/1‑150 (MQ=255)
                                                                                                                                              cTGGGCAATGCTCAGGAGACCTTGGCTGCGCTGCTCGGCGAGCATGCCGACGTGGCGGTGCTGACCGAAATCGAGCCGCGCAAGGGGCTGCATCTGCAAAATCTCTGTGAGTCGCGCCTGTGTGCGCTGTTACCGGTGGGGCACGCCTgg      >  2:1617763/1‑150 (MQ=255)
                                                                                                                                               tGGGCAATGCTCAGGAGACCTTGGCTGCGCTGCTCGGCGAGCATGCCGACGTGGCGGTGCTGACCGAAATCGAGCCGCGCAAGGGGCTGCATCTGCAAAATCTCTGTGAGTCGCGCCTGTGTGCGCTGTTACCGGTGGGGCACGCCTggg     >  2:778017/1‑150 (MQ=255)
                                                                                                                                                  gCAATGCTCAGGAGACCTTGGCTGCGCTGCTCGGCGAGCATGCCGACGTGGCGGTGCTGACCGAAATCGAGCCGCGCAAGGGGCTGCATCTGCAAAATCTCTGTGAGTCGCGCCTGTGTGCGCTGTTACCGGTGGGGCACGCCTGGGTCa  >  1:1520468/1‑150 (MQ=255)
                                                                                                                                                   cAATGCTCAGGAGACCTTGGCTGCGCTGCTCGGCGAGCATGCCGACGTGGCGGTGCTGACCGAAATCGAGCCGCGCAAGGGGCTGCATCTGCAAAATCTCTGTGAGTCgcgc                                       >  1:1937064/1‑112 (MQ=255)
                                                                                                                                                   cAATGCTCAGGAGACCTTGGCTGCGCTGCTCGGCGAGCATGCCGACGTGGCGGTGCTGACCGAAATCGAGCCGCGCAAGGGGCTGCATCTGCAAAATCTCTGTGAGTCgcgc                                       <  2:1937064/112‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
CGAAGAAGCGCAGGCCATGCTTGAGGCCAACCGCCAGTTGCTGACCGGGCGCCTGGAGGTAGCCGCCGACGGCCCGCATCGGGTCATGCCGATGCTGGCGCTGCTGCGTGAGCGGTACCCTGGTGTCACCGTCAACTTGCGCCTGGGCAATGCTCAGGAGACCTTGGCTGCGCTGCTCGGCGAGCATGCCGACGTGGCGGTGCTGACCGAAATCGAGCCGCGCAAGGGGCTGCATCTGCAAAATCTCTGTGAGTCGCGCCTGTGTGCGCTGTTACCGGTGGGGCACGCCTGGGCCA  >  NC_002947/2516907‑2517202

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