Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,039,245 A→G 9.2% intergenic (‑265/+29) PP_5290 ← / ← pyrE transposase/orotate phosphoribosyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,039,2450AG 9.2% 202.8 / 8.4 76intergenic (‑265/+29)PP_5290/pyrEtransposase/orotate phosphoribosyltransferase
Reads supporting (aligned to +/- strand):  ref base A (35/33);  new base G (5/2);  total (41/35)
Fisher's exact test for biased strand distribution p-value = 4.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.48e-01

AGGGGTGATTTCTCCTACGGTCTGTACTGCGCGAACAGTCAGAAGCGAGTTATGGCCCTCCTCCTCCCTTCAAGTCCTACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGCTGTTTCAGCTGTGGATCGTCAGCCAGGAACTCCAGCACCTGGGTCAGCGAAACAATGCTGACCACCGGGATAC  >  NC_002947/6039098‑6039392
                                                                                                                                                   |                                                                                                                                                   
aggGGTGATTTCTCCTACGGTCTGTACTGCGCGAACAGTCAGAAGCGAGTTATGGCCCTCCTCCTCCCTTCAAGTCCTACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTgaagaa                                                                                                                                                     <  2:311969/148‑1 (MQ=255)
  gggTGATTTCTCCTACGGTCTGTACTGCGCGAACAGTCAGAAGCGAGTTATGGCCCTCCTCCTCCCTTCAAGTCCTACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGt                                                                                                                                                  >  5:353024/1‑149 (MQ=255)
  gggTGATTTCTCCTACGGTCTGTACTGCGCGAACAGTCAGAAGCGAGTTATGGCCCTCCTCCTCCCTTCAAGTCCTACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGt                                                                                                                                                  >  4:198439/1‑149 (MQ=255)
   ggTGATTTCTCCTACGGTCTGTACTGCGCGAACAGTCAGAAGCGAGTTATGGCCCTCCTCCTCCCTTCAAGTCCTACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTc                                                                                                                                                 <  3:87559/149‑1 (MQ=255)
    gTGATTTCTCCTACGGTCTGTACTGCGCGAACAGTCAGAAGCGAGTTATGGCCCTCCTCCTCCCTTCAAGTCCTACCATATAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCg                                                                                                                                                >  5:95303/1‑149 (MQ=255)
            tcCTACGGTCTGTACTGCGCGAACAGTCAGAAGCGAGTTATGTCACTCCTCCTCCCTTCAAGTCCTACCATACAAGAGGGTTTACCAGCGAAAGGGCAAGATCAGTCGGCAACAAAAAATGCGACACTAGGAAGAAGTTCGCCttgta                                                                                                                                         >  5:201264/1‑145 (MQ=17)
             ccTACGGTCTGTACTGCGCGAACAGTCAGAAGCGAGTTATGGCCCTCCTCCTCCCTTCAAGTCCTACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCa                                                                                                                                       >  7:201775/1‑149 (MQ=255)
              cTACGGTCTGTACTGCGCGAACAGTCAGAAGCGAGTTATGGCCCTCCTCCTCCCTTCAAGTCCTACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCTTCTTCAGAGGGTCGCCtttttttg                                                                                                                                      <  1:115931/149‑4 (MQ=18)
              cTACGGTCTGTACTGCGCGAACAGTCAGAAGCGAGTTATGGCCCTCCTCCTCCCTTCAAGTCCTACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCttttt                                                                                                                                         <  6:166440/146‑1 (MQ=255)
              cTACGGTCTGTACTGCGCGAACAGTCAGAACCGAGTTATGGCCCTCCTCCTCCCTTCAAGTCCTACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCAt                                                                                                                                      >  2:163993/1‑149 (MQ=255)
                 cGGTCTGTACTGCGCGAACAGTCAGAAGCGAGTTATGGCCCTCCTCCTCCCTTCAAGTCCTACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTcc                                                                                                                                   >  2:63513/1‑149 (MQ=255)
                 cGGTCTGTACTGCGCGAACAGTCAGAAGCGAGTTATGGCCCTCCTCCTCCCTTCAAGTCCTACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTc                                                                                                                                    <  2:123383/148‑1 (MQ=255)
                  ggTCTGTACTGCGCGAACAGTCAGAAGCGAGTTATGGCCCTCCTCCTCCCTTCAAGTCCTACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCttttt                                                                                                                                         <  5:12673/142‑1 (MQ=255)
                  ggTCTGTACTGCGCGAACAGTCAGAAGCGAGTTATGGCCCTCCTCCTCCCTTCAAGTCCTACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCttttt                                                                                                                                         >  6:12673/1‑142 (MQ=255)
                    tCTGTACTGCGCGAACAGTCAGAAGCGAGTTATGGCCCTCCTCCTCCCTTCAAGTCCTACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTgg                                                                                                                                <  3:198439/149‑1 (MQ=255)
                     cTGTACTGCGCGAACAGTCAGAAGCGAGTTATGGCCCTCCTCCTCCCTTCAAGTCCTACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTggg                                                                                                                               >  8:227191/1‑149 (MQ=255)
                        tACTGCGCGAACAGTCAGAAGCGAGTTATGGCCCTCCTCCTCCCTTCAAGTCCTACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCa                                                                                                                            <  7:227191/149‑1 (MQ=255)
                          cTGCGCGAACAGTCAGAAGCGAGTTATGGCCCTCCTCCTCCCTTCAAGTCCTACCATATAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAg                                                                                                                           <  6:95303/148‑1 (MQ=255)
                          cTGCGCGAACAGTCAGAAGCGAGTTATGGCCCTCCTCCTCCCTTCAAGTCCTACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAg                                                                                                                           <  3:66571/148‑1 (MQ=255)
                           tGCGCGAACAGTCAGAAGCGAGTTATGGCCCTCCTCCTCCCTTCAAGTCCTACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGa                                                                                                                         >  4:124083/1‑149 (MQ=255)
                             cgcgAACAGTCAGAAGCGAGTTATGGCCCTCCTCCTCCCTTCAAGTCCTACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGAt                                                                                                                        <  6:354121/148‑1 (MQ=255)
                             cgcgAACAGTCAGAAGCGAGTTATGGCCCTCCTCCTCCCTTCAAGTCCTACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGAt                                                                                                                        <  6:140063/148‑1 (MQ=255)
                                gAACAGTCAGAAGCGAGTTATGGCCGTCCTCCTCCCTTCAAGTCCTACCATTCAAGCGGGGTTAGCCGGGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACGCTCTGAAGAAGGTCGCCTTTTTCATTCCTgg                                                                                                                                <  3:345491/137‑1 (MQ=25)
                                gAAAAGTCAGAAGCGAGTTATGGCCCTCCTCCTCCCTTCAAGTCCTACCATACAAGCGCGTTTACCCGCGAAAGCGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTCAAGAAGATCGCCTTTTTCATTCCTgg                                                                                                                                >  4:345491/1‑137 (MQ=25)
                                  aCAGTCAGAAGCGAGTTATGGCCCTCCTCCTCACTTCAAGTCCTACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATc                                                                                                                  <  8:145439/149‑1 (MQ=255)
                                    aGTCAGAAGCGAGTTATGGCCCTCCTCCTCCCTTCAAGTCCTACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATccc                                                                                                                >  6:124484/1‑149 (MQ=255)
                                    aGTCAGAAGCGAGTTATGGCCCTCCTCCTCCCTTCAAGTCCTACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCAAAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATccc                                                                                                                >  3:244145/1‑149 (MQ=255)
                                       cAGAAGCGAGTTACGGCCCTCCTCCTCCCTTCAAGTCCTACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAACAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTa                                                                                                             <  8:98164/149‑1 (MQ=255)
                                         gAACCGAGTTATGGCCCTCCTCCTCCCTTCAAGTCCTACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACt                                                                                                           <  1:163993/149‑1 (MQ=255)
                                                 ttATGGCCCTCCTCCTCCCTTCAAGTCCTACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGACGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGCTCCCGTACTGCGccccg                                                                                                   >  6:317748/1‑147 (MQ=255)
                                                 ttATGGCCCTCCTCCTCCCTTCAAGTCCTACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCgg                                                                                                   >  7:187056/1‑149 (MQ=255)
                                                       cccTCCTCCTCCCTTCAAGTCCTACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGcc                                                                                             >  3:296093/1‑149 (MQ=255)
                                                        cctcctcctCCCTTCAAGTCCTACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCt                                                                                            <  3:355458/149‑1 (MQ=255)
                                                            ctcctccCTTCAAGTCCTACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGa                                                                                         <  3:124083/148‑1 (MQ=255)
                                                            ctcctccCTTCAAGTCCTACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGa                                                                                         <  7:87056/148‑1 (MQ=255)
                                                               ctccCTTCAAGTCCTACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACa                                                                                       <  5:127175/147‑1 (MQ=255)
                                                               ctccCTTCAAGTCCTACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGc                                                                                     >  2:70979/1‑149 (MQ=255)
                                                                 cccTTCAAGTCCTACCATACAAGCGGGCTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCg                                                                                   >  2:266667/1‑149 (MQ=255)
                                                                  ccTTCAAGTCCTACCATACAAGCGGGTTTACCCGGGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCgg                                                                                  >  2:322706/1‑149 (MQ=255)
                                                                  ccTTCAAGTCCTACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGcc                                                                                    <  8:101338/147‑1 (MQ=255)
                                                                  ccTTCAAGTCCTACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCgg                                                                                  >  7:125104/1‑149 (MQ=255)
                                                                      cAAGTCCTACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAATAAGGTCGCCTTTTTCATTCGTGGGCCAGGATGAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCa                                                                                <  2:77136/147‑1 (MQ=255)
                                                                             tACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTa                                                                                                 <  5:374309/123‑1 (MQ=255)
                                                                             tACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTa                                                                                                 >  6:374309/1‑123 (MQ=255)
                                                                               ccATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGCTGttt                                                                      >  5:145834/1‑148 (MQ=255)
                                                                                 aTACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGCTGTTTCAg                                                                   >  6:376298/1‑149 (MQ=255)
                                                                                  tACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGCTGTTTCAGc                                                                  >  1:9121/1‑149 (MQ=255)
                                                                                  tACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGCTGTTTCAGc                                                                  <  8:201775/149‑1 (MQ=255)
                                                                                   aCAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGCTGTTTCAGCt                                                                 >  6:182321/1‑149 (MQ=255)
                                                                                         gggTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGCTGTTTCAGCTGTGGAt                                                           <  5:182321/149‑1 (MQ=255)
                                                                                         gggTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGCTGTTTCAGCTGTGGAt                                                           >  2:279985/1‑149 (MQ=255)
                                                                                          ggTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGCTGTTTCAGCTGTGGATc                                                          <  4:236058/149‑1 (MQ=255)
                                                                                                 ccGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTGATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCGAGACAGCCGGCAGATGCTGTTTCAGCTGGTGATCGTCAGc                                                    >  6:174567/1‑148 (MQ=255)
                                                                                                 ccGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGCTGTTTCAGCTGTGGATCGTCAGc                                                    >  8:245793/1‑148 (MQ=255)
                                                                                                  cgcgAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGCTGTTTCAGCTGTGGATCGTCAGc                                                    >  5:80470/1‑147 (MQ=255)
                                                                                                    cgAAAGGGCCAGATCAGGCGGCAAAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGCTGTTTCAGCTGTGGATCGTCAGCCAgg                                                <  4:244145/149‑1 (MQ=255)
                                                                                                         gggCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGCTGTTTCAGCTGTGGATCGTCAGCCAGGAACTc                                           <  6:353024/149‑1 (MQ=255)
                                                                                                         gggCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGCTGTTTCAGCTGTGGATCGTCAGCCAGGAACTc                                           <  2:9121/149‑1 (MQ=255)
                                                                                                            ccAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGCTGTTTCAGCTGTGGATCGTCAGCCAGGAACTCCAg                                        >  1:289514/1‑149 (MQ=255)
                                                                                                            ccAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGCTGTTTCAGCTGTGGATCGTCAGCCAGGAACTCCAg                                        >  7:11651/1‑149 (MQ=255)
                                                                                                                   aGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGCTGTTTCAGCTGTGGATCGTCAGCCAGGAACTCCAGCACCTgg                                 <  1:63513/149‑1 (MQ=255)
                                                                                                                         gaatgaaaaaGGCGACCTTCTTCAGAGGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCAGGTAGGCCTCGACAGCCGGCAGATGCTGTTTCAGCTGTGGATCGTCAGCCAGGAACTCCAGCACCTGGGTCAGc                           >  8:18229/6‑149 (MQ=38)
                                                                                                                           atgaaaaaGGCGACCTTCTTCAGAGGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGCTGTTTCAGCTGTGGATCGTCAGCCAGGAACTCCAGCACCTGGGTCAGCg                          >  4:239668/4‑148 (MQ=39)
                                                                                                                            aaaaaaaGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGAAGGCCTCGACAGCCGGCAGATGCTGTTTCAGCTGTGGATCGTCAGCCAGGAACTCCAGCACCt                                   >  8:245491/1‑138 (MQ=255)
                                                                                                                            aaaaaaaGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGAGCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGCTGTTTCAGCTGTGGATCGTCAGCCAGGAACTCCAGCACCt                                   <  7:245491/138‑1 (MQ=255)
                                                                                                                             aaaaaaGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGCTGTTTCAGCTGTGGATCGTCAGCCAGGAACTCCAGCACCTGGGTCAGCGaaa                       <  1:128404/149‑1 (MQ=255)
                                                                                                                               aaaaGGCGACCTTCTTCAGAGGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGCTGTTTCAGCTGTGGATCGTCAGCCAGGAACTCCAGCACCTGGGTCAGCGAAAc                      >  3:49212/1‑148 (MQ=255)
                                                                                                                               aaaaGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGCTGTTTCAGCTGTGGATCGTCAGCCAGGAACTCCAGCACCTGGGTCAGCGAAAc                      >  2:223557/1‑148 (MQ=255)
                                                                                                                                aaaGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGCTGTTTCAGCTGTGGATCGTCAGCCAGGAACTCCAGCACCTGGGTCAGCGAAACaa                    >  6:171685/1‑149 (MQ=255)
                                                                                                                                 aaGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGCTGTTTCAGCTGTGGATCGTCAGCCAGGAACTCCAGCACCTGGGTCAGCGAAACAAt                   <  7:189911/149‑1 (MQ=255)
                                                                                                                                  aGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGCTGTTTCAGCTGTGGATCGTCAGCCAGGAACTCCAGCACCTGGGTCAGCGAAACAATg                  <  7:245793/149‑1 (MQ=255)
                                                                                                                                   ggCGACCTTCTTCAGAGGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGCTGttt                                                                      <  2:240126/96‑1 (MQ=38)
                                                                                                                                   ggCGACCTTCTTCAGAGGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGCTGttt                                                                      >  1:240126/1‑96 (MQ=38)
                                                                                                                                   ggCGACCTTCTTCAGAGGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTACGCCACGACAGCCGGCAGATGCTGTTTCAGCTGTGGATCGTCCGCCAGGAACTCCAGCACCTGGGTCAGCCAAACAATGc                 >  4:103596/1‑149 (MQ=255)
                                                                                                                                   ggCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGCTGTTTCAGCTGTGGATCGTCAGCCAGGAACTCCCGCACCTGGGTCAGCGATACAATGc                 >  4:301304/1‑149 (MQ=255)
                                                                                                                                       aCCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGCTGTTTCAGCTGTGGATCGTCAGCCAGGAACTCCAGCACCTGGGTCAGCGAAACAATGCTGac             <  5:376298/149‑1 (MQ=255)
                                                                                                                                          ctctGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGCTGTTTCAGCTGTGGATCGTCAGCCAGGAACTCCAGCACCTGGGTCAGCGAAACAATGCTGaccac          >  7:309828/1‑149 (MQ=255)
                                                                                                                                               aagaagGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGCTGTTTCAGCTGTGGATCGTCAGCCAGGAACTCCAGCACCTGGGTCAGCGAAACAATGCTGACCACCg        <  5:171685/146‑1 (MQ=255)
                                                                                                                                               aagaagGTCGCCTTTTTCATACCTGGGCCAGGATCAGATCACGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGCTGTTTCAGCTGTGGATCGTCAGCCAGGAACTCCAGCACCTGGGTCAGCGAAACAATGCTGACCACCGGGa     >  1:25153/1‑149 (MQ=255)
                                                                                                                                                agaagGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGCTGTTTCAGCTGTGGATCGTCAGCCAGGAACTCCAGCACCTGGGTCAGCGAAACAATGCTGACCACCGGGa     >  3:325087/1‑148 (MQ=255)
                                                                                                                                                agaagGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGCTGTTTCAGCTGTGGATCGTCAGCCAGGAACTCCAGCACCTGGGTCAGCGAAACAATGCTGACCACCGGGa     >  3:330470/1‑148 (MQ=255)
                                                                                                                                                  aagGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGCTGTTTCAGCTGTGGATCGTCAGCCAGGAACTCCAGCACCTGGGTCAGCGAAACAATGCTGACCACCGGGATAc  <  6:145834/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
AGGGGTGATTTCTCCTACGGTCTGTACTGCGCGAACAGTCAGAAGCGAGTTATGGCCCTCCTCCTCCCTTCAAGTCCTACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGCTGTTTCAGCTGTGGATCGTCAGCCAGGAACTCCAGCACCTGGGTCAGCGAAACAATGCTGACCACCGGGATAC  >  NC_002947/6039098‑6039392

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