Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,508,6620TA44.2% 120.3 / 36.4 77E384D (GAA→GAT) yhfWphosphopentomutase‑related metalloenzyme superfamily protein
Reads supporting (aligned to +/- strand):  ref base T (7/36);  new base A (34/0);  total (41/36)
Fisher's exact test for biased strand distribution p-value = 2.80e-15
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCGCCATC  >  NC_000913/3508513‑3508809
                                                                                                                                                     |                                                                                                                                                   
cACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAAtt                                                                                                                                                    <  2:1010224/151‑1 (MQ=255)
        ttttAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGCGGGAAAAATCacacacc                                                                                                                                            >  1:855739/1‑150 (MQ=255)
        ttttAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCacacac                                                                                                                                             <  2:1318002/150‑1 (MQ=255)
        ttttAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCacacac                                                                                                                                             <  1:270366/150‑1 (MQ=255)
            aaaGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGCCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAg                                                                                                                                        <  1:1039871/151‑1 (MQ=255)
                ccGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGCGGGAAAAATTACACACCGGAGcgcc                                                                                                                                    >  1:1117247/1‑151 (MQ=255)
                 cGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGccccc                                                                                                                                   <  2:1087078/151‑1 (MQ=255)
                  ggAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCca                                                                                                                                  <  1:209416/151‑1 (MQ=255)
                    agagaTCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATCACACACGGTAGCCCCcac                                                                                                                                 >  1:466342/1‑150 (MQ=255)
                     gagaTCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGCGGGAAAAATCACACACGGTAGCCCCCCCAt                                                                                                                               >  1:520960/1‑151 (MQ=255)
                        aTCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATCACACACGGGAGCCCCCCCATCa                                                                                                                             >  2:362064/1‑150 (MQ=255)
                        aTCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAg                                                                                                                            <  1:795678/151‑1 (MQ=255)
                         tCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATCACACACCGGAGACCCCACATc                                                                                                                              >  2:93370/1‑148 (MQ=255)
                          cAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATTAAACACCGGAGgcccccccacaaa                                                                                                                           >  2:414275/1‑142 (MQ=255)
                           aaTATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGaaa                                                                                                                         <  1:1234594/151‑1 (MQ=255)
                            atatGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGCGGGAAAAATCACACACCGGAGCCCCCCCATCAaaaa                                                                                                                         >  1:250405/1‑150 (MQ=255)
                            atatGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAg                                                                                                                        <  1:150519/151‑1 (MQ=255)
                                         cATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGGGGGAAAAATCACAAACGGTAACCCCCCCATCAGAAAACgggggggggac                                                                                                            >  2:478630/1‑143 (MQ=255)
                                          aTGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCacac                                                                                                          <  1:667274/151‑1 (MQ=255)
                                            gggCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGCGGGAAAAATTACACACCGGAGCCCCCCCATCAGAAAACGGGGGGCGgacaacc                                                                                                        >  1:515475/1‑148 (MQ=255)
                                            gggCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCg                                                                                                        <  2:716019/151‑1 (MQ=255)
                                            gggCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCg                                                                                                        <  2:1204670/151‑1 (MQ=255)
                                             ggCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGGGGGGGAAAAATTACACACCGTAGCCCCCCCATCCGAAAGCGGGGGGCGCACACCg                                                                                                        >  2:23454/1‑150 (MQ=255)
                                               cGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGCGGGAAAAATCAAACACGGTAGCCCCCACATCAGAAAAAGTGGGGGGCACACCGAGc                                                                                                     >  2:1159654/1‑151 (MQ=255)
                                                gAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGCGGGGAAAATCAAACACCGGAGCCCCCACATCAGAAAGCGTGGTGCGGACAACCAGcc                                                                                                    >  2:431574/1‑150 (MQ=255)
                                                gAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATTACACACCGTAGCCCCCACATCAGAAAACGGGGGGCGCACACCGAGCt                                                                                                    >  1:901316/1‑151 (MQ=255)
                                                gAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATTACAAACGGTAGCCCCCACATCAGAAAGCGGGGTGCGCACACCGAGCt                                                                                                    >  1:122329/1‑151 (MQ=255)
                                                gAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATCACACACGGGAGCCCCCACATCAGAAAACGGGGGGGGCACACCCAGcc                                                                                                    >  1:265082/1‑150 (MQ=255)
                                                gAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAAATCACACACCGTAGCCCCCCCATCAGAAAGCGTGGTGCGCACACCCAGCt                                                                                                    >  1:506027/1‑151 (MQ=255)
                                                        aaCCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGCGGGAAAAATTACACACGGGAACCCCCCCATCAGAAAACGGGGGGCGCACACCGCGCTGCGTAGCg                                                                                            >  1:1040865/1‑151 (MQ=255)
                                                           ccGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATCACACACGGGAGCCCCCACATCAGAAAGCGGGGGGCGCACACCCAGCTGCGTAGCGatc                                                                                         >  1:901525/1‑151 (MQ=255)
                                                               cGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATcc                                                                                     <  2:355525/151‑1 (MQ=255)
                                                                      aaaCTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAAAAACAAACCGCAGCCCCCCCCACCGCAAAAGGGGTGGGCACACCGAGCTGCGGCGCGAGCCACCCctgctg                                                                              >  2:1358926/1‑151 (MQ=255)
                                                                        aCTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATCACACACGGGAGCCCCCACATCAGAAAAGGGGGGGGGCGCACCGAGCTGCGGAGCGATCATCCCTTGCTGAt                                                                            >  2:473385/1‑151 (MQ=255)
                                                                          tCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGCGGGAAAAATCACACACCGGAGCCCCCCCATCAGAAAACGGGGGGCGCACACCCAGCTGCGGAGCGGTCATCCCCTTCTGATa                                                                           >  1:605916/1‑150 (MQ=255)
                                                                          tCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGCGGGAAAAATCAAAAACGGGAGCCCCCCCATCAGAAAACGTGGGGCGCACACCCAGCCGCGTAGCGATCATCCCCTGCTGATa                                                                           >  1:553140/1‑150 (MQ=255)
                                                                           cATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATaaa                                                                         <  2:427125/151‑1 (MQ=255)
                                                                             tAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAAATCACACACCGTAGACCCCACATCAGAAAACGTGGGGCGCACACCCAGCCGCGTAGCGATCATCCCTTGCTGATAAAcc                                                                       >  1:784348/1‑151 (MQ=255)
                                                                              aGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAAcc                                                                       <  2:741026/150‑1 (MQ=255)
                                                                              aGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCa                                                                      <  1:1229019/151‑1 (MQ=255)
                                                                                 gtgtCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCa                                                                   <  1:584182/151‑1 (MQ=255)
                                                                                 gtgtCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCa                                                                   <  1:56260/151‑1 (MQ=255)
                                                                                 gtgtCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCa                                                                   <  2:120724/151‑1 (MQ=255)
                                                                                   gtCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCAc                                                                  <  1:694453/150‑1 (MQ=255)
                                                                                      tcCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGGGGGAAAAATCAAAAAACGGAACCCCCCCAACCGAAAACGGGGGGGGCCCACCCCGCTGCGTAGGGGTCATCCCCTGCTGATAAAACAACACTggg                                                              >  2:174303/1‑150 (MQ=255)
                                                                                       ccTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTg                                                                <  2:978412/148‑1 (MQ=255)
                                                                                            aaaCCGGAGGGAGGAGAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACt                                                        <  2:604850/151‑1 (MQ=255)
                                                                                            aaaCCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACt                                                        <  2:894185/151‑1 (MQ=255)
                                                                                            aaaCCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACt                                                        <  1:1174160/151‑1 (MQ=255)
                                                                                              aCCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTc                                                      <  1:201817/151‑1 (MQ=255)
                                                                                               ccGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAAATCACACACGGGAGCCCCCCCATCAGAAAGCGTGGGGCGCACCCCCAGCTGCGTAGCGATCATCCCTTGCTGATAAAACAACACTGGCACCACTTcc                                                     >  1:158967/1‑150 (MQ=255)
                                                                                                cGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGCGGGAAAAATTACACACGGGAGCCCCCACATCAGAAAGCGGGGGGGGCGCACCGAGCTGCGGAGCGAGCATCCCTTGCTGATAAAACAGCACTGGCACCACCTcgc                                                    >  2:1196436/1‑151 (MQ=255)
                                                                                                 ggagggagGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATCACACACCGGAGCCCCCCCATCAGAAAACGGGGGGCGCACACCCAGCTGCGTAGCGATCATCCCCTGCTGATAAAACAGCACTGGCACCACTTcgcg                                                   >  1:1005619/1‑151 (MQ=255)
                                                                                                   agggaggAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATTACAAACGGGAACCCCCCCATCAGAAAACGGGGGGGGCACACCCAGCCGCGTAGCGATCATCCCTTTCTGGTAAAACAGCACTCGGCCCACTTCGCggg                                                 >  2:620779/1‑151 (MQ=255)
                                                                                                         gaaagaaaAGAGCGACCATTTTGTGGCGGTGGCGGGCGGAAAAAATCACACACCGTATCCCCCCCATCAGAAAGCGGGGTGCGCACACCGAGCCGCGTAGCGATCATCCCCTGCTGATAAAACAGCACCGGCACCACCTCGCGGGGAtggg                                           >  1:713949/1‑150 (MQ=255)
                                                                                                          aaagaaaagaGCGACCATTTTGGGTCGGTGGCGAGCGGACAACTTCACACACGGTAGCCCCCAACTCAGAAAGCGTGGTGCGCACACCAAGCTGCGTAGCGATCATCCCTTGCTGATAACCAAAAACGGGCACCACGGCGCCGCGAtggtg                                          <  2:212125/151‑1 (MQ=255)
                                                                                                           aagaaaagaGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTAtggtgg                                         <  2:1210650/151‑1 (MQ=255)
                                                                                                            agaaaagaGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATAACACACCGGAGCCCCCCCATCAGAAAGCGGGGGGGGCACAACCCGCTGCGTAGGGGTCATCCCCTTCTGATAAAACAGCCCCGGCACCACCTCCCgggggtgggggg                                        >  2:919313/1‑142 (MQ=255)
                                                                                                                     cGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATCAAACACGGGAGCCCCCCCATCAGAAAGCGTGGGGCGCACACCCAGCTGCGGAGCGATCATCCCCTTCTGATAAAACAGCACAGGGACCCCTTCGCGGGGATGGGGGGTGGGACCa                                >  2:58938/1‑150 (MQ=255)
                                                                                                                      gACCATTTTGTGGCGGTGGCGCGCGGGAAAAATTACACACGGGAGCCCCCCCCTCAGAAAACGGGGGGGGCACACCCAGCTGCGTAGCGATCATCCCCTGCTGATAAAACAGCACCGGCACCACCTCGCGGGGATGGGGGCTGTGACCaa                               >  1:633393/1‑150 (MQ=255)
                                                                                                                      gACCATTTTGTGGCGGTGGCGCGCGGGAAAAATTACACACGGGAGCCCCCCCATCAGAAAACGTGGGGCGCACACCCAGCTGCGGAGCGATCATCCCCTTCTGATAAAACAACACTGGCACCACCTCGCGGGGATGGGGGCGGTGACCaa                               >  2:305728/1‑150 (MQ=255)
                                                                                                                         cATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGt                           <  1:310323/151‑1 (MQ=255)
                                                                                                                         cATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGt                           <  1:699387/151‑1 (MQ=255)
                                                                                                                          aTTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTc                          <  1:202484/151‑1 (MQ=255)
                                                                                                                              tgtgGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGa                       <  2:397154/150‑1 (MQ=255)
                                                                                                                                   cggtggcgCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTg                 >  1:333398/1‑151 (MQ=255)
                                                                                                                                     gtggcgCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGcc               >  2:1332812/1‑151 (MQ=255)
                                                                                                                                       ggcgCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGTGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGt             <  2:264082/151‑1 (MQ=255)
                                                                                                                                          gcgcgGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGAt          <  2:449156/151‑1 (MQ=255)
                                                                                                                                           cgcgGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATc         <  2:10758/151‑1 (MQ=255)
                                                                                                                                            gcgGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATcc        <  2:454592/151‑1 (MQ=255)
                                                                                                                                            gcgGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATcc        <  2:705527/151‑1 (MQ=255)
                                                                                                                                            gcgGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATcc        >  1:980804/1‑151 (MQ=255)
                                                                                                                                             cgGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCg       >  1:752051/1‑151 (MQ=255)
                                                                                                                                             cgGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCg       >  1:850369/1‑151 (MQ=255)
                                                                                                                                             cgGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCg       >  2:1023383/1‑151 (MQ=255)
                                                                                                                                                aaaaaTTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCGCCa    <  2:1178540/151‑1 (MQ=255)
                                                                                                                                                  aaaTTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCGCCATc  >  1:83224/1‑151 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
CACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCGCCATC  >  NC_000913/3508513‑3508809

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