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

CGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGTTGCGGGAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCG  >  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)
               aaTCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCg                                                                                                                                         >  1:923133/1‑151 (MQ=255)
                   aGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCggat                                                                                                                                      >  2:1156843/1‑148 (MQ=255)
                   aGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGttt                                                                                                                                     >  1:1256928/1‑150 (MQ=255)
                        tCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGGTTAAgtt                                                                                                                                 >  2:538194/1‑149 (MQ=255)
                           gCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGTTGCg                                                                                                                             <  2:780113/151‑1 (MQ=255)
                            cccATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGTTGCgg                                                                                                                            <  1:591475/151‑1 (MQ=255)
                                tGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTTAGGGTGCGGGGAt                                                                                                                       >  1:1274813/1‑151 (MQ=255)
                                       caccacaccaAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGTTGCGGGAATCAGCAAg                                                                                                               <  1:1270710/151‑1 (MQ=255)
                                       caccacaccaAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGTTGCGGGAATCAGCAAg                                                                                                               <  1:1007127/151‑1 (MQ=255)
                                           acaccaAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGTTGCGGGAATCAGCAAGCCCg                                                                                                           <  2:1077240/151‑1 (MQ=255)
                                               caAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGGTAAGGGTGGGGGGAATAGCaaagccccccc                                                                                                        >  2:905176/1‑142 (MQ=255)
                                                  agagATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGGTTAGGGTGCGGGGATCAGCAAACCCCCCCCttt                                                                                                    >  1:283547/1‑151 (MQ=255)
                                                    agaTACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGGTGGGGGGATCAGCAAACCCCCCCCttttt                                                                                                  >  1:862204/1‑150 (MQ=255)
                                                        aCATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTTAAGGTGCGGGGATCAGCAAAACCCCCACTTTTAcccc                                                                                              >  1:854569/1‑150 (MQ=255)
                                                               ggTCGATATGGAGGCCCCCAGCGGGGCGTGTCCCGCAAACAGCAACGGTGAAATAAACAGGAAGGTAAATCCCAGCGGTTCGGTAATCCCCCCCAGCTTGGCGGTTAGGGTTGGGGAAATCAGAAAGCCCGCCACTTTTACCCGATTTTcc                                                                                       <  2:1196055/151‑1 (MQ=255)
                                                               ggTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGAGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGATAAGGGTGCGGGGAAAAGGAAGCCCCGCCCTTTTaaaccacttttc                                                                                       >  1:919598/1‑140 (MQ=255)
                                                                gTCGACATGGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAAGGGTAAAATGAACAGGAAGATAAATTCCAGCGGTTCGGTAATCCCCACCAGCAGGGCGTTTAAGGTTGGGGAAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCg                                                                                      <  1:408879/151‑1 (MQ=255)
                                                                  cGACATTGAGGCCGCCAGCGCGGCGTGTACCGCAAACAGCAAGGGTAAAATAAACAGGAAGGTAATTCCCAGCGGTTCGGTAATCCCCCCCAGCATGGCGGTTAAGGTGGCGGGAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGt                                                                                    <  1:1204094/151‑1 (MQ=255)
                                                                    aCATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGGTTAAGGTTGGGGGAACAACAAACCCCGCACCTTTTTCCCGATTTTCCgggg                                                                                   >  1:325054/1‑149 (MQ=255)
                                                                    aCATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAAGTTGCGGGGATTAGCAAACCCCCCCCCTTTACCCGAtttttccgggc                                                                                  >  2:993973/1‑144 (MQ=255)
                                                                       tggAGGCCGCCAGCCCGGCGTGTCCCGAAAACAGAAAGGGAAAAATGAACAGGAAGGAAAATCCCGGCGTTTCGTTATTCCCCCCCACCATGGCGTTTAGGTTGGGGGAAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCg                                                                                 <  2:1131573/147‑1 (MQ=255)
                                                                        tGAGGCCGCCAGCCCGGCGTGTACCGAAAACAGCAACGGTAAAATGAACAGGAAGGTAAATCCCAGCGGTTCGGTAATCCCCACCAGCAGGGCGTTTAAGGTGGGGGGAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGc                                                                               <  1:753798/150‑1 (MQ=255)
                                                                          aGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGGTAAAGGTTGGGGGAACAGCAAACCCCCCACCTTTTCCCCATTTTTCTGGGCGGGAgg                                                                            >  2:690807/1‑150 (MQ=255)
                                                                             ccgccAGCAGGGCGTGTCCCGCAAACAGCAAGGGTGAAAGGAACAGGAAGGTAAATCCCAGCGGTTCGGTAATCCCCACCAGCAGGGCGGTTAAGGTTGCGGGAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGaa                                                                         <  1:805470/151‑1 (MQ=255)
                                                                              cgccAGCACGGCGTGTACCGAAAACAGCAACGGTAAAATGAACAGGAAGGTAAATCCCAGCGGTTCGGTAATCCCCACCAGCAGGGCGGTTAAGGTGGGGGAAATCAGAAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAg                                                                        <  1:627436/151‑1 (MQ=255)
                                                                               cccAGCACGGGGTGCACCGAAAACAGCAACGGTAAAATAAACAGGAAGGAAAATCCCGGGGTTTCGGTAATCCCCCCCAGCATGGCGGTTAAGGTGGGGGGAACCAGAAGGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGt                                                                       <  2:282843/150‑1 (MQ=255)
                                                                                ccAGCAGGGCGTGTCCCGCAAACAAAAACGGTAAAATGAACGGGAAGTAAAATTCCGGGGGTTCGGTAATCCCCACAAGCATGGCGGTTAGGGTTGGGGGAATCAGAAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTa                                                                      <  2:644250/151‑1 (MQ=255)
                                                                                  aGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAAGGTTCGGGGAACAGCAAGCCCCCCCCTTTTTCCCGCTTTTTCGGGGCgggcgggaaggac                                                                     >  2:1339581/1‑139 (MQ=255)
                                                                                   ccACGGCGTGTCCCGAAAACAGCAACGGTAAAATAAACAGAAAGATAAATCCCAGCGTTTCGTTAATCCCCACCAGCATGGCGTTTAAGGTGGCGGAAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACAt                                                                   <  1:512051/150‑1 (MQ=255)
                                                                                    ccgggggtcacccGCAAACAGCACCGGTGAAATAAAAAGGAAGTTAAATCCAGGCGGTTCGTTAATCCCCACCAGCGTGGCGGTTAGGGTGGGGGGAATGAGAAGGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATc                                                                  <  2:60556/140‑1 (MQ=255)
                                                                                    cACGGCGTCCCCCGCAAACAGCAGGGTAGAAATGAACAGGAAGTTAATTTCCAGCGTTTCGGTAATCCCCCCCAGCATGGCGTTTAGGGTGGGGGAAATCAGAAGGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATc                                                                  <  2:16672/151‑1 (MQ=255)
                                                                                       ggCGTGTCCCGCAAACAGCAACGGTAAAATGAACAGGAAGGAAAATCCCAGCGGTTCGGTAATCCCCCCCAGCAGGGCGGTTAAGGTTGCGGGAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCt                                                               <  1:15940/151‑1 (MQ=255)
                                                                                       ggCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGTTGCGGGGATCAGCAAGCCCGCCACCTTTTCCCGATTTTCCGGGGCCGCAGTGGAGTACATCGCt                                                               >  1:95756/1‑151 (MQ=255)
                                                                                        gcttgtCCCGAAAACAGCAGGGGTAAAATAAACGGGAAGGTAATTTCCGGCGTTTCGGAAATCCCCCCCAGCAGGGCGTTTAGGTTGGCGGGAATCAGAAACCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTa                                                              <  2:498750/148‑1 (MQ=255)
                                                                                          gtgtACCGAAAACAGCAAGGGTAAAATGAACAGGAAGGAAAATCCCGGGGGTTGGGAAATCCCCACCAGCAGGGGGGTTAGGGTGGGGGAAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTaaa                                                            <  1:614442/151‑1 (MQ=255)
                                                                                            gtACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGGTGCGGGAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTaa                                                             >  1:1150658/1‑148 (MQ=255)
                                                                                            gtACCGCAAACAGCAACGGTAAAATGAACAGGAAGGTAAATCCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGTTGCGGGAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAg                                                           <  1:1302270/150‑1 (MQ=255)
                                                                                             tACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGGTGCGGGAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGGGCGGCAGTGAAGTACATCGCTAAAGa                                                          >  2:682324/1‑150 (MQ=255)
                                                                                             acccGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATCCCGGCGGTTCGGTAATCCCCCCCAGCATGGCGGTTAAGGTGGGGGAAATAAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGaa                                                         <  1:365014/149‑1 (MQ=255)
                                                                                                   aaaCAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAAGTTTCGGGGATCAGCAAGCCCCCCCCTTTTTCCCCGTTTTTCGGGGCGGCCGGGAAGTACATCGCTAAAAAAAAGccc                                                   >  2:1213269/1‑151 (MQ=255)
                                                                                                   aaaCAGCAACGGAAAAATGAACAGGAAGGAAAATCCCAGCGGTTCGGAAATCCCCCCCAGCAGGGCGGTTAGGGTTGCGGGAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGccc                                                   <  2:95573/151‑1 (MQ=255)
                                                                                                     aCAGCAACGGTAAAATGAACAGGAAGGTAAATCCCAGCGGTTCGGTAATCCCCCCCAGCAGGGCGTTTAAGGTTGCGGGAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCAc                                                 <  1:267746/151‑1 (MQ=255)
                                                                                                        gCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAAGTTGCGGGAATCAGCAAGGCCCCCACCTTTACCCGATTTTTCGGTGCGGGAGTGAAGGACATCGCCAAAGAAAATCCCCCGGc                                              >  2:840727/1‑151 (MQ=255)
                                                                                                          aaCGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAAGTTGCGGGGATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGGGAAGTACATCGCTAAAGAAATGCCCACGGCAc                                            >  2:1008201/1‑151 (MQ=255)
                                                                                                           aCGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAAGGTGCGGGGATCAGCAAGCCCGCCACTTTTACCCGATTTTTCCGGGCGGGAGTGAAGTACATCGCCAAAGAAATGCCTACGGGAcc                                           >  2:152042/1‑151 (MQ=255)
                                                                                                            cGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAAGGTGCGGGGATCAGCAAGCCCGCCACCTTTAACCGGTTTTTCGGGGCGGCAGTGAAGTACATCGCCAAAAAAAAGCCCCCGGCAccc                                          >  2:385434/1‑150 (MQ=255)
                                                                                                             ggtaaaagaaaCAGAAAGGAAAATCCCAGCGGTTCGGTAATCCCCCCCAGCAGGGGGTTAAAGGTTGGGGGAACAAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCaa                                         <  1:594197/142‑1 (MQ=255)
                                                                                                              gtaaaaTGAACAGGAGGGTAAATTCCAGCGGTTCGGTAATCCCCACCAGCATGGCGGTTAAGGTTGCGGGAATCAGAAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCaaa                                        <  1:587429/148‑1 (MQ=255)
                                                                                                              gtaaaaTGAACAGGAAGGAAAATCCCAGCGGTTCGGTAATCCCCCCCAGCAGGGCGGTTAAGGTTGCGGGAACCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCaaa                                        <  1:1314089/148‑1 (MQ=255)
                                                                                                              gTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGTTGCGGGAATCAGCAAGCCCGCCACTTTTACCCGATTTTTCGGTGCGGCAGTGAAGTACATCGCTAAAAAAATGCCCACCGCACCaa                                         >  1:696656/1‑150 (MQ=255)
                                                                                                               aaaaagaaaCAGGAAGTTATTCTCCGGCGGTTCGAAATTCCCCCCCAGCGGGGCGGTTAAGGTTGCGGAAATAAGAAGGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAg                                       <  2:591288/144‑1 (MQ=255)
                                                                                                                   aggAACGGAAAGGTAAATCCCAGGGTTTGGGAAATCCCCCCCAGCTGGGCGTTTAGGGTGGGGGAAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGAttt                                   <  2:290088/149‑1 (MQ=255)
                                                                                                                    aaaaaaGGAAGGTAAATCCCAGCGGTTGGGTAATCCCCCCCAGAAGGGCGGTTAAGGTTGGGGAAACCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGAtttt                                  <  1:291197/146‑1 (MQ=255)
                                                                                                                     gAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGGTGCGGGAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGGGAAGTACATCGCTAAAAAAATGCCCACGGCACCAAAGAttttt                                 >  1:812924/1‑151 (MQ=255)
                                                                                                                     aaaCAGGAAGGAAAATTCCGGGGTTTCGATAATTCCCCCCAGCATGGCGGTAAGGGTGGCGGAAAACAGAAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGAttttt                                 <  2:516422/150‑1 (MQ=255)
                                                                                                                        aaGGAAGTAAAATCCCAGGGGTTCGGTAATCCCCCCCAGCATGGCGGTTAAGGTGGCGGGAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAg                              <  1:479922/150‑1 (MQ=255)
                                                                                                                           gAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGTTGCGGGGATCAGCCAGCCCCCCACTTTTTCCCGATTTTCCGGGGGGGCAGTGGAGTACAACGCTAAAAAAAAGCCCACCGCACCCAAAATTTTTTAGtt                            >  2:885376/1‑150 (MQ=255)
                                                                                                                              ggTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGTTTCGGGGATCAGCAAGCCCGCCACTTTTTCCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATTCCCACGGCACCCAAGAATTTTTAGTTAccc                        >  1:363959/1‑150 (MQ=255)
                                                                                                                                 aaaTCCCAGCGTTTCGTTAATCCCCCCCGGCAGGGCGGTTAGGGTGGCGGAAACCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGc                     <  2:1354125/151‑1 (MQ=255)
                                                                                                                                  aatcccggggtttCGATAATCCCCACCAGCAGGGGGTTTAGGGTGGGGGAAACCAGCAACCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCa                    <  2:1237757/140‑1 (MQ=255)
                                                                                                                                     tCCAGCGGTTCGGTAATCCCCCCCAGCAGGGGGGTTAAGGTTGCGGGAACCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAgg                  <  1:858292/150‑1 (MQ=255)
                                                                                                                                     cccAGCGGTTCGGTAATCCCCCCCAGCGTGGCGGTTAGGGTTGCGGAAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAggg                 <  1:1049080/150‑1 (MQ=255)
                                                                                                                                       cgggggTTCGTTAATCCCCCCCAGCATGGCGTTTAGGTTGGGGGAAATAAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCa               <  1:406442/147‑1 (MQ=255)
                                                                                                                                       cAGCGTTTCGTTAATCCCCACCGGCAGGGCGGTAAAGGTTGCGGAAATCAGAAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAGGATTTTTGAGTTACCGTGCAGGGCa               <  1:1081496/151‑1 (MQ=255)
                                                                                                                                           ggTTGGGTAATCCCCCCCAGCAGGGCGTTTAGGGTTGGGGAAATCAGCAACCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAAc           <  2:821916/151‑1 (MQ=255)
                                                                                                                                               cGGTAATCCCCCCCAGCGTGGCGGTTAAGGTGGGGGGAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGc        <  1:91324/150‑1 (MQ=255)
                                                                                                                                                gggtactCCCCCCAGGGTGGGGGTTAAGGTTGGTGAAAACCACAACCCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGc        <  2:95753/144‑1 (MQ=255)
                                                                                                                                                  tAATTCCCACCAGCATGGCGGTTAAGGTTGCGGGAATCAGCAAGCCCGCCCCTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCCTc    >  1:1179797/1‑151 (MQ=255)
                                                                                                                                                    aTTCCCACCAGCATGGCGGTTAAGGTTGCGGGAATCAGCAAGCCCGCCACTTTTACCCGATTTTTCGGTGCGGCAGTGAAGGACATCGCTAAAAAAAAGCCCACGGCACCAAAAATTTTTGAGTTACCGTGCAGGGCAAAACCGCCCTCCg  >  1:1109633/1‑151 (MQ=255)
                                                                                                                                                      |                                                                                                                                                       
CGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGTTGCGGGAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCG  >  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: