Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 95,665 A→C 10.8% intergenic (+57/+15) PP_0087 → / ← PP_0088 hypothetical protein/alkanal monooxygenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00294795,6650AC10.8% 280.3 / 5.8 96intergenic (+57/+15)PP_0087/PP_0088hypothetical protein/alkanal monooxygenase
Reads supporting (aligned to +/- strand):  ref base A (67/15);  new base C (2/11);  total (69/27)
Fisher's exact test for biased strand distribution p-value = 5.39e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTGGAGATCGTCAAGAAAGGCGAGCCGGCGGGCAAGGTGATTGGTAGCTGCGAGGCCGGTACCAAGGAGATTGTTTACCGCCGCGGCTAATGTGTTTGCTCATTCGCGGGCTTGCTCGCTCCCAAGTCAACCTATGGGAGCGGGCGAGCCCGCGACGAGGGTCACTCGGTCTTGAGGGCCTGCGCCATCAGCTCATATGACCGAATGCGGTCCGCATGCTCATACAGGTCACAGGTAAAGATCAGCTCATCCGCACCCGTCTGTTCCAGCAGCACTTCTACCTTGGCCTTTAC  >  NC_002947/95519‑95811
                                                                                                                                                  |                                                                                                                                                  
cTGGAGATCGTCAAGAAAGGCGAGCCGGCGGGCAAGGTGATTGGTAGCTGCGAGGCCGGTACCAAGGAGATTGTTTACCGCCGCGGCTAATGTGTTTGCTCATTCGCGGGCTTGCTCGCTCCCAAGTCAACCTATGGGAGCGGGCGAGcc                                                                                                                                                 >  1:5224737/1‑150 (MQ=255)
cTGGAGATCGTCAAGAAAGGCGAGCCGGCGGGCAAGGTGATTGGTAGCTGCGAGGCCGGTACCAAGGAGATTGTTTACCGCCGCGGCTAATGTGTTTGCTCATTCGCGGGCTTGCTCGCTCCCAAGTCAACCTATGGGAGCGGGCGAGcc                                                                                                                                                 >  1:5224767/1‑150 (MQ=255)
 tGGAGATCGTCAAGAAAGGCGAGCCGGCGGGCAAGGTGATTGGTAGCTGCGAGGCCGGTACCAAGGAGATTGTTTACCGCCGCGGCTAATGTGTTTGCTCATTCGCGGGCTTGCTCGCTCCCAAGTCAACCTATGGGAGCGGGCGAGccc                                                                                                                                                <  1:4462715/150‑1 (MQ=255)
 tGGAGATCGTCAAGAAAGGCGAGCCGGCGGGCAAGGTGATTGGTAGCTGCGAGGCCGGTACCAAGGAGATTGTTTACCGCCGCGGCTAATGTGTTTGCTCATTCGCGGGCTTGCTCGCTCCCAAGTCAACCTATGGGAGCGGGCGAGccc                                                                                                                                                <  2:1092960/150‑1 (MQ=255)
 tGGAGATCGTCAAGAAAGGCGAGCCGGCGGGCAAGGTGATTGGTAGCTGCGAGGCCGGTACCAAGGAGATTGTTTACCGCCGCGGCTAATGTGTTTGCTCATTCGCGGGCTTGCTCGCTCCCAAGTCAACCTATGGGAGCGGGCGAGccc                                                                                                                                                <  1:5419437/150‑1 (MQ=255)
  ggAGATCGTCAAGAAAGGCGAGCCGGCGGGCAAGGTGATTGGTAGCTGCGAGGCCGGTACCAAGGAGATTGTTTACCGCCGCGGCTAATGTGTTTGCTCATTCGCGGGCTTGCTCGCTCCCAAGTCAACCTATGGGAGCGTGCGAGCCcg                                                                                                                                               <  2:4998700/150‑1 (MQ=255)
  ggAGATCGTCAAGAAAGGCGAGCCGGCGGGCAAGGTGATTGGTAGCTGCGAGGCCGGTACCAAGGAGATTGTTTACCGCCGCGGCTAATGTGTTTGCTCATTCGCGGGCTTGCTCGCTCCCAAGTCAACCTATGGGAGCGGGCGAGCCcg                                                                                                                                               >  1:4985054/1‑150 (MQ=255)
  ggAGATCGTCAAGAAAGGCGAGCCGGCGGGCAAGGTGATTGGTAGCTGCGAGGCCGGTACCAAGGAGATTGTTTACCGCCGCGGCTAATGTGTTTGCTCATTCGCGGGCTTGCTCGCTCCCAAGTCAACCTATGGGAGCGGGCGAGCCcg                                                                                                                                               >  2:5177173/1‑150 (MQ=255)
  ggAGATCGTCAAGAAAGGCGAGCCGGCGGGCAAGGTGATTGGTAGCTGCGAGGCCGGTACCAAGGAGAGTGTTTACCGCCGCGGCTAATGTGTTTGCTCATTCGCGGGCTTGCTCGCTCCCAAGTCAACCTATGGGAGCGGGCGAGCCcg                                                                                                                                               <  2:5092510/150‑1 (MQ=255)
         gTCAAGAAAGGCGAGCCGGCGGGCAAGGTGATTGGTAGCTGCGAGGCCGGTACCAAGGAGATTGTTTACCGCCGCGGCTAATGTGTTTGCTCATTCGCGGGCTTGCTCGCTCCCAAGTCAACCTATGGGAGCGGGCGAGCCCGCGACGAg                                                                                                                                        >  1:1254458/1‑150 (MQ=255)
          tCAAGAAAGGCGAGCCGGCGGGCAAGGTGATTGGTAGCTGCGAGGCCGGTACCAAGGAGATTGTTTACCGCCGCGGCTAATGTGTTTGCTCATTCGCGGGCTTGCTCGCTCCCAAGTCAACCTATGGGAGCGGGCGAGCCCGCGACGAgg                                                                                                                                       >  1:635965/1‑150 (MQ=255)
               aaaGGCGAGCCGGCGGGCAAGGTGATTGGTAGCTGCGAGGCCGGTACCAAGGAGATTGTTTACCGCCGCGGCTAATGTGTTTGCTCATTCGCGGGCTTGCTCGCTCCCAAGTCAACCTATGGGAGCGGGCGAGCCCGTGACGAGGGTCAc                                                                                                                                  <  1:657301/150‑1 (MQ=255)
                aaGGCGAGCCGGCGGGCAAGGTGATTGGTAGCTGCGAGGCCGGTACCAAGGAGATTGTTTACCGCCGCGGCTAATGTGTTTGCTCATTCGCGGGCTTGCTCGCTCCCAAGTCAACCTATGGGAGCGGGCGAGCCCGCGACGAGGGTCACt                                                                                                                                 >  2:653012/1‑150 (MQ=255)
                                                tGCGAGGCCGGTACCAAGGAGATTGTTTACCGCCGCGGCTAATGTGTTTGCTCATTCGCGGGCTTGCTCGCTCCCAAGTCAACCTATGGGAGCGGGCGAGCCCGCGACGAGGGTCACTCGGTCTTGAGGGCCTGCGCCATCAGCTCatat                                                                                                 >  2:558596/1‑150 (MQ=255)
                                                  gggggCGGGTACCACGTTGTTTGTTTCCCGCCGCGGCTATTGTTTTTGCTCTTGCGGGGGATTGCGCGCACCCACGAAACCCTCGGGGCGGGGGCGCGCCCGCGACGAGGGTCACTCGGTCTTGAGGGCCTGCGCCATCAGCTCATATGa                                                                                               <  2:2234770/147‑1 (MQ=255)
                                                   gAGGCCGGTACCAAGGAGATTGTTTACCGCCGCGGCTAATGTGTTTGCTCATTCGCGGGCTTGCTCGCTCCCAAGTCAACCTATGGGAGCGGGCGAGCCCGCGACGAGGGTCACTCGGTCTTGAGGGCCTGCGCCATCAGCTCATATGAc                                                                                              >  2:344905/1‑150 (MQ=255)
                                                     ggcggttACCATGGATTTTTTTTCCCGCCGGGGCTATTGTGTTTGCTCTTGCGCGGGCTTGCTCGCTCCCCCGTCAACCTCTGGGCGGGGGGGCGCCCGCGACGAGGGTCACTCGGTCTTGAGGGCCTGCGCCATCAGCTCATATGACCg                                                                                            <  2:4448237/144‑1 (MQ=255)
                                                     ggCCGGTACCAAGGAGATTGTTTACCGCCGCGGCTAATGTGTTTGCTCATTCGCGGGCTTGCTCGCTCCCAAGTCAACCTATGGGAGCGGGCGAGCCCGCGACGAGGGTCACTCGGTCTTGAGGGCCTGCGCCATCAGCTCATAAGATCg                                                                                            >  1:1764314/1‑150 (MQ=255)
                                                               aaGGAGATTGTTTACCGCCGCGGCTAATGTGTTTGCTCATTCGCGGGCTTGCTCGCTCCCAAGTCAACCTATGGGAGCGGGCGAGCCCGCGACGAGGGTCACTCGGTCTTGAGGGCCTGCGCCATCAGCTCATATGACCGAATGCGGTcc                                                                                  >  1:1077571/1‑150 (MQ=255)
                                                               aaGGAGATTGTTTACCGCCGCGGCTAATGTGTTTGCTCATTCGCGGGCTTGCTCGCTCCCAAGTCAACCTATGGGAGCGGGCGAGCCCGCGACGAGGGTCACTCGGTCTTGAGGGCCTGCGCCATCAGCTCATATGACCGAATGCGGTcc                                                                                  >  2:3415848/1‑150 (MQ=255)
                                                                  gagaTTGTTTACCGCCGCGGCTAATGTGTTTGCTCATTCGCGGGCTTGCTCGCTCCCAAGTCAACCTATGGGAGCGGGCGAGCCCGCGACGAGGGTCACTCGGTCTTGAGGGCCTGCGCCATCAGCTCATATGACCGAATGCGGTCCGCa                                                                               >  2:2909455/1‑150 (MQ=255)
                                                                  gagaTTGTTTACCGCCGCGGCTAATGTGTTTGCTCATTCGCGGGCTTGCTCGCTCCCAAGTCAACCTATGGGAGCGGGCGAGCCCGCGACGAGGGTCACTCGGTCTTGAGGGCCTGCGCCATCAGCTCATATGACCGAATGCGGTCCGCa                                                                               >  2:1453855/1‑150 (MQ=255)
                                                                     aTTGTTTACCGCCGCGGCTAATGTGTTTGCTCATTCGCGGGCTTGCTCGCTCCCAAGTCAACCTATGGGAGCGGGCGAGCCCGCGACGAGGGTCACTCGGTCTTGAGGGCCTGCGCCATCAGCTCATATGACCGAATGCGGTCCGCATGc                                                                            >  2:5223524/1‑150 (MQ=255)
                                                                      ttGTTTACCGCCGCGGCTAATGTGTTTGCTCATTCGCGGGCTTGCTCGCTCCCAAGTCAACCTATGGGAGCGGGCGAGCCCGCGACGAGGGTCACTCGGTCTTGAGGGCCTGCGCCATCAGCTCATATGACCGAATGCGGTCCGCATGCt                                                                           >  1:4141357/1‑150 (MQ=255)
                                                                      ttGTTTACCGCCGCGGCTAATGTGTTTGCTCATTCGCGGGCTTGCTCGCTCCCAAGTCAACCTATGGGAGCGGGCGAGCCCGCGACGAGGGTCACTCGGTCTTGAGGGCCTGCGCCATCAGCTCATATGACCGAATGCGGTCCGCATGCt                                                                           >  2:3721752/1‑150 (MQ=255)
                                                                         tgtACCGCCGCGGCTAATGTGTTTGCTCATTCGCGGGCTTGCTCGCTCCCAAGTCAACCTATGGGAGCGGGCGAGCCCGCGACGAGGGTCACTCGGTCTTGAGGGCCTGCGCCATCAGCTCATATGACCGAATGCGGTCCGCATGCTCAt                                                                        >  1:2640511/3‑150 (MQ=255)
                                                                            aCCGCCGCGGCTAATGTGTTTGCTCATTCGCGGGCTTGCTCGCTCCCAAGTCAACCTATGGGAGCGGGCGAGCCCGCGACGAGGGTCACTCGGTCTTGAGGGCCTGCGCCATCAGCTCATATGACCGAATGCGGTCCGCATGCTCATACa                                                                     >  1:2666295/1‑150 (MQ=255)
                                                                                  gcgGCTAATGTGTTTGCTCATTCGCGGGCTTGCTCGCTCCCAAGTCAACCTATGGGAGCGGGCGAGCCCGCGACGAGGGTCACTCGGTCTTGAGGGCCTGCGCCATCAGCTCATATGACCGAATGCGGTCCGCATGCTCATACAGGTcac                                                               >  2:5382247/1‑150 (MQ=255)
                                                                                   gggttattgttTTTGCTCTTTGGGGGGGTTGCTCGATCCAACGTCCCCGTTTGGTTGGGGGGGAGCCCGCGACGAGGGTCACTCGGTCTTGAGGGCCTGCGCCATCAGCTCATATGACCGAATGCGGTCCGCATGCTCATACAGGTcaca                                                              <  2:391222/143‑1 (MQ=255)
                                                                                   cgGCTAATGTGTTTGCTCATTCGCGGGCTTGCTCGCTCCCAAGTCAACCTATGGGAGCGGGCGAGCCCGCGACGAGGGTCACTCGGTCTTGAGGGCCTGCGCCATCAGCTCATATGACCGAATGCGGTCCGCATGCTCATACAGGTcaca                                                              >  1:3298702/1‑150 (MQ=255)
                                                                                    ggCTATTGTTTTTGCTCTTTCGGGGGGTTGCTCGATCCAAATTCAACCTTTGGGAGCGGGCGCGCCCGCGACGAGGGTCACTCGGTCTTGAGGGCCTGCGCCATCAGCTCATATGACCGAATGCGGTCCGCATGCTCATACAGGTCACAg                                                             <  2:4021706/150‑1 (MQ=255)
                                                                                     gCTAATGTGTTTGCTCATTCGCGGGCTTGCTCGCTCCCAAGTCAACCTATGGGAGCGGGCGAGCCCGCGACGAGGGTCACTCGGTCTTGAGGGCCTGCGCCATCAGCTCATATGACCGAATGCGGTCCGCATGCTCATACAGGTCACAgg                                                            >  2:4989621/1‑150 (MQ=255)
                                                                                        aaTGTGTTTGCTCATTCGCGGGCTTGCTCGCTCCCAAGTCAACCTATGGGAGCGGGCGAGCCCGCGACGAGGGTCACTCGGTCTTGAGGGCCTGCGCCATCAGCTCATATGACCGAATGCGGTCCGCATGCTCATACAGGTCACAGGTaa                                                         >  1:550214/1‑150 (MQ=255)
                                                                                         ttttgtTTGCTCATTCGCGGGTTTGCTCGCTCCCACGACAACCTATGGGAGGGGGCGCGCCCGCGACGAGGGTCACTCGGTCTTGAGGGCCTGCGCCATCAGCTCATATGACCGAATGCGGTCCGCATGCTCATACAGGTCACAGGTaaa                                                        <  2:1789671/147‑1 (MQ=255)
                                                                                          tgttTTTGCTCTTTCGCGGGCTTGCTCGCGCCCTCGTCAACCTTGGGGCGCGGGCGCGCCCGCGACAAGGGTCACTCGGTCTTGAGGGCCTGCGCCATCAGCTCATATGACCGAATGCGGTCCGCATGCTCATACAGGTCACAGGTAAAg                                                       <  2:2629331/150‑1 (MQ=255)
                                                                                          tgtgTTGGCCCATGCGGGGGTTTGCTGGCTCCCAAGTCACCCTATGGGCGGGGGCGCGCCCGCGACAAGGGTCACTCGGTCTTGAGGGCCTGCGCCATCAGCTCATATGACCGAATGCGGTCCGCATGCTCATACAGGTCACAGGTAAAg                                                       <  2:2629327/150‑1 (MQ=255)
                                                                                             gtTTGCTCTTGCGGGGGCTGGCCCGCCCCAACGTCAACCTCGGGGCGGGGGCGCGCCCGCGACGAGGGTCACTCGGTCTTGAGGGCCTGCGCCATCAGCTCATATGACCGAATGCGGTCCGCATGCTCATACAGGTCACAGGTAAAGATc                                                    <  2:2666295/150‑1 (MQ=255)
                                                                                              tttGCTCATTCGCGGGCTTGCTCGCTCCCAAGTCAACCTATGGGAGCGGGCGAGCCCGCGACGAGGGTCACTCGGTCTTGAGGGCCTGCGCCATCAGCTCATATGACCGAATGCGGTCCGCATGCTCATACAGGTCACAGGTCAAGATCa                                                   >  1:3883333/1‑150 (MQ=255)
                                                                                              tttGCTCATTCGCGGGCTTGCTCGCTCCCAAGTCAACCTATGGGAGCGGGCGAGCCCGCGACGAGGGTCACTCGGTCTTGAGGGCCTGCGCCATCAGCTCATATGACCGAATGCGGTCCGCATGCTCATACAGGTCACAGGTAATGATCa                                                   >  1:383600/1‑150 (MQ=255)
                                                                                              tttGCTCATTCGCGGGCTTGCTCGCTCCCAAGTCAACCTATGGGAGCGGGCGAGCCCGCGACGAGGGTCACTCGGTCTTGAGGGCCTGCGCCATCAGCTCATATGACCGAATGCGGTCCGCATGCTCATACAGGTCACAGGTAAAGATCa                                                   >  1:1495949/1‑150 (MQ=255)
                                                                                                   tCATTCGCGGGCTTGCTCGCTCCCAAGTCAACCTATGGGAGCGGGCGAGCCCGCGACGAGGGTCACTCGGTCTTGAGGGCCTGCGCCATCAGCTCATATGACCGAATGCGGTCCGCATGCTCATACAGGTCACAGGTAAAGATCAGCTCa                                              >  1:1640730/1‑150 (MQ=255)
                                                                                                   tCATTCGCGGGCTTGCTCGCTCCCAAGTCAACCTATGGGAGCGGGCGAGCCCGCGACGAGGGTCACTCGGTCTTGAGGGCCTGCGCCATCAGCTCATATGACCGAATGCGGTCCGCATGCTCATACAGGTCACAGGTAAAGATCAGCTCa                                              >  2:2777406/1‑150 (MQ=255)
                                                                                                   tCATTCGCGGGCTTGCTCGCTCCCAAGTCAACCTATGGGAGCGGGCGAGCCCGCGACGAGGGTCACTCGGTCTTGAGGGCCTGCGCCATCAGCTCATATGACCGAATGCGGTCCGCATGCTCATACAGGTCACAGGTAAAGATCAGCTCa                                              >  1:337757/1‑150 (MQ=255)
                                                                                                   tCATTCGCGGGCTTGCTCGCTCCCAAGTCAACCTATGGGAGCGGGCGAGCCCGCGACGAGGGTCACTCGGTCTTGAGGGCCTGCGCCATCAGCTCATATGACCGAATGCGGTCCGCATGCTCATACAGGTCACAGGTAAAGATCAGCTCa                                              >  1:736901/1‑150 (MQ=255)
                                                                                                    cATTCGCGGGCTTGCTCGCTCCCAAGTCAACCTATGGGAGCGGGCGAGCCCGCGACGAGGGTCACTCGGTCTTGAGGGCCTGCGCCATCAGCTCATATGACCGAATGCGGTCCGCATGCTCATACAGGTCACAGGTAAAGATCAGCTCAt                                             >  2:748947/1‑150 (MQ=255)
                                                                                                    cATTCGCGGGCTTGCTCGCTCCCAAGTCAACCTATGGGAGCGGGCGAGCCCGCGACGAGGGTCACTCGGTCTTGAGGGCCTGCGCCATCAGCTCATATGACCGAATGCGGTCCGCATGCTCATACAGGTCACAGGTAAAGATCAGCTCAt                                             >  1:3249082/1‑150 (MQ=255)
                                                                                                      ttCGCGGGCTTGCTCGCTCCCAAGTCAACCTATGGGAGCGGGCGAGCCCGCGACGAGGGTCACTCGGTCTTGAGGGCCTGCGCCATCAGCTCATATGACCGAATGCGGTCCGCATGCTCATACAGGTCACAGGTAAAGATCAGCTCATcc                                           >  2:3287324/1‑150 (MQ=255)
                                                                                                           gggCTTGCTCGCTCCCAAGTCAACCTATGGGAGCGGGCGAGCCCGCGACGAGGGTCACTCGGTCTTGAGGGCCTGCGCCATCAGCTCATATGACCGAATGCGGTCCGCATGCTCATACAGGTCACAGGTAAAGATCAGCTCATCCGCAcc                                      >  2:3675473/1‑150 (MQ=255)
                                                                                                           gggCTCGCCCGCTCCCAAGTCAACCTATGGGAGCGGGCGAGCCCGCGACGAGGGTCACTCGGTGTTGAGGGCCTGCGCCATCAGCTCATATGACCGAATGAGGTCCGCATGCTCATACAGGTCACAGGTAAAGATCAGCTCATCCGCAcc                                      >  1:111382/1‑150 (MQ=255)
                                                                                                               ttGCTCGCTCCCAAGTCAACCTATGGGAGCGGGCGAGCCCGCGACGAGGGTCACTCGGTCTTGAGGGCCTGCGCCATCAGCTCATATGACCGAATGCGGTCCGCATGCTCATACAGGTCACAGGTAAAGATCAGCTCATCCGCACCCGTc                                  >  1:2601440/1‑150 (MQ=255)
                                                                                                               ttGCTCGCTCCCAAGTCAACCTATGGGAGCGGGCGAGCCCGCGACGAGGGTCACTCGGTCTTGAGGGCCTGCGCCATCAGCTCATATGACCGAATGCGGTCCGCATGCTCATACAGGTCACAGGTAAAGATCAGCTCATCCGCACCCGTc                                  >  2:2834971/1‑150 (MQ=255)
                                                                                                               ttGCTCGCTCCCAAGTCAACCTATGGGAGCGGGCGAGCCCGCGACGAGGGTCACTCGGTCTTGAGGGCCTGCGCCATCAGCTCATATGACCGAATGCGGTCCGCATGCTCATACAGGTCACAGGTAAAGATCAGCTCATCCGCACCCGTc                                  >  2:2834954/1‑150 (MQ=255)
                                                                                                                 gctcgctcCCAAGTCAACCTATGGGAGCGGGCGAGCCCGCGACGAGGGTCACTCGGTCTTGAGGGCCTGCGCCATCAGCTCATATGACCGAATGCGGTCCGCATGCTCATACAGGTCACAGGTAAAGATCAGCTCATCCGCACCCGTCTg                                >  2:5146302/1‑150 (MQ=255)
                                                                                                                 gctcgctcCCAAGTCAACCTATGGGAGCGGGCGAGCCCGCGACGAGGGTCACTCGGTCTTGAGGGCCTGCGCCATCAGCTCATATGACCGAATGCGGTCCGCATGCTCATACAGGTCACAGGTAAAGATCAGCTCATCCGCACCCGTCTg                                >  2:4383395/1‑150 (MQ=255)
                                                                                                                 gctcgctcCCAAGTCAACCTATGGGAGCGGGCGAGCCCGCGACGAGGGTCACTCGGTCTTGAGGGCCTGCGCCATCAGCTCATATGACCGAATGCGGTCCGCATGCTCATACAGGTCACAGGTAAAGATCAGCTCATCCGCACCCGTCTg                                >  2:879413/1‑150 (MQ=255)
                                                                                                                   tcgctcCCAAGTCAACCTATGGGAGCGGGCGAGCCCGCGACGAGGGTCACTCGGTCTTGAGGGCCTGCGCCATCAGCTCATATGACCGAATGCGGTCCGCATGCTCATACAGGTCACAg                                                             >  2:4146499/1‑119 (MQ=255)
                                                                                                                   tcgctcCCAAGTCAACCTATGGGAGCGGGCGAGCCCGCGACGAGGGTCACTCGGTCTTGAGGGCCTGCGCCATCAGCTCATATGACCGAATGCGGTCCGCATGCTCATACAGGTCACAg                                                             <  1:4146499/119‑1 (MQ=255)
                                                                                                                   tcgctcCCAAGTCAACCTATGGGAGCGGGCGAGCCCGCGACGAGGGTCACTCGGTCTTGAGGGCCTGCGCCATCAGCTCATATGACCGAATGCGGTCCGCATGCTCATACAGGTCACAGGTACAGATCAGCTCATCCGCACCCGTCTGtt                              >  2:2837187/1‑150 (MQ=255)
                                                                                                                   tcgctcCCAAGTCAACCTATGGGAGCGGGCGAGCCCGCGACGAGGGTCACTCGGTCTTGAGGGCCTGCGCCATCAGCTCATATGACCGAATGCGGTCCGCATGCTCATACAGGTCACAGGTAAAGATCAGCTCATCCGCACCCGTCTGtt                              >  1:2886884/1‑150 (MQ=255)
                                                                                                                   tcgctcCCAAGTCAACCTATGGGAGCGGGCGAGCCCGCGACGAGGGTCACTCGGTCTTGAGGGCCTGCGCCATCAGCTCATATGACCGAATGCGGTCCGCATGCTCATACAGGTCACAGGTAAAGATCAGCTCATCCGCACCCGTCTGtt                              >  2:5265814/1‑150 (MQ=255)
                                                                                                                   tcgctcCCAAGTCAACCTATGGGAGCGGGCGAGCCCGCGACGAGGGTCACTCGGTCTTGAGGGCCTGCGCCATCAGCTCATATGACCGAATGCGGTCCGCATGCTCATACAGGTCACAGGTAAAGATCAGCTCATCCGCACCCGTCTGtt                              >  1:1229314/1‑150 (MQ=255)
                                                                                                                   tcgctcCCAAGTCAACCTATGGGAGCGGGCGAGCCCGCGACGAGGGTCACTCGGTCTTGAGGGCCTGCGCCATCAGCTCATATGACCGAATGCGGTCCGCATGCTCATACAGGTCACAGGTAAAGATCAGCTCATCCGCACCCGTCTGtt                              >  2:2539192/1‑150 (MQ=255)
                                                                                                                   tcgctcCCAAGTCAACCTATGGGAGCGGGCGAGCCCGCGACGAGGGTCACTCGGTCTTGAGGGCCTGCGCCATCAGCTCATATGACCGAATGCGGTCCGCATGCTCATACAGGTCACAGGTAAAGATCAGCTCATCCGCACCCGTCTGtt                              >  2:3617638/1‑150 (MQ=255)
                                                                                                                   tcgctcCCAAGTCAACCTATGGGAGCGGGCGAGCCCGCGACGAGGGTCACTCGGTCTTGAGGGCCTGCGCCATCAGCTCATATGACCGAATGCGGTCCGCATGCTCATACAGGTCACAGGTAAAGATCAGCTCATCCGCACCCGTCTGtt                              >  2:3617549/1‑150 (MQ=255)
                                                                                                                   tcgctcCCAAGTCAACCTATGGGAGCGGGCGAGCCCGCGACGAGGGTCACTCGGTCTTGAGGGCCTGCGCCATCAGCTCATATGACCGAATGCGGTCCGCATGCTCATACAGGTCACAGGTAAAGATCAGCTCATCCGCACCCGTCTGtt                              >  1:3993495/1‑150 (MQ=255)
                                                                                                                   tcgctcCCAAGTCAACCTATGGGAGCGGGCGAGCCCGCGACGAGGGTCACTCGGTCTTGAGGGCCTGCGCCATCAGCTCATATGACCGAATGCGGTCCGCATGCTCATACAGGTCACAGGTAAAGATCAGCTCATCCGCACCCGTCTGtt                              >  1:887033/1‑150 (MQ=255)
                                                                                                                   tcgctcCCAAGTCAACCTATGGGAGCGGGCGAGCCCGCGACGAGGGTCACTCGGTCTTGAGGGCCTGCGCCATCAGCTCATATGACCGAATGCGGTCCGCATGCTCATACAGGTCACAGGTAAAGATCAGCTCATCCGCACCCGTCTGtt                              >  1:282207/1‑150 (MQ=255)
                                                                                                                   tcgctcCCAAGTCAACCTATGGGAGCGGGCGAGCCCGCGACGAGGGTCACTCGGTCTTGAGGGCCTGCGCCATCAGCTCATATGACCGAATGCGGTCCGCATGCTCATACAGGTCACAGGTAAAGATCAGCTCATCCGCACCCGTCTGtt                              >  1:887127/1‑150 (MQ=255)
                                                                                                                   tcgctcCCAAGTCAACCTATGGGAGCGGGCGAGCCCGCGACGAGGGTCACTCGGTCTTGAGGGCCTGCGCCATCAGCTCATATGACCGAATGCGGTCCGCATGCTCATACAGGTCACAGGTAAAGATCAGCTCATCCGCACCCGTCTGtt                              >  2:1470161/1‑150 (MQ=255)
                                                                                                                   tcgctcCCAAGTCAACCTATGGGAGCGGGCGAGCCCGCGACGAGGGTCACTCGGTCTTGAGGGCCTGCGCCATCAGCTCATATGACCGAATGCGGTCCGCATGCTCATACAGGTCACAGGTAAAGATCAGCTCATCCGCACCCGTCTGtt                              >  1:1422212/1‑150 (MQ=255)
                                                                                                                      ctcCCACGACAACCTTGGGGTGCGGGCGCGCCCGCGACGAGGGTCACTCGGTCTTGAGGGCCTGCGCCATCAGCTCATATGACCGAATGCGGTCCGCATGCTCATACAGGTCACAGGTAAAGATCAGCTCATCCGCACCCGTCTGTTCca                           <  2:4765269/150‑1 (MQ=255)
                                                                                                                      ctcCCAAGTCAACCTTTGGGAGCGGGCGAGCCCGCGACGAGGGTCACTCGGTCTTGAGGGCCTGCGCCATCAGCTCATATGACCGAATGCGGTCCGCATGCTCATACAGGTCACAGGTAAAGATCAGCTCATCCGCACCCGTCTGTTCca                           <  1:2834954/150‑1 (MQ=255)
                                                                                                                      ctcCCAAGTCAACCTTTGGGAGCGGGCGAGCCCGCGACGAGGGTCACTCGGTCTTGAGGGCCTGCGCCATCAGCTCATATGACCGAATGCGGTCCGCATGCTCATACAGGTCACAGGTAAAGATCAGCTCATCCGCACCCGTCTGTTCca                           <  1:2834971/150‑1 (MQ=255)
                                                                                                                        cccAAGTCAACCTATGGGAGCGGGCGAGCCCGCGACGAGGGTCACTCGGTCTTGAGGGCCTGCGCCATCAGCTCATATGACCGAATGCGGTCCGCATGCTCATACAGGTCACAGGTAAAGATCAGCTCATCCGCACCCGTCTGTTCcagc                         >  1:1660407/1‑150 (MQ=255)
                                                                                                                         ccAAGTCAACCTATGGGAGCGGGCGAGCCCGCGACGAGGGTCACTCGGTCTTGAGGGCCTGCGCCATCAGCTCATATGACCGAATGCGGTCCGCATGCTCATACAGGTCACAGGTAAAGATCAGCTCATCCGCACCCGTCTGTTCcagca                        >  1:1814512/1‑150 (MQ=255)
                                                                                                                         ccAAGTCAACCTATGGGAGCGGGCGAGCCCGCGACGAGGGTCACTCGGTCTTGAGGGCCTGCGCCATCAGCTCATATGACCGAATGCGGTCCGCATGCTCATACAGGTCACAGGTAAAGATCAGCTCATCCGCACCCGTCTGTTCcagca                        >  1:3619997/1‑150 (MQ=255)
                                                                                                                         ccAAGTCAACCTATGGGAGCGGGCGAGCCCGCGACGAGGGTCACTCGGTCTTGAGGGCCTGCGCCATCAGCTCATATGACCGAATGCGGTCCGCATGCTCATACAGGTCACAGGTAAAGATCAGCTCATCCGCACCCGTCTGTTCcagca                        >  2:557493/1‑150 (MQ=255)
                                                                                                                           aaGTCAACCTATGGGAGCGGGCGAGCCCGCGACGAGGGTCACTCGGTCTTGAGGGCCTGCGCCATCAGCTCATATGACCGAATGCGGTCCGCATGCTCATACAGGTCACAGGTAAAGATCAGCTCATCCGCACCCGTCTGTTCCagcagc                      >  2:234530/1‑150 (MQ=255)
                                                                                                                            acaaaaCCTTGGGGTGCGGCCTCGCCCGCGACGAGGGTCACTCGGTCTTGAGGGCCTGCGCCATCAGCTCATATGACCGAATGCGGTCCGCATGCTCATACAGGTCACAGGTAAAGATCAGCTCATCCGCACCCGTCTGTTCCAgcagca                     <  2:2923048/146‑1 (MQ=255)
                                                                                                                             gTCAACCTATGGGAGCGGGCGAGCCCGCGACGAGGGTCACTCGGTCTTGAGGGCCTGCGCCATCAGCTCATATGACCGAATGCGGTCCGCATGCTCATACAGGTCACAGGTAAAGATCAGCTCATCCGCACCCGTCTGTTCCAGCAGCAc                    <  2:2886884/150‑1 (MQ=255)
                                                                                                                              tCAACCTATGGGAGCGGGCGAGCCCGCGACGAGGGTCACTCGGTCTTGAGGGCCTGCGCCATCa                                                                                                         <  2:921914/64‑1 (MQ=255)
                                                                                                                              tCAACCTATGGGAGCGGGCGAGCCCGCGACGAGGGTCACTCGGTCTTGAGGGCCTGCGCCATCa                                                                                                         >  1:921914/1‑64 (MQ=255)
                                                                                                                                 aCCTATGGGAGCGGGCGAGCCCGCGACGAGGGTCACTCGGTCTTGAGGGCCTGCGCCATCAGCTCATATGACCGAATGCGGTCCGCATGCTCATACAGGTCACAGGTAAAGATCAGCTCATCCGCACCCGTCTGTTCCAGCAGCACTTCt                >  2:3755935/1‑150 (MQ=255)
                                                                                                                                  ccTATGGGAGCGGGCGAGCCCGCGACGAGGGTCACTCGGTCTTGAGGGCCTGCGCCATCAGCTCATATGACCGAATGCGGTCCGCATGCTCATACAGGTCACAGGTAAAGATCAGCTCATCCGCACCCGTCTGTTCCAGCAGCACTTCTa               >  1:3332693/1‑150 (MQ=255)
                                                                                                                                     ttGGGAGCGAGCAAGCCCGCGACGAGGGTCACTCGGTCTTGAGGGCCTGCGCCATCAGCTCATATGACCGAATGCGGTCCGCATGCTCATACAGGTCACAGGTAAAGATCAGCTCATCCGCACCCGTCTGTTCCAGCAGCACTTCTACCt            >  2:3908903/2‑150 (MQ=255)
                                                                                                                                      ggggggggggCGCCCCCGCGACGAGGGTCACTCGGTCTTGAGGGCCTGCGCCATCAGCTCCTATGACCGAATGCGGTCCGCATGCTCATACAGGTCACAGGTAAAGATCAGCTCATCCGCACCCGTCTGTTCCAGCAGCACTTCTACCtt           <  1:1080406/143‑1 (MQ=255)
                                                                                                                                       gggtgggggCGTCCCCGCGACGAGGGTCACTCGGTCTTGAGGGCCTGCGCCATCAGCTCATATGACCGAATGCGGTCCGCATGCTCATACAGGTCACAGGTAAAGATCAGCTCATCCGCACCCGTCTGTTCCAGCAGCACTTCTACCTTg          <  2:3151757/144‑1 (MQ=255)
                                                                                                                                        ggAGCGGGCGAGCCCGCGACGAGGGTCACTCGGTCTTGAGGGCCTGCGCCATCAGCTCATATGACCGAATGCGGTCCGCATGCTCATACAGGTCACAGGTAAAGATCAg                                                  >  2:4431353/1‑109 (MQ=255)
                                                                                                                                        ggAGCGGGCGAGCCCGCGACGAGGGTCACTCGGTCTTGAGGGCCTGCGCCATCAGCTCATATGACCGAATGCGGTCCGCATGCTCATACAGGTCACAGGTAAAGATCAg                                                  <  1:4431353/109‑1 (MQ=255)
                                                                                                                                         tcgCGGGCTCGCCCGCGACGAGGGTCACTCGGTCTTGAGGGCCTGCGCCATCAGCTCATATGACCGAATGCGGTCCGCATGCTCATACAGGTCACAGGTAAAGATCAGCTCATCCGCACCCGTCTGTTCCAGCAGCACTTCTACCTTGGc        >  2:4389538/3‑150 (MQ=255)
                                                                                                                                         tcgCGGGCTCGCCCGCGACGAGGGTCACTCGGTCTTGAGGGCCTGCGCCATCAGCTCATATGACCGAATGCGGTCCGCATGCTCATACAGGTCACAGGTAAAGATCAGCTCATCCGCACCCGTCTGTTCCAGCAGCACTTCTACCTTGGc        >  2:3996723/3‑150 (MQ=255)
                                                                                                                                         gtgggggCGCGCCCGCGACGAGGGTCACTCGGTCTTGAGGGCCTGCGCCATCAGCTCATATGACCGAATGCGGTCCGCATGCTCATACAGGTCACAGGTAAAGATCAGCTCATCCGCACCCGTCTGTTCCAGCAGCACTTCTACCTTGGc        <  2:1254458/146‑1 (MQ=255)
                                                                                                                                            cGGGCGAGCCCGCGACGAGGGTCACTCGGTCTTGAGGGCCTGCGCCATCAGCTCATATGACCGAATGCGGTCCGCATGCTCATACAGGTCACAGGTAAAGATCAGCTCATCCGCACCCGTCTGTTCCAGCAGCACTTCTACCTTGGCCtt     >  2:1610682/1‑150 (MQ=255)
                                                                                                                                             gggCGAGCCCGCGACGAGGGTCACTCGGTCTTGAGGGCCTGCGCCATCAGCTCATATGACCGAATGCGGTCCGCATGCTCATACAGGTCACAGGTAAAGATCAGCTCATCCGCACCCGTCTGTTCCAGCAGCACTTCTACCTTGGCCttt    <  1:4989621/150‑1 (MQ=255)
                                                                                                                                               gCGAGCCCGCGACGAGGGTCACTCGGTCTTGAGGGCCTGCGCCATCAGCTCATATGACCGAATGCGGTCCGCATGCTCATACAGGTCACAGGTAAAGATCAGCTCATCCGCACCCGTCTGTTCCAGCAGCACTTCTACCTTGGCCTTTAc  <  1:5146302/150‑1 (MQ=255)
                                                                                                                                               gCGAGCCCGCGACGAGGGTCACTCGGTCTTGAGGGCCTGCGCCATCAGCTCATATGACCGAATGCGGTCCGCATGCTCATACAGGTCACAGGTAAAGATCAGCTCATCCGCACCCGTCTGTTCCAGCAGCACTTCTACCTTGGCCTTTAc  >  1:3197840/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
CTGGAGATCGTCAAGAAAGGCGAGCCGGCGGGCAAGGTGATTGGTAGCTGCGAGGCCGGTACCAAGGAGATTGTTTACCGCCGCGGCTAATGTGTTTGCTCATTCGCGGGCTTGCTCGCTCCCAAGTCAACCTATGGGAGCGGGCGAGCCCGCGACGAGGGTCACTCGGTCTTGAGGGCCTGCGCCATCAGCTCATATGACCGAATGCGGTCCGCATGCTCATACAGGTCACAGGTAAAGATCAGCTCATCCGCACCCGTCTGTTCCAGCAGCACTTCTACCTTGGCCTTTAC  >  NC_002947/95519‑95811

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