Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,504,077 T→A 5.4% 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,0770TA 5.4% 294.1 / 6.3 93intergenic (+35/‑228)rpsI/petA30S ribosomal protein S9/ubiquinol‑cytochrome c reductase iron‑sulfur subunit
Reads supporting (aligned to +/- strand):  ref base T (37/51);  new base A (2/3);  total (39/54)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.38e-01

CCTGCGTGGCGCTCTGCGTCGTGCTGGCTACGTCACCCGCGACGCTCGTGAAGTTGAGCGTAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAGCAAAGCGTGGGTTGTCGCCCTTGCCGGGGTAATCACCTTGTCAGTGTGTGGATTTGTTTCTTACCATGACG  >  NC_002947/1503928‑1504225
                                                                                                                                                     |                                                                                                                                                     
ccTGCGTGGCGCTCTGCGTCGTGCTGGCTACGTCACCCGCGACGCTCGTGAAGTTGAGCGTAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCtt                                                                                                                                                       <  1:1346553/150‑1 (MQ=255)
ccTGCGTGGCGCTCTGCGTCGTGCTGGCTACGTCACCCGCGACGCTCGTGAAGTTGAGCGTAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCtt                                                                                                                                                       <  2:1732925/150‑1 (MQ=255)
        gcgcTCTGCGTCGTGCTGGCTACGTCACCCGCGACGCTCGTGAAGTTGAGCGTAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCgg                                                                                                                                               <  2:2023779/150‑1 (MQ=255)
         cgcTCTGCGTCGTGCTGGCTACGTCACCCGCGACGCTCGTGAAGTTGAGCGTAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCAAGGAACCggg                                                                                                                                              <  2:1935736/150‑1 (MQ=255)
          gcTCTGCGTCGTGCTGGCTACGTCACCCGCGACGCTCGTGAAGTTGAGCGTAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGc                                                                                                                                             <  1:2048816/150‑1 (MQ=255)
          gcTCTGCGTCGTGCTGGCTACGTCACCCGCGACGCTCGTGAAGTTGAGCGTAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGc                                                                                                                                             <  1:1911479/150‑1 (MQ=255)
           ctctGCGTCGTGCTGGCTACGTCACCCGCGACGCTCGTGAAGTTGAGCGTAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCt                                                                                                                                           <  1:1993193/150‑1 (MQ=255)
             ctGCGTCGTGCTGGCTACGTCACCCGCGACGCTCGTGAAGTTGAGCGTAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCttt                                                                                                                                         <  1:724803/150‑1 (MQ=255)
                  tcgtGCTGGCTACGTCACCCGCGACGCTCGTGAAGTTGAGCGTAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGAACAGTCCGTACTACAACGATTGATGCTTCGGCTGGAGTCGAGAAAAGCCCGGTTCCTTGGAACCGGGCtttttttt                                                                                                                                    <  1:656861/150‑1 (MQ=255)
                      gCTGGCTACGTCACCCGCGACGCTCGTGAAGTTGAGCGTAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGc                                                                                                                                <  1:1635090/150‑1 (MQ=255)
                            tACGTCACCCGCGACGCTCGTGAAGTTGAGCGTAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGaa                                                                                                                          >  2:1423176/1‑150 (MQ=255)
                             aCGTCACCCGCGACGCTCGTGAAGTTGAGCGTAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTCCCTTGGAACCGGGCTTTTTTTTTTAGCCTTGaa                                                                                                                          <  1:512880/150‑1 (MQ=255)
                                tCACCCGCGACGCTCGTGAAGTTGAGCGTAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTaa                                                                                                                      <  2:514324/150‑1 (MQ=255)
                                tCACCCGCGACGCTCGTGAAGTTGAGCGTAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTaa                                                                                                                      <  2:2124783/150‑1 (MQ=255)
                                   cccGCGACGCTCGTGAAGTTGAGCGTAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCt                                                                                                                   >  2:2030615/1‑150 (MQ=255)
                                          cGCTCGTGAAGTTGAGCGTAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTcagca                                                                                                            <  1:2259097/150‑1 (MQ=255)
                                              cGTGAAGTTGAGCGTAAGAAAGTGGGTCTGCATAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGGCAATCTCCCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGcagtca                                                                                                        <  1:419749/150‑1 (MQ=255)
                                               gTGAAGTTGAGCGTAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGaa                                                                                                                          <  1:1751164/131‑1 (MQ=255)
                                               gTGAAGTTGAGCGTAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGaa                                                                                                                          >  2:1751164/1‑131 (MQ=255)
                                                tGAAGTTGAGCGTAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAAc                                                                                                                     >  1:266819/1‑135 (MQ=255)
                                                tGAAGTTGAGCGTAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAAc                                                                                                                     <  2:266819/135‑1 (MQ=255)
                                                 gAAGTTGAGCGTAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCt                                                                                                                              <  2:2167533/125‑1 (MQ=255)
                                                 gAAGTTGAGCGTAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCt                                                                                                                              >  1:2167533/1‑125 (MQ=255)
                                                    gTTGAGCGTAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTg                                                                                                                            <  1:1499224/124‑1 (MQ=255)
                                                    gTTGAGCGTAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTg                                                                                                                            >  2:1499224/1‑124 (MQ=255)
                                                     ttGAGCGTAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGc                                                                                                                                <  2:1890238/119‑1 (MQ=255)
                                                     ttGAGCGTAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGc                                                                                                                                >  1:1890238/1‑119 (MQ=255)
                                                       gAGCGTAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCt                                                                                               >  2:1025385/1‑150 (MQ=255)
                                                        aGCGTAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCAAGGAACCggg                                                                                                                                              <  1:1971019/103‑1 (MQ=255)
                                                        aGCGTAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCAAGGAACCggg                                                                                                                                              >  2:1971019/1‑103 (MQ=255)
                                                          cGTAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGt                                                                                            >  1:354399/1‑150 (MQ=255)
                                                          cGTAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGt                                                                                            >  2:1751833/1‑150 (MQ=255)
                                                            tAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCGGGGCGAAAAAGCACGCATCCTTGAAACCGGGCTGTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCg                                                                                          >  1:1374281/1‑150 (MQ=255)
                                                              aGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCa                                                                                        <  2:2087998/150‑1 (MQ=255)
                                                              aGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCa                                                                                        >  1:2087998/1‑150 (MQ=255)
                                                              aGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTCGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCa                                                                                        <  1:1945659/150‑1 (MQ=255)
                                                               gAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCAt                                                                                       >  1:1502139/1‑150 (MQ=255)
                                                                    tGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGACCCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACg                                                                                  <  1:623239/150‑1 (MQ=255)
                                                                        tCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGa                                                                              <  1:1423176/150‑1 (MQ=255)
                                                                            cGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCaa                                                                          >  1:928417/1‑150 (MQ=255)
                                                                            cGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCaa                                                                          <  1:1036673/150‑1 (MQ=255)
                                                                                      gTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAGCAAAGCGTg                                                                <  2:1502139/150‑1 (MQ=255)
                                                                                          gCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCgg                                                                                                                                               <  1:775211/68‑1 (MQ=255)
                                                                                          gCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCgg                                                                                                                                               >  2:775211/1‑68 (MQ=255)
                                                                                          gCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCt                                                                                                                   >  1:646589/1‑95 (MQ=255)
                                                                                          gCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCt                                                                                                                   <  2:646589/95‑1 (MQ=255)
                                                                                               cTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGcagtcagt                                                                                                      <  2:2098399/103‑1 (MQ=255)
                                                                                               cTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCGTTAACTAACCTGTCAGcagtcagt                                                                                                      >  1:2098399/1‑103 (MQ=255)
                                                                                                tCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCgg                                                                                                                                               >  1:1018912/1‑62 (MQ=255)
                                                                                                tCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCgg                                                                                                                                               <  2:1018912/62‑1 (MQ=255)
                                                                                                tCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAGCAAAGCGTGGGTTGTCGcc                                                      >  1:700643/1‑150 (MQ=255)
                                                                                                tCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAGCAAAGCGTGGGTTGTCGcc                                                      >  1:244664/1‑150 (MQ=255)
                                                                                                  aGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAGCAAAGCGTGGGTTGTCGCCCt                                                    >  2:710151/1‑150 (MQ=255)
                                                                                                  aGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAGCAAAGCGTGGGTTGTCGCCCt                                                    >  1:1804305/1‑150 (MQ=255)
                                                                                                   actgccgaaGCGGTATTAACGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAGCAAAGCGTGGGTTGTCGccc                                                     <  1:956904/141‑1 (MQ=255)
                                                                                                       tCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCg                                                                                          <  2:2098031/107‑1 (MQ=255)
                                                                                                       tCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCg                                                                                          >  1:2098031/1‑107 (MQ=255)
                                                                                                       tCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAGCAAAGCGTGGGTTGTCGCCCTTGCCg                                               >  1:154475/1‑150 (MQ=255)
                                                                                                         cAAGCGTTAATCCCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACAAGGATCTCTATCAAAAAACCAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAGCAAAGCGTGGGTTGTCGCCCTTGCCggg                                             <  1:550335/150‑1 (MQ=255)
                                                                                                         cAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGACCCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAGCAAAGCGTGGGTTGTCGCCCTTGCCggg                                             <  1:1640880/150‑1 (MQ=255)
                                                                                                         cAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAGCAAAGCGTGGGTTGTCGCCCTTGCCggg                                             >  2:484627/1‑150 (MQ=255)
                                                                                                         cAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAGCAAAGCGTGGGTTGTCGCCCTTGCCggg                                             >  1:1279381/1‑150 (MQ=255)
                                                                                                          aaGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAGCAAAGCGTGGGTTGTCGCCCTTGCCgggg                                            <  2:244664/150‑1 (MQ=255)
                                                                                                          aaGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAGCAAAGCGTGGGTTGTCGCCCTTGCCgggg                                            <  2:2197842/150‑1 (MQ=255)
                                                                                                              gTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTTTCGCATAGACGTAGATCAAGCAAAGCGTGGGTTGTCGCCCTTGCCGGGGt                                           <  2:1418401/147‑1 (MQ=255)
                                                                                                              gTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAGCAAAGCGTGGGTTGTCGCCCTTGCCGGGGt                                           >  1:1418401/1‑147 (MQ=255)
                                                                                                               ttAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGTCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAGCAAAGCGTGGGTTGTc                                                         >  1:2112498/1‑132 (MQ=255)
                                                                                                               ttAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGTCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAGCAAAGCGTGGGTTGTc                                                         <  2:2112498/132‑1 (MQ=255)
                                                                                                               ttAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAGCAAAGCGTGGGTTGTCGCCCTTGCCGGGGTAATc                                       <  1:470389/150‑1 (MQ=255)
                                                                                                                tAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTca                                                                                                               >  2:1842525/1‑77 (MQ=255)
                                                                                                                tAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTca                                                                                                               <  1:1842525/77‑1 (MQ=255)
                                                                                                                 aaTACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAGCAAAGCGTGGGTTGTCGCCCTTGCCGGGGTAATCAc                                     >  2:375992/1‑150 (MQ=255)
                                                                                                                    agctcTTGGGCAGTCACAAAAAGCCCCGGTCCTTGGGTCCGGGCTTTTTCTTTAGCCTTGAACGAACCTGTCAGCTGTTAGTGGCCGTTTGTCGCTTAGACGTATATCAAGCAAAGCCTGGGTTGTCGCCCTTGCCGGGGTAATCACCtt                                  <  2:6846/146‑1 (MQ=255)
                                                                                                                    aCCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAGCAAAGCGTGGGTTGTCGCCCTTGCCGGGGTAATCACCtt                                  <  1:365632/150‑1 (MQ=255)
                                                                                                                       gCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCa                                                                                        >  1:1793215/1‑93 (MQ=255)
                                                                                                                       gCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCa                                                                                        <  2:1793215/93‑1 (MQ=255)
                                                                                                                       gCTTCGGCAGTCGAAAAAAGCCCGGTTCATTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAGCAAAGCGTGGGTTGTCGCCCTTGCCGGGGTAATCACCTTGTc                               <  2:1088393/150‑1 (MQ=255)
                                                                                                                            ggCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAATCAAAGCGTGGGTTGTCTCCCTTGCCGGGGTAATCACCTTGTCAgtgt                          <  1:1025385/150‑1 (MQ=255)
                                                                                                                            ggCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAGCAAAGCGTGGGTTGTCGCCCTTGCCGGGGTAATCACCTTGTCAgtgt                          <  1:1751833/150‑1 (MQ=255)
                                                                                                                             gCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAGCAAAGCGTGGGTTGTCGCCCTTGCCGGGGt                                           >  2:1966599/1‑132 (MQ=255)
                                                                                                                             gCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAGCAAAGCGTGGGTTGTCGCCCTTGCCGGGGt                                           <  1:1966599/132‑1 (MQ=255)
                                                                                                                              cAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCATTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAGCAAAGCGTGGGTTGTCGCCCTTGCCGGGGTAATCACCTTGTCAgtgtgt                        <  2:1279381/150‑1 (MQ=255)
                                                                                                                                   aaaaaaaGCCCGGTTCCAAGGAACCGGGCTTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAGCAAAGCGTGGGTTGTc                                                         >  1:1855608/2‑113 (MQ=38)
                                                                                                                                   aaaaaaaGCCCGGTTCCAAGGAACCGGGCTTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAGCAAAGCGTGGGTTGTc                                                         <  2:1855608/112‑1 (MQ=38)
                                                                                                                                      aaaaGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAGCAAAGCGTGGGTTGTCGCCCTTGCCGGGGTAATCACCTTGTCAGTGTGTGGATTTGt                <  1:1334976/150‑1 (MQ=255)
                                                                                                                                       atctcccGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAGCAAAGCGTGGGTTGTCGCCCTTGCCGGGGTAATCACCTTGTCAGTGTGTGGATTTGtt               <  1:504656/146‑1 (MQ=255)
                                                                                                                                       aaaGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAGCAAAGCGTGGGTTGTCGccc                                                     <  2:2130559/112‑1 (MQ=255)
                                                                                                                                       aaaGCCCGGTTCCTTGGAACCGGGCTGTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAGCAAAGCGTGGGTTGTCGccc                                                     >  1:2130559/1‑112 (MQ=255)
                                                                                                                                        aaGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTTGCCTTGAACTAATCTGTTATCAGTTAGTGGCCTGTTTTCGGATAGACGTAGATTAAGCAAAGCGTGGGTTGTCGCCCTTGCCGGGGTAATCACCTTGTCAGTGTTTGGTTTTGttt              >  2:2048066/1‑150 (MQ=255)
                                                                                                                                        aaGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAGCAAAGCGTGGGTTGTCGCCCTTGCCGGGGTAATCACCTTGTCAGTGTGTGGATTTGttt              <  2:700643/150‑1 (MQ=255)
                                                                                                                                         aGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAGCAAAGCGTGGGTTGTCGCCCTTGCCGGGGTAATCACCTTGTCAGTGTGTGGATTTGTTTc             <  1:375992/150‑1 (MQ=255)
                                                                                                                                            ccGGTTCCTTGGAACCGGGCTTTTTTTTTTTTTTTTTACTAAACTATCAGCCGTCAGTGGTCGCTTGTCGTGTAGACGTAGATAAAGCAAAGCGTCGCTTGTCGCCCTTGCCGAGGTAATCATCTTGTTAGTGTGTGGatttatttatgt          >  2:504656/1‑146 (MQ=255)
                                                                                                                                            ccGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAGCAAAGCGt                                                                 >  1:2135316/1‑95 (MQ=255)
                                                                                                                                            ccGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAGCAAAGCGt                                                                 <  2:2135316/95‑1 (MQ=255)
                                                                                                                                            ccGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAGCAAAGCGTGGGTTGTCGCCCTTGCCGGGGTAATCACCTTGTCAGTGTGTGGATTTGTTTCTTa          >  1:1403492/1‑150 (MQ=255)
                                                                                                                                                    ttGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAGCAAAGCGTGGGTTGTCGCCCTTGCCGGGGTAATCACCTTGTCAGTGTGTGGATTTGTTTCTTACCATGACg  >  1:27247/1‑150 (MQ=255)
                                                                                                                                                    ttGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAGCAAAGCGTGGGTTGTCGCCCTTGCCGGGGTAATCACCTTGTCAGTGTGTGGATTTGTTTATTAGCacgacg  >  2:1376492/1‑150 (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: