Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 403,280 G→T 6.2% intergenic (‑80/+89) 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,2800GT 6.2% 773.6 / 23.5 246intergenic (‑80/+89)msrA/PP_0337peptide methionine sulfoxide reductase MsrA/sensory box protein
Reads supporting (aligned to +/- strand):  ref base G (91/129);  new base T (12/3);  total (103/143)
Fisher's exact test for biased strand distribution p-value = 5.49e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 251 total aligned reads displayed.

GGCCAGGCAGGGCCTGTTCCGCAGTGGGCATGACGTTTTTGTTCACCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCCTTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCAGGCCACCGGGCCACGCGGATAGCGCTTCAGCCGCTCGGCCAGTTCGCGCCCGGGAATC  >  NC_002947/403131‑403429
                                                                                                                                                     |                                                                                                                                                     
  ccAGGCAGGGCCTGTTCCGCAGTGGGCATGACGTTTTTGTTCACCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCCTTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGtt                                                                                                                                                     >  1:559887/1‑150 (MQ=255)
   caggcaggGCCTGTTCCGCAGTGGGCATGACTTTTTTGTTCACCAGGATTTCCGAACGCAGGACCATGGCGGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCCTTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGttt                                                                                                                                                    <  2:3972350/150‑1 (MQ=255)
   caggcaggGCCTGTTCCGCAGTGGGCATGACGTTTTTGTTCACCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCCTTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGttt                                                                                                                                                    <  1:5929063/150‑1 (MQ=255)
    aggcaggGCCTGTTCCGCAGTGGGCATGACGTTTTTGTTCACCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCCTTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTc                                                                                                                                                   >  2:5906689/1‑150 (MQ=255)
      gcaggGCCTGTTCCGCAGTGGGCATGACGTTTTTGTTCACCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCCTTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCaa                                                                                                                                                 >  1:3730291/1‑150 (MQ=255)
      gcaggGCCTGTTCCGCAGTGGGCATGACGTTTTTGTTCACCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCCTTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCaa                                                                                                                                                 >  1:3548214/1‑150 (MQ=255)
      gcaggGCCTGTTCCGCAGTGGGCATGACGTTTTTGTTCACCAGGATTTCCGAACGCAGGACATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCCTTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCaaa                                                                                                                                                <  2:2214834/150‑1 (MQ=255)
       caggGCCTGTTCCGCAGTGGGCATGACGTTTTTGTTCACCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCCTTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCaaa                                                                                                                                                <  2:5238321/150‑1 (MQ=255)
       agccaagttctaagcaacaccctgctctgcttcacctgaagcgggaaacgccgaccccatccggcgtttcccgcttctttttgtctgtgtgaaatgcgccggcctcttcgcgggctTGCCCGCTCCCACAGGTACAGCACAGTTCTcaag                                                                                                                                                >  2:1761492‑M1/117‑149 (MQ=255)
        aggGCCTGTTCCGCAGTGGGCATGACGTTTTTGTTCACCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGCTTATTCAGTGGGGTAGGGGTTTGGGCCTTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAg                                                                                                                                               <  2:4384141/150‑1 (MQ=255)
        aggGCCTGTTCCGCAGTGGGAATGACGTTTTTGTTCACCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGCTTATTCAGTGGGGTAGGGGTTTGGGCCTTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAg                                                                                                                                               <  2:6569512/150‑1 (MQ=255)
         gggCCTGTTCCGCAGTGGGCATGACGTTTTTGTTCACCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCCTTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAgg                                                                                                                                              >  2:3264121/1‑150 (MQ=255)
          ggCCTGTTCCGCAGTGGGCATGACGTTTTTGTTCACCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCCTTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGt                                                                                                                                             >  1:2247726/1‑150 (MQ=255)
              tGTTCCGCAGTGTGCATGACGTTTTTTTTCACCAGGATTTCCGAACGCTGGACCATGGCCGTTCCTGTGTGTTTTTTCAGTGTGGTAGTGGTTTGGGCCTTCGCGGGCATGCCCGCTCCCTCAGGTACAGCACAGGTTTCAAAGGTtgtg                                                                                                                                         <  1:4095101/150‑1 (MQ=255)
              tGTTCCGCAGTGGGCATGACGTTTTTGTTCACCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCCTTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTtgtg                                                                                                                                         <  1:1116043/150‑1 (MQ=255)
              tGTTCCGCAGTGGGCATGACGTTTTTGTTCACCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCCTTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTtgtg                                                                                                                                         <  1:2514879/150‑1 (MQ=255)
              tGTTCCGCAGTGGGCATGACGTTTTTGTTCACCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCCTTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTtgtg                                                                                                                                         >  2:6384829/1‑150 (MQ=255)
               gTTCCGCAGTGGGCATGACGTTTTTGTTCACCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCCTTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTgtgg                                                                                                                                        <  1:5989291/150‑1 (MQ=255)
                ttCCGCAGTGGGCATGACGTTTTTGTTCACCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCCTTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTgtggt                                                                                                                                       <  2:1581155/150‑1 (MQ=255)
                ttCCGCAGTGGGCATGACGTTTTTGTTCACCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGGTTAGTCAGTGGGGTAGGGGTTTGGGCCTTCGCGGGCATGCCCGCTCCCACAGGGCCAGCACAGGTTTCAAAGGGTgtggt                                                                                                                                       >  2:3952083/1‑150 (MQ=255)
                 tCCGCAGTGGGCATGACGTTTTTGTTCACCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCCTTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTgtggtg                                                                                                                                      <  1:1572188/150‑1 (MQ=255)
                   cGCAGTGGGCATGACGTTTTTGTTCACCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCCTTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGtt                                                                                                                                    >  2:5420835/1‑150 (MQ=255)
                    gCAGTGGGCATGACGTTTTTGTTCACCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCCTTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTc                                                                                                                                   <  1:3264121/150‑1 (MQ=255)
                    gCAGTGGGCATGACGTTTTTGTTCACCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCCTTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTc                                                                                                                                   <  2:1674983/150‑1 (MQ=255)
                    gCAGTGGGCATGACGTTTTTGTTCACCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCCTTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTc                                                                                                                                   <  2:1940708/150‑1 (MQ=255)
                     cAGTGGGCATGACGTTTTTGTTCACCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCCTTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTcc                                                                                                                                  <  2:102180/150‑1 (MQ=255)
                        tGGGCATGACGTTTTTGTTCACCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCCTTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCtgt                                                                                                                               <  2:2216748/150‑1 (MQ=255)
                        tGGGCATGACGTTTTTGTTCACCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCCTTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCtgt                                                                                                                               >  1:3386144/1‑150 (MQ=255)
                         gggCATGACGTTTTTGTTCACCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCCTTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCtgtg                                                                                                                              <  2:2003394/150‑1 (MQ=255)
                          ggCATGACGTTTTTGTTCACCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCCTTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTgg                                                                                                                             <  1:4554211/150‑1 (MQ=255)
                          ggCATGACGTTTTTGTTCACCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCCTTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTgg                                                                                                                             <  2:6843890/150‑1 (MQ=255)
                            cATGACGTTTTTGTTCACCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCCTTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGa                                                                                                                           >  1:4899568/1‑150 (MQ=255)
                            cATGACGTTTTTGTTCACCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCCTTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGa                                                                                                                           <  2:313620/150‑1 (MQ=255)
                            cATGACGTTTTTGTTCACCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCCTTCGCGGGCATGCCCGCTCCCACAGGAACACCACAACCTTTGAAACCTGTGCTGTACCTGTGGGa                                                                                                                           <  2:6761685/150‑1 (MQ=21)
                             tgctctgcttcacctgatgcgggaaacgccgaccccatccggcgtttcccgcttctttttgtctgtgtgaaatgcgccggccccttcgcgggctTGCCCGCTCCCACAGGTACAGCACAGTTCTCAAGGTCTGTGGTGATACTGTGGGAg                                                                                                                          >  2:1188202‑M1/95‑150 (MQ=255)
                             ggctctgcttcacctgaagcgggaaacgccgaccccatccggcgtttcccgcttctttttgtctgtgtgaaatgcgccggccccttcgcgggctTGCCCGCTCCCACAGGTACAGCACAGTTCTCAAGGTCTGTGGTGATACTGTGGGAg                                                                                                                          >  1:5322409‑M1/95‑150 (MQ=255)
                              tGACGTTTTTGTTCACCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCCTTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGc                                                                                                                         <  2:1480348/150‑1 (MQ=255)
                              tGACGTTTTTGTTCACCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCCTTCGCGGGCACGCCCGCTCCCACAGGAACACCACAACCTTTGAAACCTGTGCTGTACCTGTGGGAGc                                                                                                                         <  1:1888749/150‑1 (MQ=18)
                                tctgcttcacctgaagcgggaaccgccgaccccatccggcgtttcccgcttctttttgtctgtgtgaaatgcgccggccccttcgcgggctTGCCCGCTCCCACAGGTACAGCACAGTTCTCAAGGTCTGTGGTGATACTGTGGGAGCgg                                                                                                                       >  2:4810769‑M1/92‑150 (MQ=255)
                                tctgcttcacctgaagcgggaaacgccgaccccatccggcgtttcccgcttctttttgtctgtgtgaaatgcgacggcctcttcgcgggctTGCCCGCTCCCACAGGTACAGCACAGTTCTCAAGGTCTGTGGTGATACTGTGGGAGCgg                                                                                                                       >  1:795198‑M1/92‑150 (MQ=255)
                                  gTTTTTGTTCACCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCCTTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGc                                                                                                                     >  1:3502762/1‑150 (MQ=255)
                                  gTTTTTGTTCACCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCCTTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGc                                                                                                                     <  1:2997830/150‑1 (MQ=255)
                                  gTTTTTGTTCACCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCCTTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGc                                                                                                                     <  2:6723680/150‑1 (MQ=255)
                                  gTTTTTGTTCACCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCCTTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGc                                                                                                                     >  1:90343/1‑150 (MQ=255)
                                    ttttGTTCACCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCCTTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGt                                                                                                                   <  1:5124048/150‑1 (MQ=255)
                                    ttttGTTCACCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCCTTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGt                                                                                                                   <  1:4726146/150‑1 (MQ=255)
                                    ttttGTTCACCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTAGTCAGGGGGGTAGGGGTTTGGGCCTTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGt                                                                                                                   <  2:6514968/150‑1 (MQ=255)
                                     cctGTTCACCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCCTTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTg                                                                                                                  >  2:5100471/3‑150 (MQ=255)
                                          tCACCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCCTTCGCGGGCATGCCCGCTCCCACAGGTACAGCACCGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCcgc                                                                                                             >  2:5884376/1‑150 (MQ=255)
                                            aCCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCCTTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGa                                                                                                           >  1:2797317/1‑150 (MQ=255)
                                             ccAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCCTTCGCGGGCATGCCCGCTCCCACAGGTACAGCACCGGTTTCAAAGGGGGTGGTGTTCCAGTGGGAGCGGGCGTGCCCGCGaa                                                                                                          >  1:4446428/1‑150 (MQ=255)
                                             ccAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCCTTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGaa                                                                                                          >  2:4992354/1‑150 (MQ=255)
                                             ccAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCCTTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGaa                                                                                                          >  2:3042822/1‑150 (MQ=255)
                                              cAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCCTTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGaaa                                                                                                         >  2:5692836/1‑150 (MQ=255)
                                              cAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCCTTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGaaa                                                                                                         >  2:2655641/1‑150 (MQ=255)
                                              cAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCCTTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGaaa                                                                                                         >  2:5264065/1‑150 (MQ=255)
                                              cAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCCTTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGaaa                                                                                                         >  2:3795660/1‑150 (MQ=255)
                                              cAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCCTTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGGGGGAGCGGGCGGGCCCGCGaaa                                                                                                         >  1:1231175/1‑150 (MQ=255)
                                              cAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCCTTCGCGGGCAAGCCCCCACCAACAGGTACAGCACAGGTTTCAATGGGTGTGGTGTTACTGTGGGTGCGGGCGTGCCCGCGaaa                                                                                                         >  1:805006/1‑150 (MQ=255)
                                                   tttCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCCTTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGcc                                                                                                    >  1:3951156/1‑150 (MQ=255)
                                                       cGAACGCAGGACCTTGGCCGTTCCTGTGTGTTTATTCTGTGGGGTTGGGGTTTGGGCCTTGGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGGGTTCCTGTGGGTGCGGGCGTGCCCGCGAAAGGGCCCGaa                                                                                                <  1:5567868/150‑1 (MQ=38)
                                                        gAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGGGGGGTAGGGGTTTGGGCCTTCGCGGGCATGCCCGCTCCCACAGGGACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGCAt                                                                                               >  1:4626787/1‑150 (MQ=255)
                                                                   cATTGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCCTTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCATGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAg                                                                                    >  2:2352599/1‑150 (MQ=255)
                                                                   cATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCCTTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAg                                                                                    >  1:5637874/1‑150 (MQ=255)
                                                                    aTGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCCTTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCtata                                                                                     <  2:5536586/148‑1 (MQ=255)
                                                                                 tgtgTTTATTCAGTGGGGTAGGGGTTTGGGCCTTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTAc                                                                      >  2:308832/1‑150 (MQ=255)
                                                                                    ttttATTCTGTGGGGTCGGGGTTTGGGCCTTGGGGGGCTGGCCCGCTCCCACTGGTACAGCACTGTTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCa                                                                   <  1:5243727/149‑1 (MQ=38)
                                                                                    gtTTATTCAGTGGGGTAGGGGTTTGGGCCTTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCa                                                                   >  1:2620981/1‑150 (MQ=255)
                                                                                     tttATTCAGTGGGGTAGGGGTTTGGGCCTTCGCGGGCAAGCCCGCACCCACAGGTACAGCACAGGTTTCAAAGGGTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCCAAAGGGCCCCAATGGCTTCGATAGTGCCCGAGGATTACGCAg                                                                  >  1:198534/1‑150 (MQ=38)
                                                                                          tCAGTGGGGTAGGGGTTTGGGCCTTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTGGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTc                                                             <  2:3555472/150‑1 (MQ=255)
                                                                                           cAGTGGGGTAGGGGTTTGGGCCTTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCACAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCa                                                            >  2:1433746/1‑150 (MQ=255)
                                                                                           cAGTGGGGTAGGGGTTTGGGCCTCCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAACGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCa                                                            <  1:6184448/150‑1 (MQ=255)
                                                                                            aGTGGGGTCGGGGTTTGGGCCTTCGCGGGCATGCCCGCTCCCACAGGTACAGCACTGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCAg                                                           <  1:2177324/150‑1 (MQ=255)
                                                                                              tGGGGTAGGGGTTTGGGCCTTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCAGGc                                                         >  2:3648019/1‑150 (MQ=255)
                                                                                                      gggTTTGGGCCTTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCAGGCCACCGGGc                                                 >  2:2739901/1‑150 (MQ=255)
                                                                                                      gggTTTGGGCCTTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCAGGCCACCGGGc                                                 <  2:1103639/150‑1 (MQ=255)
                                                                                                      gggTTTGGGCCTTCGCGGGCATGCCCGCTCCCACAGGGACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGTGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCAGGCCACCGGGc                                                 <  2:56353/150‑1 (MQ=255)
                                                                                                      gggTTTGGGCCTTCGCGGGCATGCCCGCACCCACAGGTACAGCACAGGTTTAAAAGGTTGTGGTGTTACTGTGGGAGCGGGCGTGCCCGCGACAGGGCCCGAATGGGTTCTATAGTGCCCGAGGGTTACGCAGCGGTCAGGCCACCGGGc                                                 >  1:1683188/1‑150 (MQ=38)
                                                                                                                cTTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCAGGCCACCGGGCCACGCGGATa                                       >  1:103269/1‑150 (MQ=255)
                                                                                                                      gggCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGGTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGACAGGGCCCGAATGGCTGCTCTAGTGCCCGAGGAGTACGCAGCTGTCAGGCCACCGGGCCACGCGGATAGCGCtg                                 >  2:3874169/1‑149 (MQ=38)
                                                                                                                       ggCATGCCCGCTCCCCCAGGTACAGCACCGGTGTCAAAGGTGGGGGTGTTCCTGTGGGTGCGGGCGGGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCAGGCCACCTGGCCACGCGGATAGCGCTTc                                <  2:4587196/150‑1 (MQ=255)
                                                                                                                        gCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCAGGCCACCGGGCCACGCGGATAGCGCTTCa                               >  1:2323429/1‑150 (MQ=255)
                                                                                                                          aTGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCAGGCCACCGGGCCACGCGGATAGCGCTTCAGc                             <  2:6334744/150‑1 (MQ=255)
                                                                                                                           TGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTCTAGTGCCCGAGGATTACGCAGCTGTCAGGCCACCGGGCCACGCGGATAGCGCTTCAGcc                            <  2:3151650‑M1/150‑1 (MQ=255)
                                                                                                                            GCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCAGGCCACCGGGCCACGCGGCTCGCGCTTCAGCCg                           <  1:3058566‑M1/150‑1 (MQ=255)
                                                                                                                            GCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCAGGCCACCGGGCCACGCGGATAGCGCTTCAGCCg                           <  1:6652535‑M1/150‑1 (MQ=255)
                                                                                                                              CCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCAGGCCACCGGGCCACGCGGATAGCGCTTCAGCCGCt                         <  2:2391274‑M1/150‑1 (MQ=255)
                                                                                                                               CGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCAGGCCACCGGGCCACGCGGATAGCGCTTCAGCCGCTc                        >  1:5506405‑M1/1‑150 (MQ=255)
                                                                                                                                  cCCCACCGGTACAGCACCGGTTTCAACGGTTGTGGTGTTCGTGTGGGCGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCAGGCCACCGGGCCACGCGGATAGCGCTTCAGCCGCTCGGc                     >  1:6384829‑M1/2‑150 (MQ=255)
                                                                                                                                  TCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGCAAGGGCCCGAATGGCTTCTAGAGTGCCCGAGGGTTACGCCGCTGTCAGGCCCCCGGGCCACGCGGCTAGCGCTGCAGCCGCTCGGc                     <  2:4736984‑M1/150‑1 (MQ=255)
                                                                                                                                    CCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCAGGCCACCGGGCCACGCGGATAGCGCTTCAGCCGCTCGGCCa                   <  1:5727626‑M1/150‑1 (MQ=255)
                                                                                                                                    CCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCAGGCCACCGGGCCACGCGGATAGCGCTTCAGCCGCTCGGCCa                   <  2:724109‑M1/150‑1 (MQ=255)
                                                                                                                                    CCACAGGTACAGCACAGGTTTCAAAGGGTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCAGGCCACCGGGCCACGCGGATAGCGCTTCAGCCGCTCGGCCa                   >  2:957917‑M1/1‑150 (MQ=255)
                                                                                                                                     CACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCAGGCCACCGGGCCACGCGGATAGCGCTTCAGCCGCTCGGCCAg                  <  2:840664‑M1/150‑1 (MQ=255)
                                                                                                                                      ACAGGAACAGCACAGCTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCAGGCCACCGGGCCACGCGGATAGCGCTTCAGCCGCTCGGCCAGt                 <  1:1358731‑M1/150‑1 (MQ=255)
                                                                                                                                        cGGTACAGCACAGGTTTCAAAGGTCGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCAGGCCACCGGGCCACGCGGATAGCGCTTCAGCCGCTCGGCCAGTTc               >  2:2085577‑M1/2‑150 (MQ=255)
                                                                                                                                          gTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCAGGCCACCGGGCCACGCGGATAGCGCTTCAGCCGCTCGGCCAGTTcgc             <  1:4993176/150‑1 (MQ=255)
                                                                                                                                              aGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCAGGCCACCGGGCCACGCGGATAGCGCTTCAGCCGCTCGGCCAGTTCGCGccc         >  2:3243910/1‑150 (MQ=255)
                                                                                                                                                cacaGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCAGGCCACCGGGCCACGCGGATAGCGCTTCAGCCGCTCGGCCAGTTCGCGCCCgg       >  1:583173/1‑150 (MQ=255)
                                                                                                                                                    ggTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCAGGCCACCGGGCCACGCGGATAGCGCTTCAGCCGCTCGGCCAGTTCGCGCCCGGGAAt   >  2:887062/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GGCCAGGCAGGGCCTGTTCCGCAGTGGGCATGACGTTTTTGTTCACCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCCTTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCAGGCCACCGGGCCACGCGGATAGCGCTTCAGCCGCTCGGCCAGTTCGCGCCCGGGAATC  >  NC_002947/403131‑403429

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: