Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 403,299 T→A 5.3% intergenic (‑99/+70) msrA ← / ← PP_0337 peptide methionine sulfoxide reductase MsrA/sensory box protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947403,2990TA 5.3% 493.3 / 13.5 150intergenic (‑99/+70)msrA/PP_0337peptide methionine sulfoxide reductase MsrA/sensory box protein
Reads supporting (aligned to +/- strand):  ref base T (86/56);  new base A (3/5);  total (89/61)
Fisher's exact test for biased strand distribution p-value = 2.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.09e-01

Only 100 of 154 total aligned reads displayed.

CGCAGTGGGCATGACGTTTTTGTTCACCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCAGGCCACCGGGCCACGCGGATAGCGCTTCAGCCGCTCGGCCAGTTCGCGCCCGGGAATCGGCCGGTCGAACAGGTAGC  >  NC_002947/403150‑403448
                                                                                                                                                       |                                                                                                                                                     
cGCAGTGGGCATGACGTTTTTGTTCACCAGGATTTCCGAACGCAGGACCATGTCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGtt                                                                                                                                                       <  2:4355207/150‑1 (MQ=255)
cGCAGTGGGCATGACGTTTTTGTTCACCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGtt                                                                                                                                                       >  2:4396676/1‑150 (MQ=255)
 gCAGTGGGCATGACGTTTTTGTTCACCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTc                                                                                                                                                      <  1:2733150/150‑1 (MQ=255)
    gTGGGCATGACGTTTTTGTTCACCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCtg                                                                                                                                                   <  1:906243/150‑1 (MQ=255)
     tGGGCATGACGTTTTTGTTCACCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCtgt                                                                                                                                                  <  2:2260368/150‑1 (MQ=255)
     tGGGCATGACGTTTTTGTTCACCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCtgt                                                                                                                                                  <  2:2645965/150‑1 (MQ=255)
           tGACGTTTTTGTTCACCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCCCTTTCGCGGGCACGCCCGCTCCCACAGGAACACCACAACCTTTGAAACCTGTGCTGTACCTGTGGGa                                                                                                                                              <  1:2470481/150‑1 (MQ=17)
            gACGTTTTTGTTCACCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCCCTTTCGCGGGCACGCCCGCTCCCACAGGAACACCACAACCTTTGAAACCTGTGCTGTACCTGTGGGAg                                                                                                                                             <  2:313878/150‑1 (MQ=17)
            gACGTTTTTGTTCACCAGGATTTCCGAACGAAGGACCATGGCAGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCACAGGTACAGCGCACGTTTCAAAGGTTGTGGTGTTACTGTGGGTGCg                                                                                                                                           >  1:3450973/1‑150 (MQ=39)
             aCGTTTTTGTTCACCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCgg                                                                                                                                          <  2:1255609/150‑1 (MQ=255)
             aCGTTTTTGTTCACCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCgg                                                                                                                                          <  2:1255647/150‑1 (MQ=255)
              cttttttGTACACCTGTATTTCGGAACGCAGTACCTTGGCCGTTCCTATGTGTTTATGCAGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCggg                                                                                                                                         <  1:347522/148‑1 (MQ=37)
               gTTTTTGTTCACCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGc                                                                                                                                        <  2:3672769/150‑1 (MQ=255)
               gTTTTTGTTCACCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGc                                                                                                                                        >  2:3573576/1‑150 (MQ=255)
                   ttGTTCACCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGc                                                                                                                                    <  1:2299515/150‑1 (MQ=255)
                   ttGTTCACCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGc                                                                                                                                    <  1:1938178/150‑1 (MQ=255)
                     gTTCACCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCCTGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGccc                                                                                                                                  >  1:2042111/1‑150 (MQ=255)
                      ttCACCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTGTTCAGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTGTCAAAGGTTGTGGTGTTCCTGTGGGCGCGGGCGTGCCcg                                                                                                                                 >  2:3740494/1‑150 (MQ=255)
                         aCCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGa                                                                                                                              >  2:180876/1‑150 (MQ=255)
                         aCCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCGGTGGGAGCGGGCGTGCCCGCGa                                                                                                                              >  2:2469071/1‑150 (MQ=255)
                           cAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGaaa                                                                                                                            >  2:2757016/1‑150 (MQ=255)
                            aGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCATTCCCGCTCCCACAGGTACAGCACCGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAg                                                                                                                           >  1:1999375/1‑150 (MQ=255)
                            aGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAg                                                                                                                           >  1:46593/1‑150 (MQ=255)
                            aGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAg                                                                                                                           >  1:632825/1‑150 (MQ=255)
                                     gAACGCAGGCCCATGGCCGGTCCTCTGTGTTTATGCAGTGGGGGAGGGGTTGGGGCC‑‑GTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGGTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAAt                                                                                                                  <  2:1421942/150‑1 (MQ=37)
                                      aaCGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATg                                                                                                                 >  2:1072617/1‑150 (MQ=255)
                                      aaCGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATg                                                                                                                 >  1:3260434/1‑150 (MQ=255)
                                       aCGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATgg                                                                                                                <  2:4277359/150‑1 (MQ=255)
                                         gCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCt                                                                                                              <  2:700216/150‑1 (MQ=255)
                                           aGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTc                                                                                                            >  2:3014782/1‑150 (MQ=255)
                                           aGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGGTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTc                                                                                                            <  2:3867152/150‑1 (MQ=255)
                                             gACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCATTCCCGCTCCCAAAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCCCGAAAGGGCCCCAATGGCttaaa                                                                                                          >  1:3204815/1‑147 (MQ=39)
                                             gACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCta                                                                                                          <  2:3702320/150‑1 (MQ=255)
                                             gACCATGGCCGTTCCTCTGTGTTTATTCAGGGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCta                                                                                                          >  2:121283/1‑150 (MQ=255)
                                               ccATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTGGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCtata                                                                                                        <  2:1360249/150‑1 (MQ=255)
                                                cATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAg                                                                                                       >  2:1704407/1‑150 (MQ=255)
                                                          ccTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGa                                                                                             <  2:3012036/150‑1 (MQ=255)
                                                            tctGTGTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGAtt                                                                                           >  2:2458197/1‑150 (MQ=255)
                                                            tctGTGTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGAtt                                                                                           >  2:2591727/1‑150 (MQ=255)
                                                             ctGTGTTTATTCTGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTa                                                                                          <  2:1231213/150‑1 (MQ=255)
                                                              tgtgTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTAc                                                                                         >  1:3375111/1‑150 (MQ=255)
                                                              tgtgTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTAc                                                                                         >  2:1566091/1‑150 (MQ=255)
                                                                   ttATTCAGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGGGTTCCTGGGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCGGCTATAGTGCCCGAGGATTACGCAGc                                                                                    >  2:2627891/1‑150 (MQ=255)
                                                                         cggggggTCGGGGTTGGGGCC‑‑TTCGGGGGCATGCCCGCTCCCCCAGGTACAGCACCGGTTCCAATGGTTGTGGTGTTCGTGTGGGCGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCAg                                                                              <  1:3941727/147‑1 (MQ=255)
                                                                         aGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCAg                                                                              >  1:2906667/1‑150 (MQ=255)
                                                                         aGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCAg                                                                              >  1:4054373/1‑150 (MQ=255)
                                                                         aGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCAg                                                                              >  2:323403/1‑150 (MQ=255)
                                                                         aGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCAAGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGATCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTGCTATAGAGCCCGAGGATTACGCAGCTGTCAg                                                                              >  2:2046365/1‑150 (MQ=255)
                                                                          gTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGGTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCGATAGTGCCCGAGGATTACGCAGCTGTCAgg                                                                             >  1:2191844/1‑150 (MQ=255)
                                                                           tGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCACAGGGACCGCCCAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGCACGGGCCCGAATGGCGTCTATAGTGCCCGTGGATTACGCAGCGGTCAGGc                                                                            >  1:4420213/1‑150 (MQ=38)
                                                                               gTAGGGGTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCACAGGGACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATATTGCCCGAGGATTACGCAGCTGTCAGGCCAcc                                                                        >  1:998369/1‑150 (MQ=255)
                                                                               gTAGGGGTTGGGGCC‑‑TTCGCGGGCAGGCCCGCGCCCCCCGGGACAGCCCAGGGTGCAACGGTTGTGGGGGTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATCGGGCCCGAGGATTTCGCCGCTGTCAGGCCAcc                                                                        <  2:2869957/150‑1 (MQ=255)
                                                                                   gggTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCAGGCCACCGGGc                                                                    >  1:1026398/1‑150 (MQ=255)
                                                                                    ggTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCAGGCCACCGGGcc                                                                   <  2:1874691/150‑1 (MQ=255)
                                                                                     gTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCAGGCCACCGGGCCa                                                                  <  2:682348/150‑1 (MQ=255)
                                                                                            cc‑‑TTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTACTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCACGAGGATTACGCAGCTGTCAGGCCACCGGGCCACGCGGAt                                                           >  1:2575366/1‑150 (MQ=255)
                                                                                               ctcggggggCTTGCCCGCCCCCCCAGGTACAGCACAGGTTACACCGGTTGTGGTGTTCCTGTGGGAGGGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCAGGCCACCGGGCCACGCGGATa                                                          <  1:3803702/144‑1 (MQ=255)
                                                                                                ggagcgGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCGGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCAGGCCACCGGGCCACGCGGATAg                                                         >  2:3686181/4‑150 (MQ=255)
                                                                                                    cgGGCAGGCCCGCTCCCACAGGTACCGCACCGGTTCCAAAGGGTGGGGTGTTCCTGTGGGCGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGGGCCCGAGGATTACGCAGCTGTCAGGCCACCGGGCCACGCGGATAGCGCt                                                     <  2:101311/150‑1 (MQ=255)
                                                                                                        ctggCCCGCGCCCACTGTTACAGCACTGGTTTCACCGGTTGTGGTGTTCCTGTGGGTGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCAGGCCACCGGGCCACGCGGATAGCGCTTCAg                                                 <  1:1217942/147‑1 (MQ=255)
                                                                                                         aTGCCCGCTCCCCCAGGTACAGCACAGGTTCCAACGGTTGTGGGGTTCCTGTGGGAGCGGGGGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCAGGCCACCGGGCCACGCGGATAGCGCTTCAGc                                                <  2:910465/150‑1 (MQ=255)
                                                                                                         aTGCCCGCTCCCACAGGTACAGCACAGGTTTCACAGGGTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCGTCTATAGTGCCCGAGGATTACGCAGCTGTCAGGCCACCGGGCCACGCGGATAGCGCTTCAGc                                                >  2:1698191/1‑150 (MQ=255)
                                                                                                         aTGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCAGGCCACCGGGCCACGCGGATAGCGCTTCAGc                                                >  2:2660838/1‑150 (MQ=255)
                                                                                                           gCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCAGGCCACCGGGCCACGCGGATAGCGCTTCAGCCg                                              >  1:1008504/1‑150 (MQ=255)
                                                                                                             ccGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCAGGCCACCGGGCCACGCGGATAGCGCTTCAGCCGCt                                            >  1:925988/1‑150 (MQ=255)
                                                                                                             ccGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCAGGCCACCGGGCCACGCGGATAGCGCTTCAGCCGCt                                            >  1:2141951/1‑150 (MQ=255)
                                                                                                             ccGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCAGGCCACCGGGCCACGCGGATAGCGCTTCAGCCGCt                                            >  1:867643/1‑150 (MQ=255)
                                                                                                              cGCTCCCACAGGTACAGCACAGGTTTCAACGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCAGGCCACCGGGCCACGCGGATAGCGCTTCAGCCGCTc                                           <  1:2732257/150‑1 (MQ=255)
                                                                                                               gCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCAGGCCACCGGGCCACGCGGATAGCGCTTCAGCCGCTCg                                          >  2:1658374/1‑150 (MQ=255)
                                                                                                               gCTCCCACAGGAACACCACAACCTTTGAAACCTGTGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCAGGCCACCGGGCCACGCGGATAGCGCTTCAGCCGCTCg                                          >  2:3213966/1‑150 (MQ=21)
                                                                                                                cTCCCACAGGAACACCACAACCTTTGAAACCTGTGCTGTACCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCAGGCCACCGGGCCACGCGGATAGCGCTTCAGCCGCTCgg                                         >  1:3388624/1‑150 (MQ=18)
                                                                                                                   ccACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCAGGCCACCGGGCCACGCGGATAGCGCTTCAGCCGCTCGGCCa                                      >  1:3096991/1‑150 (MQ=255)
                                                                                                                     acaGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCAGGCCACCGGGCCACGCGGATAGCGCTTCAGCCGCTCGGCCAGt                                    >  2:2484623/1‑150 (MQ=255)
                                                                                                                     acaGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCAGGCCACCGGGCCACGCGGATAGCGCTTCAGCCGCTCGGCCAGt                                    >  2:3144529/1‑150 (MQ=255)
                                                                                                                        ggTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCAGGCCACCGGGCCACGCGGATAGCGCTTCAGCCGCTCGGCCAGTTcg                                 >  1:1343480/1‑150 (MQ=255)
                                                                                                                          tACAGCCCTGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCAGGCCACCGGGCCACGCGGATAGCGCTTCAGCCGCTCGGCCAGTTcgcg                               <  1:2307657/150‑1 (MQ=255)
                                                                                                                          tACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCAGGCCACCGGGCCACGCGGATAGCGCTTCAGCCGCTCGGCCAGTTcgcg                               >  2:3536498/1‑150 (MQ=255)
                                                                                                                              gCACAGGTTTCAATGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCAGGCCACCGGGCCACGCGGATAGCGCTTCAGCCGCTCGGCCAGTTCGCGCCCg                           <  1:3107506/150‑1 (MQ=255)
                                                                                                                              gCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCAGGCCACCGGGCCACGCGGATAGCGCTTCAGCCGCTCGGCCAGTTCGCGCCCg                           <  2:2721780/150‑1 (MQ=255)
                                                                                                                                acaGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCAGGCCACCGGGCCACGCGGATAGCGCTTCAGCCGCTCGGCCAGTTCGCGCCCggg                         >  2:3787123/1‑150 (MQ=255)
                                                                                                                                acaGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCAGGCCACCGGGCCACGCGGATAGCGCTTCAGCCGCTCGGCCAGTTCGCGCCCggg                         >  1:1673241/1‑150 (MQ=255)
                                                                                                                                  cggTTTCAACGGGTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCAGGCCACCGGGCCACGCGGATAGCGCTTCAGCCGCTCGGCCAGTTCGCGCCCGGGaa                       <  1:2227649/149‑1 (MQ=255)
                                                                                                                                    gTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCAGGCCACCGGGCCACGCGGATAGCGCTTCAGCCGCTCGGCCAGTTc                                  >  1:2051527/1‑137 (MQ=255)
                                                                                                                                    gTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCAGGCCACCGGGCCACGCGGATAGCGCTTCAGCCGCTCGGCCAGATCGCGCCCGGGAATc                     >  2:2292303/1‑150 (MQ=255)
                                                                                                                                      ttCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCAGGCCACCGGGCCACGCGGATAGCGCTTCAGCCGCTCGGCCAGTTCGCGCCCGGGAATCgg                   <  1:1566091/150‑1 (MQ=255)
                                                                                                                                       tCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCAGGCCACCGGGCCACGCGGATAGCGCTTCAGCCGCTCGGCCAGTTCGCGCCCGGGAATCGGc                  >  1:1654150/1‑150 (MQ=255)
                                                                                                                                       tCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCAGGCCACCGGGCCACGCGGATAGCGCTTCAGCCGCTCGGCCAGTTCGCGCCCGGGAATCGGc                  >  1:294841/1‑150 (MQ=255)
                                                                                                                                       tCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCAGGCCACCGGGCCACGCGGATAGCGCTTCAGCCGCTCGGCCAGTTCGCGCCCGGGAATCGGc                  >  2:2536283/1‑150 (MQ=255)
                                                                                                                                         aaaGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCAGGCCACCGGGCCACGCGGATAGCGCTTCAGCCGCTCGGCCAGTTCGCGCCCGGGAATCGGCCg                >  1:180463/1‑150 (MQ=255)
                                                                                                                                         aaaGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCAGGCCACCGGGCCACGCGGATAGCGCTTCAGCCGCTCGGCCAGTTCGCGCCCGGGAATCGGCCg                >  1:4175459/1‑150 (MQ=255)
                                                                                                                                            ggTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCAGGCCACCGGGCCACGCGGATAGCGCTTCAGCCGCTCGGCCAGTTCGCGCCCGGGAATCGGCCGGTc             <  1:4349116/150‑1 (MQ=255)
                                                                                                                                             gTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCAGGCCACCGGGCCACGCGGATAGCGCTTCAGCCGCTCGGCCAGTTCGCGCCCGGGAATCGGCCGGTCg            >  2:3726053/1‑150 (MQ=255)
                                                                                                                                             gTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCAGGCCACCGGGCCACGCGGATAGCGCTTCAGCCGCTCGGCCAGTTCGCGCCCGGGAATCGGCCGGTCg            >  1:3259704/1‑150 (MQ=255)
                                                                                                                                                gtggtgTTCCTGGGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCAGGCCACCGGGCCACGCGGATAGCGCTTCAGCCGCTCGGCCAGTTCGCGCCCGGGAATCGGCCGGTCGAAc         <  1:2639130/150‑1 (MQ=255)
                                                                                                                                                   gtgtTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCAGGCCACCGGGCCACGCGGATAGCGCTTCAGCCGCTCGGCCAGTTCGCGCCCGGGAATCGGCCGGTCGAACAgg      <  1:2469071/150‑1 (MQ=255)
                                                                                                                                                   gtgtTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCAGGCCACCGGGCCACGCGGATAGCGCTTCAGCCGCTCGGCCAGTTCGCGCCCGGGAATCGGCCGGTCGAACAgg      <  2:4132892/150‑1 (MQ=255)
                                                                                                                                                    tgtTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCAGGCCACCGGGCCACGCGGATAGCGCTTCAGCCGCTCGGCCAGTTCGCGCCCGGGAATCGGCCGGTCGAACAGGt     >  2:2187886/1‑150 (MQ=255)
                                                                                                                                                     gtTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCAGGCCACCGGGCCACGCGGATAGCGCTTCAGCCGCTCGGCCAGTTCGCGCCCGGGAATCGGCCGGTCGAACAGGTa    >  1:2900935/1‑150 (MQ=255)
                                                                                                                                                      ttCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCAGGCCACCGGGCCACGCGGATAGCGCTTCAGCCGCTCGGCCAGTTCGCGCCCGGGAATCGGCCGGTCGAACAGGTAg   >  2:2083430/1‑150 (MQ=255)
                                                                                                                                                       tCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCGTCGATAGTGCCGGAGGAGTACGCCGCTGTCTGGCCACCGGGCCACGCGGATAGCGCGTCAGCCGCTCGGCCCGTTCGCGCCCGGGTCGCGGCCGGGCGAACAGGTAGc  >  1:3233917/1‑150 (MQ=255)
                                                                                                                                                       |                                                                                                                                                     
CGCAGTGGGCATGACGTTTTTGTTCACCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCAGGCCACCGGGCCACGCGGATAGCGCTTCAGCCGCTCGGCCAGTTCGCGCCCGGGAATCGGCCGGTCGAACAGGTAGC  >  NC_002947/403150‑403448

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