Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,504,068 T→A 9.1% 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,0680TA 9.1% 210.3 / 9.0 76intergenic (+35/‑228)rpsI/petA30S ribosomal protein S9/ubiquinol‑cytochrome c reductase iron‑sulfur subunit
Reads supporting (aligned to +/- strand):  ref base T (35/34);  new base A (3/4);  total (38/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.18e-01

CTGCGTGGCGCTCTGCGTCGTGCTGGCTACGTCACCCGCGACGCTCGTGAAGTTGAGCGTAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAGCAAAGCGTGGGTTGTCGCCCTTGCCGGGGTAATCACCTTGTCAGTGTGTGGATTTGTTTCTTACCATGA  >  NC_002947/1503920‑1504214
                                                                                                                                                    |                                                                                                                                                   
cTGCGTGGCGCTCTGCGTCGTGCTGGCTACGTCACCCGCGACGCTCGTGAAGTTGAGCGTAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCtt                                                                                                                                                     <  2:325662/149‑1 (MQ=255)
  gCGTGGCGCTCTGCGTCGTGCTGGCTACGTCACCCGCGACGCTCGTGAAGTTGAGCGTAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTg                                                                                                                                                    >  3:334512/1‑148 (MQ=255)
   cGTGGCGCTCTGCGTCGTGCTGGCTACGTCACCCGCGACGCTCGTGAAGTTGAGCGTAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGACAGTCGAAAAAAGCCCGGTTCCTTg                                                                                                                                                    <  4:303750/147‑1 (MQ=255)
      ggCGCGCTGCGTCGAGCTGGCTACGTCACCCGCGACGCTCGTGAACTTGAGCGTAAGAAAGTGCGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAACCGTTAATACCGCTTCGGCAGTCCAAAAAAGCCCGGTTCCTTGGAAcc                                                                                                                                               <  8:384066/149‑1 (MQ=255)
        cgcTCTGCGTCGTGCTGGCTACGTCACCCGCGACGCTCGTGAAGTTGAGCGTAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCgg                                                                                                                                             <  4:342639/149‑1 (MQ=255)
        cgcTCTGCGTCGTGCTGGCTACGTCACCCGCGACGCTCGTGAAGTTGAGCGTAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCAAGGAACCgg                                                                                                                                             <  2:89623/149‑1 (MQ=255)
         gcTCTGCGTCGAGCTGGCTACGTCACACGCGACGCTCGTGAAGTTGAGCGTAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCTTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCAAGGAACCggg                                                                                                                                            <  1:396823/149‑1 (MQ=255)
            ctGCGTCGTGCTGGCTACGTCACCCGCGACGCTCGTGAAGTTGAGCGTAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCtt                                                                                                                                        >  1:248943/1‑149 (MQ=255)
             tGCGTCGTGCTGGCTACGTCACCCGCGACGCTCGTGAAGTTGAGCGTAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCAAGGAACCGGGCt                                                                                                                                         <  7:73845/147‑1 (MQ=255)
               cgtcgtGCTGGCTACGTCACCCGCGACGCTCGTGAAGTTGAGCGTAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCtttt                                                                                                                                      <  8:55184/148‑1 (MQ=255)
               cgtcgtGCTGGCTACGTCACCCGCGACGCTCGTGAAGTTGAGCGTAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCttt                                                                                                                                       <  4:352423/147‑1 (MQ=255)
               cgtcgtGCTGGCTACGTCACCCGCGACGCTCGTGAAGTTGAGCGTAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCttt                                                                                                                                       <  2:133055/147‑1 (MQ=255)
                gtcgtGCTGGCTACGTCACCCGCGACGCTCGTGAAGTTGAGCGTAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTACTACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCtttttt                                                                                                                                    <  1:99503/149‑1 (MQ=255)
                   gtGCTGGCTACGTCACCCGCGACGCTCGTGAAGTTGAGCGTAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGGTCCTTGGAACCGGGCttttttttt                                                                                                                                 <  5:138240/149‑1 (MQ=255)
                     gCTGGCTACGTCACCCGCGACGCTCGTGAAGTTGAGCGTAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTa                                                                                                                                >  8:368109/1‑148 (MQ=255)
                     gCTGGCTACGTCACCCGCGACGCTCGTGAAGTTGAGCGTAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGCTCCTTGGAACCGGGCTTTTTTTTTa                                                                                                                                <  7:368109/148‑1 (MQ=255)
                      cTGGCTACGTCACCCGCGACGCTCGTGAAGTTGTGCGTAAGAAAGTGGGTCTGCGTAAAGCGCGTTAGCGTCCTCAGTACTCCAAGCGTTATTACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGc                                                                                                                              <  5:155296/149‑1 (MQ=255)
                      cTGGCGACGTCACCCGCGACGCTCGTGAAGTTGAGCGTAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGc                                                                                                                              <  7:83902/149‑1 (MQ=255)
                         gCTACGTCACCCGCGACGCTCGTGAAGTTGAGCGTAAGAAAGGGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGcc                                                                                                                             <  4:50738/147‑1 (MQ=255)
                                 aCCCGCGACGCTCGTGAAGTTGAGCGTAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAAc                                                                                                                   <  7:95598/149‑1 (MQ=255)
                                   ccGCGACGCTCGTGAAGTTGAGCGTAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTAGCCTTGAACTAACCTg                                                                                                                <  2:273228/149‑1 (MQ=255)
                                        tcGCTAGTGAAGTTGAGCGTAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTGCTCAGTACTCAAATCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACGGGGCTTTTTTTTTAGCCTTGAACTAACCTGTcag                                                                                                            <  5:13408/148‑1 (MQ=255)
                                          gCTCGTGAAGTTGAGCGTAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTcag                                                                                                            >  8:306244/1‑147 (MQ=255)
                                               tGAAGTTGAGCGTAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGcagtcag                                                                                                     >  5:10870/1‑149 (MQ=255)
                                                gAAGTTGAGCGTAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGcagtcagt                                                                                                    >  6:312765/1‑149 (MQ=255)
                                                     tGAGCGTAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCg                                                                                               >  3:310848/1‑149 (MQ=255)
                                                         cGTAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACAGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCGTGAACTAACCTGTCAGCAGTAAGTGGCCGCTTg                                                                                           >  5:321142/1‑149 (MQ=255)
                                                         cGTAAGAAAGTGGGTCTGCGTAAAGCGCGAAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAATAAAGCCCGGTACCTTGGAACCGGGCGTTTTTTTTAGCCTAGAACTAGCCGGTCAGCAGTCAGAGGCCGCTTg                                                                                           >  4:61747/1‑149 (MQ=255)
                                                           tAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTc                                                                                         >  5:181113/1‑149 (MQ=255)
                                                           tAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTc                                                                                         <  7:306244/149‑1 (MQ=255)
                                                            aaGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTc                                                                                         <  4:334512/148‑1 (MQ=255)
                                                            aaGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCg                                                                                        >  2:186761/1‑149 (MQ=255)
                                                                    gggTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCCAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACg                                                                                >  8:352122/1‑149 (MQ=255)
                                                                     ggTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGTAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACg                                                                                <  2:271526/148‑1 (MQ=255)
                                                                       tCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGc                                                                                              >  4:361288/1‑132 (MQ=255)
                                                                       tCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGc                                                                                              <  3:361288/132‑1 (MQ=255)
                                                                         tGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGa                                                                            >  6:167348/1‑148 (MQ=255)
                                                                          gCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGAt                                                                           >  8:327392/1‑148 (MQ=255)
                                                                           cGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGAt                                                                           <  4:310848/147‑1 (MQ=255)
                                                                                aGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAGCa                                                                     >  2:171937/1‑148 (MQ=255)
                                                                                     gTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAGCAAAGCGt                                                               <  1:11601/149‑1 (MQ=255)
                                                                                           gTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAGCAAAGCGTGGGTTg                                                         >  8:183934/1‑149 (MQ=255)
                                                                                                   tACTCCAAGCGTTAATACCGCTTCGGCAGTCCAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAGCAAAGCGTGGGTTGTCGCCCtt                                                 <  7:352122/149‑1 (MQ=255)
                                                                                                    aCTCCAATCGTTAATACCGCATCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAGCaaa                                                                   <  8:373878/130‑1 (MQ=255)
                                                                                                     cTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCGTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATTGACGTAGATCAAGCAAAGCGTGGGTTGTTGCCCTTGc                                               >  2:14683/1‑149 (MQ=255)
                                                                                                     cTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTAGTTGGAACCGGGCTTTTTTTTTAGCCATGAACAAACCTGTCAGCAGTCAGTGGCCCCTTGAAGCATAGACGTACATCAAGCAAACCGTCGGTTGTAGcccatt                                                >  4:226420/1‑145 (MQ=255)
                                                                                                      tCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAGCAAAGCGTGGGTTGTCGCCCTTGcc                                              >  5:344099/1‑149 (MQ=255)
                                                                                                      ccgaaGCGGTATTAACGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTa                                                                              <  3:316750/115‑1 (MQ=255)
                                                                                                      ccgaaGCGGTATTAACGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTa                                                                              >  4:316750/4‑118 (MQ=255)
                                                                                                        gaaGCGGTATTAACGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCAGTCACCAGTCAGTGGCCGCTTGTAGCATAGACGTAGAACAAGCAAAGCGTCGGTTGTCGCCCTTGACgg                                            >  6:39139/2‑149 (MQ=255)
                                                                                                           gagTTAATACCGCTTCGGCATTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAGCAAAGCGTGGGTTGTCGCCCTTGCCGGGGt                                         <  6:321142/147‑1 (MQ=255)
                                                                                                           gCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAGCAAAGCGTGGGTTGTCGCCCTTGCCgggg                                          >  7:58954/1‑148 (MQ=255)
                                                                                                                aaTACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAGCAAAGCGTGGGTTGTCGCCCTTGCCGGGGTAATc                                     <  1:186761/148‑1 (MQ=255)
                                                                                                                    ccGCTTCGGCAGTCTAAAAAAGCCAGGGTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAGCAAAGCGTGGGTTGTCGCCCTTGCCGGGGTAATCACCtt                                <  3:61747/149‑1 (MQ=255)
                                                                                                                    acGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGTAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAGCAAAGCGTGGGTTGTCGCCCTTGCCGGGGTAATCACCtt                                <  5:39139/148‑1 (MQ=255)
                                                                                                                     cGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAGCAAAGCGTGGGTTGTCGCCCTTGCCGGGGTAATCACCTTg                               >  6:352158/1‑149 (MQ=255)
                                                                                                                      gCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAGCAAAGCGTGGGTTGTCGCCCTTGCCGGGGTAATCACCTTg                               >  2:42173/1‑148 (MQ=255)
                                                                                                                      gCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAGCAAAGCGTGGGTTGTCGCCCTTGCCGGGGTAATCACCTTGt                              >  1:75544/1‑149 (MQ=255)
                                                                                                                       cTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAGCAAAGCGTGGGTTGTCGCCCTTGCCGGGGTAATCACCTTGt                              >  2:231064/1‑148 (MQ=255)
                                                                                                                       cTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTATTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAGCAAAGCGTGGGTTGTCGCCCTTGCCGGGGTAATCACCTTGTc                             <  5:352158/149‑1 (MQ=255)
                                                                                                                         tCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAGCAAAGCGTGGGTTGTCGCCCTTGCCgggg                                          <  2:392049/134‑1 (MQ=255)
                                                                                                                         tCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAGCAAAGCGTGGGTTGTCGCCCTTGCCgggg                                          >  1:392049/1‑134 (MQ=255)
                                                                                                                          cGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAGCAAAGCGTgg                                                             <  8:109853/114‑1 (MQ=255)
                                                                                                                           ggCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAGCAAAGCGTGGGTTGTCGCCCTTGCCGGGGTAATCACCTTGTCAgt                          <  5:312765/148‑1 (MQ=255)
                                                                                                                           ggCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCATGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAGCAAAGCGTGGGTTGTCGCCCTTGCCGGGGTAATCACCTTGTCAgtg                         >  2:206606/1‑149 (MQ=255)
                                                                                                                             ctgTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTAAGCAGTCAGTGGGCGCTTGTCGCATAGACGTAGATCAAGCAAAGCGTGGGTTGTCGCCCTTGCCGGGGTAATGACCTTGTCAgtgt                        <  3:163283/146‑1 (MQ=255)
                                                                                                                              aGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAGCAAAGCGTGGGTTGTCGCCCTTGCCGGGGTAATCACCTTGTCAgtgtgt                      >  7:303888/1‑149 (MQ=255)
                                                                                                                               gTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAGCAAAGCGTGGGTTGTCGCCCTTGCCGGGGTAATCACCTTGTCAGtgtgtg                     <  5:167348/149‑1 (MQ=255)
                                                                                                                                tCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGACCCATAGACGTAGATCAAGCAAAGCTTGGGTTGTCGCCCTTGCAGGGGTAATCAAATTGTAAgtct                        >  4:163283/1‑143 (MQ=255)
                                                                                                                                       aaGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAGCAAAGCGTGGGTTGTCGCCCTTGCCGGGGTAATCACCTTGTCAGTGTGTGGATTTGtt             >  3:166254/1‑149 (MQ=255)
                                                                                                                                        aGCCCGGTTCCAAGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAGCAAAGCGTGGGTTGTCGCCCTTGCCgggg                                          >  6:370026/1‑119 (MQ=255)
                                                                                                                                        aGCCCGGTTCCAAGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAGCAAAGCGTGGGTTGTCGCCCTTGCCgggg                                          <  5:370026/119‑1 (MQ=255)
                                                                                                                                        aGCCCGGTTCCAAGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAGCAAAGCGTGGGTTGTCGCCCTTGCCGGGGTAATCACCTTGTCAGTGTGTGGATTTGttt            >  7:83061/1‑149 (MQ=255)
                                                                                                                                         gCCCGGTTCCAAGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAGCAAAGCGTGGGTTGTCGCCCTTGCCGGGGTAATCACCTTGTCAGTGTGTGGATTTGTTTc           >  7:374821/1‑149 (MQ=255)
                                                                                                                                          ccccGTTCCTTGGAACCGGGCTTTTATTTTAGCCTTTAACTAACCTGTAAGCAGTCACTGCACGCTTGTCGCATATACAAAAATCAATAAAAGCGAGGTATCTCGCACTTGCAGGGGTAAACACCTTGTCAGTGTGTGGATTTGtttaa          >  4:270103/1‑147 (MQ=255)
                                                                                                                                          cccGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAGCAAAGCGTGGGTTGTCGCCCTTGCCGGGGTAATCACCTTGTCAGTGTGTGGATTTGTTTCt          >  5:177370/1‑149 (MQ=255)
                                                                                                                                           ccGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAGCAAAGCGTGGGTTGTCGCCCTTGCCGGGGTAATCACCTTGTCAGTGTGTGGATTTGTTTCt          >  3:120622/1‑148 (MQ=255)
                                                                                                                                                   ttGGAACCGGGCTTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAGCAAAGCGTGGGTTGTCGCCCTTGCCGGGGTAATCACCTTGTCAGTGTGTGGATTTGTTTCTTACCATGa  >  1:174588/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
CTGCGTGGCGCTCTGCGTCGTGCTGGCTACGTCACCCGCGACGCTCGTGAAGTTGAGCGTAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAGCAAAGCGTGGGTTGTCGCCCTTGCCGGGGTAATCACCTTGTCAGTGTGTGGATTTGTTTCTTACCATGA  >  NC_002947/1503920‑1504214

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