Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,862,6380TC44.0% 115.3 / 14.9 75pseudogene (966/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 (32/10);  new base C (0/33);  total (32/43)
Fisher's exact test for biased strand distribution p-value = 1.25e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GGTTAAGGTTGCGGG‑AATCAGCAA‑GCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCG  >  NC_000913/3862488‑3862786
                                                                                                                                                      |                                                                                                                                                       
cGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAAtt                                                                                                                                                         <  1:536697/151‑1 (MQ=255)
  gTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTcc                                                                                                                                                       >  1:63958/1‑151 (MQ=255)
     aaaCCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCccac                                                                                                                                                    >  2:890165/1‑151 (MQ=255)
       aCCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCccacc                                                                                                                                                   >  1:1238029/1‑150 (MQ=255)
         cTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGc                                                                                                                                                >  2:735537/1‑151 (MQ=255)
           aGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCAt                                                                                                                                              <  1:1169760/151‑1 (MQ=255)
           aGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCAt                                                                                                                                              <  2:1087141/151‑1 (MQ=255)
            gTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATg                                                                                                                                             <  2:1286628/151‑1 (MQ=255)
               aaTCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑g                                                                                                                                         >  1:923133/1‑151 (MQ=255)
                   aGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑ggat                                                                                                                                      >  2:1156843/1‑148 (MQ=255)
                   aGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GGttt                                                                                                                                     >  1:1256928/1‑150 (MQ=255)
                        tCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGGTTAAgtt                                                                                                                                 >  2:538194/1‑149 (MQ=255)
                           gCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GGTTAAGGTTGCg                                                                                                                             <  2:780113/151‑1 (MQ=255)
                            cccATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GGTTAAGGTTGCgg                                                                                                                            <  1:591475/151‑1 (MQ=255)
                                tGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GGTTTAGGGTGCGGG‑GAt                                                                                                                       >  1:1274813/1‑151 (MQ=255)
                                       caccacaccaAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GGTTAAGGTTGCGGG‑AATCAGCAA‑g                                                                                                               <  1:1270710/151‑1 (MQ=255)
                                       caccacaccaAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GGTTAAGGTTGCGGG‑AATCAGCAA‑g                                                                                                               <  1:1007127/151‑1 (MQ=255)
                                           acaccaAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GGTTAAGGTTGCGGG‑AATCAGCAA‑GCCCg                                                                                                           <  2:1077240/151‑1 (MQ=255)
                                               caAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GGGTAAGGGTGGGGG‑GAATAGCaaagccccccc                                                                                                        >  2:905176/1‑142 (MQ=255)
                                                  agagATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GGGTTAGGGTGCGGG‑GATCAGCAA‑ACCCCCCCCttt                                                                                                    >  1:283547/1‑151 (MQ=255)
                                                    agaTACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GGTTAAGGGTGGGGG‑GATCAGCAA‑ACCCCCCCCttttt                                                                                                  >  1:862204/1‑150 (MQ=255)
                                                        aCATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GGTTTAAGGTGCGGG‑GATCAGCAA‑AACCCCCACTTTTAcccc                                                                                              >  1:854569/1‑150 (MQ=255)
                                                               ggTCGATATGGAGGCCCCCAGCGGGGCGTGTCCCGCAAACAGCAACGGTGAAATAAACAGGAAGGTAAATCCCAGCGGTTCGGTAATCCCCCCCAGCTTGGC‑GGTTAGGGTTGGGGA‑AATCAGAAA‑GCCCGCCACTTTTACCCGATTTTcc                                                                                       <  2:1196055/151‑1 (MQ=255)
                                                               ggTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGAGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GGATAAGGGTGCGGG‑GAAAAGGAA‑GCCCCGCCCTTTTaaaccacttttc                                                                                       >  1:919598/1‑140 (MQ=255)
                                                                gTCGACATGGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAAGGGTAAAATGAACAGGAAGATAAATTCCAGCGGTTCGGTAATCCCCACCAGCAGGGC‑GTTTAAGGTTGGGGA‑AATCAGCAA‑GCCCGCCACTTTTACCCGATTTTCCg                                                                                      <  1:408879/151‑1 (MQ=255)
                                                                  cGACATTGAGGCCGCCAGCGCGGCGTGTACCGCAAACAGCAAGGGTAAAATAAACAGGAAGGTAATTCCCAGCGGTTCGGTAATCCCCCCCAGCATGGC‑GGTTAAGGTGGCGGG‑AATCAGCAA‑GCCCGCCACTTTTACCCGATTTTCCGGt                                                                                    <  1:1204094/151‑1 (MQ=255)
                                                                    aCATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGGTTAAGGTTGGGGG‑AACAACAAA‑CCCCGCACCTTTTTCCCGATTTTCCgggg                                                                                   >  1:325054/1‑149 (MQ=255)
                                                                    aCATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GGTTAAAGTTGCGGG‑GATTAGCAA‑ACCCCCCCCCTTTACCCGAtttttccgggc                                                                                  >  2:993973/1‑144 (MQ=255)
                                                                       tggAGGCCGCCAGCCCGGCGTGTCCCGAAAACAGAAAGGGAAAAATGAACAGGAAGGAAAATCCCGGCGTTTCGTTATTCCCCCCCACCATGGC‑‑GTTTAGGTTGGGGGAAATCAGCAA‑GCCCGCCACTTTTACCCGATTTTCCGGTGCg                                                                                 <  2:1131573/147‑1 (MQ=255)
                                                                        tGAGGCCGCCAGCCCGGCGTGTACCGAAAACAGCAACGGTAAAATGAACAGGAAGGTAAATCCCAGCGGTTCGGTAATCCCCACCAGCAGGGC‑GTTTAAGGTGGGGGG‑AATCAGCAA‑GCCCGCCACTTTTACCCGATTTTCCGGTGCGGc                                                                               <  1:753798/150‑1 (MQ=255)
                                                                          aGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGGTAAAGGTTGGGGG‑AA‑CAGCAA‑ACCCCCCACCTTTTCCCCATTTTTCTGGGCGGGAgg                                                                            >  2:690807/1‑150 (MQ=255)
                                                                             ccgccAGCAGGGCGTGTCCCGCAAACAGCAAGGGTGAAAGGAACAGGAAGGTAAATCCCAGCGGTTCGGTAATCCCCACCAGCAGGGC‑GGTTAAGGTTGCGGG‑AATCAGCAA‑GCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGaa                                                                         <  1:805470/151‑1 (MQ=255)
                                                                              cgccAGCACGGCGTGTACCGAAAACAGCAACGGTAAAATGAACAGGAAGGTAAATCCCAGCGGTTCGGTAATCCCCACCAGCAGGGC‑GGTTAAGGTGGGGGA‑AATCAGAAA‑GCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAg                                                                        <  1:627436/151‑1 (MQ=255)
                                                                               cccAGCACGGGGTGCACCGAAAACAGCAACGGTAAAATAAACAGGAAGGAAAATCCCGGGGTTTCGGTAATCCCCCCCAGCATGGC‑GGTTAAGGTGGGGGG‑AACCAGAAG‑GCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGt                                                                       <  2:282843/150‑1 (MQ=255)
                                                                                ccAGCAGGGCGTGTCCCGCAAACAAAAACGGTAAAATGAACGGGAAGTAAAATTCCGGGGGTTCGGTAATCCCCACAAGCATGGC‑GGTTAGGGTTGGGGG‑AATCAGAAA‑GCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTa                                                                      <  2:644250/151‑1 (MQ=255)
                                                                                  aGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GGTTAAAGGTTCGGG‑GAACAGCAA‑GCCCCCCCCTTTTTCCCGCTTTTTCGGGGCgggcgggaaggac                                                                     >  2:1339581/1‑139 (MQ=255)
                                                                                   ccACGGCGTGTCCCGAAAACAGCAACGGTAAAATAAACAGAAAGATAAATCCCAGCGTTTCGTTAATCCCCACCAGCATGGC‑GTTTAAGGTGGCGGA‑AATCAGCAA‑GCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACAt                                                                   <  1:512051/150‑1 (MQ=255)
                                                                                    ccgggggtcacccGCAAACAGCACCGGTGAAATAAAAAGGAAGTTAAATCCAGGCGGTTCGTTAATCCCCACCAGCGTGGC‑GGTTAGGGTGGGGGG‑AATGAGAAG‑GCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATc                                                                  <  2:60556/140‑1 (MQ=255)
                                                                                    cACGGCGTCCCCCGCAAACAGCAGGGTAGAAATGAACAGGAAGTTAATTTCCAGCGTTTCGGTAATCCCCCCCAGCATGGC‑GTTTAGGGTGGGGGA‑AATCAGAAG‑GCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATc                                                                  <  2:16672/151‑1 (MQ=255)
                                                                                       ggCGTGTCCCGCAAACAGCAACGGTAAAATGAACAGGAAGGAAAATCCCAGCGGTTCGGTAATCCCCCCCAGCAGGGC‑GGTTAAGGTTGCGGG‑AATCAGCAA‑GCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCt                                                               <  1:15940/151‑1 (MQ=255)
                                                                                       ggCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GGTTAAGGTTGCGGG‑GATCAGCAA‑GCCCGCCACCTTTTCCCGATTTTCCGGGGCCGCAGTGGAGTACATCGCt                                                               >  1:95756/1‑151 (MQ=255)
                                                                                        gcttgtCCCGAAAACAGCAGGGGTAAAATAAACGGGAAGGTAATTTCCGGCGTTTCGGAAATCCCCCCCAGCAGGGC‑GTTTAGGTTGGCGGG‑AATCAGAAA‑CCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTa                                                              <  2:498750/148‑1 (MQ=255)
                                                                                          gtgtACCGAAAACAGCAAGGGTAAAATGAACAGGAAGGAAAATCCCGGGGGTTGGGAAATCCCCACCAGCAGGGG‑GGTTAGGGTGGGGGA‑AATCAGCAA‑GCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTaaa                                                            <  1:614442/151‑1 (MQ=255)
                                                                                            gtACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GGTTAAGGGTGCGGG‑AATCAGCAA‑GCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTaa                                                             >  1:1150658/1‑148 (MQ=255)
                                                                                            gtACCGCAAACAGCAACGGTAAAATGAACAGGAAGGTAAATCCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GGTTAAGGTTGCGGG‑AATCAGCAA‑GCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAg                                                           <  1:1302270/150‑1 (MQ=255)
                                                                                             tACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GGTTAAGGGTGCGGG‑AATCAGCAA‑GCCCGCCACTTTTACCCGATTTTCCGGGGCGGCAGTGAAGTACATCGCTAAAGa                                                          >  2:682324/1‑150 (MQ=255)
                                                                                             acccGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATCCCGGCGGTTCGGTAATCCCCCCCAGCATGGC‑GGTTAAGGTGGGGGA‑AATAAGCAA‑GCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGaa                                                         <  1:365014/149‑1 (MQ=255)
                                                                                                   aaaCAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GGTTAAAGTTTCGGG‑GATCAGCAA‑GCCCCCCCCTTTTTCCCCGTTTTTCGGGGCGGCCGGGAAGTACATCGCTAAAAAAAAGccc                                                   >  2:1213269/1‑151 (MQ=255)
                                                                                                   aaaCAGCAACGGAAAAATGAACAGGAAGGAAAATCCCAGCGGTTCGGAAATCCCCCCCAGCAGGGC‑GGTTAGGGTTGCGGG‑AATCAGCAA‑GCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGccc                                                   <  2:95573/151‑1 (MQ=255)
                                                                                                     aCAGCAACGGTAAAATGAACAGGAAGGTAAATCCCAGCGGTTCGGTAATCCCCCCCAGCAGGGC‑GTTTAAGGTTGCGGG‑AATCAGCAA‑GCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCAc                                                 <  1:267746/151‑1 (MQ=255)
                                                                                                        gCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GGTTAAAGTTGCGGG‑AATCAGCAA‑GGCCCCCACCTTTACCCGATTTTTCGGTGCGGGAGTGAAGGACATCGCCAAAGAAAATCCCCCGGc                                              >  2:840727/1‑151 (MQ=255)
                                                                                                          aaCGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GGTTAAAGTTGCGGG‑GATCAGCAA‑GCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGGGAAGTACATCGCTAAAGAAATGCCCACGGCAc                                            >  2:1008201/1‑151 (MQ=255)
                                                                                                           aCGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GGTTAAAGGTGCGGG‑GATCAGCAA‑GCCCGCCACTTTTACCCGATTTTTCCGGGCGGGAGTGAAGTACATCGCCAAAGAAATGCCTACGGGAcc                                           >  2:152042/1‑151 (MQ=255)
                                                                                                            cGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GGTTAAAGGTGCGGG‑GATCAGCAA‑GCCCGCCACCTTTAACCGGTTTTTCGGGGCGGCAGTGAAGTACATCGCCAAAAAAAAGCCCCCGGCAccc                                          >  2:385434/1‑150 (MQ=255)
                                                                                                             ggtaaaagaaaCAGAAAGGAAAATCCCAGCGGTTCGGTAATCCCCCCCAGCAGGGG‑GTTAAAGGTTGGGGG‑AACAAGCAA‑GCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCaa                                         <  1:594197/142‑1 (MQ=255)
                                                                                                              gtaaaaTGAACAGGAGGGTAAATTCCAGCGGTTCGGTAATCCCCACCAGCATGGC‑GGTTAAGGTTGCGGG‑AATCAGAAA‑GCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCaaa                                        <  1:587429/148‑1 (MQ=255)
                                                                                                              gtaaaaTGAACAGGAAGGAAAATCCCAGCGGTTCGGTAATCCCCCCCAGCAGGGC‑GGTTAAGGTTGCGGG‑AACCAGCAA‑GCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCaaa                                        <  1:1314089/148‑1 (MQ=255)
                                                                                                              gTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GGTTAAGGTTGCGGG‑AATCAGCAA‑GCCCGCCACTTTTACCCGATTTTTCGGTGCGGCAGTGAAGTACATCGCTAAAAAAATGCCCACCGCACCaa                                         >  1:696656/1‑150 (MQ=255)
                                                                                                               aaaaagaaaCAGGAAGTTATTCTCCGGCGGTTCGAAATTCCCCCCCAGCGGGGC‑GGTTAAGGTTGCGGA‑AATAAGAAG‑GCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAg                                       <  2:591288/144‑1 (MQ=255)
                                                                                                                   aggAACGGAAAGGTAAATCCCAGGGTTTGGGAAATCCCCCCCAGCTGGGC‑GTTTAGGGTGGGGGA‑AATCAGCAA‑GCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGAttt                                   <  2:290088/149‑1 (MQ=255)
                                                                                                                    aaaaaaGGAAGGTAAATCCCAGCGGTTGGGTAATCCCCCCCAGAAGGGC‑GGTTAAGGTTGGGGA‑AACCAGCAA‑GCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGAtttt                                  <  1:291197/146‑1 (MQ=255)
                                                                                                                     gAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GGTTAAGGGTGCGGG‑AATCAGCAA‑GCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGGGAAGTACATCGCTAAAAAAATGCCCACGGCACCAAAGAttttt                                 >  1:812924/1‑151 (MQ=255)
                                                                                                                     aaaCAGGAAGGAAAATTCCGGGGTTTCGATAATTCCCCCCAGCATGGC‑GGTAAGGGTGGCGGA‑AAACAGAAA‑GCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGAttttt                                 <  2:516422/150‑1 (MQ=255)
                                                                                                                        aaGGAAGTAAAATCCCAGGGGTTCGGTAATCCCCCCCAGCATGGC‑GGTTAAGGTGGCGGG‑AATCAGCAA‑GCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAg                              <  1:479922/150‑1 (MQ=255)
                                                                                                                           gAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GGTTAAGGTTGCGGG‑GATCAGCCA‑GCCCCCCACTTTTTCCCGATTTTCCGGGGGGGCAGTGGAGTACAACGCTAAAAAAAAGCCCACCGCACCCAAAATTTTTTAGtt                            >  2:885376/1‑150 (MQ=255)
                                                                                                                              ggTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GGTTAAGGTTTCGGG‑GATCAGCAA‑GCCCGCCACTTTTTCCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATTCCCACGGCACCCAAGAATTTTTAGTTAccc                        >  1:363959/1‑150 (MQ=255)
                                                                                                                                 aaaTCCCAGCGTTTCGTTAATCCCCCCCGGCAGGGC‑GGTTAGGGTGGCGGA‑AACCAGCAA‑GCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGc                     <  2:1354125/151‑1 (MQ=255)
                                                                                                                                  aatcccggggtttCGATAATCCCCACCAGCAGGGG‑GTTTAGGGTGGGGGA‑AACCAGCAA‑CCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCa                    <  2:1237757/140‑1 (MQ=255)
                                                                                                                                     tCCAGCGGTTCGGTAATCCCCCCCAGCAGGGG‑GGTTAAGGTTGCGGG‑AACCAGCAA‑GCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAgg                  <  1:858292/150‑1 (MQ=255)
                                                                                                                                     cccAGCGGTTCGGTAATCCCCCCCAGCGTGGC‑GGTTAGGGTTGCGGA‑AATCAGCAA‑GCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAggg                 <  1:1049080/150‑1 (MQ=255)
                                                                                                                                       cgggggTTCGTTAATCCCCCCCAGCATGGC‑‑GTTTAGGTTGGGGGAAATAAGCAA‑GCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCa               <  1:406442/147‑1 (MQ=255)
                                                                                                                                       cAGCGTTTCGTTAATCCCCACCGGCAGGGC‑GGTAAAGGTTGCGGA‑AATCAGAAA‑GCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAGGATTTTTGAGTTACCGTGCAGGGCa               <  1:1081496/151‑1 (MQ=255)
                                                                                                                                           ggTTGGGTAATCCCCCCCAGCAGGGC‑GTTTAGGGTTGGGGA‑AATCAGCAA‑CCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAAc           <  2:821916/151‑1 (MQ=255)
                                                                                                                                               cGGTAATCCCCCCCAGCGTGGC‑GGTTAAGGTGGGGGG‑AATCAGCAA‑GCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGc        <  1:91324/150‑1 (MQ=255)
                                                                                                                                                gggtactCCCCCCAGGGTGGG‑GGTTAAGGTTGGTGA‑AAACCACAACCCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGc        <  2:95753/144‑1 (MQ=255)
                                                                                                                                                  tAATTCCCACCAGCATGGC‑GGTTAAGGTTGCGGG‑AATCAGCAA‑GCCCGCCCCTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCCTc    >  1:1179797/1‑151 (MQ=255)
                                                                                                                                                    aTTCCCACCAGCATGGC‑GGTTAAGGTTGCGGG‑AATCAGCAA‑GCCCGCCACTTTTACCCGATTTTTCGGTGCGGCAGTGAAGGACATCGCTAAAAAAAAGCCCACGGCACCAAAAATTTTTGAGTTACCGTGCAGGGCAAAACCGCCCTCCg  >  1:1109633/1‑151 (MQ=255)
                                                                                                                                                      |                                                                                                                                                       
CGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GGTTAAGGTTGCGGG‑AATCAGCAA‑GCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCG  >  NC_000913/3862488‑3862786

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