Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,862,6610TG40.0% 49.2 / 16.8 56pseudogene (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 (10/23);  new base G (20/2);  total (31/25)
Fisher's exact test for biased strand distribution p-value = 7.78e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGTTGCGGGAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCT  >  NC_000913/3862513‑3862807
                                                                                                                                                       |                                                                                                                                                      
ccGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCCGTTCCGTAATTCCCACCAGCATGGCGGTTAAGGtt                                                                                                                                                       <  2:515570/150‑1 (MQ=255)
  gCCCATGGTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGGACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCAACAGCATGGCGGGTTAAGGtgcg                                                                                                                                                     >  2:330836/1‑148 (MQ=255)
     cATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGGTGggggg                                                                                                                                                >  1:524431/1‑150 (MQ=255)
         ttGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGGTGCggg                                                                                                                                                 >  1:292656/1‑146 (MQ=255)
         ttGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGGAATTCCCACCAGCATGGCGGATCAAGTTTGGGGaaac                                                                                                                                             >  2:181987/1‑148 (MQ=255)
           gCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGTTGCGGGAATCAGc                                                                                                                                         <  2:327315/151‑1 (MQ=255)
            cccACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGGTAAAGGTGCggggaaca                                                                                                                                            >  2:397941/1‑143 (MQ=255)
             ccaccaCACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGTTGCGGGAATCAGCaa                                                                                                                                       <  2:1191532/151‑1 (MQ=255)
                ccaccccaAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGTTGCGGGAATCAGCAAGcc                                                                                                                                    <  2:962054/151‑1 (MQ=255)
                 cacaccaAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGTTGCGGGAATCAGCAAGccc                                                                                                                                   <  1:458097/151‑1 (MQ=255)
                            gaTACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGGTGGGGGGAACAGCAAAGCCGCCACcttt                                                                                                                        >  1:908671/1‑149 (MQ=255)
                               aCATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGGTAAGGGTGCGGGGAACAGCAAACCCCCCCCtttttacc                                                                                                                    >  2:1061053/1‑146 (MQ=255)
                                            cATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGGTAAGGGTGCGGGGATCAACCAGGCCGCCACCTTTTACCCATTTTTCGGGGCg                                                                                                      >  2:458466/1‑151 (MQ=255)
                                             atggAGGCCGCCAGCCGGGCGTGACCCGCAAAAAGCAACGGTAAAATAAAAAGGAAGGTAAATCCCAGCGGTTCGTTAATCCCCCCCAGCATGGCGTTTAGGGTGGCGGGAATAAGCAAGCCCGCCACTTTTACCCGATTTTCCg                                                                                                           <  2:292813/142‑1 (MQ=255)
                                              ttGAGGCCGCCAGCACGGGGTGTACCGCAAAAAGCAAGGGAAAAATGAAAAGGAAGTAAATTCCCGGCGGTTCGTTAATCCCCACCAGCAGGGCGTTTAAGGTTGCGGGAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGt                                                                                                         <  1:853587/146‑1 (MQ=255)
                                               tGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAAGGTTCGGGAATCAGCAAACCCCCCCCCTTTTCCCCCTTTTTCGGGGCg                                                                                                      >  2:373652/1‑148 (MQ=255)
                                               tGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGGTGAGGGTGCGGGGATAAAGCAACCCGCCACCTTTAACAGCTTTTTTGGTGCgg                                                                                                     >  2:279650/1‑149 (MQ=255)
                                                 aGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGATTAGGGTGCGGGGAACCCCCAACCCCCCCCCTTTTCCCCGTTTTTCgggggggca                                                                                                    >  2:1245872/1‑143 (MQ=255)
                                                    ccgccAGCGGGGGGTTTACCGAAAACAGCAAGGGTAAAATAAACAGGAAGTAAATTTCCGGCGGTTGGATAATCCCCACCAGCAGGGGGTTTAGGGTGGGGGGAATCAGCAGCCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGa                                                                                               <  2:310957/150‑1 (MQ=255)
                                                    ccgccAGCACGGCGTGTCCCGAAAACAGCAACGGTGAAATGAACAGAAAGGTAAATTCCGGCGGTTCGGTAATCCCCACCAGCATGGCGGTTAAGGTGGGGGAAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTg                                                                                                <  1:752279/149‑1 (MQ=255)
                                                     cgccAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTTAGGGTGCGGGGATCAGCAAACCCGCCCCCTTTTACCCATTTTTCGGGGCGGGAGTg                                                                                                >  1:1090874/1‑148 (MQ=255)
                                                       ccAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAAGGTGCGGGGAACAGCAAACCCCCCCCTTTTTCCCGATTTTTCGGGGCGGCAGGGaa                                                                                              >  1:373718/1‑148 (MQ=255)
                                                              ggCGTGTACCCAAAACAACACGGGTAAAAAAAGCAGAAAGGAAAATCCCGGCGTTTCGTTAATCCCCACCAGCAGGGCGTTTAGGGTTGCGGAAATAAGCAACCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACAt                                                                                        <  2:339164/147‑1 (MQ=255)
                                                                  tgtACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGCAATTCCCACCAGCATGGCGGTTAAGATATCCGGAATCAGCAAACCCGCCCCCCTTACCCGCGTTTTCGTTGCgccggcgaa                                                                                              >  2:1109231/1‑129 (MQ=255)
                                                                      cgccAAACAGAAACGGTAAAATGAACAGGAGGTAAAATCCAGGGGGTTCGGTATTTCCCACGAGCATGGCGTTTAGGGTGGGGGAAATCACCAGGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAAt                                                                            <  1:180604/148‑1 (MQ=255)
                                                                       cGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGTTGCGGGAATCAGCAAACCCGCCCCTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAaaaaaa                                                                             >  1:756192/1‑149 (MQ=255)
                                                                        gCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGGTGGGGGAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGGGGGGCAGGGAAGTAAATCGCTAAAGaaa                                                                             >  2:525861/1‑148 (MQ=255)
                                                                          aaaCAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGGTGCGGGGATCAGCAAGCCCGCCACCTTTACCCGATTTTTCGGTGCGGCAGGGAAGGACATCGCCAAAGAAATAccc                                                                        >  2:152670/1‑151 (MQ=255)
                                                                          aaaCAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGGTGAGGGAAACAGCAAACCCGCCACCTTTACCCGATTTTCCGGTGCGGGAGTGAAGTTAATCGCCAaaaaaaaaccc                                                                        >  2:1239704/1‑146 (MQ=255)
                                                                            aCAGCAAGGGAAAAATAAAAAGGAAGGTAAATTCCGGGGTTGCGGTAATCCCCGCCAGCATGGCGTTTAGGGTTGGGGAAAACAGAAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCAc                                                                      <  1:1250134/151‑1 (MQ=255)
                                                                              aGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGGTGCGGGAATCAGCAAGCCCGCCCCCTTTACCCGAATTTTCGGGTCGGGAGTGAAGGAAATCCCTAAAAAAATAccccc                                                                      >  2:577019/1‑147 (MQ=255)
                                                                                 aagggTGAAATGAACGGGAAGGTAAATCCCAGCGGTTCGGTAATCCCCCCCAGCATGGCGGTTAAGTTTGCGGGAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCAc                                                                      <  1:1168770/143‑1 (MQ=255)
                                                                                 aaCGGTAAAATGAACAGGAAGGTAAATTCCGGCGGTTCGGTAATCCCCCCCAGCATGGCGGTTAGGGTTGGGGGAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCAc                                                                 <  1:1060477/151‑1 (MQ=255)
                                                                                    ggTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGGTGCGGGAATCAACAAGCCCGCCACCTTTTCCCCATTTTCCCGTGCGGCAGGGGAGGACATCGCTAAAGAAATGCCCACGGCAcc                                                                >  1:684871/1‑149 (MQ=255)
                                                                                    ggTAAAAAGAAAAGAAAGGTAAATTCCAGCGTTTGGGTAATCCCCCCCAGCATGGCGGTTAAGGTGGGGGGAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCAc                                                                 <  1:297867/148‑1 (MQ=255)
                                                                                          ataaaCAGGAAGGAAAATTCCAGCGTTCCGGTAATCCCCACCAGCATGGCGGTAAAGGTTGGGGAAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGAttt                                                        <  1:208735/148‑1 (MQ=255)
                                                                                              aCAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGTTGCGGGAATCAGCAAACCCCCCACTTTTACCCGATTTTTCGGGGCGGCCGTGAAGTAAATCGCTAAAAAAAAGCCCACGGCCCCCAAGAttttt                                                      >  1:1105495/1‑149 (MQ=255)
                                                                                              aCAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGGTGCGGGAATCAGCAAGCCCGCCCCTTTTACCCGATTTTCCGGGGGGGCAGTGAAGTAGATCGCTAAAGAAAAGCCCACGGCACCAAAGAttttt                                                      >  1:451020/1‑149 (MQ=255)
                                                                                                     gatAAATCCCAGCGGTTGGGTAATCCCCACCAGCGTGGCGTTTAGGGGGGGGGAAATAAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTAcc                                              <  2:463821/148‑1 (MQ=255)
                                                                                                         aaTTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTTAAGATGGGGGGATCTGCCAGCCCGCCACATTTACCCGATTTTTCGGGGCGGCAGGGAAGGACCTCGCTCAAAAAATGCCCACGGCAACAAAGATTTTTTAGTTAcc                                              >  2:551345/1‑146 (MQ=255)
                                                                                                         aaTTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGTTGCGGGGAGCCGCAAGCCCGCCCCCTTTACCCGATTTTCCGGTGCGGCGGTGAAGTACATCGCACAAAAAATACCCAAGGCACCAAAGATTTTTGAGTTACCCTg                                           >  1:1180869/1‑149 (MQ=255)
                                                                                                         aaTTCCAGCGGTTCGGTAATCCCCCCCAGCTTGGCGTTTAAGGTTGGGGGAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCa                                         <  1:340323/151‑1 (MQ=255)
                                                                                                         aaTTCCACCGGTGGGTTAATCCCCCCGAGCATGGCGGTAAAGGGGGCGGGAAACAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCa                                         <  2:738340/151‑1 (MQ=255)
                                                                                                           ttCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGGTTCGGGGATCAGCAAGCCCGCCACCTTTTCCCGATTTTCCGGTGCGGGAGTGAAGGACAACGCTCAAAAAATGCCCACGGCACCAAAGAATTTTTAGGTTACGTGCAgg                                       >  2:684781/1‑151 (MQ=255)
                                                                                                            tCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGTTGCGGGAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGGGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAggg                                      >  2:712136/1‑151 (MQ=255)
                                                                                                            tCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGGTGCGGGAATCAGCAAGCCCGCCCCTTTTTCCCCCTTTTTCGGGGCGGCGGTGAAGTACATCGCTAAAAAAATACCCACCGCCCCAAAGATTTTTTAGTTACCCTGCAggg                                      >  1:592698/1‑151 (MQ=255)
                                                                                                                 accGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGTTGCGGGGATCAGCAAACCCGCCCCTTTTACCCCATTTTTCGGTGCGGCAGTGGAGTACATCGCCAAAGAAAAGCCCACGGCCCCCAAGATTTTTGAGTTACCCCGCAggg                                      >  1:354752/3‑148 (MQ=255)
                                                                                                                    ggTTCGGTAATTCCCACCAGCATGGCGGTTAAGGGTGCGGGGATCAGCAAGCCCGCCCCTTTTACCCGATTTTTCGGGGCGGCAGTGAAGGACATCGCTAAAAAAATGCCCACGGCACCCAAGATTTTTGGGTTACCGTGCAGGGCaaa                                  >  1:750269/1‑149 (MQ=255)
                                                                                                                      ttCGGTAATTCCCACCAGCATGGCGGTTAAAGTTGCGGGAATCAGCAAACCCGCCCCCTTTACCCCATTTTCCGGTGCGGCCGTGAAGTACATCGCTAAAGAAATACCCACGGCACCAAAGATTTTTGAGTTACCGTGCCGGGCaaaa                                 >  2:1239143/1‑148 (MQ=255)
                                                                                                                       tCGGTAATTCCCACCAGCATGGCGGTTAAGGTTGCGGGGATCAGCAAACCCGCCACTTTTTACCGAATTTTCGGTGGGGCAGTGGAGTACATCCGTAAAAAAATGCCCACGGGCCCCAAGATTTTTTAGTTACCCTGCAGGGCAAAAACGc                             >  1:1062945/1‑151 (MQ=255)
                                                                                                                       tCGGTAATTCCCACCAGCATGGCGGTTAAGGTTGCGGGAATCAGCAAACCCGCCACTTTTACCCCATTTTCCGGTGCGGCAGTGAAGTACATCGCCCAAAAAATACCCACCGCACCAAAGATTTTTGAGGTACCGTGCAGGGCAAAACCGc                             >  2:684356/1‑151 (MQ=255)
                                                                                                                                ccccccAGCTTGGCGTTTAAGGTTGCGGGAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCtt                          <  1:581950/145‑1 (MQ=255)
                                                                                                                                ccccccAGCAGGGCGTTTAAGGTGGCGGAAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTc                         <  1:1245181/146‑1 (MQ=255)
                                                                                                                                ccccccAGCAGGGCGGTTAAGGTGGCGGGAAACAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCgg                      <  1:98228/149‑1 (MQ=255)
                                                                                                                                ccccccAGAATGGAGATTAGGGTTGGGGGAAGCAACAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTc                         <  1:330667/146‑1 (MQ=255)
                                                                                                                                 cccccaGCTTGGGGGTTAAGGTTGGGGGAACCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGaa                   <  1:488696/148‑1 (MQ=255)
                                                                                                                                           ggCGGTTAGGGTGGCGGAAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGAttt          <  1:1185677/150‑1 (MQ=255)
                                                                                                                                                   aaGTTTGCGGAAATCAGAAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCt  <  1:181886/151‑1 (MQ=255)
                                                                                                                                                       |                                                                                                                                                      
CCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGTTGCGGGAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCT  >  NC_000913/3862513‑3862807

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