Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,403,895 C→G 5.3% intergenic (+123/‑119) aroQ‑III → / → PP_3004 type II 3‑dehydroquinate dehydratase/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,403,8950CG 5.3% 709.0 / 21.7 208intergenic (+123/‑119)aroQ‑III/PP_3004type II 3‑dehydroquinate dehydratase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (83/114);  new base G (6/5);  total (89/119)
Fisher's exact test for biased strand distribution p-value = 5.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.81e-01

Only 100 of 209 total aligned reads displayed.

GCGTGCGGCGATGAGCATGGTCTGAGCAGGGCGCAGGGGTCGGATCTCCGAGCGCGTCTACCCGCGAAGCAGGTGGCGCGGGGCATGCGGCGACGGCGCCTGAAGTGACCGTAGGTGTGGTGGATGTCGTGCTGCCTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTCCGCCTTCAACAAGCCACCCGGTTCGTTGCCTCCAACCCCACGGACCGGGGCTCTCATGACCTGGGAGGTCGAATGGACGTCGCTACAACGATCAGTCAACAAG  >  NC_002947/3403748‑3404044
                                                                                                                                                   |                                                                                                                                                     
gcgtgcgGCGATGAGCATGGTCTGAGCAGGGCGCAGGGGTCGGATCTCCGAGCGCGTCTACCCGCGAAGCAGGTGGCGCGGGGCATGCGGCGACGGCGCCTGAAGTGACCGTAGGTGTGGTGGATGTCGTGCTGCCTACCCATCCAGCat                                                                                                                                                     >  2:933977/1‑150 (MQ=255)
 cgtgcgGCGATGAGCATGGTCTGAGCAGGGCGCAGGGGTCGTATCTCCGATCGCGTCTACCCGCGAAGCAGGTGGCGCGGGGCATGCTGCGACGGCGCCTGAAGTGACCGTAGGTGTGGTGGATGTCGTGCTGCCTACCCATCCAGCatg                                                                                                                                                    <  2:375296/150‑1 (MQ=255)
    gcggcgATGAGCATGGTCTGAGCAGGGCGCAGGGGTCGGATCTCCGAGCGCGTCTACCCGCGAAGCAGGTGGCGCGGGGCATGCGGCGACGGCGCCTGAAGTGACCGTAGGTGTGGTGGATGTCGTGCTGCCTACCCATCCAGCatggat                                                                                                                                                 <  1:1642173/150‑1 (MQ=255)
    gcggcgATGAGCATGGTCTGAGCAGGGCGCAGGGGTCGGATCTCCGAGCGCGTCTACCCGCGAAGCAGGTGGCGCGGGGCATGCGGCGACGGCGCCTGAAGTGACCGTAGGTGTGGTGGATGTCGTGCTGCCTACCCATCCAGCatggat                                                                                                                                                 >  1:2496202/1‑150 (MQ=255)
        cgATGAGCATGGTCTGAGCAGGGCGCAGGGGTCGGATCTCCGAGCGCGTCTACCCGCGAAGCAGGTGGCGCGGGGCATGCGGCGACGGCGCCTGAAGTGACCGTAGGTGTGGTGGATGTCGTGCTGCCTACCCATCCAGCATGGATGGaa                                                                                                                                             >  1:3324759/1‑150 (MQ=255)
        cgATGAGCATGGTCTGAGCAGGGCGCAGGGGTCGGATCTCCGAGCGCGTCTACCCGCGAAGCAGGTGGCGCGGGGCATGCGGCGACGGCGCCTGAAGTGACCGTAGGTGTGGTGGATGTCGTGCTGCCTACCCATCCAGCATGGATGGaa                                                                                                                                             <  1:738442/150‑1 (MQ=255)
          aTGAGCATGGTCTGAGCAGGGCGCAGGGGTCGGATCTCCGAGCGCGTCTACCCGCGAAGCAGGTGGCGCGGGGCATGCGGCGACGGCGCCTGAAGTGACCGTAGGTGTGGTGGATGTCGTGCTGCCTACCCATCCATGctggatgggtag                                                                                                                                           <  1:3427618/150‑5 (MQ=255)
             aGCATGGTCTGAGCAGGGCGCAGGGGTCGGATCTCCGAGCGCGTCTACCCGCGAAGCAGGTGGCGCGGGGCATGCGGCGACGGCGCCTGAAGTGACCGTAGGTGTGGTGGATGTCGTGCTGCCTACCCATCCAGCATGGATGGAAAGCCa                                                                                                                                        <  2:5187177/150‑1 (MQ=255)
              gCATGGTCTGAGCAGGGCGCAGGGGTCGGATCTCCGGGCGCGTCTACCCGCGAAGCAGGTGGCGCGGGGCATGCGGCGACGGCGCCTGAAGTGACCGTAGGTGTGGTGGATGTCGTGCTGCCTACCCATCCAGCATGGATGGAAAGCCAg                                                                                                                                       >  2:1930894/1‑150 (MQ=255)
              gCATGGTCTGAGCAGGGCGCAGGGGTCGGATCTCCGAGCGCGTCTACCCGCGAAGCAGGTGGCGCGGGGCATGCGGCGACGGCGCCTGAAGTGACCGTAGGTGTGGTGGATGTCGTGCTGCCTACCCATCCAGCATGGATGGAAAGCCAg                                                                                                                                       >  1:5009093/1‑150 (MQ=255)
               cATGGTCTGAGCAGGGCGCAGGGGTCGGATCTCCGAGCGCGTCTACCCGCGAAGCAGGTGGCGCGGGGCATGCGGCGACGGCGCCTGAAGTGACCGTAGGTGTGGTGGATGTCGTGCTGCCTACCCATCCAGCATGGATGGAAAGCCAGt                                                                                                                                      <  1:1567486/150‑1 (MQ=255)
                    tCTGAGCAGGGCGCAGGGGTCGGATCTCCGAGCGCGTCTACCCGCGAAGCAGGTGGCGCGGGGCATGCGGCGACGGCGCCTGAAGTGACCGTAGGTGTGGTGGATGTCGTGCTGCCTACCCATCCAGCATGGATGGAAAGCCAGTAGCGa                                                                                                                                 <  2:1372792/150‑1 (MQ=255)
                      tGAGCAGGGCGCAGGGGTCGGATCTCCGAGCGCGTCTACCCGCGAAGCAGGTGGCGCGGGGCATGCGGCGACGGCGCCTGAAGTGACCGTAGGTGTGGTGGATGTCGTGCTGCCTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACt                                                                                                                               >  2:3473185/1‑150 (MQ=255)
                       gAGCAGGGCGCAGGGGTCGGATCTCCGAGCGCGTCTACCCGCGAAGCAGGTGGCGCGGGGCATGCGGCGACGGCGCCTGAAGTGACCGTAGGTGTGGTGGATGTCGTGCTGCCTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACtt                                                                                                                              <  1:5605632/150‑1 (MQ=255)
                         gCAGGGCGCAGGGGTCGGATCTCCGAGCGCGTCTACCCGCGAAGCAGGTGGCGCGGGGCATGCGGCGACGGCGCCTGAACTGACCGTAGGTGTGGTGGATGTCGTGCTGCCTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTcc                                                                                                                            <  2:1673422/150‑1 (MQ=255)
                           aGGGCGCAGGGGTCGGATCTCCGAGCGCGTCTACCCGCGAAGCAGGTGGCGCGGGGCATGCGGCGACGGCGCCTGAAGTGACCGTAGGTGTGGTGGATGTCGTGCTGCCTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGt                                                                                                                          <  2:1889638/150‑1 (MQ=255)
                                gcAGGGGTCGGATCTCCGAGCGCGTCTACCCGCGAAGCAGGTGGCGCGGGGCATGCGGCGACGGCGCCTGAAGTGACCGTAGGTGTGGTGGATGTCGTGCTGCCTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTg                                                                                                                     >  1:4480215/1‑150 (MQ=255)
                                  aGGGGTCGGATCTCCGAGCGCGTCTACCCGCGAAGCAGGTGGCGCGGGGCATGCGGCGACGGCGCCTGAAGTGACCGTAGGTGTGGTGGATGTCGTGCTGCCTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGct                                                                                                                   >  1:5474215/1‑149 (MQ=255)
                                  aGGGGTCGGATCTCCGAGCGCGTCTACCCGCGAAGCAGGTGGCGCGGGGCATGCGGCGACGGCGCCTGAAGTGACCGTAGGTGTGGTGGATGTCGTGCTGCCTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCa                                                                                                                   >  2:3169797/1‑150 (MQ=255)
                                  aGGGGTCGGATCTCCGAGCGCGTCTACCCGCGAAGCAGGTGGCGCGGGGCATGCGGCGACGGCGCCTGAAGTGACCGTAGGTGTGGTGGATGTCGTGCTGCCTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCa                                                                                                                   >  1:1849325/1‑150 (MQ=255)
                                    gggTCGGATCTCCGAGCGCGTCTACCCGCGAAGCAGGTGGCGCGGGGCATGCGGCGACGGCGCCTGAAGTGACCGTAGGTGTGGTGGATGTCGTGCTGCCTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATa                                                                                                                 <  1:1736820/150‑1 (MQ=255)
                                         ggATCTCCGAGCGCGTCTACCCGCGAAGCAGGTGGCGCGGGGCATGCGGCGACGGCGCCTGAAGTGACCGTAGGTGTGGTGGATGTCGTGCTGCCTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACg                                                                                                            >  2:1824242/1‑150 (MQ=255)
                                          gATCTCCGAGCGCGTCTACCCGCGAAGCAGGTGGCGCGGGGCATGCGGCGACGGCGCCTGAAGTGACCGTAGGTGTGGTGGATGTCGTGCTGCCTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACgg                                                                                                           <  1:5583023/150‑1 (MQ=255)
                                          gATCTCCGAGCGCGTCTACCCGCGAAGCAGGTGGCGCGGGGCATGCGGCGACGGCGCCTGAAGTGACCGTAGGTGTGGTGGATGTCGTGCTGCCTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACgg                                                                                                           <  1:3879245/150‑1 (MQ=255)
                                          gATCTCCGAGCGCGTCTACCCGCGAAGCAGGTGGCGCGGGGCATGCGGCGACGGCGCCTGAAGTGACCGTAGGTGTGGTGGATGTCGTGCTGCCTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACgg                                                                                                           <  1:2331974/150‑1 (MQ=255)
                                          gATCTCCGAGCGCGTCTACCCGCGAAGCAGGTGGCGCGGGGCATGCGGCGACGGCGCCTGAAGTGACCGTAGGTGTGGTGGATGTCGTGCTGCCTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACgg                                                                                                           <  1:1062105/150‑1 (MQ=255)
                                           aTCTCCGAGCGCGTCTACCCGCGAAGCAGGTGGCGCGGGGCATGCGGCGACGGCGCCTGAAGTGACCGTAGGTGTGGTGGATGTCGTGCTGCCTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGt                                                                                                          >  1:2970315/1‑150 (MQ=255)
                                            tctcCGAGCGCGTCTACCCGCGAAGCAGGTGGCGCGGGGCATGCGGCGACGGCGCCTGAAGTGACCGTAGGTGTGGTGGATGTCGTGCTGCCTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTc                                                                                                         <  1:2554389/150‑1 (MQ=255)
                                              tcCGAGCGCGTCTACCCGCGAAGCAGGTGGCGCGGGGCATGCGGCGACGGCGCCTGAAGTGACCGTAGGTGTGGTGGATGTCGTGCTGCCTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTCCg                                                                                                       <  2:5131990/150‑1 (MQ=255)
                                               ccGAGCGCGTCTACCCGCGAAGCAGGTGGCGCGGGGCATGCGGCGACGGCGCCTGAAGTGACCGTAGGTGTGGTGGATGTCGTGCTGCCTACCCATCCATGctggatgggt                                                                                                                                             <  2:5124978/111‑3 (MQ=255)
                                               ccGAGCGCGTCTACCCGCGAAGCAGGTGGCGCGGGGCATGCGGCGACGGCGCCTGAAGTGACCGTAGGTGTGGTGGATGTCGTGCTGCCTACCCATCCATGctggatgggt                                                                                                                                             >  1:5124978/1‑109 (MQ=255)
                                                cGAGCGCGTCTACCCGCGAAGCAGGTGGCGCGGGGCATGCGGCGACGGCGCCTGAAGTGACCGTAGGTGTGGTGGATGTCGTGCTGCCTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGt                                                                                                                          <  2:5537842/129‑1 (MQ=255)
                                                   gcgcGTCTACCCGCGAAGCAGGTGGCGCGGGGCATGCGGCGACGGCGCCTGAAGTGACCGTAGGTGTGGTGGATGTCGTGCTGCCTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTCCGCCTTc                                                                                                  <  2:2776528/150‑1 (MQ=255)
                                                   gcgcGTCTACCCGCGAAGCAGCTGGCGCGGGGCATGCGGCGACGGCGCCTGAAGTGACCGTAGGTGTGGAGGATGTCGAGCTGCCTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTCCGCCTTc                                                                                                  <  2:4759614/150‑1 (MQ=255)
                                                      cgTCTACCCGCGAAGCAGGTGGCGCGGGGCATGCGGCGACGGCGCCTGAAGTGACCGTAGGTGTGGTGGATGTCGTGCTGCCTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTCCGCCTTcaac                                                                                               <  2:4993286/150‑1 (MQ=255)
                                                        tCTACCCGCGAAGCAGGTGGCGCGGGGCATGCGGCGACGGCGCCTGAAGTGACCGTAGGTGTGGTGGATGTCGTGCTGCCTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTCCGCCTTcaacaa                                                                                             <  2:4084152/150‑1 (MQ=255)
                                                        tCTACCCGCGAAGCAGGTGGCGCGGGGCATGCGGCGACGGCGCCTGAAGTGACCGTAGGTGTGGTGGATGTCGTGCTGCCTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTCCGCCTTcaacaa                                                                                             <  1:1982525/150‑1 (MQ=255)
                                                        tCTACCCGCGAAGCAGGTGGCGCGGGGCATGCGGCGACGGCGCCTGAAGTGACCGTAGGTGTGGTGGATGTCGTGCTGCCTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTCCGCCTTcaacaa                                                                                             <  1:2999118/150‑1 (MQ=255)
                                                          tACCCGCGAAGCAGGTGGCGCGGGGCATGCGGCGACGGCGCCTGAAGTGACCGTAGGTGTGGTGGATGTCGTGCTGCCTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTCCGCCTTCAACAAGc                                                                                           >  2:4828322/1‑150 (MQ=255)
                                                           aCCCGCGACGCAGGTGGCGCGGGGCATGCGGCGACGGCGCCTGAAGTGACCGTCGGTGGGGTGGATGTCGTGCTGCCTACCCATCCCGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATCGTACGGTCCGCCTGCAACACGcc                                                                                          >  2:3436973/1‑150 (MQ=255)
                                                           aCCCGCGAAGCAGGTGGCGCGGGGCATGCGGCGACGGCGCCTGAAGTGACCGTAGGTGTGGTGGATGTCGTGCTGCCTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTCCGCCTTCAACAAGc                                                                                           <  1:4322124/149‑1 (MQ=255)
                                                             ccGCGAAGCAGGTGGCGCGGGGCATGCGGCGACGGCGCCTGAAGTGACCGTAGGTGTGGTGGATGTCGTGCTGCCTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTCCGCCTTCAACAAGCCAc                                                                                        >  2:4245130/1‑150 (MQ=255)
                                                             ccGCGAAGCAGGTGGCGCGGGGCATGCGGCGACGGCGCCTGAAGTGACCGTAGGTGTGGTGGATGTCGTGCTGCCTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTCCGCCTTCAACAAGCCAc                                                                                        >  2:1174799/1‑150 (MQ=255)
                                                             ccGCGAAGCAGGTGGCGCGGGGCATGCGGCGACGGCGCCTGAAGTGACCGTAGGTGTGGTGGATGTCGTGCTGCCTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTCCGCCTTCAACAAGCCAc                                                                                        >  2:1027205/1‑150 (MQ=255)
                                                             ccGCGAAGCAGGTGGCGCGGGGCATGCGGCGACGGCGCCTGAAGTGACCGTAGGTGTGGTGGATGTCGTGCTGCCTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCAGTGACTGCATAGTACGGTCCGCCTTCAACAAGCCAc                                                                                        >  1:5305567/1‑150 (MQ=255)
                                                              cgcgAAGCAGGTGGCGCGGGGCATGCGGCGACGGCGCCTGAAGTGACCGTAGGTGTGGTGGATGTCGTGCTGCCTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTCCGCCTTCAACAAGCCAcc                                                                                       >  2:3258758/1‑150 (MQ=255)
                                                               gcgAAGCAGGTGGCGCGGGGCATGCGGCGACGGCGCCTGAAGTGACCGTAGGTGTGGTGGATGTCGTGCTGCCTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTCCGCCTTCAACAAGCCAccc                                                                                      >  2:4930101/1‑150 (MQ=255)
                                                               gcgAAGCAGGTGGCGCGGGGCATGCGGCGACGGCGCCTGAAGTGACCGTAGGTGTGGTGGATGTCGTGCTGCCTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTCCGCCTTCAACAAGCCAccc                                                                                      >  2:1491498/1‑150 (MQ=255)
                                                               gcgAAGCAGGTGGCGCGGGGCATGCGGCGACGGCGCCTGAAGTGACCGTAGGTGTGGTGGATGTCGTGCTGCCTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTCCGCCTTCAACAAGCCAccc                                                                                      <  2:2684061/150‑1 (MQ=255)
                                                                  aaGCAGGTGGCGCGGGGCATGCGGCGACGGCGCCTGAAGTGACCGTAGGTGTGGTGGATGTCGTGCTGCCTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTCCGCCTTCAACAAGCCACCCGGt                                                                                   <  1:5030997/150‑1 (MQ=255)
                                                                  aaGCAGGTGGCGCGGGGCATGCGGCGACGGCGCCTGAAGTGACCGTAGGTGTGGTGGATGTCGTGCTGCCTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTCCGCCTTCAACAAGCCACCCGGt                                                                                   <  1:4443942/150‑1 (MQ=255)
                                                                    gCAGGTGGCGCGGGGCATGCGGCGACGGCGCCTGAAGTGACCGTAGGTGTGGTGGATGTCGTGCTGCCTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTCCGCCTTCAACAAGCCACCCGGTTc                                                                                 <  1:5171427/150‑1 (MQ=255)
                                                                      aGGTGGCGCGGGGCATGCGGCGACGGCGCCTGAAGTGACCGTAGGTGTGGTGGATGTCGTGCTGCCTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTCCGCCTTCAACAAGCCACCCGGTTCGt                                                                               >  1:1068968/1‑150 (MQ=255)
                                                                         tGGCGCGGGGCATGCGGCGACGGCGCCTGAAGTGACCGTAGGTGTGGTGGATGTCGTGCTGCCTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTCCGCCTTCAACAAGCCACCCGGTTCGTTGc                                                                            <  2:227849/150‑1 (MQ=255)
                                                                          ggCGCGGGGCATGCGGCGACGGCGCCTGAAGTGACCGTAGGTGTGGTGGATGTCGTGCTGCCTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTCCGCCTTCAACAAGCCACCCGGTTCGTTGcc                                                                           <  2:4741141/150‑1 (MQ=255)
                                                                            cgcgGGGCATGCGGCGACGGCGCCTGAAGTGACCGTAGGTGTGGTGGATGTCGTGCTGCCTACCCATCCAGCATGGATGGAAAGCCAGTAGCGAATTCCGTGACTGCATAGTACGGTCCGCCTTCAACAAGCCACCCGGTTCGTTGCCTc                                                                         >  1:2673628/1‑150 (MQ=255)
                                                                             gcgGGGCATGCGGCGACGGCGCCTGAAGTGACCGTAGGTGTGGTGGATGTCGTGCTGCCTACCCATCCATGctggatgggta                                                                                                                                            <  1:207621/82‑4 (MQ=255)
                                                                             gcgGGGCATGCGGCGACGGCGCCTGAAGTGACCGTAGGTGTGGTGGATGTCGTGCTGCCTACCCATCCATGctggatgggta                                                                                                                                            >  2:207621/1‑79 (MQ=255)
                                                                               ggggCATGCGGCGACGGCGCCTGAAGTGACCGTAGGTGTGGTGGATGTCGTGCTGCCTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTCCGCCTTCAACAAGCCACCCGGTTCGTTGCCTCCaa                                                                      <  1:2887862/150‑1 (MQ=255)
                                                                                  gCATGCGGCGACGGCGCCTGAAGTGACCGTAGGTGTGGTGGATGTCGTGCTGCCTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTCCGCCTTCAACAAGCCACCCGGTTCGTTGCCTCCAAccc                                                                   >  2:948094/1‑150 (MQ=255)
                                                                                  gCATGCGGCGACGGCGCCTGAAGTGACCGTAGGTGTGGTGGATGTCGTGCTGCCTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTCCGCCTTCAACAAGCCACCCGGTTCGTTGCCTCCAAccc                                                                   <  1:3034868/150‑1 (MQ=255)
                                                                                    aTGCGGCGACGGCGCCTGAAGTGACCGTAGGTGTGGTGGATGTCGTGCTGCCTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTCCGCCTTCAACAAGCCACCCGGTTCGTTGCCTCCAACCCCa                                                                 <  2:1553412/150‑1 (MQ=255)
                                                                                      gcggcgACGGCGCCTGAAGTGACCGTAGGTGTGGTGGATGTCGTGCTGCCTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTCCGCCTTCAACAAGCCACCCGGTTCGTTGCCTCCAACCCCACg                                                               <  2:3658402/150‑1 (MQ=255)
                                                                                           gACGGCGCCTGAAGTGACCGTAGGTGTGGTGGATGTCGTGCTGCCTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTCCGCCTTCAACAAGCCACCCGGTTCGTTGCCTCCAACCCCACGGACCg                                                          <  1:3258758/150‑1 (MQ=255)
                                                                                             cGGCGCCTGAAGTGACCGTAGGTGTGGTGGATGTCGTGCTGCCTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTCCGCCTTCAACAAGCCACCCGGTTCGTTGCCTCCAACCCCACGGACCggg                                                        <  1:5527839/150‑1 (MQ=255)
                                                                                              ggCGCCTGAAGTGACCGTAGGTGTGGTGGATGTCGTGCTGCCTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTCCGCCTTCAACAAGCCACCCGGTTCGTTGCCTCCAACCCCACGGACCgggg                                                       >  1:3442570/1‑150 (MQ=255)
                                                                                              ggCGCCTGAAGTGACCGTAGGTGTGGTGGATGTCGTGCTGCCTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTCCGCCTTCAACAAGCCACCCGGTTCGTTGCCTCCAACCCCACGGACCgggg                                                       <  1:3905015/150‑1 (MQ=255)
                                                                                                cgcCTGAAGTGACCGTAGGTGTGGTGGATGTCGTGCTGCCTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTCCGCCTTCAACAAGCCACCCGGTTCGTTGCCTCCAACCCCACGGACCGGGGct                                                     >  2:3008807/1‑150 (MQ=255)
                                                                                                cgcCTGAAGTGACCGTAGGTGTGGTGGATGTCGTGCTGCCTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTCCGCCTTCAACAAGCCACCCGGTTCGTTGCCTCCAACCCCACGGACCGGGGct                                                     >  2:3825596/1‑150 (MQ=255)
                                                                                                    tGAAGTGACCGTAGGTGTGGTGGATGTCGTGCTGCCTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTCCGCCTTCAACAAGCCACCCGGTTCGTTGCCTCCAACCCCACGGACCGGGGCTCTCa                                                 <  1:1114653/150‑1 (MQ=255)
                                                                                                    tGAAGTGACCGTAGGTGTGGTGGATGTCGTGCTGCCTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTCCGCCTTCAACAAGCCACCCGGTTCGTTGCCTCCAACCCCACGGACCGGGGCTCTCa                                                 <  1:3527129/150‑1 (MQ=255)
                                                                                                     gAAGTGACCGTAGGTGTGGTGGATGTCGTGCTGCCTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTCCGCCTTCAACAAGCCACCCGGTTCGTTGCCTCCAACCCCACGGACCGGGGCTCTCAt                                                <  2:2183884/150‑1 (MQ=255)
                                                                                                     gAAGTGACCGTAGGTGTGGTGGATGTCGTGCTGCCTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTCCGCCTTCAACAAGCCACCCGGTTCGTTGCCTCCAACCCCACGGACCGGGGCTCTCAt                                                >  1:3798675/1‑150 (MQ=255)
                                                                                                       aGTGACCGTAGGTGTGGTGGATGTCGTGCTGCCTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTCCGCCTTCAACAAGCCACCCGGTTCGTTGCCTCCAACCCCACGGACCGGGGCTCTCATGa                                              <  1:3732432/150‑1 (MQ=255)
                                                                                                        gTGACCGTAGGTGTGGTGGATGTCGTGCTGCCTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTCCGCCTTCAACAAGCCACCCGGTTCGTTGCCTCCAACCCCACGGACCGGGGCTCTCATGAc                                             <  2:847516/150‑1 (MQ=255)
                                                                                                          gACCGTAGGTGTGGTGGATGTCGTGCTGCCTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTCCGCCTTCAACAAGCCACCCGGTTCGTTGCCTCCAACCCCACGGACCGGGGCTCTCATGACCt                                           <  1:4505053/150‑1 (MQ=255)
                                                                                                           aCCGTAGGTGTGGTGGATGTCGTGCTGCCTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTCCGCCTTCAACAAGCCACCCGGTTCGTTGCCTCCAACCCCACGGACCGGGGCTCTCATGACCTg                                          <  2:3966037/150‑1 (MQ=255)
                                                                                                            ccGTAGGTGTGGTGGATGTCGTGCTGCCTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTCCGCCTTCAACAAGCCACCCGGTTCGTTGCCTCCAACCCCACGGACCGGGGCTCTCATGACCTgg                                         <  2:5289983/150‑1 (MQ=255)
                                                                                                            ccGTAGGTGTGGTGGATGTCGTGCTGCCTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTCCGCCTTCAACAAGCCACCCGGTTCGTTGCCTCCAACCCCACGGACCGGGGCTCTCATGACCTgg                                         <  2:5043910/150‑1 (MQ=255)
                                                                                                              gTAGGTGTGGTGGATGTCGTGCTGCCTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTCCGCCTTCAACAAGCCACCCGGTTCGTTGCCTCCAACCCCACGGACCGGGGCTCTCATGACCTGGGa                                       <  2:364685/150‑1 (MQ=255)
                                                                                                                aGGTGTGGTGGATGTCGTGCTGCCTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTCCGCCTTCAACAAGCCACCCGGTTCGTTGCCTCCAACCCCACGGACCGGGGCTCTCATGACCTGGGAgg                                     <  2:1429873/150‑1 (MQ=255)
                                                                                                                  gtgtGGTGGATGTCGTGCTGCCTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTCCGCCTTCAACAAGCCACCCGGTTCGTTGCCTCCAACCCCACGGACCGGGGCTCTCATGACCTGGGAGGTc                                   >  1:721452/1‑150 (MQ=255)
                                                                                                                   tgtgGTGGATGTCGTGCTGCCTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTCCGCCTTCAACAAGCCACCCGGTACGTTGCCTCCAACCCCACGGACCGGGGCTCTCATGACCTGGGAGGTCg                                  >  2:2500302/1‑150 (MQ=255)
                                                                                                                        tggATGTCGTGCTGCCTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTCCGCCTTCAACAAGCCACCCGGTTCGTTGCCTCCAACCCCACGGACCGGGGCTCTCATGACCTGGGAGGTCGAATgg                             >  2:4337011/1‑150 (MQ=255)
                                                                                                                         ggATGTCGTGCTGCCTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTCCGCCTTCAACAAGCCACCCGGTTCGTTGCCTCCAACCCCACGGACCGGGGCTCTCATGACCTGGGAGGTCGAATGGa                            <  1:1655004/150‑1 (MQ=255)
                                                                                                                          gATGTCGTGCTGCCTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTCCGCCTTCAACAAGCCACCCGGTTCGTTGCCTCCAACCCCACGGACCGGGGCTCTCATGACCTGGGAGGTCGAATGGAc                           >  2:695353/1‑150 (MQ=255)
                                                                                                                               cGTGCTGCCTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTCCGCCTTCAACAAGCCACCCGGTTCGTTGCCTCCAACCCCACGGACCGGGGCTCTCATGACCTGGGAGGTCGAATGGACGTCGc                      <  2:2970315/150‑1 (MQ=255)
                                                                                                                                gTGCTGCCTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTCCGCCTTCAACAAGCCACCCGGTTCGTTGCCTCCAACCCCACGGACCGGGGCTCTCATGACCGGGGAGGTCGAATGGACGTCGCt                     >  2:3556129/1‑150 (MQ=255)
                                                                                                                                   ctgcCTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTCCGCCTTCAACAAGCCACCCGGTTCGTTGCCTCCAACCCCACGGACCGGGGCTCTCATGACCTGGGAGGTCGAATGGACGTCGCTACa                  >  2:1201779/1‑150 (MQ=255)
                                                                                                                                    tgcCTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTCCGCCTTCAACAAGCCACCCGGTTCGTTGCCTCCAACCCCACGGAGCGGGGCTCTCATGACCTGGGAGGTCGAATGGACGTCGCTACaa                 >  2:359314/1‑150 (MQ=255)
                                                                                                                                     ggctttccatccatGCTGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTCCGCCTTCAACAAGCCACCCGGTTCGTTGCCTCCAACCCCACGGACCGGGGCTCTCATGACCTGGGAGGTCGAATGGACGTCGCTACAAc                >  1:2981033/7‑150 (MQ=255)
                                                                                                                                     gcCTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTCCGCCTTCAACAAGCCACCCGGTTCGTTGCCTCCAACCCCACGGACCGGGGCTCTCATGACCTGGGAGGTCGAAt                               <  1:4148685/135‑1 (MQ=255)
                                                                                                                                     gcCTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTCCGCCTTCAACAAGCCACCCGGTTCGTTGCCTCCAACCCCACGGACCGGGGCTCTCATGACCTGGGAGGTCGAATGGACGTCGCTACAAc                >  2:1046161/1‑150 (MQ=255)
                                                                                                                                     gcCTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACGGCATAGTACGGTCCGCCTTCAACAAGCCACCCGGTTCGTTGCCTCCAACCCCACGGACCGGGGCTCTCATGACCTGGGAGGTCGAAt                               >  2:4148685/1‑135 (MQ=255)
                                                                                                                                        tttccatccatGCTGGATGGGAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTCCGCCTTCAACAAGCCACCCGGTTCGTTGCCTCCAACCCCACGGACCGGGGCTCTCATGACCTGGGAGGTCGAATGGACGTCGCTACAACGAt             >  2:2317500/4‑150 (MQ=255)
                                                                                                                                        tACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTCCGCCTTCAACAAGCCACCCGGTTCGTTGCCTCCAACCCCACGGACCGGGGCTCTCATGACCTGGGAGGTCGAATGGACGTCGCTACAACGAt             >  2:1301117/1‑150 (MQ=255)
                                                                                                                                          cccATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTCCGCCTTCAACAAGCCACCCGGTTCGTTGCCTCCAACCCCACGGACCGGGGCTCTCATGACCTGGGAGGTCGAATGGACGTCGCTACAACGATCa           <  1:958572/150‑1 (MQ=255)
                                                                                                                                             aTCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTCCGCCTTCAACAAGCCACCCGGTTCGTTGCCTCCAACCCCACGGACCGGGGCTCTCATGACCTGGGAGGTCGAATGGACGTCGCTACAACGATCAGTc        >  2:644451/1‑150 (MQ=255)
                                                                                                                                                cAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTCCGCCTTCAACAAGCCACCCGGTTCGTTGCCTCCAACCCCACGGACCGGGGCTCTCATGACCTGGGAGGTCGAATGGACGTCGCTACAACGATCAGTcaac     >  2:4409590/1‑150 (MQ=255)
                                                                                                                                                cAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTCCGCCTTCAACAAGCCACCCGGTTCGTTGCCTCCAACCCCACGGACCGGGGCTCTCATGACCTGGGAGGTCGAATGGACGTCGCTACAACGATCAGTcaac     >  1:4666762/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
GCGTGCGGCGATGAGCATGGTCTGAGCAGGGCGCAGGGGTCGGATCTCCGAGCGCGTCTACCCGCGAAGCAGGTGGCGCGGGGCATGCGGCGACGGCGCCTGAAGTGACCGTAGGTGTGGTGGATGTCGTGCTGCCTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTCCGCCTTCAACAAGCCACCCGGTTCGTTGCCTCCAACCCCACGGACCGGGGCTCTCATGACCTGGGAGGTCGAATGGACGTCGCTACAACGATCAGTCAACAAG  >  NC_002947/3403748‑3404044

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