Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,504,077 T→A 15.8% intergenic (+35/‑228) rpsI → / → petA 30S ribosomal protein S9/ubiquinol‑cytochrome c reductase iron‑sulfur subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,504,0770TA15.8% 196.7 / 29.6 76intergenic (+35/‑228)rpsI/petA30S ribosomal protein S9/ubiquinol‑cytochrome c reductase iron‑sulfur subunit
Reads supporting (aligned to +/- strand):  ref base T (30/34);  new base A (7/5);  total (37/39)
Fisher's exact test for biased strand distribution p-value = 5.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.50e-01

CCTGCGTGGCGCTCTGCGTCGTGCTGGCTACGTCACCCGCGACGCTCGTGAAGTTGAGCGTAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAGCAAAGCGTGGGTTGTCGCCCTTGCCGGGGTAATCACCTTGTCAGTGTGTGGATTTGTTTCTTACCATGACG  >  NC_002947/1503928‑1504225
                                                                                                                                                     |                                                                                                                                                     
ccTGCGTGGCGCTCTGCGTCGTGCTGGCTACGTCACCCGCGACGCTCGTGAAGTTGAGCGTAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCtt                                                                                                                                                       >  1:1325798/1‑150 (MQ=255)
 cTGCGTGGCGCTCTGCGTCGTGCTGGCTACGTCACCCGCGACGCTCGTGAAGTTGAGCGTAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTg                                                                                                                                                      <  1:1139711/150‑1 (MQ=255)
        gcgcTCTGCGTCGTGCTGGCTACGTCACCCGCGACGCTCGTGAAGTTGAGCGTAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCgg                                                                                                                                               <  1:1055137/150‑1 (MQ=255)
         cgcTCTGCGTCGTGCTGGCTACGTCACCCGCGACGCTCGTGAAGTTGAGCGTAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCggg                                                                                                                                              <  2:1977995/150‑1 (MQ=255)
         cgcTCTGCGTCGTGCTGGCTACGTCACCCGCGACGCTCGTGAAGTTGAGCGTAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCggg                                                                                                                                              <  2:1115940/150‑1 (MQ=255)
            tctGCGTCGTGCTGGCTACGTCACCCGCGACGCTCGTGAAGTTGAGCGTAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCtt                                                                                                                                           >  1:567860/1‑150 (MQ=255)
              tGCGTCGTGCTGGCTACGTCACCCGCGACGCTCGTGAAGTTGAGCGTAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCtttt                                                                                                                                         <  1:2197242/150‑1 (MQ=255)
                 gtcgtGCTGGCTACGTCACCCGCGACGCTCGTGAAGTTGAGCGTAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCttttttt                                                                                                                                      <  2:1980692/150‑1 (MQ=255)
                    gtGCTGGCTACGTCACCCGCGACGCTCGTGAAGTTGAGCGTAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTa                                                                                                                                   <  1:214325/150‑1 (MQ=255)
                         ggCTACGTCACCCGCGACGCTCGTGAAGTTGAGCGTAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCtt                                                                                                                             <  1:750798/150‑1 (MQ=255)
                         ggCTACGTCACCCGCGACGCTCGTGAAGTTGAGCGTAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCtt                                                                                                                             <  1:1689629/150‑1 (MQ=255)
                            tACGTCACCCGCGACGCTCGTGAAGTTGAGCGTAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCtt                                                                                                                                                       >  1:1085953/1‑122 (MQ=255)
                            tACGTCACCCGCGACGCTCGTGAAGTTGAGCGTAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCtt                                                                                                                                                       <  2:1085953/122‑1 (MQ=255)
                               gTCACCCGCGACCCTCGTGAAGTTGAGCGGAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGGCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTa                                                                                                                       >  2:730875/1‑150 (MQ=255)
                                 cACCCGCGACGCTCGTGAAGTTGAGCGTAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAAc                                                                                                                     <  2:1325798/150‑1 (MQ=255)
                                 cACCCGCGACGCTCGTGAAGTTGAGCGTAAGAAAGTGCGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCgg                                                                                                                                               >  2:1748442/1‑125 (MQ=255)
                                 cACCCGCGACGCTCGTGAAGTTGAGCGTAAGAAAGTGCGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCgg                                                                                                                                               <  1:1748442/125‑1 (MQ=255)
                                      gcgACGCTCGTGAAGTTGAGCGTAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTc                                                                                                                <  1:2002432/150‑1 (MQ=255)
                                               gTGAAGTTGAGCGTAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTTGCCCTTGAACTAACCTTTCAGaggtca                                                                                                        >  2:1356601/1‑150 (MQ=255)
                                                   aGTTGAGCGTAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGc                                                                                                   <  2:1914506/150‑1 (MQ=255)
                                                     ttGAGCGTAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCg                                                                                                 >  2:1063568/1‑150 (MQ=255)
                                                          cGTAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGt                                                                                            >  2:1216275/1‑150 (MQ=255)
                                                            tAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCg                                                                                          >  2:849597/1‑150 (MQ=255)
                                                            aaaGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGATTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCAGCTTGTCg                                                                                          <  2:567860/149‑1 (MQ=255)
                                                             aaGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTGCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGc                                                                                         >  1:514987/1‑150 (MQ=255)
                                                               gAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCAt                                                                                       >  2:2108093/1‑150 (MQ=255)
                                                                  aGTGGGTCTGCGTACAGCGCGTAAGCGGCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGa                                                                                    <  1:1063568/150‑1 (MQ=255)
                                                                     gggTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGt                                                                                 <  1:928555/150‑1 (MQ=255)
                                                                       gTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCAAGGAACCggg                                                                                                                                              <  2:1918755/88‑1 (MQ=255)
                                                                       gTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTCCGGCAGTCGAAAAAAGCCCGGTTCCAAGGAACCggg                                                                                                                                              >  1:1918755/1‑88 (MQ=38)
                                                                        tCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGa                                                                              <  1:770346/150‑1 (MQ=255)
                                                                        tCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCAAGGAACCGGGCttttttc                                                                                                                                      >  1:494779/1‑94 (MQ=255)
                                                                        tCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCAAGGAACCGGGCttttttc                                                                                                                                      <  2:494779/95‑2 (MQ=255)
                                                                           gCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTcagca                                                                                                            >  1:519206/1‑117 (MQ=255)
                                                                           gCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTcagca                                                                                                            <  2:519206/117‑1 (MQ=255)
                                                                            cGTAAAGCGCGTAAGCGTCCTCAGCACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAg                                                                         <  1:1682253/150‑1 (MQ=255)
                                                                                          gCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCAt                                                                                       <  1:380738/123‑1 (MQ=255)
                                                                                          gCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCAt                                                                                       >  2:380738/1‑123 (MQ=255)
                                                                                            gTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAGCAAAGCGTGGGTTGt                                                          >  2:946873/1‑150 (MQ=255)
                                                                                               cTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGa                                                                              >  2:467590/1‑127 (MQ=255)
                                                                                               cTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGa                                                                              <  1:467590/127‑1 (MQ=255)
                                                                                               cTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCAAGGAACCGGGCtttttt                                                                                                                                       <  1:1179911/71‑1 (MQ=25)
                                                                                               cTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCAAGGAACCGGGCtttttt                                                                                                                                       >  2:1179911/1‑71 (MQ=25)
                                                                                                         cAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTg                                                                                             <  2:2288461/102‑1 (MQ=255)
                                                                                                         cAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTg                                                                                             >  1:2288461/1‑102 (MQ=255)
                                                                                                         cAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAGCAAAGCGTGGGTTGTCGCCCTTGCCggg                                             >  2:1373577/1‑150 (MQ=255)
                                                                                                         cAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAGCAAAGCGTGGGTTGTCGCCCTTGCCggg                                             >  2:1555439/1‑150 (MQ=255)
                                                                                                         cAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAGCAAAGCGTGGGTTGTCGCCCTTGCCggg                                             >  1:1121219/1‑150 (MQ=255)
                                                                                                         cAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAGCAAAGCGTGGGTTGTCGCCCTTGCCggg                                             <  1:849597/150‑1 (MQ=255)
                                                                                                         cAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAATGGCCGCTTGTCGCATAGACGTAGATCAAGCAAAGCGTGGGTTGTCGCCCTTGCCggg                                             >  2:991663/1‑150 (MQ=255)
                                                                                                         cAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGACCGGTTCCTTGGAACCGAGAAATTTTTTTAGCCTTGAACTAACCTGTCAGCTGTCAGTTGCCGATTGTTGCATAGAAGTAGATCAAGCAAAGCGTGGGTTGTCGCCTTTGCCggg                                             >  1:418087/1‑150 (MQ=255)
                                                                                                          aaGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAGCAAAGCGTGGGTTGTCGCCCTTGCCgggg                                            >  2:1725970/1‑150 (MQ=255)
                                                                                                          aaGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCAAGGAACCg                                                                                                                                                >  2:769878/1‑51 (MQ=25)
                                                                                                          aaGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCAAGGAACCg                                                                                                                                                <  1:769878/51‑1 (MQ=25)
                                                                                                               ttAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAGCAAAGCGTGGGTTGTCGCCCTTGCCGGGGTAATc                                       <  1:1373577/150‑1 (MQ=255)
                                                                                                                  aTACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAGCAAAGCGTGGGTTGTCGCCCTTGCCGGGGTAATCAcc                                    >  2:1517921/1‑150 (MQ=255)
                                                                                                                  aTACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAGCAAAGCGTGGGTTGTCGCCCTTGCCGGGGTAATCAcc                                    <  1:946873/150‑1 (MQ=255)
                                                                                                                   tACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAGCAAAGCGTGGGTTGTCGCCCTTGCCGGGGTAATCACCt                                   >  1:1021874/1‑150 (MQ=255)
                                                                                                                    aCCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTTAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAACAAAGCGTGGGTTGTCGCCCTTGCCGGGGTAATCACCtt                                  >  2:1308178/1‑150 (MQ=255)
                                                                                                                     ccGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAGCAAAGCGTGGGTTGTc                                                         >  2:1324211/1‑126 (MQ=255)
                                                                                                                     ccGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAGCAAAGCGTGGGTTGTc                                                         <  1:1324211/126‑1 (MQ=255)
                                                                                                                      cGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAGCAAAGCGTGGGTTGTCGCCCTTGCCGGGGTAATCACCTTGt                                <  1:1725970/150‑1 (MQ=255)
                                                                                                                              gctaaaaaaaaaGCCCGGTTCCAAGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAGCAAAGCGTGGGTTGTCGCCCTTGCCGGGGTAATCACCTTGTCAgtgtgt                        <  2:208487/144‑1 (MQ=255)
                                                                                                                               aggcGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAGCCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAGCAAAGCGTGGGTTGTCGCCCTTGCCGGGGTAATCACCTTGTCAGtgtgtg                       >  1:83127/4‑150 (MQ=255)
                                                                                                                               aGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAGCAAAGCGTGGGTTGTCGCCCTTGCCGGGGTAATCACCTTGTCAGtgtgtg                       >  1:83357/1‑150 (MQ=255)
                                                                                                                                 tCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTa                                                                                                                                   >  2:125185/1‑41 (MQ=255)
                                                                                                                                 tCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTa                                                                                                                                   <  1:125185/41‑1 (MQ=255)
                                                                                                                                  cGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAg                                                                         >  2:541597/1‑97 (MQ=255)
                                                                                                                                  cGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAg                                                                         <  1:541597/97‑1 (MQ=255)
                                                                                                                                  cGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAGCAAAGCGTGGGTTGTCGCCCTTGCCGGGGTAATCACCTTGTCAGTGTGTGGAt                    <  2:318236/150‑1 (MQ=255)
                                                                                                                                        aaGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAGCAAAGCGTGGGTTGTCGCCCTTGCCGGGGTAATCACCTTGTCAGTGTGTGGATTTGttt              <  2:514987/150‑1 (MQ=255)
                                                                                                                                         aGCCCGGTTCCAAGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAGCAAAGCGTGGGTTGTCGCCCTTGCCGGGGTAATCACCTTGTCAGTGTGTGGATTTGTTTc             >  1:1390110/1‑150 (MQ=255)
                                                                                                                                           cccGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAGCAAAGCGTGGGTTGTCGCCCTTGCCGGGGTAATCACCTTGTCAGTGTGTGGATTTGTTTCtt           >  1:688100/1‑150 (MQ=255)
                                                                                                                                           cccGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAATGGCCGCTTGTCGCATAGACGTAGATCAAGCAAAGCGTGGGTTGTCGCCCTTGCCGGGGTAATCACCTTGTCAGTGTGTGGATTTGTTTCtt           <  1:991663/150‑1 (MQ=255)
                                                                                                                                           cccGGTTCCAAGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAGCAAAGCGTGGGTTGTCGCCCTTGCCGGGGTAATCACATTGTCAGTGTGTGGATTTGTTTCtt           >  1:1645772/1‑150 (MQ=255)
                                                                                                                                            ccGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAGCAAAGCGTGGGTTGTCGCCCTTGCCGGGGTAATCACCTTGTCAGTGTGTGGATTTGTTTCTTa          <  2:1345536/150‑1 (MQ=255)
                                                                                                                                               gTTCCAAGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCGAGCACAGCGTGGATTGTCGCCCTTGCCGGGGTAAGTACCTTGTTAGTGTTTGGGTGTGTTTGTTACCa       >  1:277119/1‑150 (MQ=255)
                                                                                                                                                ttCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGa                                                                              >  2:1608173/1‑78 (MQ=255)
                                                                                                                                                ttCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGa                                                                              <  1:1608173/78‑1 (MQ=255)
                                                                                                                                                    ttGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAGCAAAGCGTGGGTTGTCGCCCTTGCCGGGGTAATCACCTTGTCAGTGTGTGGATTTGTTTCTTACCATGACg  <  1:1555439/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CCTGCGTGGCGCTCTGCGTCGTGCTGGCTACGTCACCCGCGACGCTCGTGAAGTTGAGCGTAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAGCAAAGCGTGGGTTGTCGCCCTTGCCGGGGTAATCACCTTGTCAGTGTGTGGATTTGTTTCTTACCATGACG  >  NC_002947/1503928‑1504225

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