Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 403,295 G→C 5.6% intergenic (‑95/+74) 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,2950GC 5.6% 482.2 / 15.1 142intergenic (‑95/+74)msrA/PP_0337peptide methionine sulfoxide reductase MsrA/sensory box protein
Reads supporting (aligned to +/- strand):  ref base G (82/52);  new base C (3/5);  total (85/57)
Fisher's exact test for biased strand distribution p-value = 2.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.03e-01

Only 100 of 150 total aligned reads displayed.

GTTCCGCAGTGGGCATGACGTTTTTGTTCACCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCAGGCCACCGGGCCACGCGGATAGCGCTTCAGCCGCTCGGCCAGTTCGCGCCCGGGAATCGGCCGGTCGAACAGG  >  NC_002947/403146‑403444
                                                                                                                                                       |                                                                                                                                                     
gTTCCGCAGTGGGCATGACGTTTTTGTTCACCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTgtgg                                                                                                                                                       >  2:2011716/1‑150 (MQ=255)
    cGCAGTGGGCATGACGTTTTTGTTCACCAGGATTTCCGAACGCAGGACCATGTCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGtt                                                                                                                                                   <  2:4355207/150‑1 (MQ=255)
        gTGGGCATGACGTTTTTGTTCACCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCtg                                                                                                                                               <  1:906243/150‑1 (MQ=255)
         tGGGCATGACGTTTTTGTTCACCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCtgt                                                                                                                                              <  2:2645965/150‑1 (MQ=255)
         tGGGCATGACGTTTTTGTTCACCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCtgt                                                                                                                                              <  2:2260368/150‑1 (MQ=255)
          gggCATGACGTTTTTGTTCACCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCtgtg                                                                                                                                             <  1:312262/150‑1 (MQ=255)
                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)
                 aCGTTTTTGTTCACCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCACAGGAACACCACAACCTTTGAAACCTGTGCTGTACCTGTGGGAGCgg                                                                                                                                      <  1:894617/150‑1 (MQ=21)
                  cttttttGTACACCTGTATTTCGGAACGCAGTACCTTGGCCGTTCCTATGTGTTTATGCAGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCggg                                                                                                                                     <  1:347522/148‑1 (MQ=37)
                   gTTTTTGTTCACCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGc                                                                                                                                    >  2:3573576/1‑150 (MQ=255)
                   gTTTTTGTTCACCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGc                                                                                                                                    <  2:3672769/150‑1 (MQ=255)
                       ttGTTCACCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGc                                                                                                                                <  1:2299515/150‑1 (MQ=255)
                       ttGTTCACCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGc                                                                                                                                <  1:1938178/150‑1 (MQ=255)
                          ttCACCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCcg                                                                                                                             >  2:2227649/1‑150 (MQ=255)
                             aCCAGGTTTTCCGAACGCTGGACCTGGGCGTTTCCTCTGTGTTTTTTCTGGGGGGTTGGGGTTTGGGCC‑‑TTGGGGGGCATGCCCGCTCCCACAGGTACAGCACCGGTTTCAATGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGa                                                                                                                          <  1:499510/150‑1 (MQ=25)
                             aCCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGa                                                                                                                          >  2:180876/1‑150 (MQ=255)
                               ctggATTTCCGAACGCTGGACCATGGCCGTTCCTATGTGTTTATTCAGTGGGTTAGGGGTTTGGGGC‑‑TTCGCGGGCTTGCCCGCTCCCACATGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGaaa                                                                                                                        <  1:2543641/148‑1 (MQ=37)
                               cAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTCGGGGTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGaaa                                                                                                                        >  2:4349116/1‑150 (MQ=255)
                               cAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGaaa                                                                                                                        >  2:2342141/1‑150 (MQ=255)
                               cAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGaaa                                                                                                                        >  1:2721780/1‑150 (MQ=255)
                               cAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGaaa                                                                                                                        >  2:2757016/1‑150 (MQ=255)
                                aGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCATTCCCGCTCCCACAGGTACAGCACCGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAg                                                                                                                       >  1:1999375/1‑150 (MQ=255)
                                aGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAg                                                                                                                       >  1:632825/1‑150 (MQ=255)
                                         gAACGCAGGCCCATGGCCGGTCCTCTGTGTTTATGCAGTGGGGGAGGGGTTGGGGCC‑‑GTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGGTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAAt                                                                                                              <  2:1421942/150‑1 (MQ=37)
                                         gAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCCTGCCCCCGCCCACACGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGGGCGGGCGTGCCCGCGAAAGGGCCCGCAt                                                                                                              >  1:2264354/1‑150 (MQ=39)
                                          aaCGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATg                                                                                                             >  1:3260434/1‑150 (MQ=255)
                                          aaCGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATg                                                                                                             >  2:1072617/1‑150 (MQ=255)
                                               aGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTc                                                                                                        >  2:3014782/1‑150 (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)
                                                              ccTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTACAACGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGa                                                                                         <  1:286244/150‑1 (MQ=255)
                                                                tctGTGTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGAtt                                                                                       >  2:2591727/1‑150 (MQ=255)
                                                                tctGTGTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGAtt                                                                                       >  2:2458197/1‑150 (MQ=255)
                                                                 ctGTGTTTATTCTGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTa                                                                                      <  2:1231213/150‑1 (MQ=255)
                                                                 ctGTGTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAAAGGCTTCTATAGTGCCCGAGGATTa                                                                                      >  1:2412939/1‑150 (MQ=255)
                                                                  tgtgTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTAc                                                                                     >  1:3375111/1‑150 (MQ=255)
                                                                  tgtgTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTAc                                                                                     >  2:1566091/1‑150 (MQ=255)
                                                                   tttttttttCAGTGGGGTTGGGGTTTGGGCC‑‑TTCGCGGGAATGCCCGCTCCCCCTGGTACAGCACAGGTTTCAATGGTTGTGGTGTTCGTGTGGGAGGGGGCTTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACg                                                                                    <  1:980168/147‑1 (MQ=255)
                                                                   gtgtTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACg                                                                                    >  1:4023084/1‑150 (MQ=255)
                                                                             cggggggTCGGGGTTGGGGCC‑‑TTCGGGGGCATGCCCGCTCCCCCAGGTACAGCACCGGTTCCAATGGTTGTGGTGTTCGTGTGGGCGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCAg                                                                          <  1:3941727/147‑1 (MQ=255)
                                                                             aGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCAg                                                                          >  1:2906667/1‑150 (MQ=255)
                                                                             aGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCAg                                                                          >  2:323403/1‑150 (MQ=255)
                                                                             aGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCAAGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGATCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTGCTATAGAGCCCGAGGATTACGCAGCTGTCAg                                                                          >  2:2046365/1‑150 (MQ=255)
                                                                              gTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCAgg                                                                         >  2:4048404/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)
                                                                                            tGGGCC‑‑TTCGCGGGCAGCCCCGCTCCCCCAGGTACAGCACTGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCAGGCCACCGGGCCAcgc                                                           <  1:3088941/150‑1 (MQ=255)
                                                                                                tgtgggagcgGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCAGGCCACCGGGCCACGCgg                                                         >  2:2471616/8‑150 (MQ=255)
                                                                                                cc‑‑TTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTACTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCACGAGGATTACGCAGCTGTCAGGCCACCGGGCCACGCGGAt                                                       >  1:2575366/1‑150 (MQ=255)
                                                                                                   ctcggggggCTTGCCCGCCCCCCCAGGTACAGCACAGGTTACACCGGTTGTGGTGTTCCTGTGGGAGGGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCAGGCCACCGGGCCACGCGGATa                                                      <  1:3803702/144‑1 (MQ=255)
                                                                                                    ggagcgGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCGGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCAGGCCACCGGGCCACGCGGATAg                                                     >  2:3686181/4‑150 (MQ=255)
                                                                                                     gagcgGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCAGGCCACCGGGCCACGCGGATAgc                                                    >  2:2705237/3‑150 (MQ=255)
                                                                                                        cgGGCTTGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTGGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCAGGCCACCGGGCCACGCGGATAGCGCt                                                 <  2:101339/150‑1 (MQ=255)
                                                                                                           gCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCAGGCCACCGGGCCACGCGGATAGCGCTTCa                                              >  2:4403245/1‑150 (MQ=255)
                                                                                                            ctggCCCGCGCCCACTGTTACAGCACTGGTTTCACCGGTTGTGGTGTTCCTGTGGGTGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCAGGCCACCGGGCCACGCGGATAGCGCTTCAg                                             <  1:1217942/147‑1 (MQ=255)
                                                                                                             aTGCCCGCTCCCCCAGGTACAGCACAGGTTCCAACGGTTGTGGGGTTCCTGTGGGAGCGGGGGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCAGGCCACCGGGCCACGCGGATAGCGCTTCAGc                                            <  2:910465/150‑1 (MQ=255)
                                                                                                             aTGCCCGCTCCCACAGGTACAGCACAGGTTTCACAGGGTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCGTCTATAGTGCCCGAGGATTACGCAGCTGTCAGGCCACCGGGCCACGCGGATAGCGCTTCAGc                                            >  2:1698191/1‑150 (MQ=255)
                                                                                                               gCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCAGGCCACCGGGCCACGCGGATAGCGCTTCAGCCg                                          >  1:1008504/1‑150 (MQ=255)
                                                                                                                  cGCTCCCACAGGTACAGCACAGGTTTCAACGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCAGGCCACCGGGCCACGCGGATAGCGCTTCAGCCGCTc                                       <  1:2732257/150‑1 (MQ=255)
                                                                                                                   gCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCAGGCCACCGGGCCACGCGGATAGCGCTTCAGCCGCTCg                                      >  2:1658374/1‑150 (MQ=255)
                                                                                                                   gCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCAGGCCACCGGGCCACGCGGATAGCGCTTCAGCCGCTCg                                      <  2:254235/150‑1 (MQ=255)
                                                                                                                    cTCCCACAGGAACACCACAACCTTTGAAACCTGTGCTGTACCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCAGGCCACCGGGCCACGCGGATAGCGCTTCAGCCGCTCgg                                     >  1:3388624/1‑150 (MQ=18)
                                                                                                                      cccACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCAGGCCACCGGGCCACGCGGATAGCGCTTCAGCCGCTCGGcc                                   >  2:797295/1‑150 (MQ=255)
                                                                                                                       ccACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCAGGCCACCGGGCCACGCGGATAGCGCTTCAGCCGCTCGGCCa                                  <  1:4403245/150‑1 (MQ=255)
                                                                                                                       ccACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCAGGCCACCGGGCCACGCGGATAGCGCTTCAGCCGCTCGGCCa                                  >  1:3096991/1‑150 (MQ=255)
                                                                                                                        cacaGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCAGGCCACCGGGCCACGCGGATAGCGCTTCAGCCGATCGGCCAg                                 >  1:4127390/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)
                                                                                                                              tACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCAGGCCACCGGGCCACGCGGATAGCGCTTCAGCCGCTCGGCCAGTTcgcg                           >  2:3536498/1‑150 (MQ=255)
                                                                                                                                  gCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCAGGCCACCGGGCCACGCGGATAGCGCTTCAGCCGCTCGGCCAGTTCGCGCCCg                       >  2:2628235/1‑150 (MQ=255)
                                                                                                                                  gCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCAGGCCACCGGGCCACGCGGATAGCGCTTCAGCCGCTCGGCCAGTTCGCGCCCg                       <  2:2721780/150‑1 (MQ=255)
                                                                                                                                    acaGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCAGGCCACCGGGCCACGCGGATAGCGCTTCAGCCGCTCGGCCAGTTCGCGCCCggg                     >  1:1673241/1‑150 (MQ=255)
                                                                                                                                    acaGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCAGGCCACCGGGCCACGCGGATAGCGCTTCAGCCGCTCGGCCAGTTCGCGCCCggg                     >  2:3787123/1‑150 (MQ=255)
                                                                                                                                     caGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCAGGCCACCGGGCCACGCGGATAGCGCTTCAGCCGCTCGGCCAGTTCGCGCCCGGGa                    >  2:376090/1‑150 (MQ=255)
                                                                                                                                        gTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCAGGCCACCGGGCCACGCGGATAGCGCTTCAGCCGCTCGGCCAGTTc                              >  1:2051527/1‑137 (MQ=255)
                                                                                                                                          ttCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCAGGCCACCGGGCCACGCGGATAGCGCTTCAGCCGCTCGGCCAGTTCGCGCCCGGGAATCgg               <  1:1566091/150‑1 (MQ=255)
                                                                                                                                           tCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCAGGCCACCGGGCCACGCGGATAGCGCTTCAGCCGCTCGGCCAGTTCGCGCCCGGGAATCGGc              >  1:1654150/1‑150 (MQ=255)
                                                                                                                                             aaaGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCAGGCCACCGGGCCACGCGGATAGCGCTTCAGCCGCTCGGCCAGTTCGCGCCCGGGAATCGGCCg            >  1:180463/1‑150 (MQ=255)
                                                                                                                                             aaaGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCAGGCCACCGGGCCACGCGGATAGCGCTTCAGCCGCTCGGCCAGTTCGCGCCCGGGAATCGGCCg            >  1:2117591/1‑150 (MQ=255)
                                                                                                                                             aaaGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCAGGCCACCGGGCCACGCGGATAGCGCTTCAGCCGCTCGGCCAGTTCGCGCCCGGGAATCGGCCg            >  1:4175459/1‑150 (MQ=255)
                                                                                                                                                ggTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCAGGCCACCGGGCCACGCGGATAGCGCTTCAGCCGCTCGGCCAGTTCGCGCCCGGGAATCGGCCGGTc         <  2:3173502/150‑1 (MQ=255)
                                                                                                                                                 gTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCAGGCCACCGGGCCACGCGGATAGCGCTTCAGCCGCTCGGCCAGTTCGCGCCCGGGAATCGGCCGGTCg        >  2:3726053/1‑150 (MQ=255)
                                                                                                                                                 gTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCAGGCCACCGGGCCACGCGGATAGCGCTTCAGCCGCTCGGCCAGTTCGCGCCCGGGAATCGGCCGGTCg        >  1:3259704/1‑150 (MQ=255)
                                                                                                                                                    gtggtgTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCAGGCCACCGGGCCACGCGGATAGCGCTTCAGCCGCTCGGCCAGTTCGCGCCCGGGAATCGGCCGGTCGAAc     >  1:2146244/1‑150 (MQ=255)
                                                                                                                                                      ggtgTTCCTGTGGGCGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCAGGCCACCGGGCCACGCGGATAGCGCTTCAGCCGCTCGGCCAGTTCGCGCCCGGGAATCGGCCGGTCGAACAg   <  1:3573576/150‑1 (MQ=255)
                                                                                                                                                       gtgtTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCAGGCCACCGGGCCACGCGGATAGCGCTTCAGCCGCTCGGCCAGTTCGCGCCCGGGAATCGGCCGGTCGAACAgg  <  2:4132892/150‑1 (MQ=255)
                                                                                                                                                       gtgtTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCAGGCCACCGGGCCACGCGGATAGCGCTTCAGCCGCTCGGCCAGTTCGCGCCCGGGAATCGGCCGGTCGAACAgg  <  1:2469071/150‑1 (MQ=255)
                                                                                                                                                       |                                                                                                                                                     
GTTCCGCAGTGGGCATGACGTTTTTGTTCACCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCAGGCCACCGGGCCACGCGGATAGCGCTTCAGCCGCTCGGCCAGTTCGCGCCCGGGAATCGGCCGGTCGAACAGG  >  NC_002947/403146‑403444

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: