Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,901,439 G→A 8.8% intergenic (+78/‑130) PP_4304 → / → sbp‑I VIC family cation transporter/sulfate ABC transporter substrate‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,901,4390GA 8.8% 580.8 / 35.2 212intergenic (+78/‑130)PP_4304/sbp‑IVIC family cation transporter/sulfate ABC transporter substrate‑binding protein
Reads supporting (aligned to +/- strand):  ref base G (106/76);  new base A (6/12);  total (116/96)
Fisher's exact test for biased strand distribution p-value = 4.91e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.63e-01

Only 100 of 216 total aligned reads displayed.

CACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAAT‑CCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCAACACCCTTCAAGAACAAGGACGCAACGTGAAAAAACTCTTCACCGC  >  NC_002947/4901290‑4901588
                                                                                                                                                         |                                                                                                                                                     
cacaAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAAT‑CCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCggg                                                                                                                                                       <  1:1334292/150‑1 (MQ=255)
 ggggcagattgccagcgactataccgggtgggtttatacgggttaaagaacaaaaaatgctttggttatgcttctgcctgtGCCGGCCTCTTCGCGGGTAAAT‑CCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGt                                                                                                                                                      >  2:5621521‑M1/82‑150 (MQ=255)
 acaAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAAT‑CCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGt                                                                                                                                                      >  2:1932340‑M2/1‑149 (MQ=255)
  caAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAAT‑CCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGtt                                                                                                                                                     >  2:3558310‑M2/1‑148 (MQ=255)
    gcagattgccagcgactataacgggtgggtttatacgggttaaagaacaaaaaatgctttggttatgcttctgcctgtGCCGGCCTCTTCGCGGGTAAAT‑CCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTa                                                                                                                                                   >  1:3800032‑M1/79‑150 (MQ=255)
     gaagaCCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAAT‑CCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGtttac                                                                                                                                                  >  2:5356861‑M2/1‑145 (MQ=255)
     gaagaCCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAAT‑CCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGtttac                                                                                                                                                  >  2:1241079‑M2/1‑145 (MQ=255)
     gaagaCCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCGCGTCGCGGGTAAAT‑CCGCTCCTACAGGCGCGGACCCGGCGCGGGCC‑‑‑CCCGGAGGAGCGGGtttac                                                                                                                                                  <  1:393964‑M2/150‑6 (MQ=255)
       agaCCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATTGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAAT‑CCGCTCCTACCGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGtttaccc                                                                                                                                                >  2:2443295‑M2/1‑143 (MQ=255)
       agaCCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAAT‑CCGCTCCTACAGGCGCGGACCGTGCGCGGGCC‑‑‑CCCGGAGGAGCGGGtttaccc                                                                                                                                                >  2:5266905‑M2/1‑143 (MQ=255)
         aCCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAAT‑CCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGtttacccgc                                                                                                                                              >  1:4750169‑M2/1‑141 (MQ=255)
           ccACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCGCGTCGCGGGTAAAT‑CCCCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGGAGGAGCGGGtttacccccgc                                                                                                                                            <  2:2518981‑M2/150‑12 (MQ=255)
            ccagcgactataacgggtgggtttatacgggttaaagaacaaaaaatgctttggttatgcttctgcctgtGCTGGCCTCTTCGCGGGTAAAT‑CCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGaa                                                                                                                                           >  1:3728324‑M1/71‑150 (MQ=255)
             aCGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCGTCGCGGGTAAAT‑ACCCTCCTACACGCGCCGGCCCTGCGCGGTGC‑‑‑CCCGTGGGGGCGGGtttaaccccgaag                                                                                                                                          <  2:1394378‑M2/150‑14 (MQ=255)
              cGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAAT‑CCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGtttacccgcgaaga                                                                                                                                         >  2:3908094‑M2/1‑136 (MQ=255)
               gAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAAT‑CCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCTTAGGAGCGGGtttacccgcgaagag                                                                                                                                        >  1:4362503‑M2/1‑135 (MQ=255)
                   aGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAAT‑CCGCACCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGtttacccgcgaagaggccg                                                                                                                                    <  2:1642503‑M2/150‑20 (MQ=255)
                      gggCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCGTCGCCGGTAAAT‑ACGCTCCTACAGGCGCGGCCCGTGCGCGGGCC‑‑‑CCCGGGGGAGCGGGTTTGACCGCCAAGAGGCCGGca                                                                                                                                 >  1:393451/1‑150 (MQ=255)
                          ggcgTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAAT‑CCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAgg                                                                                                                             <  2:3633964/150‑1 (MQ=255)
                                   ttctCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTTAAT‑CCGCGCCTACAGGCGCGGACCGGGCGCGGTCC‑‑‑CCCGTTGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCAt                                                                                                                    >  1:3722574/1‑150 (MQ=255)
                                   ttctCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAAT‑CCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCAt                                                                                                                    <  2:4232285/150‑1 (MQ=255)
                                   ttctCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCGTCGCGGGTAAAT‑CCGCGCCTACAGGCGCGGCCCGTCCGCGGTCC‑‑‑CCCGGAGGAGCGGGTTTTACCGCGACGCGGCCGGCACAGGCAGAAGCag                                                                                                                    >  2:3475896/1‑149 (MQ=255)
                                         cTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAAT‑CCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCaa                                                                                                              <  2:2797689/150‑1 (MQ=255)
                                                aaTGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAAT‑CCGCTCCTACAGGCGCGGACCGGGCGCGGTCC‑‑‑CCCGGAGGCGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCCGAAGCATAACCAAAGCAttt                                                                                                       >  2:4312003/1‑150 (MQ=255)
                                                   gggTTGTTTCCCAAAGCGTGATAGCGCCTGCGCGGGCCTCTTCGCGGGTAAAT‑CCGCTCCTACAGGCGCGTACGGTGCGCGGTCC‑‑‑CCCTTTGGCGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCAtttttt                                                                                                    <  2:2364286/150‑1 (MQ=38)
                                                   gggTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAAT‑CCGCTCCTCCAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCTTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCAtttttt                                                                                                    <  2:5168943/150‑1 (MQ=255)
                                                    ggTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAAT‑CCCCTCCTACAGGCGCGGCCCGTGCGCGGGCC‑‑‑CCCCTAGGAGCGGGTTTACCCGCCAAGAGGCCGGCACAGGCAGAAGCATAAACAAAGCATTTTTTg                                                                                                   >  1:5070762/1‑150 (MQ=38)
                                                    ggTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAAT‑ACCCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCAatttttt                                                                                                   >  1:5071034/1‑149 (MQ=255)
                                                     gttgttTCCCAAAGCGTGATAGCGCCTGCGCCGGCCAATTCGCGGGTAAAT‑ACGCTACTAAAGGCGCGGACCCTGCGCGGTCC‑‑‑CCCCTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAAACAAAGCATTTTTTGt                                                                                                  >  2:776701/1‑150 (MQ=38)
                                                           tCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAAT‑CCGCTCCTACAGGCGCGGACCGTGCGCGGGCC‑‑‑CCCGGAGGCGCGGGTTTACCCGCGAAGAGGCCGGCCCCGGCAGAAGCATAACCCAAGCATTTTTGGTTCTTTa                                                                                            >  2:3872170/1‑150 (MQ=38)
                                                             ccAAAGCGTGATACCGCCTGCGCCGGCCTCTTCGCGGGTAAAT‑CCGCTCCTACAGGCGCGGACCGTGCGCGGACC‑‑‑CCCGTAGGAGCGTTTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAAc                                                                                          <  2:2079380/150‑1 (MQ=255)
                                                                      gATTGCGCCTGCGCCGGCCTCTTCGCGGGTAAAT‑CCGCTCCTCCCGGCGCGGACCGTGCGCGGCCC‑‑‑CCCGTTGGAGGGGTTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATaaa                                                                                 <  2:1446989/150‑1 (MQ=255)
                                                                           cgcgTGCGCGGGCCTCTT‑GGGGGTAAATCCCGCCCCTACCGGCGGGGACGGTGCGGGGCCC‑‑‑CCCGTCGGAGGGTTTTTACCCGCTAAGCGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAacccac                                                                            <  1:2414019/150‑1 (MQ=255)
                                                                               tGCGCCGGCCTCTTCGCGGGTAAAT‑CCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGtttac                                                                                                                                                  >  1:2161661‑M2/1‑71 (MQ=255)
                                                                                         TCTTCGGGGGTAAAT‑CCGCACCTACCGGCGCGGACCGTGCGCGGCCC‑‑‑CCCGTAGGAGCGGTTTAACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCt                                                              >  2:2354642‑M1/1‑150 (MQ=255)
                                                                                         TCTTCGCGGGTAAAT‑CCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGc                                                               >  2:5236745‑M1/1‑149 (MQ=255)
                                                                                         TCTTCGCGGGTAAAT‑CCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGc                                                               <  1:5237134‑M1/149‑1 (MQ=255)
                                                                                         TCTTCGCGGGTAAAT‑CCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCt                                                              <  1:913926‑M1/150‑1 (MQ=255)
                                                                                           TTCGCGGGTAAAC‑CCGCTCCTACGGG‑‑‑GGACCGCGCACGGTCCGCGCCTGTAGGAGCGGATTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTgg                                                            <  1:102031‑M1/150‑1 (MQ=255)
                                                                                           TTCGCGGGTAAAC‑CCGCTCCTACGGG‑‑‑GGACCGCGCACGGTCC‑‑‑GCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTATTTGTTCTTTAACCCGTATAAACCCACCCGGTATAGTCGCAggcta                                                         <  2:407765‑M1/150‑3 (MQ=255)
                                                                                            TCGCGGGTAAAC‑CCGCTCCTACGGG‑‑‑GGACCGCGCACGGTCCGCGCCTGTAGGAGCGGATTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGc                                                           <  2:3226270‑M1/150‑1 (MQ=255)
                                                                                               CGGGTAAAC‑CCGCTCCTACGGG‑‑‑GGACCGCGCACGGTCCGCGCCTGTAGGAGCGGATTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAAt                                                        <  2:4358658‑M1/150‑1 (MQ=255)
                                                                                                 GGTAAAT‑CCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCt                                                      <  1:3211262‑M1/150‑1 (MQ=255)
                                                                                                 GGTAAAC‑CCGCTCCTACGGG‑‑‑GGACCGCGCACGGTCCGCGCCTGTAGGAGCGGATTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCt                                                      <  1:5441932‑M1/150‑1 (MQ=255)
                                                                                                   TAAAT‑CCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCATTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGc                                                    <  2:4593551‑M1/150‑1 (MQ=255)
                                                                                                    AAAT‑CCGCTCCTACCGGCGCGGACGGTGCGCGGCCC‑‑‑CCCGTTGGTGGGGTTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGcc                                                   >  1:3356963‑M1/1‑150 (MQ=255)
                                                                                                    AAAT‑CCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGcc                                                   <  1:3628743‑M1/150‑1 (MQ=255)
                                                                                                     AAT‑CCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTCTAGTCGCTGGCACTCTGccc                                                  <  2:5017562‑M1/150‑1 (MQ=255)
                                                                                                       t‑CCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTCTAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCa                                                >  2:3955130/1‑150 (MQ=255)
                                                                                                       t‑CCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAGCCCACCCGTTATAGTCGCTGGCAATCTGCCCCa                                                >  1:3833705/1‑150 (MQ=255)
                                                                                                       t‑CCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCa                                                >  2:570172/1‑150 (MQ=255)
                                                                                                       t‑CCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCa                                                >  1:4916142/1‑150 (MQ=255)
                                                                                                       t‑CCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCa                                                >  2:2636957/1‑150 (MQ=255)
                                                                                                       t‑CCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCa                                                >  1:4540227/1‑150 (MQ=255)
                                                                                                       t‑CCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCa                                                >  1:2834005/1‑150 (MQ=255)
                                                                                                       t‑CCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTGATAGTCGCTGGCAATCTGCCCCa                                                >  2:3661247/1‑150 (MQ=255)
                                                                                                       t‑CCGCGCCTGTAGGAGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCATATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCa                                                >  2:5077666/1‑150 (MQ=38)
                                                                                                         ccGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTATACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAATCCCACCCGTTATAGTCGCTGGCAATCTGCCCCaa                                               >  1:4937029/1‑150 (MQ=255)
                                                                                                           gCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGtttacccgcgaagaggccggcgcaggcgctatcacgctttgggaaacaacccattcccacagcgagaacagaacgccgccccctgctcgtgggtcttcttgtggca                                             <  2:2000872‑M2/150‑107 (MQ=255)
                                                                                                           gCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCAaca                                             >  2:1159291/1‑150 (MQ=255)
                                                                                                               cTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCAACACCCt                                         >  2:2709190/1‑150 (MQ=255)
                                                                                                                tACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGtttacccgcgaagaggccggcgcaggcgctatcacgctttgggaaacaacccattcccacagcgagaacagaacgccgccccctgctcgtgggtcttcttgtggcaggtcg                                        <  1:303240‑M2/150‑112 (MQ=255)
                                                                                                                tACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGtttacccgcgaagaggccggcgcaggcgctatcacgctttgggaaacaacccattcccacagcgagaacagaacgccgccccctgctcgtgggtcttcttgtggcaggtcg                                        <  2:2380474‑M2/150‑112 (MQ=255)
                                                                                                                tACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGtttacccgcgaagaggccggcgcaggcgctatcacgctttgggaaacaacccattcccacagcgagaacagaacgccgccccctgctcgtgggtcttcttgtggcaggtcg                                        <  2:33942‑M2/150‑112 (MQ=255)
                                                                                                                tACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCACTCTGCCCCAACACCCtt                                        >  2:3992040/1‑150 (MQ=255)
                                                                                                                tACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCAACACCCtt                                        >  1:2608382/1‑150 (MQ=255)
                                                                                                                tACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGAtttacccgcgaagaggccggcgcaggcgctatcacgctttgggaaacaacccattc                                                                                               <  1:1127267‑M2/95‑57 (MQ=255)
                                                                                                                  cAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGtttacccgcgaagaggccggcgcaggcgctatcacgctttgggaaacaacccattcccacagcgagaacagaacgccgccccctgctcgtgggtcttcttgtggcaggtcggg                                      <  1:5180307‑M2/150‑114 (MQ=255)
                                                                                                                  cAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGAtttacccgcgaagaggccggcgcaggcgctatcacgctttgggaaacaacccattcccacagcgagaacagaacgccgccccctgctcgtgggtcttcttgtggcaggtcggg                                      <  2:3362179‑M2/150‑114 (MQ=255)
                                                                                                                  cAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGAtttacccgcgaagaggccggcgcaggcgctatcacgctttg                                                                                                              >  2:1487984‑M2/1‑37 (MQ=255)
                                                                                                                   cggCGGGGCCCGGGCGCGGCCC‑‑‑CCCGGTGGCGGGGGTTTACCCGCTAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCAACACCCTTCaa                                     <  1:693818/149‑1 (MQ=255)
                                                                                                                   aGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGAtttacccgcgaagaggccggcgcaggcgctatcacgctttgggaaacaacccattcccacagcgagaacagaacgccgccccctgctcgtgggtcttcttgtggcaggtcgggc                                     <  2:4111090‑M2/150‑115 (MQ=255)
                                                                                                                     gcgcGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCAACACCCTTCAAGa                                   <  2:1356712/150‑1 (MQ=255)
                                                                                                                       gcgGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCAACACCCTTCAAGAAc                                 >  2:3014165/1‑150 (MQ=255)
                                                                                                                        gggaCCGCGCACGGTCC‑‑‑CCCGTAGGAGCGGAtttacccgcgaagaggccggcgcaggcgctatcacgctttgggaaacaacccattcccacagcgagaacagaacgccgccccctgctcgtgggtcttcttgtggcaggtcgggcaatca                                <  1:181878‑M2/149‑120 (MQ=255)
                                                                                                                         gtaCCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGAATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGATGGCAATCTGCCCCAAAACCCTTCAAGAACaa                               <  2:457155/148‑1 (MQ=255)
                                                                                                                         ggACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAACCATATATTGTTCTTTAACCCGTATAAACCCACCCCTTATAGTCGCTGGCAATCTGCCCCAACACCCTTcaacaacaa                               <  1:4012487/150‑1 (MQ=38)
                                                                                                                           aCCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGACGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCAACACCCTTCAAGAACAAgg                             >  2:2867403/1‑150 (MQ=255)
                                                                                                                             ggtggggggTCC‑‑‑CCGGTTGGTGGGGGTTTACCGGCTAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCAACACCCTTCAAGAACAAGGAc                           <  1:2518981/143‑1 (MQ=255)
                                                                                                                               tGCGGGGTCC‑‑‑CCCGTTGGAGGGGTTTTACCGGCGAAGAGGCCGGCACAGTCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAATCCCACGCGTTATAGTCGCTGGCAATCTGCCCCAACACCCTTCAAGAACAAGGACGc                         <  2:2941155/150‑1 (MQ=255)
                                                                                                                               tGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCAACACCCTTCAAGAACAAGGACGc                         >  1:2944815/1‑150 (MQ=255)
                                                                                                                                 cgcgGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCAACACCCTTCAAGAACAAGGACGCaa                       >  2:614173/1‑150 (MQ=255)
                                                                                                                                 cgcgGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCAACACCCTTCAAGAACAAGGACGCaa                       >  1:1591454/1‑150 (MQ=255)
                                                                                                                                  gcgGTCC‑‑‑CCCGTAGGAGCGGGtttacccgcgaagaggccggcgcaggcgctatcacgctttgggaaacaacccattcccacagcgagaacagaacgccgccccctgctcgtgggtcttcttgtggcaggtcgggcaatcatgctgcagct                      <  1:1503100‑M2/150‑130 (MQ=255)
                                                                                                                                     gTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCAACACCCTTCAAGAACAAGGACGCAACGTg                   >  1:804513/1‑150 (MQ=255)
                                                                                                                                        c‑‑‑cccGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCAACACCCTTCAAGAACAAGGACGCAACGTGaaa                >  1:3427671/1‑150 (MQ=255)
                                                                                                                                            cccGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCAACACCCTTCAAGAACAAGGACGCAACGTGaaaa               >  2:574088/1‑150 (MQ=255)
                                                                                                                                            cccGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCAACACCCTTCAAGAACAAGGACGCAACGTGaaaa               >  2:5275105/1‑150 (MQ=255)
                                                                                                                                              cGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCAACACCCTTCAAGAACAAGGACGCAACGTGaaaaaa             <  2:393964/150‑1 (MQ=255)
                                                                                                                                              cGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCAACACCCTTCAAGAACAAGGACGCAACGTGaaaaaa             <  1:2321727/150‑1 (MQ=255)
                                                                                                                                                 aGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCACAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCAACACCCTTCAAGAACAAGGACGCAACGTGAAAAAActc          <  2:2293446/150‑1 (MQ=255)
                                                                                                                                                 aGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCAACACCCTTCAAGAACAAGGACGCAACGTGAAAAAActc          <  1:5126434/150‑1 (MQ=255)
                                                                                                                                                   gAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCAACACCCTTCAAGAACAAGGACGCAACGTGAAAAAACTCtt        >  2:3162945/1‑150 (MQ=255)
                                                                                                                                                   gAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCAACACCCTTCAAGAACAAGGACGCAACGTGAAAAAACTCtt        >  2:1943567/1‑150 (MQ=255)
                                                                                                                                                     gCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCAACACCCTTCAAGAACAAGGACGCAACGTGAAAAAACTCTTCa      >  1:331428/1‑150 (MQ=255)
                                                                                                                                                     gCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCAACACCCTTCAAGAACAAGGACGCAACGTGAAAAAACTCTTCa      >  1:1105788/1‑150 (MQ=255)
                                                                                                                                                     gCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCAACACCCTTCAAGAACAAGGACGCAACGTGAAAAAACTCTTCa      >  1:3438321/1‑150 (MQ=255)
                                                                                                                                                     gCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCAACACCCTTCAAGAACAAGGACGCAACGTGAAAAAACTCTTCa      >  1:2083974/1‑150 (MQ=255)
                                                                                                                                                         gtttacccgcgaagaggccggcgcaggcgctatcacgctttgggaaacaacccattcccacagcgagaacagaacgccgccccctgctcgtgggtcttcttgtggcaggtcgggcaatcatgctgcagctgctccccacgcatggcattg  <  2:4686146‑M2/150‑150 (MQ=255)
                                                                                                                                                         gTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCAACACCCTTCAAGAACAAGGACGCAACGTGAAAAAACTCTTCACCGc  >  1:3044623/1‑150 (MQ=255)
                                                                                                                                                         |                                                                                                                                                     
CACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAAT‑CCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCAACACCCTTCAAGAACAAGGACGCAACGTGAAAAAACTCTTCACCGC  >  NC_002947/4901290‑4901588

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: ‑