Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,508,6620TA40.6% 98.8 / 22.5 64E384D (GAA→GAT) yhfWphosphopentomutase‑related metalloenzyme superfamily protein
Reads supporting (aligned to +/- strand):  ref base T (10/28);  new base A (26/0);  total (36/28)
Fisher's exact test for biased strand distribution p-value = 4.44e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCGCCATC  >  NC_000913/3508512‑3508809
                                                                                                                                                      |                                                                                                                                                   
ccACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAAt                                                                                                                                                     <  1:571436/151‑1 (MQ=255)
  aCGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTc                                                                                                                                                   <  2:916617/151‑1 (MQ=255)
     ttgttTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATca                                                                                                                                                  >  2:683668/1‑149 (MQ=255)
      tgttTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGCAAAAtt                                                                                                                                                    >  2:123628/1‑146 (MQ=255)
        tttttAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATtacaca                                                                                                                                              >  2:354632/1‑150 (MQ=255)
        tttttAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGGCGCGGGCGGGAAAAATcaaaa                                                                                                                                               >  2:309492/1‑146 (MQ=255)
           ttAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACgg                                                                                                                                           <  1:767005/150‑1 (MQ=255)
              aaGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGGGGAAAAATTCACACACCGTa                                                                                                                                         >  2:894954/1‑149 (MQ=255)
                gCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGCGGGAAAAATCAAACACGGGAGccc                                                                                                                                     >  1:527361/1‑151 (MQ=255)
                  cGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAAATCACACACAGTAGccccc                                                                                                                                   >  2:806126/1‑151 (MQ=255)
                   ggAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGCGGGAAAAATCACACACCGTAGacccc                                                                                                                                   >  2:800079/1‑150 (MQ=255)
                     agagaTCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAAATCAAACACCGGAGccccccc                                                                                                                                 >  2:67407/1‑148 (MQ=255)
                        gaTCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCa                                                                                                                             <  1:184360/151‑1 (MQ=255)
                               atGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGc                                                                                                                       <  2:874211/150‑1 (MQ=255)
                                 gCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGTGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCgtg                                                                                                                    <  1:623432/151‑1 (MQ=255)
                                      tGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGAGCGGGGAAAAAATCACAAAAAGTAAGCCCCACCTCAAAAAGCgt                                                                                                                     >  2:33965/1‑145 (MQ=255)
                                         tCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGCGGGAAAATTCACACACGGTAGCCCCCACATCAGAAAACGTGGGGCGcac                                                                                                            >  2:912233/1‑151 (MQ=255)
                                         tCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCCCGGAAAAAATGACACACCGTAGGCCCCACACCAGAAAgggaggcg                                                                                                                 >  1:151029/1‑139 (MQ=255)
                                          cATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGCGGGAAAAATCACAAACGGTAGCCCCCACATCAGAAAGCGGGGTGGGcac                                                                                                            >  2:131979/1‑150 (MQ=255)
                                          cATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGcac                                                                                                            <  2:616311/150‑1 (MQ=255)
                                           aTGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAATTCACACACGGGAGCCCCCCCATCAGAAAGCGTGGGGCGCacac                                                                                                          >  2:233122/1‑151 (MQ=255)
                                                cGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGc                                                                                                     <  1:472916/151‑1 (MQ=255)
                                                 gAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATCACAAACGGGAGCCCCCCCCTTAGAAAGCGGGGTGCGGACACCGAGc                                                                                                     >  1:280590/1‑150 (MQ=255)
                                                 gAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCt                                                                                                    <  1:954275/151‑1 (MQ=255)
                                                   gggTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGc                                                                                                  <  1:360278/151‑1 (MQ=255)
                                                    ggTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCg                                                                                                 <  1:792131/151‑1 (MQ=255)
                                                     gtgtAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGt                                                                                                <  2:483390/151‑1 (MQ=255)
                                                     gtgtAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGt                                                                                                <  2:977948/151‑1 (MQ=255)
                                                       gtAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTa                                                                                               <  1:172218/150‑1 (MQ=255)
                                                           cccGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAAATCACACACGGTAGCCCCCACATCAGAAAACGTGGGGCGCACACCCAGCTGCGTAGCg                                                                                            >  2:257946/1‑149 (MQ=255)
                                                                cGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATcc                                                                                     <  2:302533/151‑1 (MQ=255)
                                                                   aTCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATCACACACGGGAGCCCCCACATCAGAAAGCGTGGGGGGCACACCCAGCTGCGTAGCGATCATCCCtt                                                                                  >  2:858540/1‑151 (MQ=255)
                                                                   aTCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCtt                                                                                  <  2:780590/151‑1 (MQ=255)
                                                                     cAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTg                                                                                 <  1:172906/150‑1 (MQ=255)
                                                                      aaaaCTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATAAAACACCGTAGCCCCCACACCACAAAAAGGGGGGCGCACACCCAGCGGCGTAGCGATCATCCCCTGc                                                                                >  1:777692/1‑150 (MQ=255)
                                                                       aaaCTCATAGGGTGTCTCATGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGc                                                                                <  1:228348/149‑1 (MQ=255)
                                                                             aTAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATCACACACGGGAGCCCCCCCATCAGAAAACGGGGGGCGCACACCCAGCTGCGTAGCGATCATCCCCTGCTGATaaa                                                                         >  2:336328/1‑150 (MQ=255)
                                                                                       tcCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATCACACACCGGAACCCCCACATCAGGAAGCGTGGGGGGCACACCGAGCTGCGGAGCGAGCATCCCCTGCTGATAAAACAGCAcatg                                                               >  1:114673/1‑147 (MQ=255)
                                                                                           gCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGcacca                                                          <  1:968648/151‑1 (MQ=255)
                                                                                           gCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGcacc                                                           <  2:352110/150‑1 (MQ=255)
                                                                                           gCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGcacc                                                           <  2:795632/150‑1 (MQ=255)
                                                                                            cAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAAATTACACACGGGAGCCCCCCCATCAGAAAACGTGGGGCGCACACCCAGCTGCGTAGCGATCATCCCTTGCTGGTAAACCAGCACTGGgaccac                                                         >  1:968583/1‑151 (MQ=255)
                                                                                                 cGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAATTCACACACGGGAGCCCCCCCCTCAGAAAGCGGGGGGCGCACACCCAGCTGCGTAGCGATCATCCCCTTCTGATAAACCAACACCGGCACCACCTcg                                                     >  2:924802/1‑150 (MQ=255)
                                                                                                 cGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATTACACACGGGAGCCCCCCCATCAGAAAACGGGGGGCGCACACCGAGCTGCGTAGCGATCATCCCTTTCTGATAAAACAGCACTGGCACCACTTcgc                                                    >  2:738225/1‑151 (MQ=255)
                                                                                                 cGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATCACACACGGTAGCCCCCCCACCAGAAAGCGGGGTGCGCACACCCAGCTGCGGAGCGATCATCCCTTGCTGATAAAACAGCACTGGcaccacct                                                       >  2:281169/1‑146 (MQ=255)
                                                                                                  ggagggagGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATCACACACGGGAGCCCCCCCATCAGAAAGCGGGGGGCGCACACCCAGCTGCGGAGCGATCATCCCCTGCTGATAAACCAGCACTCGCACCACCTcg                                                     >  1:51599/1‑149 (MQ=255)
                                                                                                    agggaggAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATTAAACACGGTAGCCCCCCCATCAGAAAACGTGGGGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAAACAGCACTGGCACCACTTCGCggg                                                 >  2:382218/1‑151 (MQ=255)
                                                                                                    agggaggAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAAATTAAACACGGTAGCCCCCCCATCAGAAAGCGTGGTGCGCACACCCAGCTGCGGAGCGATCATCCCCTGCTGATAAACCAACACCGGCACCACCTCGCggg                                                 >  2:134145/1‑151 (MQ=255)
                                                                                                          gaaagaaaAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTAtggt                                           <  1:33965/151‑1 (MQ=255)
                                                                                                            aagaaaagaGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATTACACACCGGAACCCCCCCATCAGAAAGCGTGGGGCGCACACCCAGCTGCGGAGCGATCATCCCCTTCTGGTAAACCAACACCGGCACCACTTCGCGGGGAtggtgg                                         >  2:469760/1‑151 (MQ=255)
                                                                                                            aagaaaagaGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAAATCACAAACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCCAGCCGCGTAGCGATCATCCCTTGCTGATAAACCAGCACAGGCACCACTTTTCGGGTAtgggg                                          >  2:824664/1‑148 (MQ=255)
                                                                                                             agaaaagaGCGACCATTTTGTGGCGGTGGCGCGAGGAAAAATTCACACACCGTAGCCCCCACAATAAAAAGTGTGGCGCGCACCGCCGGCGTCGGTGCGGTAAACCCTTGTTGATAAAACAGCACTGGTGCCACTTTTCGGGTCtggttg                                         >  2:457522/1‑148 (MQ=255)
                                                                                                                aaagaGCGACCATTTTGTGGCGGGGGGGGGGGGGAAAAATAACACACAGGAACCCCCCCATCACAAAGCGTGGGGGGGACAACCGGCTGCGGAGCGATCATCCCCTGCCGGTAAAACAACACCGGCCCCCCTTCGCGGGGATGGGGGGtg                                      >  1:871580/1‑150 (MQ=255)
                                                                                                                 aagaGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCtgtg                                    <  2:226121/151‑1 (MQ=255)
                                                                                                                  agagCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGa                                   <  1:436191/151‑1 (MQ=255)
                                                                                                                          cATTTTGTGGCGGTGGCGCGCGGGAAAAATCACACACCGTAGCCCCCCCCTCAGAAAACGGGGGGCGCCCACCCAGCTGCGTAGCGGTCATCCCCTGCTGATAAACCAACACCGGCACCACCTCGCGGGGATGGTGGGTGTGACCCATgg                            >  2:515537/1‑150 (MQ=255)
                                                                                                                           aTTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTc                          <  1:235736/151‑1 (MQ=255)
                                                                                                                                gtggcggtggcgCGCGGAAAAATTCAC‑C‑CGGTAGCCCCCCCATCAGAAAGCGTGGGGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAAACAGCACTGGCACCACTTCGCGGGGATGGTGGCTGTGACCAATTGTCGGATCg                    >  2:553/1‑150 (MQ=255)
                                                                                                                                  ggcggtggcgCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCg                    <  1:736451/150‑1 (MQ=255)
                                                                                                                                    cggtggcgCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGtt                  >  2:345017/1‑150 (MQ=255)
                                                                                                                                     ggtggcgCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGc                <  1:246182/151‑1 (MQ=255)
                                                                                                                                      gtggcgCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGcc               <  1:232980/151‑1 (MQ=255)
                                                                                                                                              cgGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCg       >  2:788386/1‑151 (MQ=255)
                                                                                                                                                 aaaaaTTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCCCTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCGCCa    <  2:524653/151‑1 (MQ=255)
                                                                                                                                                  aaaaTTCACACACGGTAGCCCCCACATCAGAAAGCGTGGCGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGAATGGTGGCTGTGACCAATGGCCGGATCGTTGCCGTGATCCGCCAt   >  2:293975/1‑151 (MQ=255)
                                                                                                                                                   aaaTTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCGCCATc  >  2:669314/1‑151 (MQ=255)
                                                                                                                                                      |                                                                                                                                                   
CCACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCGCCATC  >  NC_000913/3508512‑3508809

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