Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,862,6610TG40.0% 38.9 / 27.8 67pseudogene (943/1617 nt)glvCpseudogene, arbutin specific enzyme IIC component of PTS;enzyme; Transport of small molecules: Carbohydrates, organic acids, alcohols; PTS system, arbutin‑like IIB component; PTS system, arbutin‑like IIC component
Reads supporting (aligned to +/- strand):  ref base T (15/24);  new base G (24/2);  total (41/26)
Fisher's exact test for biased strand distribution p-value = 9.34e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.20e-01

CCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGTTGCGGGAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGC  >  NC_000913/3862513‑3862803
                                                                                                                                                     |                                                                                                                                                     
ccGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGTTg                                                                                                                                                     <  2:65432/151‑1 (MQ=255)
      aTGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTTAGGGTGgggg                                                                                                                                                >  1:304337/1‑149 (MQ=255)
      aTGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGTTGCGGGaa                                                                                                                                             <  1:904400/151‑1 (MQ=255)
      aTGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGTTGCGGGaa                                                                                                                                             <  2:804427/151‑1 (MQ=255)
           gCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGTTGCGGGAATCAGc                                                                                                                                       <  2:223634/151‑1 (MQ=255)
                  acaccaAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGTTGTGGGAATCAGCAAGgcc                                                                                                                                 >  2:69452/1‑150 (MQ=255)
                     ccaAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGGTAAAGGTTCGGGGAGCAGAAAGAGCGcc                                                                                                                             >  2:85790/1‑150 (MQ=255)
                     ccaAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGATCGGTAATTCCCACCAGCATGGCGGTTTAGGGTGCGGGGATCAACAAGGCCTcc                                                                                                                             >  2:12024/1‑150 (MQ=255)
                       aaaGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGGAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTTAAGGTGCGGGAAACAGCAAAgccgcccc                                                                                                                           >  2:626041/1‑148 (MQ=255)
                          gagaTACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTTAGGGTGCGGGGATCAGCAAACCCCCCACcttt                                                                                                                      >  1:776671/1‑151 (MQ=255)
                             aTACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAAGTTGGGGGGATCAGCAAACCCCCCCCcttt                                                                                                                      >  1:45783/1‑148 (MQ=255)
                                   tACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGGAAATTCCAGCGGGTCGGTAATTCCCACCAGCATGGCGGTTAAAGGTGCGGGGAGGAACAAACCCCCCCCtttttctccagtt                                                                                                             >  2:1037453/1‑142 (MQ=255)
                                      ggacGACATGGGGCCCCCCAGCCGGGCGTGCCCCGAAAACAGCAACGGAAAAAAAAAAAGGAAGGAAAATCCCGGGGGTTCGGTACTCCCCCCCAGCTGGGGGGTTAGGTTGGGGGAAATAAGAAAGCCCGCCACTTTTACCCGAtt                                                                                                             <  2:976772/144‑1 (MQ=255)
                                      ggTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAAGGTGCGGGGATCAGCAAACCCCCCCCTTTTTCCCGCtttt                                                                                                           >  1:863402/1‑149 (MQ=255)
                                       gTCGACATGGAGGCCGCCAGCGCGGCGTGTCCCACAAACAGCAACGGTAAAATGAACAGGAAGGTAATTCCCAGCGGTTCGGTAATTCCCACCAGCATGGCGTTTAGGTTGGCGGGAAACAGCAAGCCCGCCACTTTTACCCGATTTTCCg                                                                                                        <  1:1171068/151‑1 (MQ=255)
                                        tCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGGTTCGGGGATTAGGAAAACCGCCACCTTTTACCCGTTTTcc                                                                                                         >  1:707100/1‑149 (MQ=255)
                                         cGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTTAAGGTGCGGGGAACAACAAACCCCCCCCTTTTTCCCCAtttt                                                                                                           >  1:566746/1‑146 (MQ=255)
                                          gACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGGTAAAGATTGGGGGAGCAGCAAACCCCCCCCCTTTTcccccgtttttggt                                                                                                      >  2:926714/1‑139 (MQ=255)
                                           acttggggcccgccAACACGGCGTTCCCCCCAAACAACAACGTGGAAAGGAACAGGAAGGTAAACCCCAGGGTTCCGGTAATCCCCCCCAGCTTGGCGGTTAGGTTGGGGGAAATAAGAACGCCCGCCACTTTTACCCGATTTTCCGGt                                                                                                      <  2:90360/140‑1 (MQ=255)
                                            cATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAAGGTGCGGGGATCACCAAACCCGCCACTTTTACCCCCTTTTCCggcg                                                                                                     >  1:1058488/1‑147 (MQ=255)
                                                   gccgccAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGGTTCGGGGATCAGCAAACCCGCCCCTTTTTCCCGATTTTCCCGGGCGGGAggg                                                                                             >  1:226635/1‑148 (MQ=255)
                                                     cgccAGCAGGGGGTGTACCGCAAACAAAAACGATAAAAGGAACAGGAAGAAAAATTCCAGGGTTTCGATATTCCCCACCAGCAGGGCGGTTAAGGTTGCGGAATTCAGAAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGaa                                                                                           <  1:1071373/150‑1 (MQ=255)
                                                         aGCACGGCGTTAACCGAAAACAGAAGCGGTGAAATAAACAGGAAGGAAAATTCCAGCGGTTGGGAAATCCCCCCCGTCAGGGCGTTTAAGGTGGCGGGAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTa                                                                                        <  1:1279273/149‑1 (MQ=255)
                                                         aGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAAGGTTCGGGGATCAGCAAACCCCCCCCTTTTTCCCCATTTTTCCGTGCgg                                                                                                  >  1:608603/1‑139 (MQ=255)
                                                            aCGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGGAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGGTTAGGGTGGGGGGAAACACAAAACCCCCCCTTTTTCCCCCTTTTTCGGGGGGGccggggagaac                                                                                       >  1:151502/1‑137 (MQ=255)
                                                              ggCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGATAAGGTTGTGGGAATTAGCCAGCCCGCCCCTTTTATCCGATTTTTCCGTTCGGCAGGGAAgga                                                                                        >  1:416442/1‑142 (MQ=255)
                                                              ggCGTGTACCGAAACCAGCAAGGGTGAAATGAACGGGAAGGAAAATCCCAGGGGTTCGTTATTCCCCCCCAGCGGGGCGGTTAGGGTGGCGGGAATCAGAAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCg                                                                                   <  2:622405/149‑1 (MQ=255)
                                                               gCGTGTACCGAAAACAGCAACGGAGAAATAAACAGAAAGGTAAATTCCGGGGGTTCGGTAATCCCCACCAGCATGGCGGTTAGGTTGGCGGAAATCAGAAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTa                                                                                <  1:255709/151‑1 (MQ=255)
                                                                cGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGGTGGGGGAAATAGCAAACCCCCCACTTTTTTCCGCTTTTTCGGGGCGGCGGGGGAGTACCACGCTa                                                                                >  1:529298/1‑150 (MQ=255)
                                                                 gtgacccGAAAACAGCACCGGTAAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATCCCCCCCAGTGGGGCGGTTAAGGTTGGGGGAACCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTaaa                                                                              <  1:384441/146‑1 (MQ=255)
                                                                    tACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGGTGCGGGGATCAGCAAGACCGCCCCTTTTTCCCGATTTTTCGGGGCGGGAGTGGAGGAAATCGCTaaa                                                                              >  2:343285/1‑148 (MQ=255)
                                                                      ccccAAACAGCAACGGGGAAGAAAGAAGGAAGTAAAATCCCAGCGGTTGTTAAACCCCCCCCGGAGGGTCGTTTAGGGGGGGGGAAACAAGCACGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGaa                                                                           <  2:1060169/146‑1 (MQ=255)
                                                                         cAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAAGGTGCGGGGATCAGCAAACCCCCCCCTTTTACCCGATTTTCCGGTGCGGCAGTGAAATACATCGCTAAAGAAAAGcc                                                                      >  2:696331/1‑151 (MQ=255)
                                                                           accagcagGGGAGAAAGGAAAGGGAAGAAAAGTCCCAGGGGTTCGTTATTCCCCCCCAGCAGGGGGTTGGGGGTGGGGAAAACCACCACGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGcc                                                                      <  2:331321/147‑1 (MQ=255)
                                                                           aaCAGCAACGGAAAAAGAAAGAGGAAGGAAAATTCCAGGGGTTCGTTAATCCCCACCAGCAGGGCGTTTAGGGTGGGGGAAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCa                                                                    <  2:776643/151‑1 (MQ=255)
                                                                            aCAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGGTTCGGGGATCAGCAAGCCCGCCACTTTTACCCGATTTTTCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATACCCAc                                                                   >  1:1251632/1‑151 (MQ=255)
                                                                                  aCGGTGAAATAAACGGGAAGGTAAATTCCGGCGTTTCGTTAATCCCCCCCGGCATGGCGGTTAGGGTTGGGGGAATCAGAAGGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGc                                                                <  2:137656/148‑1 (MQ=255)
                                                                                   cGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGTTGCGGGAAACAGCAAACCCGCCACTTTTAACCCATTTTCCCGGGCGGCAGGGGAGTAAAACGCTAAAAAAAAACCCACCGCAcc                                                            >  1:48815/1‑150 (MQ=255)
                                                                                    ggaaaaaTGAAGGGAAAGGTAACTCCCGGCGGTGCGGTAATCCCCACCAGCGGGGCGGTGAGGGTGGGGGAAACAAGCAACCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCa                                                           <  2:1209162/146‑1 (MQ=255)
                                                                                    ggTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGTTGCGGGGATCAGCAAACCCGCCCCTTTTACCCGATTTTCCGGGGCGGCAGGGAAGGACAACGCTTAAGAAAAGCCCACGGCAcc                                                            >  1:269675/1‑149 (MQ=255)
                                                                                    ggTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAAGTTGCGGGGATCAGCCAGCCCGCCCCCTTTTACCTAATTTTCCGTGCGGGCGGGAAGGACATCGCaaaaaaaaatccccagg                                                                 >  2:1200109/1‑136 (MQ=255)
                                                                                        aaaggaaaaggaAGTTAATCTCCAGCGGTTCGTAAATCCCCCCCAGATTGGGGGTAAAGGTGGGGGGAAAAAGAAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGa                                                       <  2:1006857/150‑1 (MQ=255)
                                                                                         aataaaagggAAGATAAATCCCCGGGGTTCGGAAATCCCCCCCAGCGTGGCGGTTAGGTTTGCGGAAATCAGCAACCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGAtt                                                     <  2:451738/143‑1 (MQ=255)
                                                                                         aaggacagggAGGGTAAATTCAGGCGGTTGGGTAACCCCCCCCGGCATGGCGGTAGAGGTGGGGGAAATCAGCAACCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGAtt                                                     <  2:637260/143‑1 (MQ=255)
                                                                                         aaTGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGGTAAAGTTGGGGGGATCAGGAAGCCCCCCCCTTTTTCCCGATTTTTCGGGGGGGGCGGGGAGTACATCGCTAAAGAAAAACCCACCGCACCAAAGaat                                                     >  2:420945/1‑149 (MQ=255)
                                                                                         aaTGAAAGGGAAGTTAAATCCCAGCGGTTCGGTATTCCCCCCCAGCAGGGCGTTAAAGGTGGCGGAAATAAGAAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGAtt                                                     <  1:800034/151‑1 (MQ=255)
                                                                                          aTGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAAGATTGGGGGAGCAGGAAGCCCGCCCCTTTTTCCCCGTTTTTCGGGGCGGCAGGGGAGAAAATCGCCAAAGAAAAACCCACCGCCCCCAAGAatt                                                    >  2:1063286/1‑151 (MQ=255)
                                                                                           tGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAAGGTGCGGGAATCAGCAAGCCCGCCACTTTTATCCGATTTTTCCGTGCGGCCGAGAAGTACATCGCTCAAGAAATACCCACCGCACCCAAAAtttt                                                   >  1:1151637/1‑151 (MQ=255)
                                                                                                aGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGGTGCGGGGATCAGCAAGCCCGCCACTTTTACCCCATTTTTCCGGGCGGCAGGGAAGAACATCGCTAAAGAAAAACCCACCGCACCAAAAAttttt                                                  >  1:544760/1‑147 (MQ=255)
                                                                                                aGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGGTGCGGGAAGCAGCAAACCCCCCCCTTTTTACCCATTTTCCCGTGCGGCAGGGAAAGACAACGCTAAAGAAATGCCCACGGCACCAAAAAttttttggt                                              >  1:62370/1‑147 (MQ=255)
                                                                                                        aaaTCCCAGGGGTTCGGTAATCCCCCCCAGCAGGGGGTTTAAGGTTGGGGGAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGt                                        <  1:518150/149‑1 (MQ=255)
                                                                                                         aaTTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTTAAGGTGCGGGGATCAGCAAGGCCCCCACTTTTACCCCGTTTTTCCGGGGGGGAGTGAAGGACATCGCTTAAAAAAAACCCACCGCACCCAAAAAttttttgttaacctgc                                       >  2:659775/1‑139 (MQ=255)
                                                                                                         aaTTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGTTGCGGGGATCAGCAAGCCCGCCACTTTTTCCCGATTTTCCGGTGGGGCAGTGAAGGACATCGCTAAAGAAATGCCCACCGCACCAAAGATTTTTTAGTTACCGt                                        >  1:765428/1‑148 (MQ=255)
                                                                                                         aaTTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGGTGCGGGAATAAGCAAGCCCCCCACTTTTATCCCATTTTCCCGGGCGGGCGTGGAGTAAATCGCTAAAGAAATGCCCACCGCACCCAAAATTTTTTAGTTACCgg                                        >  2:492924/1‑147 (MQ=255)
                                                                                                           ttCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGGTGCGGGAATCAGCAAACCCGCCACCTTTTCCCCATTTTCCGGTGCGGCCGGGAAGTAAATCGCTAAAGAAAAGCCCACGGCACCAAAGATTTCTTAGATACCGTGCCgg                                   >  2:1034315/1‑151 (MQ=255)
                                                                                                              cAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGTTGCGGGAATCAGCAAGCCCGCCACTTTTACCCCATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAAAGCCCACGGCACCCAAGATTTTTGAGTTACCGTGCAggg                                  >  1:834886/1‑149 (MQ=255)
                                                                                                                gCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAAGTTGCGGGGATCAGCAAACCCCCCCCTTTTTCCCGATTTTTCGGTGCGGCAGGGAAGGACATCGCCAAAGAAAAGCCCACGGCCCCCAAGATTTTTTAGTTACCGTGCAggg                                  >  2:617927/1‑147 (MQ=255)
                                                                                                                 cGGTTCGGTAATTCCCACCAGCATGGCGGTTAAAGGTGCGGGGATCAGCCAGCCCCCCACCTTTACCCGCTTTTCCGGGGGGGGAGTGGAGTACAACGCTAAAGAAAAGCCCACCGCACCAAAGAATTTTTAGTTTCCGTGCAgggggaa                              >  2:909522/1‑146 (MQ=255)
                                                                                                                  ggTTCGGTAATTCCCACCAGCATGGCGGTTAAGGTTGCGGGAATCAACAAACCCGCCCCCTTTTCCCCATTTTTCGGGGCGGGAGGGGAGGACCTCGCTAAAAAAAAGCCCACGGCACCAAAAAATTTTTAGGTAACGTGCAGGGGaaaa                             >  2:582215/1‑150 (MQ=255)
                                                                                                                    ttCGGTAATCCCCACCAGCAGGGGGTTTAGGGTTGGGGGATACAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAAcc                           <  1:926752/150‑1 (MQ=255)
                                                                                                                       ggTAATCCCCCCCAGAAGGGGGTTTAGGTTGGGGGAAACAAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCg                          <  2:812286/148‑1 (MQ=255)
                                                                                                                        gTAATTCCCACCAGCATGGCGGTTAAGGTTGCGGGAATAAGCAAGCCCGCCCCTTTTACCCGATTTTTCGGTGCGGCAGTGAAGTACATCGCCAAAGAAAAACCCACGGGCCCAAAGATTTTTTAGTTACCGTGCAGGGCAAAAACGc                         >  1:684348/1‑148 (MQ=255)
                                                                                                                           aTTCCCACCAGCATGGCGGTTTAGGTTGCGGGGATCAGCAAGCCCCCCACTTTTCCCGATTTTCCGGTGCGGGAGTGAAGGACATCGCTTAAGAAATGCCCACCGCACCCAAGATTTTTTAGGTAACGGGCAGGGCAAAAACGCCCTcc                    >  1:720082/1‑149 (MQ=255)
                                                                                                                           aTTCCCACCAGCATGGCGGTTAAGGTTGCGGGAATCAACAAGCCCCCCACTTTTACCCGCATTTTCGGTGGGGCAGTGAAGTAAATCGCTAAAGAAATGCCCACGGCACCCAAAATTTTTGTGTTACCGCGCAGGGCAAAAACGCCtttc                    >  2:709151/1‑148 (MQ=255)
                                                                                                                                  ccaGCATGGCGGTTAAGGGTGCGGGGATCAGCAAGCCCGCCACTTTTTCCCGATTTTCCGGGGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCCAAGATTTTTTAGTTACCGTGCCGGGCAAAACCGCCTTCCGGGaa               >  2:1046934/1‑148 (MQ=255)
                                                                                                                                         ggCGGTAAGGGTGGAGGGAATGGGGAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTc     <  1:645759/151‑1 (MQ=255)
                                                                                                                                             ggTTAGGGTTGCGGGAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGc  <  1:423088/151‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGTTGCGGGAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGC  >  NC_000913/3862513‑3862803

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