Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 403,270 T→A 5.8% intergenic (‑70/+99) 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,2700TA 5.8% 378.1 / 10.5 121intergenic (‑70/+99)msrA/PP_0337peptide methionine sulfoxide reductase MsrA/sensory box protein
Reads supporting (aligned to +/- strand):  ref base T (47/66);  new base A (2/5);  total (50/71)
Fisher's exact test for biased strand distribution p-value = 6.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.07e-01

Only 100 of 124 total aligned reads displayed.

GGCGTGTCGCGGCCAGGCAGGGCCTGTTCCGCAGTGGGCATGACGTTTTTGTTCACCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCCTTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCAGGCCACCGGGCCACGCGGATAGCGCTTCAGCCGCTCGGCCAGTTCGCG  >  NC_002947/403121‑403419
                                                                                                                                                      |                                                                                                                                                     
ggCGTGTCGCGGCCAGGCAGGGCCTGTTCCGCAGTGGGCATGACGTTTTTGTTCACCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCCTTCGCGGGCATGCCCGCTCCCACAGGt                                                                                                                                                       <  1:1163452/150‑1 (MQ=255)
   gtgtCGCGGCCAGGCAGGGCCTGTTCCGCAGTGGGCATGACGTTTTTGTTCACCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCCTTCGCGGGCATGCCCGCTCCCACAGGTACa                                                                                                                                                    <  1:4520372/150‑1 (MQ=255)
        gcgGCCAGGCAGGGCCTGTTCCGCAGTGGGCATGACGTTTTTGTTCACCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCCTTCGCGGGCATGCCCGCTCCCACAGGTACAGcaca                                                                                                                                               <  2:1563861/150‑1 (MQ=255)
        gcgGCCAGGCAGGGCCTGTTCCGCAGTGGGCATGACGTTTTTGTTCACCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCCTTCGCGGGCATGCCCGCTCCCACAGGTACAGcaca                                                                                                                                               <  1:2970266/150‑1 (MQ=255)
          ggCCAGGCAGGGCCTGTTCCGCAGTGGGCATGACGTTTTTGTTCACCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCCTTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAgg                                                                                                                                             >  1:4615643/1‑150 (MQ=255)
             caggcaggGCCTGTTCCGCAGTGGGCCTGACGTTGTTGTACACCAGGATCTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCCTTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGttt                                                                                                                                          <  1:3159294/150‑1 (MQ=255)
             caggcaggGCCTGTTCCGCAGTGGGCATGACGTTTTTGTTCACCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCCTTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGttt                                                                                                                                          <  2:2135913/150‑1 (MQ=255)
              aggcaggGCCTGTTCCGCAGTGGGCATGACGTTTTTGTTCACCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCTTTCGCGGGCACGCCCGCTCCCACAGGAACACcacaaccttt                                                                                                                                         <  2:4905068/150‑7 (MQ=38)
              aggcaggGCCTGTTCCGCAGTGGGCATGACGTTTTTGTTCACCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCTTTCGCGGGCACGCCCGCTCCCACAGGAACACcacaaccttt                                                                                                                                         <  2:4905012/150‑7 (MQ=38)
              aggcaggGCCTGTTCCGCAGTGGGCATGACGTTTTTGTTCACCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCCTTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTc                                                                                                                                         <  1:697917/150‑1 (MQ=255)
              aggcaggGCCTGTTCCGCAGTGGGCATGACGTTTTTGTTCACCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCCTTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTc                                                                                                                                         >  2:2971881/1‑150 (MQ=255)
                gcaggGCCTGTTCCGCAGTGGGCATGACGTTTTTGTTCACCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCCTTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCaa                                                                                                                                       >  1:2052989/1‑150 (MQ=255)
                 caggGCCTGTTCCGCAGTGGGCATGACGTTTTTGTTCACCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCCTTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCaaa                                                                                                                                      <  1:1425456/150‑1 (MQ=255)
                     gCCTGTTCCGCAGTGGGCATGACGTTTTTGTTCACCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCCTTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGtt                                                                                                                                  <  1:3305422/150‑1 (MQ=255)
                     gCCTGTTCCGCAGTGGGCATGACGTTTTTGTTCACCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCCTTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGtt                                                                                                                                  <  1:251247/150‑1 (MQ=255)
                     gCCTGTTCCGCAGTGGGCATGACGTTTTTGTTCACCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCCTTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGtn                                                                                                                                  <  1:11321/150‑2 (MQ=255)
                       cTGTTCCGCAGTGGGCATGACGTTTTTGTTCACCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCCTTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTtgt                                                                                                                                <  2:5011216/150‑1 (MQ=255)
                       cTGTTCCGCAGTGGGCATGACGTTTTTGTTCACCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCCTTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTtgt                                                                                                                                >  2:567189/1‑150 (MQ=255)
                        tGTTCCGCAGTGGGCATGACGTTTTTGTTCACCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCCTTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTtgtg                                                                                                                               >  2:2459561/1‑150 (MQ=255)
                        tGTTCCGCAGTGGGCATGACGTTTTTGTTCACCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCCTTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTtgtg                                                                                                                               <  2:68541/150‑1 (MQ=255)
                        tGTTCCGCAGTGGGCATGACGTTTTTGTTCACCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCCTTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTtgtg                                                                                                                               <  2:4033545/150‑1 (MQ=255)
                        tGTTCCGCAGTGGGCATGACGTTTTTGTTCACCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCCTTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTtgtg                                                                                                                               >  1:5105872/1‑150 (MQ=255)
                        tGTTCCGCAGTGGGCATGACGTTTTTGTTCACCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCCTTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTtgtg                                                                                                                               <  2:3476445/150‑1 (MQ=255)
                        tGTTCCGCAGTGGGCATGACGTTTTTGTTCACCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCCTTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTtgtg                                                                                                                               <  1:2020264/150‑1 (MQ=255)
                        tGTTCCGCAGTGGGCATGACGTTTTTGTTCACCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGGGGGGGAGGGGTTTGGGCCTTCGCGGGCATGCCCGCTCCCACAGGCACAGCACAGGTTTCAAAGGTtgtg                                                                                                                               >  1:1798784/1‑150 (MQ=255)
                        gggtCCGCAGTGGGCATGACGTTTTTGTTCACCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCCTTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTtgtg                                                                                                                               <  2:3408716/147‑1 (MQ=255)
                         gTTCCGCAGTGGGCATGACGTTTTTGTTCACCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCCTTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTgtgg                                                                                                                              <  1:3878386/150‑1 (MQ=255)
                          ttCCGCAGTGGTCATGTCGTTTTTGTTCACCAGGATTTCCGTACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCCTTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTgtggt                                                                                                                             <  1:4179469/150‑1 (MQ=255)
                          ttCCGCAGTGGGCATGACGTTTTTGTTCACCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCCTTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTgtggt                                                                                                                             <  1:4179400/150‑1 (MQ=255)
                          gtCCGCAGTGGGCATGACGTTTTTGTTCACCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCCTTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTgtggt                                                                                                                             <  2:209737/149‑1 (MQ=255)
                                  tGGGCATGACGTTTTTGTTCACCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCCTTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCtgt                                                                                                                     <  2:3596341/150‑1 (MQ=255)
                                   gggCATGACGTTTTTGTTCACCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCCTTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCtgtg                                                                                                                    <  2:1135515/150‑1 (MQ=255)
                                   gggCATGACGTTTTTGTTCACCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCCTTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCtgtg                                                                                                                    <  2:1135357/150‑1 (MQ=255)
                                    ggCATGTCGTTTTTGTTCTCCAGGATTTCCGATCGCAGGACCATGGCCGTTCCTCTGTGTTGTTTCTGTGGGGTAGGTGTTTGGGCCTACGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAACGGTTGTGGTGTTCCTGTgg                                                                                                                   <  2:5107918/150‑1 (MQ=38)
                                    ggCATGACGTTTTTGTTCACCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCCCTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTgg                                                                                                                   <  1:3877577/150‑1 (MQ=255)
                                     gCATGACGTTTTTGTTCACCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCCTTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTAAAAGGTTGTGGTGTTACTGTggg                                                                                                                  >  1:1805765/1‑150 (MQ=255)
                                       aTGACGTTTTTGTTCACCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGCCCTTTCGCGGGCACGCCCGCTCCCACAGGAACACCACAACCTTTGAAACCTGTGCTGTACCTGTGGGa                                                                                                                 <  1:3461064/150‑1 (MQ=17)
                                        tGACGTTTTTGTTCACCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCTGTGGGGTAGGGGTTTGGGCCTTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGc                                                                                                               <  1:2844026/150‑1 (MQ=255)
                                        tGACGTTTTTGTTCACCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCCTTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGc                                                                                                               >  1:3491600/1‑150 (MQ=255)
                                          aCGTTTTTGTTCACCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCCTTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCgg                                                                                                             <  2:3807465/150‑1 (MQ=255)
                                          aCGTTTTTGTTCACCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCCTTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCgg                                                                                                             >  2:2091467/1‑150 (MQ=255)
                                            gTTTTTGTTCACCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCCTTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGc                                                                                                           <  1:2024029/150‑1 (MQ=255)
                                                   ttCACCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCCTTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCcg                                                                                                    >  1:412686/1‑150 (MQ=255)
                                                      aCCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCCTTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGa                                                                                                 >  1:4362916/1‑150 (MQ=255)
                                                       ccAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCCTTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGaa                                                                                                >  2:465703/1‑150 (MQ=255)
                                                       ccAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCCTTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGaa                                                                                                >  1:753027/1‑150 (MQ=255)
                                                        cAGGCTTTCCGACCGCCGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTGGGGCCTTCGCGGGCATGCCCGCTCCCACCGGTACAGCACAGGTTTCAACGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGaaa                                                                                               <  1:2415936/150‑1 (MQ=38)
                                                        cAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCCTTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTc                                                                                                                         >  2:3775634/1‑124 (MQ=255)
                                                        cAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCCTTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTc                                                                                                                         <  1:3775634/124‑1 (MQ=255)
                                                        cAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCCTTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGaaa                                                                                               >  2:4401213/1‑150 (MQ=255)
                                                        cAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCCTTCGCGGGCATGCCCGCACCCACAGGTACAGCACAGGTTTCAAAGGGTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGaaa                                                                                               >  2:743254/1‑150 (MQ=255)
                                                        cAGGATGTCCGAACGCAGGACCCTGGCCGTCCCTCTGTGTTTCGTCAGGGGGGTCGGGGGTGGGGCCTTCGCGGGCATGCCCGCTCCCACCGTTCCCGCACCGGTTTCAAAGGTTGTGGTGGGCCTGTGGGCGCGGGCGTGCCCGCGaac                                                                                               <  2:3014020/150‑2 (MQ=25)
                                                        cAGGAGTTCCGAACGCAGGACCATGGCCGGTCCTCTGGGTTTATTCAGTGGGGTAGGGGTTGGGGCCTTCGCGGGCAGGCCCGCTCCCACAGGTACAGCACAGGTTGCAAAGGTTGTGGTGTTCCTGGGGGAGCGGGCGTGCCCGCGaaa                                                                                               <  2:3013684/150‑1 (MQ=38)
                                                         aGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCCTTCGCGGGCACGCCCGCTCCCACAGGAACACCACAACCTTTGAAACCTGTGCTGTACCTGTGGGAGCGGGCATGCCCGCGaagg                                                                                              <  1:4080404/150‑3 (MQ=14)
                                                            aTTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCCTTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAATGGTTTTGGTTTTCCTGTGGGAGCGGGCGTGCCCGCGAAAgggg                                                                                           >  1:2494162/1‑149 (MQ=255)
                                                                  gAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCCTTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAAt                                                                                     >  2:1672416/1‑150 (MQ=255)
                                                                  gAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCCTTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCGGTGGGAGCGGGCGGGCCCGCGACAGGGCCCGAAt                                                                                     >  1:103848/1‑150 (MQ=255)
                                                                  gAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCCTTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTGTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAAt                                                                                     >  1:2344206/1‑150 (MQ=255)
                                                                     cGCAGGACCATGGCCGTGCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCCTTCGCGGGCATGCCCGCTCCCACCGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGc                                                                                  <  2:53868/150‑1 (MQ=255)
                                                                        aGGACCTTGGCCGTTCCTCTGTGTTTATGCAGTGGGGTAGGGGTTTGGGCCTTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTc                                                                               <  2:2472693/150‑1 (MQ=255)
                                                                        aGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCCTTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCATGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATggagtg                                                                               >  2:449000/1‑146 (MQ=255)
                                                                            ccATGGCCGTTCCTCTGTTTTTATTCAGTGGGGTAGGGGTTTGGGCCTTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCtata                                                                           <  1:5044095/150‑1 (MQ=255)
                                                                            ccATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTTGGGGTTTGGGCCTTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCtata                                                                           <  2:2762936/150‑1 (MQ=255)
                                                                              aTGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCCTTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGt                                                                         >  1:1598853/1‑150 (MQ=255)
                                                                                ggCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCCTTCGCGGGCAAGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGc                                                                       >  2:2757063/1‑150 (MQ=255)
                                                                                               tttATTCAGTGGGGTAGGGGTTTGGGCCTTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAg                                                        >  1:3344783/1‑150 (MQ=255)
                                                                                               tttATTCAGTGGGGTAGGGGTTTGGGCCTTCGCGGGCATGCCCGCACCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTACTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAg                                                        >  1:4754188/1‑150 (MQ=255)
                                                                                                     ccgTGGGGTTGGGGTTTGGCCCTTCGGGGGCATGCCCGCTCCCACAGGTACAGCACCGTTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGGGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCa                                                  <  1:2569745/148‑1 (MQ=37)
                                                                                                     ccgTGGGGTTGGGGTTGGGGCCTTCGCGGGCTGGCCCGCTCCCCCCGGTACAGCCCAGGTTTCAATGGTTGTGGTGTTCCTGTGGGCGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCa                                                  <  1:674669/148‑1 (MQ=255)
                                                                                                     cAGTGGGGTAGGGGTTTGGGCCTTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCa                                                  <  2:3204802/150‑1 (MQ=255)
                                                                                                            gTAGGGGTTTGGGCCTTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGGTGGGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCAGGCCAcc                                           >  2:3968977/1‑150 (MQ=255)
                                                                                                               ggggTTTGGGCCTTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCAGGCCACCggg                                        >  2:4505098/1‑150 (MQ=255)
                                                                                                                 ggTTTGGGCCTTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCAGGCCACCGGGcc                                      <  2:4090676/150‑1 (MQ=255)
                                                                                                                   ttggggCCTGCGCGGGCAGGCCCGCTCCCACCGGTACAGCACCGGTTCCAACGGTTGTGGGGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCATCTGTCAGGCCACCGGGCCAc                                    <  1:291526/147‑1 (MQ=38)
                                                                                                                     tGGGCCGTCGCGGGCAGGCCCGCTCCCCCAGGTACAGCACAGGGTTCACCGGTGGTGGGGTTCCTGGGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCAGGCCACCGGGCCAcgc                                  <  2:1113550/150‑1 (MQ=38)
                                                                                                                        gCCTTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCAGGCCACCGGGCCACGCggt                               >  2:2533472/1‑149 (MQ=255)
                                                                                                                         ccTTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAATGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCAGGCCACCGGGCCACGCGGAt                              >  1:4336316/1‑150 (MQ=255)
                                                                                                                          cTTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAACGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCAGGCCACCGGGCCACGCGGATa                             <  1:2459561/150‑1 (MQ=255)
                                                                                                                          cTTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCAGGCCACCGGGCCACGCGGATa                             <  1:3306605/150‑1 (MQ=255)
                                                                                                                          cTTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCAGGCCACCGGGCCACGCGGATa                             <  1:2971881/150‑1 (MQ=255)
                                                                                                                            ttCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCtgtg                                                                                                                    <  2:614150/62‑1 (MQ=255)
                                                                                                                               gcgGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCAGGCCACCGGGCCACGCGGATAgcgc                         <  2:2052989/150‑1 (MQ=255)
                                                                                                                                  ggccggCCCGCTCCCACCGGTCCAGCACAGGTTTCAATGGTTGTGGTGTTCCGGTGGGTGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCAGGCCACCGGGCCACGCGGATAGCGCTTc                      <  2:1312502/145‑1 (MQ=255)
                                                                                                                                  ggCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTGGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCAGGCCACCGGGCCACGCGGATAGCGCTTc                      <  2:4615834/150‑1 (MQ=255)
                                                                                                                                  ggCATGCCCGCGCCCACAGGTACAGCACAGGTTCCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCAGGCCACCGGGCCACGCGGATAGCGCTTc                      <  1:3878210/150‑1 (MQ=255)
                                                                                                                                         ccGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCAGGCCACCGGGCCACGCGGATAGCGCTTCAGCCGCt               >  1:2472381/1‑150 (MQ=255)
                                                                                                                                         ccGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCAGGCCACCGGGCCACGCGGATAGCGCTTCAGCCGCt               >  1:4821872/1‑150 (MQ=255)
                                                                                                                                           gCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCAGGCCACCGGGCCACGCGGATAGCGCTTCAGCCGCTCg             <  2:3772175/150‑1 (MQ=255)
                                                                                                                                           gCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCAGGCCACCGGGCCACGCGGATAGCGCTTCAGCCGCTCg             >  2:673306/1‑150 (MQ=255)
                                                                                                                                           gCTCCCACAGGAACACCACAACCTTTGAAACCTGTGCTGTACCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCAGGCCACCGGGCCACGCGGATAGCGCTTCAGCCGCTCg             >  1:1554055/1‑150 (MQ=18)
                                                                                                                                            cTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCAGGCCACCGGGCCACGCGGATAGCGCTTCAGCCGCTCgg            >  1:1448855/1‑150 (MQ=255)
                                                                                                                                               ccACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCAGGCCACCGGGCCACGCGGATAGCGCTTCAGCCGCTCGGCCa         <  2:412686/150‑1 (MQ=255)
                                                                                                                                                 acaGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTCTAGTGCCCGAGGATTAGGCAGCTGTCAGGCCACCGGGCCACGCGGATAGCGCTTCAGCCGCTCGGCCAGt       >  2:2699521/1‑150 (MQ=255)
                                                                                                                                                 acaGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCAGGCCACCGGTCCACGCGGATAGCGCTTCAGCCGCTCGGCCAGt       >  2:3845838/1‑150 (MQ=255)
                                                                                                                                                 acaGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCAGGCCACCGGGCCACGCGGATAGCGCTTCAGCCGCTCGGCCAGt       >  2:2153896/1‑150 (MQ=255)
                                                                                                                                                 acaGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCAGGCCACCGGGCCACGCGGATAGCGCTTCAGCCGCTCGGCCAGt       >  2:3444225/1‑150 (MQ=255)
                                                                                                                                                   aGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCAGGCCACCGGGCCACGCGGATAGCGCTTCAGCCGCTCGGCCAGTTc     >  2:4551750/1‑150 (MQ=255)
                                                                                                                                                     gTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGACCGAGGATTACGCAGCTGTCAGGCCACCGGGCCACGCGGATAGCGCTTCAGCCGCTCGGCCAGTTcgc   >  1:712612/1‑150 (MQ=255)
                                                                                                                                                      tACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCAGGCCACCGGGCCACGCGGATAGCGCTTCAGCCGCTCGGCCAGTTcgcg  >  1:2070536/1‑150 (MQ=255)
                                                                                                                                                      tACAGAACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGACAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCAGGCCACCGGGCCACGCGGATAGCGCTTCAGCCGCTCGGCCAGTTcgcg  >  2:3273960/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
GGCGTGTCGCGGCCAGGCAGGGCCTGTTCCGCAGTGGGCATGACGTTTTTGTTCACCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCCTTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCAGGCCACCGGGCCACGCGGATAGCGCTTCAGCCGCTCGGCCAGTTCGCG  >  NC_002947/403121‑403419

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: