Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,530,544 A→C 5.6% H161P (CAC→CCC) ‡ mnmA → tRNA 2‑thiouridine(34) synthase MnmA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,530,5440AC 5.6% 998.3 / 33.6 288H161P (CAC→CCC) ‡mnmAtRNA 2‑thiouridine(34) synthase MnmA
Reads supporting (aligned to +/- strand):  ref base A (140/132);  new base C (6/10);  total (146/142)
Fisher's exact test for biased strand distribution p-value = 3.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.35e-01

Only 100 of 289 total aligned reads displayed.

AATCAAGTTCAAGGCGTTCCTCGACTACGCCCTGTCGCTGGGTGCCGACCTGATCGCCACCGGCCACTATGTGCGCCGCCGCGATACTGGCGCGCTCACCGAACTGCTCAAGGGCCTGGACCCGAACAAGGACCAGAGCTACTTCCTGCACGCCGTCGGCGGCCAGGAAATCGCCCGCACACTGTTCCCGGTCGGCGAACTGGAAAAGCCCGAAGTACGCGCCATCGCCGAAAAACACGGCCTGGCCACCGCCAAGAAGAAAGATTCCACCGGCATCTGCTTCATTGGCGAGCGCCGC  >  NC_002947/4530395‑4530692
                                                                                                                                                     |                                                                                                                                                    
  tCAAGTTCAAGGCGTTCCTCGACTACGCCCTGTCGCTGGGTGCCGACCTGATCGCCACCGGCCACTATGTGCGCCGCCGCGATACTGGCGCGCTCACCGAACTGCTCAAGGGCCTGGACCCGAACAAGGACCAGAGCTACTTCCTGCACg                                                                                                                                                    >  2:3371110/1‑150 (MQ=255)
  tCAAGTTCAAGGCGTTCCTCGACTACGCCCTGTCGCTGGGTGCCGACCTGATCGCCACCGGCCACTATGTGCGCCGCCGCGATACTGGCGCGCTCACCGAACTGCTCAAGGGCCTGGACCCGAACAAGGACCAGAGCTACTTCCTGCACg                                                                                                                                                    >  1:2622422/1‑150 (MQ=255)
  tCAAGTTCAAGGCGTTCCTCGACTACGCCCTGTCGCTGGGTGCCGACCTGATCGCCACCGGCCACTATGTGCGCCGCCGCGATACTGGCGCGCTCACCGAACTGCTCAAGGGCCTGGACCCGAACAAGGACCAGAGCTACTTCCTGCACg                                                                                                                                                    >  1:763966/1‑150 (MQ=255)
  tCAAGTTCAAGGCGTTCCTCGACTACGCCCTGTCGCTGGGTGCCGACCTGATCGCCACCGGCCACTATGTGCGCCGCCGCGATACTGGCGCGCTCACCGAACAGCTCAACGGCCAGGACCCGAACAAGGACCAGAGCTACTTCCTGCACg                                                                                                                                                    >  2:1204215/1‑150 (MQ=255)
    aaGGTCAAGGCGTTCCTCGACTACGCCCTGTCGCTGGGTGCCGACCTGATCGCCACCGGCCACTATGTGCGCCGCCGCGATACTGGCGCGCTCACCGAACTGCTCAAGGGCCTGGACCCGAACAAGGACCAGAGCTACTTCCTGCACGcc                                                                                                                                                  <  2:5437269/150‑1 (MQ=255)
        tCAAGGCGTTCCTCGACTACGCCCTGTCGCTGGGTGCCGACCTGATCGCCACCGGCCACTATGTGCGCCGCCGCGATACTGGCGCGCTCACCGAACTGCTCAAGGGCCTGGACCCGAACAAGGACCAGAGCTACTTCCTGCACGCCGTcg                                                                                                                                              >  1:5119378/1‑150 (MQ=255)
         cAAGGCGTTCCTCGACTACGCCCTGTCGCTGGGTGCCGACCTGATCGCCACCGGCCACTATGTGCGCCGCCGCGATACTGGCGCGCTCACCGAACTGCTCAAGGGCCTGGACCCGAACAAGGACCAGAGCTACTTCCTGCACGCCGTcgg                                                                                                                                             >  1:4892801/1‑150 (MQ=255)
          aaGGCGTTCCTCGACTACGCCCTGTCGCTGGGTGCCGACCTGATCGCCACCGGCCACTATGTGCGCCGCCGCGATACTGGCGCGCTCACCGAACTGCTCAAGGGCCTGGACCCGAACAAGGACCAGAGCTACTTCCTGCACGCCGTcggc                                                                                                                                            <  2:5476251/150‑1 (MQ=255)
           aGGCGTTCCTCGACTACGCCCTGTCGCTGGGTGCCGACCTGATCGCCACCGGCCACTATGTGCGCCGCCGCGATACTGGCGCGCTCACCGAACTGCTCAAGGGCCTGGACCCGAACAAGGACCAGAGCTACTTCCTGCACGCCGTcggcg                                                                                                                                           <  1:907611/150‑1 (MQ=255)
           aGGCGTTCCTCGACTACGCCCTGTCGCTGGGTGCCGACCTGATCGCCACCGGCCACTATGTGCGCCGCCGCGATACTGGCGCGCTCACCGAACTGCTCAAGGGCCTGGACCCGAACAAGGACCAGAGCTACTTCCTGCACGCCGTcggcg                                                                                                                                           <  2:4974480/150‑1 (MQ=255)
               gTTCCTCGACTACGCCCTGTCGCTGGGTGCCGACCTGATCGCCACCGGCCACTATGTGCGCCGCCGCGATACTGGCGCGCTCACCGAACTGCTCAAGGGCCTGGACCCGAACAAGGACCAGAGCTACTTCCTGCACGCCGTCGGCGGCCa                                                                                                                                       <  2:503883/150‑1 (MQ=255)
                 tCCTCGACTACGCCCTGTCGCTGGGTGCCGACCTGATCGCCACCGGCCACTATGTGCGCCGCCGCGATACTGGCGCGCTCACCGAACTGCTCAAGGGCCTGGACCCGAACAAGGACCAGAGCTACTTCCTGGCCGCCGACGGCGTGCAgg                                                                                                                                     <  2:4924366/150‑1 (MQ=255)
                    tCGACTACGCCCTGTCGCTGGGTGCCGACCTGATCGCCACCGGCCACTATGTGCGCCGCCGCGATACTGGCGCGCTCACCGAACTGCTCAAGGGCCTGGACCCGAACAAGGACCAGAGCTACTTCCTGGCCGCCGACGGCGTGCAGGaag                                                                                                                                  <  1:2766638/150‑2 (MQ=255)
                      gacgacGCCCTGTCGCTGGGTGCCGACCTGATCGCCACCGGCCACTATGTGCGCCGCCGCGATACTGGCGCGCTCACCGAACTGCTCAAGGGCCTGGACCCGAACAAGGACCAGAGCTACTTCCTGCACGCCGTCGGCGGCCAGGAAATc                                                                                                                                <  2:1986345/150‑1 (MQ=255)
                      gACTACGCCCTGTCGCTGGGTGCCGACCTGATCGCCACCGGCCACTATGTGCGCCGCCGCGATACTGGCGCGCTCACCGAACTGCTCAAGGGCCTGGACCCGAACAAGGACCAGAGCTACTTCCTGGCCGCCGACGGCGTGCAGGaagta                                                                                                                                <  1:447671/150‑4 (MQ=255)
                      gACTACGCCCTGTCGCTGGGTGCCGACCTGATCGCCACCGGCCACTATGTGCGCCGCCGCGATACTGGCGCGCTCACCGAACTGCTCAAGGGCCTGGACCCGAACAAGGACCAGAGCTACTTCCTGCACGCCGTCGGCGGCCAGGAAATc                                                                                                                                >  1:2727560/1‑150 (MQ=255)
                      gACTACGCCCTGTCGCTGGGTGCCGACCTGATCGCCACCGGCCACTATGTGCGCCGCCGCGATACTGGCGCGCTCACCGAACTGCTCAAGGGCCTGGACCCGAACAAGGACCAGAGCTACTTCCTGCACGCCGTCGGCGGCCAGGAAATc                                                                                                                                <  1:2182465/150‑1 (MQ=255)
                      gACTACGCCCTGTCGCTGGGTGCCGACCTGATCGCCACCGGCCACTATGTGCGCCGCCGCGATACTGGCGCGCTCACCGAACTGCTCAAGGGCCTGGACCCGAACAAGGACCAGAGCTACTTCCTGCACGCCGTCGGCGGCCAGGAAATc                                                                                                                                >  2:312454/1‑150 (MQ=255)
                      gACTACGCCCTGTCGCTGGGTGCCGACCTGATCGCCACCGGCCACGATGTGCGCCGCCGCGATACGGGCGCGCTCACCGAACTGCTCAAGGGCCTGGACCCGAACAAGGACCAGAGCTACTTCCTGCACGCCGTCGGCGGCCAGGAAATc                                                                                                                                <  1:3596391/150‑1 (MQ=255)
                       aCTACGCCCTGTCGCTGGGTGCCGACCTGATCGCCACCGGCCACTATGTGCGCCGCCGCGATACTGGCGCGCTCACCGAACTGCTCAAGGGCCTGGACCCGAACAAGGACCAGAGCTACTTCCTGGCCGCCGACGGCGTGCAGGaagtag                                                                                                                               <  1:1918479/150‑5 (MQ=255)
                        cTACGCCCTGTCGCTGGGTGCCGACCTGATCGCCACCGGCCACTATGTGCGCCGCCGCGATACTGGCGCGCTCACCGAACTGCTCAAGGGCCTGGACCCGAACAAGGACCAGAGCTACTTCCTGCACGCCGTCGGCGGCCAGGAAATCGc                                                                                                                              <  1:2915141/150‑1 (MQ=255)
                         tACGCCCTGTCGCTGGGTGCCGACCTGATCGCCACCGGCCACTATGTGCGCCGCCGCGATACTGGCGCGCTCACCGAACTGCTCAAGGGCCTGGACCCGAACAAGGACCAGAGCTACTTCCTGCACGCCGTCGGCGGCCAGGAAATCGcc                                                                                                                             <  2:4622332/150‑1 (MQ=255)
                           cGCCCTGTCGCTGGGTGCCGACCTGATCGCCACCGGCCACTATGTGCGCCGCCGCGATACTGGCGCGCTCACCGAACTGCTCAAGGGCCTGGACCCGAACAAGGACCAGAGCTACTTCCTGCACGCCGTCGGCGGCCAGGAAATCGCCCg                                                                                                                           >  2:4792792/1‑150 (MQ=255)
                            gCCCTGTCGCTGGGTGCCGACCTGATCGCCACCGGCCACTATGTGCGCCGCCGCGATACTGGCGCGCTCACCGAACTGCTCAAGGGCCTGGACCCGAACAAGGACCAGAGCTACTTCCTGCACGCCGTCGGCGGCCAGGAAATCGCCCGc                                                                                                                          >  1:4712937/1‑150 (MQ=255)
                                  tCGCTGGGTGCCGACCTGATCGCCACCGGCCACTATGTGCGCCGCCGCGATACTGGCGCGCTCACCGAACTGCTCAAGGGCCTGTACCCGAACAAGGACCAGAGCTACTTCCTGCACGCCGTCGGCGGCCAGGAAATCGCCCGCACACTg                                                                                                                    <  2:4057071/150‑1 (MQ=255)
                                  tCGCTGGGTGCCGACCTGATCGCCACCGGCCACTATGTGCGCCGCCGCGATACTGGCGCGCTCACCGAACTGCTCAAGGGCCTGGACCCGAACAAGGACCAGAGCTACTTCCTGCACGCCGTCGGCGGCCAGGAAATCGCCCGCACACTg                                                                                                                    >  1:838984/1‑150 (MQ=255)
                                     cTGGGTGCCGACCTGATCGCCACCGGCCACTATGTGCGCCGCCGCGATACTGGCGCGCTCACCGAACTGCTCAAGGGCCTGGACCCGAACAAGGACCAGAGCTACTTCCTGCACGCCGTCGGCGGCCAGGAAATCGCCCGCACACTGTTc                                                                                                                 >  1:4977778/1‑150 (MQ=255)
                                      tGGGTGCCGACCTGATCGCCACCGGCCACTATGTGCGCCGCCGCGATACTGGCGCGCTCACCGAACTGCTCAAGGGCCTGGACCCGAACAAGGACCAGAGCTACTTCCTGCACGCCGTCGGCGGCCAGGAAATCGCCCGCACACTGTTcc                                                                                                                >  2:1406297/1‑150 (MQ=255)
                                        ggTGCCGACCTGATCGCCACCGGCCACTATGTGCGCCGCCGCGATACTGGCGCGCTCACCGAACTGCTCAAGGGCCTGGACCCGAACAAGGACCAGAGCTACTTCCTGCACGCCGTCGGCGGCCAGGAAATCGCCCGCACACTGTTCCCg                                                                                                              >  1:3653644/1‑150 (MQ=255)
                                         gTGCCGACCTGATCGCCACCGGCCACTATGTGCGCCGCCGCGATACTGGCGCGCTCACCGAACTGCTCAAGGGCCTGGACCCGAACAAGGACCAGAGCTACTTCCTGCACGCCGTCGGCGGCCAGGAAATCGCCCGCACACTGTTCCCgg                                                                                                             >  2:4660402/1‑150 (MQ=255)
                                            ccGACCTGATCGCCACCGGCCACTATGTGCGCCGCCGCGATACTGGCGCGCTCACCGAACTGCTCAAGGGCCTGGACCCGAACAAGGACCAGAGCTACTTCCTGCACGCCGTCGGCGGCCAGGAAATCGCCCGCACACTGTTCCCGGTCg                                                                                                          <  2:4417683/150‑1 (MQ=255)
                                              gACCTGATCGCCACCGGCCACTATGTGCGCCGCCGCGATACTGGCGCGCTCACCGAACTGCTCAAGGGCCTGGACCCGAACAAGGACCAGAGCTACTTCCTGCACGCCGTCGGCGGCCAGGAAATCGCCCGCACACTGTTCCCGGTCGGc                                                                                                        <  2:3447482/150‑1 (MQ=255)
                                                  tGATCGCCACCGGCCACTATGTGCGCCGCCGCGATACTGGCGCGCTCACCGAACTGCTCAAGGGCCTGGACCCGAACAAGGACCAGAGCTACTTCCTGCACGCCGTCGGCGGCCAGGAAATCGCCCGCACACTGTTCCCGGTCGGCGAAc                                                                                                    <  1:2470706/150‑1 (MQ=255)
                                                  tGATCGCCACCGGCCACTATGTGCGCCGCCGCGATACTGGCGCGCTCACCGAACTGCTCAAGGGCCTGGACCCGAACAAGGACCAGAGCTACTTCCTGCACGCCGTCGGCGGCCAGGAAATCGCCCGCACACTGTTCCCGGTCGGCGAAc                                                                                                    <  2:2289714/150‑1 (MQ=255)
                                                  tGATCGCCACCGGCCACTATGTGCGCCGCCGCGATACTGGCGCGCTCACCGAACTGCTCAAGGGCCTGGACCCGAACAAGGACCAGAGCTACTTCCTGCACGCCGTCGGCGGCCAGGAAATCGCCCGCACACTGTTCCCGGTCGGCGAAc                                                                                                    <  2:62514/150‑1 (MQ=255)
                                                   gATCGCCACCGGCCACTATGTGCGCCGCCGCGATACTGGCGCGCTCACCGAACTGCTCAAGGGCCTGGACCCGAACAAGGACCAGAGCTACTTCCTGCACGCCGTCGGCGGCCAGGAAATCGCCCGCACACTGTTCCCGGTCGGCGAACt                                                                                                   >  1:2377384/1‑150 (MQ=255)
                                                   gATCGCCACCGGCCACTATGTGCGCCGCCGCGATACTGGCGCGCTCACCGAACTGCTCAAGGGCCTGGACCCGAACAAGGACCAGAGCTACTTCCTGCACGCCGTCGGCGGCCAGGAAATCGCCCGCACACTGTTCCCGGTCGGCGAACt                                                                                                   <  1:2253827/150‑1 (MQ=255)
                                                    aTCGCCACCGGCCACTATGTGCGCCGCCGCGATACTGGCGCGCTCACCGAACTGCTCAAGGGCCTGGACCCGAACAAGGACCAGAGCTACTTCCTGCACGCCGTCGGCGGCCAGGAAATCGCCCGCACACTGTTCCCGGTCGGCGAACTg                                                                                                  <  1:3963697/150‑1 (MQ=255)
                                                     tCGCCACCGGCCACTATGTGCGCCGCCGCGATACTGGCGCGCTCACCGAACTGCTCAAGGGCCTGGACCCGAACAAGGACCAGAGCTACTTCCTGCACGCCGTCGGCGGCCAGGAAATCGCCCGCACACTGTTCCCGGTCGGCGAACTgg                                                                                                 <  1:3408283/150‑1 (MQ=255)
                                                       gCCACCGGCCACTATGTGCGCCGCCGCGATACTGGCGCGCTCACCGAACTGCTCAAGGGCCTGGACCCGAACAAGGACCAGAGCTACTTCCTGCACGCCGTCGGCGGCCAGGAAATCGCCCGCACACTGTTCCCGGTCGGCGAACTGGaa                                                                                               <  1:4541878/150‑1 (MQ=255)
                                                        ccACCGGCCACTATGTGCGCCGCCGCGATACTGGCGCGCTCACCGAACTGCTCAAGGGCCTGGACCCGAACAAGGACCAGAGCTACTTCCTGCACGCCGTCGGCGGCCAGGAAATCGCCCGCACACTGTTCCCGGTCGGCGAACTGGaaa                                                                                              >  1:3762087/1‑150 (MQ=255)
                                                        ccACCGGCCACTATGTGCGCCGCCGCGATACTGGCGCGCTCACCGAACTGCTCAAGGGCCTGGACCCGAACAAGGACCAGAGCTACTTCCTGCACGCCGTCGGCGGCCAGGAAATCGCCCGCACACTGTTCCCGGTCGGCGAACTGGaaa                                                                                              >  2:3134414/1‑150 (MQ=255)
                                                        ccACCGGCCACTATGTGCGCCGCCGCGATACTGGCGCGCTCACCGAACTGCTCAAGGGCCTGGACCCGAACAAGGACCAGAGCTACTTCCTGCACGCCGTCGGCGGCCAGGAAATCGCCCGCACACTGTTCCCGGTCGGCGAACTGGaaa                                                                                              >  2:2937536/1‑150 (MQ=255)
                                                         cACCGGCCACTATGTGCGCCGCCGCGATACTGGCGCGCTCACCGAACTGCTCAAGGGCCTGGACCCGAACAAGGACCAGAGCTACTTCCTGCACGCCGTCGGCGGCCAGGAAATCGCCCGCACACTGTTCCCGGTCGGCGAACTGGaaaa                                                                                             >  2:3985022/1‑150 (MQ=255)
                                                         cACCGGCCACTATGTGCGCCGCCGCGATACTGGCGCGCTCACCGAACTGCTCAAGGGCCTGGACCCGAACAAGGACCAGAGCTACTTCCTGCACGCCGTCGGCGGCCAGGAAATCGCCCGCACACTGTTCCCGGTCGGCGAACTGGaaaa                                                                                             >  1:2077613/1‑150 (MQ=255)
                                                            cGGCCACTATGTGCGCCGCCGCGATACTGGCGCGCTCACCGAACTGCTCAAGGGCCTGGACCCGAACAAGGACCAGAGCTACTTCCTGCACGCCGTCGGCGGCCAGGAAATCGCCCGCACACTGTTCCCGGTCGGCGAACTGGAAAAGcc                                                                                          <  2:1472616/150‑1 (MQ=255)
                                                             ggCCACTATGTGCGCCGCCGCGATACTGGCGCGCTCACCGAACTGCTCAAGGGCCTGGACCCGAACAAGGACCAGAGCTACTTCCTGCACGCCGTCGGCGGCCAGGAAATCGCCCGCACACTGTTCCCGGTCGGCGAACTGGAAAAGccc                                                                                         <  1:1718534/150‑1 (MQ=255)
                                                             ggCCACTATGTGCGCCGCCGCGATACTGGCGCGCTCACCGAACTGCTCAAGGGCCTGGACCCGAACAAGGACCAGAGCTACTTCCTGCACGCCGTCGGCGGCCAGGAAATCGCCCGCACACTGTTCCCGGTCGGCGAACTGGAAAAGccc                                                                                         >  2:3743467/1‑150 (MQ=255)
                                                               ccACTATGTGCGCCGCCGCGATACTGGCGCGCTCACCGAACTGCTCAAGGGCCTGGACCCGAACAAGGACCAGAGCTACTTCCTGCACGCCGTCGGCGGCCAGGAAATCGCCCGCACACTGTTCCCGGTCGGCGAACTGGAAAAGCCCGa                                                                                       >  1:234135/1‑150 (MQ=255)
                                                                cACTATGTGCGCCGCCGCGATACTGGCGCGCTCACCGAACTGCTCAAGGGCCTGGACCCGAACAAGGACCAGAGCTACTTCCTGCACGCCGTCGGCGGCCAGGAAATCGCCCGCACACTGTTCCCGGTCGGCGAACTGGAAAAGCCCGaa                                                                                      <  2:187736/150‑1 (MQ=255)
                                                                cACTATGTGCGCCGCCGCGATACTGGCGCGCTCACCGAACTGCTCAAGGGCCTGGACCCGAACAAGGACCAGAGCTACTTCCTGCACGCCGTCGGCGGCCAGGAAATCGCCCGCACACTGTTCCCGGTCGGCGAACTGGAAAAGCCCGaa                                                                                      <  1:3460410/150‑1 (MQ=255)
                                                                 aCTATGTGCGCCGCCGCGATACTGGCGCGCTCACCGAACTGCTCAAGGGCCTGGACCCGAACAAGGACCAGAGCTACTTCCTGCACGCCGTCGGCGGCCAGGAAATCGCCCGCACACTGTTCCCGGTCGGCGAACTGGAAAAGCCCGAAg                                                                                     >  2:2131412/1‑150 (MQ=255)
                                                                        gcgcCGCCGCGATACTGGCGCGCTCACCGAACTGCTCAAGGGCCTGGACCCGAACAAGGACCAGAGCTACTTCCTGCACGCCGTCGGCGGCCAGGAAATCGCCCGCACACTGTTCCCGGTCGGCGAACTGGAAAAGCCCGAAGTACgcgc                                                                              >  1:1770615/1‑150 (MQ=255)
                                                                             gccgcGATACTGGCGCGCTCACCGAACTGCTCAAGGGCCTGGACCCGAACAAGGACCAGAGCTACTTCCTGCACGCCGTCGGCGGCCAGGAAATCGCCCGCACACTGTTCCCGGTCGGCGAACTGGAAAAGCCCGAAGTACGCGCCATCg                                                                         >  1:2869820/1‑150 (MQ=255)
                                                                               cgcgATACTGGCGCGCTCACCGAACTGCTCAAGGGCCTGGACCCGAACAAGGACCAGAGCTACTTCCTGCACGCCGTCGGCGGCCAGGAAATCGCCCGCACACTGTTCCCGGTCGGCGAACTGGAAAAGCCCGAAGTACGCGCCATCGcc                                                                       <  1:4241330/150‑1 (MQ=255)
                                                                                   aTACTGGCGCGCTCACCGAACTGCTCAAGGGCCTGGACCCGAACAAGGACCAGAGCTACTTCCTGCACGCCGTCGGCGGCCAGGAAATCGCCCGCACACTGTTCCCGGTCGGCGAACTGGAAAAGCCCGAAGTACGCGCCATCGCCGaaa                                                                   <  2:1087392/150‑1 (MQ=255)
                                                                                       tGGCGCGCTCACCGAACTGCTCAAGGGCCTGGACCCGAACAAGGACCAGAGCTACTTCCTGCACGCCGTCGGCGGCCAGGAAATCGCCCGCACACTGTTCCCGGTCGGCGAACTGGAAAAGCCCGAAGTACGCGCCATCGCCGAAAAaca                                                               <  1:3179782/150‑1 (MQ=255)
                                                                                       tGGCGCGCTCACCGAACTGCTCAAGGGCCTGGACCCGAACAAGGACCAGAGCTACTTCCTGCACGCCGTCGGCGGCCAGGAAATCGCCCGCACACTGTTCCCGGTCGGCGAACTGGAAAAGCCCGAAGTACGCGCCATCGCCGAAAAaca                                                               >  2:2510730/1‑150 (MQ=255)
                                                                                         gcgcgcTCACCGAACTGCTCAAGGGCCTGGACCCGAACAAGGACCAGAGCTACTTCCTGCACGCCGTCGGCGGCCAGGATATCGCCCGCACACTGTTCCCGGTCGGCGAACTGGAAAAGCCCGAAGTACGCGCCATCGCCGAAAAACACg                                                             <  2:4497834/150‑1 (MQ=255)
                                                                                         gcgcgcTCACCGAACTGCTCAAGGGCCTGGACCCGAACAAGGACCAGAGCTACTTCCTGCACGCCGTCGGCGGCCAGGAAATCGCCCGCACACTGTTCCCGGTCGGCGAACTGGAAAAGCCCGAAGTACGCGCCATCGCCGAAAAACACg                                                             <  1:5607049/150‑1 (MQ=255)
                                                                                           gcgcTCACCGAACTGCTCAAGGGCCTGGACCCGAACAAGGACCAGAGCTACTTCCTGCACGCCGTCGGCGGCCAGGAAATCGCCCGCACACTGTTCCCGGTCGGCGAACTGGAAAAGCCCGAAGTACGCGCCATCGCCGAAAAACACGGc                                                           >  2:4815750/1‑150 (MQ=255)
                                                                                             gcTCACCGAACTGCTCAAGGGCCTGGACCCGAACAAGGACCAGAGCTACTTCCTGCACGCCGTCGGCGGCCAGGAAATCGCCCGCACACTGTTCCCGGTCGGCGAACTGGAAAAGCCCGAAGTACGCGCCATCGCCGAAAAACACGGCCt                                                         <  1:1487084/150‑1 (MQ=255)
                                                                                               tCACCGAACTGCTCAAGGGCCTGGACCCGAACAAGGACCAGAGCTACTTCCTGCACGCCGTCGGCGGCCAGGAAATCGCCCGCACACTGTTCCCGGTCGGCGAACTGGAAAAGCCCGAAGTACGCGCCATCGCCGAAAAACACGGCCTgg                                                       >  2:5166347/1‑150 (MQ=255)
                                                                                               tCACCGAACTGCTCAAGGGCCTGGACCCGAACAAGGACCAGAGCTACTTCCTGCACGCCGTCGGCGGCCAGGAAATCGCCCGCACACTGTTCCCGGTCGGCGAACTGGAAAAGCCCGAAGTACGCGCCATCGCCGAAAAACACGGCCTgg                                                       >  1:3238804/1‑150 (MQ=255)
                                                                                                cACCGAACTGCTCAAGGGCCTGGACCCGAACAAGGACCAGAGCTACTTCCTGCACGCCGTCGGCGGCCAGGAAATCGCCCGCACACTGTTCCCGGTCGGCGAACTGGAAAAGCCCGAAGTACGCGCCATCGCCGAAAAACACGGCCTGGc                                                      <  1:3353764/150‑1 (MQ=255)
                                                                                                 aCCGAACTGCTCAAGGGCCTGGACCCGAACAAGGACCAGAGCTACTTCCTGCACGCCGTCGGCGGCCAGGAAATCGCCCGCACACTGTTCCCGGTCGGCGAACTGGAAAAGCCCGAAGTACGCGCCATCGCCGAAAAACACGGCCTGGcc                                                     <  2:3045993/150‑1 (MQ=255)
                                                                                                 aCCGAACTGCTCAAGGGCCTGGACCCGAACAAGGACCAGAGCTACTTCCTGCACGCCGTCGGCGGCCAGGAAATCGCCCGCACACTGTTCCCGGTCGGCGAACTGGAAAAGCCCGAAGTACGCGCCATCGCCGAAAAACACGGCCTGGcc                                                     >  1:3148818/1‑150 (MQ=255)
                                                                                                 aCCGAACTGCTCAAGGGCCTGGACCCGAACAAGGACCAGAGCTACTTCCTGCACGCCGTCGGCGGCCAGGAAATCGCCCGCACACTGTTCCCGGTCGGCGAACTGGAAAAGCCCGAAGTACGCGCCATCGCCGAAAAACACGGCCTGGcc                                                     <  1:2294375/150‑1 (MQ=255)
                                                                                                    gAACTGCTCAAGGGCCTGGACCCGAACAAGGACCAGAGCTACTTCCTGCACGCCGTCGGCGGCCAGGAAATCGCCCGCACACTGTTCCCGGTCGGCGAACTGGAAAAGCCCGAAGTACGCGCCATCGCCGAAAAACACGGCCTGGCCAcc                                                  >  1:3926464/1‑150 (MQ=255)
                                                                                                        tGCTCAAGGGCCTGGACCCGAACAAGGACCAGAGCTACTTCCTGCACGCCGTCGGCGGCCAGGAAATCGCCCGCACACTGTTCCCGGTCGGCGAACTGGAAAAGCCCGAAGTACGCGCCATCGCCGAAAAACACGGCCTGGCCACCGCCa                                              >  1:714148/1‑150 (MQ=255)
                                                                                                         gCTCAAGGGCCTGGACCCGAACAAGGACCAGAGCTACTTCCTGCACGCCGTCGGCGGCCAGGAAATCGCCCGCACACTGTTCCCGGTCGGCGAACTGGAAAAGCCCGAAGTACGCGCCATCGCCGAAAAACACGGCCTGGCCACCGCCaa                                             <  2:1277270/150‑1 (MQ=255)
                                                                                                          cTCAAGGGCCTGGACCCGAACAAGGACCAGAGCTACTTCCTGCACGCCGTCGGCGGCCAGGAAATCGCCCGCACACTGTTCCCGGTCGGCGAACTGGAAAAGCCCGAAGTACGCGCCATCGCCGAAAAACACGGCCTGGCCACCGCCaag                                            >  1:22214/1‑150 (MQ=255)
                                                                                                          cTCAAGGGCCTGGACCCGAACAAGGACCAGAGCTACTTCCTGCACGCCGTCGGCGGCCAGGAAATCGCCCGCACACTGTTCCCGGTCGGCGAACTGGAAAAGCCCGAAGTACGCGCCATCGCCGAAAAACACGGCCTGGCCACCGCCaa                                             >  1:1344963/1‑149 (MQ=255)
                                                                                                              aGGGCCTGGACCCGAACAAGGACCAGAGCTACTTCCTGCACGCCGTCGGCGGCCAGGAAATCGCCCGCACACTGTTCCCGGTCGGCGAACTGGAAAAGCCCGAAGTACGCGCCATCGCCGAAAAACACGGCCTGGCCACCGCCAagaaga                                        <  1:1423073/150‑1 (MQ=255)
                                                                                                              aGGGCCTGGACCCGAACAAGGACCAGAGCTACTTCCTGCACGCCGTCGGCGGCCAGGAAATCGCCCGCACACTGTTCCCGGTCGGCGAACTGGAAAAGCCCGAAGTACGCGCCATCGCCGAAAAACACGGCCTGGCCACCGCCAagaaga                                        <  2:74826/150‑1 (MQ=255)
                                                                                                               gggCCTGGACCCGAACAAGGACCAGAGCTACTTCCTGCACGCCGTCGGCGGCCAGGAAATCGCCCGCACACTGTTCCCGGTCGGc                                                                                                        <  2:2753487/85‑1 (MQ=255)
                                                                                                               gggCCTGGACCCGAACAAGGACCAGAGCTACTTCCTGCACGCCGTCGGCGGCCAGGAAATCGCCCGCACACTGTTCCCGGTCGGCGAACTGGAAAAGCCCGAAGTACGCGCCATCGCCGAAAAACACGGCCTGGCCACCGCCAAgaagaa                                       <  2:4563175/150‑1 (MQ=255)
                                                                                                                ggCCTGGACCCGAACAAGGACCAGAGCTACTTCCTGCACGCCGTCGGCGGCCAGGAAATCGCCCGCACACTGTTCCCGGTCGGCGAACTGGAAAAGCCCGAAGTACGCGCCATCGCCGAAAAACACGGCCTGGCCACCGCCAAGaagaaa                                      <  2:838984/150‑1 (MQ=255)
                                                                                                                ggCCTGGACCCGAACAAGGACCAGAGCTACTTCCTGCACGCCGTCGGCGGCCAGGAAATCGCCCGCACACTGTTCCCGGTCGGCGAACTGGAAAAGCCCGAAGTACGCGCCATCGCCGAAAAACACGGCCTGGCCACCGCCAAGaagaaa                                      <  2:3148818/150‑1 (MQ=255)
                                                                                                                ggCCTGGACCCGAACAAGGACCAGAGCTACTTCCTGCACGCCGTCGGCGGCCAGGAAATCGCCCGCACACTGTTCCCGGTCGGCGAACTGGAAAAGCCCGAAGTACGCGCCATCGCCGAAAAACACGGCCTGGCCACCGCCAAGaagaaa                                      <  1:5530283/150‑1 (MQ=255)
                                                                                                                ggCCTGGACCCGAACAAGGACCAGAGCTACTTCCTGCACGCCGTCGGCGGCCAGGAAATCGCCCGCACACTGTTCCCGGTCGGCGAACTGGAAAAGCCCGAAGTACGCGCCATCGCCGAAAAACACGGCCTGGCCACCGCCAAGaagaaa                                      <  2:3611453/150‑1 (MQ=255)
                                                                                                                 gCCTGGACCCGAACAAGGACCAGAGCTACTTCCTGCACGCCGTCGGCGGCCAGGAAATCGCCCGCACACTGTTCCCGGTCGGCGAACTGGAAAAGCCCGAAGTACGCGCCATCGCCGAAAAACACGGCCTGGCCACCGCCAAGaagaaag                                     <  1:105106/150‑1 (MQ=255)
                                                                                                                    tGGACCCGAACAAGGACCAGAGCTACTTCCTGCACGCCGTCGGCGGCCAGGAAATCGCCCGCACACTGTTCCCGGTCGGCGAACTGGAAAAGCCCGAAGTACGCGCCATCGCCGAAAAACACGGCCTGGCCACCGCCAAGAAGAAAGAtt                                  <  1:4411963/150‑1 (MQ=255)
                                                                                                                     ggACCCGAACAAGGACCAGAGCTACTTCCTGCACGCCGTCGGCGGCCAGGAAATCGCCCGCACACTGTTCCCGGTCGGCGAACTGGAAAAGCCCGAAGTACGCGCCATCGCCGAAAAACACGGCCTGGCCACCGCCAAGAAGAAAGATTc                                 >  2:3851976/1‑150 (MQ=255)
                                                                                                                        cccGAACAAGGACCAGAGCTACTTCCTGCACGCCGTCGGCGGCCAGGAAATCGCCCGCACACTGTTCCCGGTCGGCGAACTGGAAAAGCCCGAAGTACGCGCCATCGCCGAAAAACACGGCCTGGCCACCGCCAAGAAGAAAGATTCCAc                              >  2:2371740/1‑150 (MQ=255)
                                                                                                                         ccGAACAAGGACCAGAGCTACTTCCTGCACGCCGTCGGCGGCCAGGAAATCGCCCGCACACTGTTCCCGGTCGGCGAACTGGAAAAGCCCGAAGTACGCGCCATCGCCGAAAAACACGGCCTGGCCACCGCCAAGAAGAAAGATTCCAcc                             >  1:5087969/1‑150 (MQ=255)
                                                                                                                           gAACAAGGACCAGAGCTACTTCCTGCACGCCGTCGGCGGCCAGGAAATCGCCCGCACACTGTTCCCGGTCGGCGAACTGGAAAAGCCCGAAGTACGCGCCATCGCCGAAAAACACGGCCTGGCCACCGCCAAGAAGAAAGATTCCACCgg                           >  2:968604/1‑150 (MQ=255)
                                                                                                                            aaCAAGGACCAGAGCTACTTCCTGCACGCCGTCGGCGGCCAGGAAATCGCCCGCACACTGTTCCCGGTCGGCGAACTGGAAAAGCCCGAAGTACGCGCCATCGCCGAAAAACACGGCCTGGCCACCGCCAAGAAGAAAGATTCCACCGGc                          >  1:4097618/1‑150 (MQ=255)
                                                                                                                             aCAAGGACCAGAGCTACTTCCTGCACGCCGTCGGCGGCCAGGAAATCGCCCGCACACTGTTCCCGGTCGGCGAACTGGAAAAGCCCGAAGTACGCGCCATCGCCGAAAAACACGGCCTGGCCACCGCCAAGAAGAAAGATTCCACCGGCa                         >  1:546870/1‑150 (MQ=255)
                                                                                                                                aGGACCAGAGCTACTTCCTGCACGCCGTCGGCGGCCAGGAAATCGCCCGCACACTGTTCCCGGTCGGCGAACTTGAAAAGCCCGAAGTACGCGCCATCGCCGAAAAACACGGCCTGGCCACCGCCAAGAAGAAAGATTCCACCGGCATCt                      <  1:5322612/150‑1 (MQ=255)
                                                                                                                                aGGACCAGAGCTACTTCCTGCACGCCGTCGGCGGCCAGGAAATCGCCCGCACACTGTTCCCGGTCGGCGAACTGGAAAAGCCCGAAGTACGCGCCATCGCCGAAAAACACGGCCTGGCCACCGCCAAGAAGAAAGATTCCACCGGCATCt                      >  2:2245289/1‑150 (MQ=255)
                                                                                                                                 ggACCAGAGCTACTTCCTGCACGCCGTCGGCGGCCAGGAAATCGCCCGCACACTGTTCCCGGTCGGCGAACTGGAAAAGCCCGAAGTACGCGCCATCGCCGAAAAACACGGCCTGGCCACCGCCAAGAAGAAAGATTCCACCGGCATCTg                     >  2:3230948/1‑150 (MQ=255)
                                                                                                                                 cgACCAGAGCTACTTCCTGCACGCCGTCGGCGGCCAGGAAATCGCCCGCACACTGTTCCCGGTCGGCGAACTGGAAAAGCCCGAAGTACGCGCCATCGCCGAAAAACACGGCCTGGCCACCGCCAAGAAGAAAGATTCCACCGGCATCTg                     <  2:4254613/149‑1 (MQ=255)
                                                                                                                                    ccAGAGCTACTTCCTGCACGCCGTCGGCGGCCAGGAAATCGCCCGCACACTGTTCCCGGTCGGCGAACTGGAAAAGCCCGAAGTACGCGCCATCGCCGAAAAACACGGCCTGGCCACCGCCAAGAAGAAAGATTCCACCGGCATCTGCtt                  <  2:2800841/150‑1 (MQ=255)
                                                                                                                                    ccAGAGCTACTTCCTGCACGCCGTCGGCGGCCAGGAAATCGCCCGCACACTGTTCCCGGTCGGCGAACTGGAAAAGCCCGAAGTACGCGCCATCGCCGAAAAACACGGCCTGGCCACCGCCAAGAAGAAAGATTCCACCGGCATCTGCtt                  <  1:1525053/150‑1 (MQ=255)
                                                                                                                                        agCTACTTCCTGCACGCCGTCGGCGGCCAGGAAATCGCCCGCACACTGTTCCCGGTCGGCGAACTGGAAAAGCCCGAAGTACGCGCCATCGCCGAAAAACACGGCCTGGCCACCGCCAAGAAGAAAGATTCCACCGGCATCTGCTTCAtt              <  1:1154764/150‑1 (MQ=255)
                                                                                                                                           tACTTCCTGCACGCCGTCGGCGGCCAGGAAATCGCCCGCACACTGTTCCCGGTCGGCGAACTGGAAAAGCCCGAAGTACGCGCCATCGCCGAAAAACACGGCCTGGCCACCGCCAAGAAGAAAGATTCCACCGGCATCTGCTTCATTGGc           <  1:2071143/150‑1 (MQ=255)
                                                                                                                                             tttCCTGGCCGCCGACGGCGTGCAGGAAATCGCCCGCACACTGTTCCCGGTCGGCGAACTGGAAATGCCCGAAGTACGCGCCATCGCCGAAAAACACGGCCTGGCCACCGCCAAGAAGAAAGATTCCACCGGCATCTGCTTCATTGGCGa         >  2:2066983/2‑150 (MQ=255)
                                                                                                                                                  tGCACGCCGTCGGCGGCCAGGAAATCGCCCGCACACTGTTCCCGGTCGGCGAACTGGAAAAGCCCGAAGTACGCGCCATCGCCGAAAAACACGGCCTGGCCACCGCCAAGAAGAAAGATTCCACCGGCATCTGCTTCATTGGCGAGcgcc    >  2:1797976/1‑150 (MQ=255)
                                                                                                                                                  tGCACGCCGTCGGCGGCCAGGAAATCGCCCGCACACTGTGCCCGGTCGGCGAACTGGAAAAGCCCGAAGTACGCGCCATCGCCGAAAAACACGGCCTGGCCACCGCCAAGAAGAAAGATTCCACCGGCATCTGCTTCATTGGCGAGcgcc    <  1:3292590/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
AATCAAGTTCAAGGCGTTCCTCGACTACGCCCTGTCGCTGGGTGCCGACCTGATCGCCACCGGCCACTATGTGCGCCGCCGCGATACTGGCGCGCTCACCGAACTGCTCAAGGGCCTGGACCCGAACAAGGACCAGAGCTACTTCCTGCACGCCGTCGGCGGCCAGGAAATCGCCCGCACACTGTTCCCGGTCGGCGAACTGGAAAAGCCCGAAGTACGCGCCATCGCCGAAAAACACGGCCTGGCCACCGCCAAGAAGAAAGATTCCACCGGCATCTGCTTCATTGGCGAGCGCCGC  >  NC_002947/4530395‑4530692

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