Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,508,6620TA43.8% 116.0 / 37.7 81E384D (GAA→GATyhfWphosphopentomutase‑related metalloenzyme superfamily protein
Reads supporting (aligned to +/- strand):  ref base T (9/36);  new base A (35/0);  total (44/37)
Fisher's exact test for biased strand distribution p-value = 1.24e-14
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCGCCATCAC  >  NC_000913/3508517‑3508811
                                                                                                                                                   |                                                                                                                                                     
ttgttTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTcac                                                                                                                                                   <  2:1161268/150‑1 (MQ=255)
  gttTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTcacac                                                                                                                                                 <  2:578483/150‑1 (MQ=255)
  gttTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGCGGCGGGCGGGAAAAtttaaac                                                                                                                                                 >  2:855615/1‑145 (MQ=255)
      ttAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGGGCGGCAAAAATaacaaa                                                                                                                                                >  1:503920/1‑145 (MQ=255)
      ttAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGCGGGAAAAATTacacacag                                                                                                                                             >  2:1190874/1‑148 (MQ=255)
      ttAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGCGGGAAAAATTACAAACgg                                                                                                                                             >  2:103190/1‑150 (MQ=255)
      ttAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGAGGCGGGCAGGAAAAATtacaca                                                                                                                                                >  1:356010/1‑147 (MQ=255)
         aaGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATcaaa                                                                                                                                                  >  2:1206299/1‑140 (MQ=255)
             cGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGCGGGAAAAATCACACACGGGAGccccc                                                                                                                                     >  2:1151980/1‑151 (MQ=255)
             cGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGccccc                                                                                                                                     <  2:354840/151‑1 (MQ=255)
             cGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGccccc                                                                                                                                     <  2:298873/151‑1 (MQ=255)
               gagagaTCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGCGGGAAAATTCAAACAAAGGACCCCCcac                                                                                                                                   >  1:1205479/1‑151 (MQ=255)
                 gagaTCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGGGCGGGAAAAATTACACACGGGAGccccccc                                                                                                                                   >  2:736745/1‑147 (MQ=255)
                 gagaTCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACAt                                                                                                                                 <  2:315069/151‑1 (MQ=255)
                   gaTCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGCGGGAAAAATCACAAAAGGCAGgcccccc                                                                                                                                   >  1:1008133/1‑145 (MQ=255)
                       aaTATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATAACAAACGGGAGCCCCCCCATCAGaa                                                                                                                            >  1:872479/1‑150 (MQ=255)
                         tatGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGGGGGGCGCGGGAAAAATCACACCCGGTTGCCACCACATCAGaa                                                                                                                            >  2:21667/1‑148 (MQ=255)
                           tGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGGCCGGGGAAAATCAAACAAACTACCCCCCCCATCaccaa                                                                                                                           >  1:836412/1‑143 (MQ=255)
                              aGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGCGGGAAAAAATAAACAAGGTAGACCCCCCCTCAGaaaaagcggg                                                                                                                    >  1:100759/1‑144 (MQ=255)
                              aGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCgtggt                                                                                                                    <  1:1201674/151‑1 (MQ=255)
                                aTGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTgc                                                                                                                  <  2:1082202/151‑1 (MQ=255)
                                    tCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGCGGGAAAAATCACACACGGTAGCCCCCCCAACAGAAAACGTGGGgcg                                                                                                                 >  2:100209/1‑148 (MQ=255)
                                     cATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGcaca                                                                                                             <  1:303997/151‑1 (MQ=255)
                                          gCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATTACAAAAGGGAGCCCCCCCAACAGAAAgaggggggggcccac                                                                                                            >  1:1108920/1‑133 (MQ=255)
                                          gCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAg                                                                                                        <  1:659663/151‑1 (MQ=255)
                                            gAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGCAAAGAGCGACCATTTTGTGGCGGCCGGGCGCGGAAAAAATCACACAAGGGAGCACCCCCCTCACAAAACGTGGGGGGCCCAccc                                                                                                          >  1:880635/1‑146 (MQ=255)
                                            gAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATTACACACGGGAGCCCCCCCATCAGAAAGCGTGGGGCGCACAcc                                                                                                           >  2:454303/1‑146 (MQ=255)
                                               ggTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGCGGGAAAAATCACACACGGGAGGCCCCACATCAGAAAAAGTGGGGCGCCCACCGAGCt                                                                                                      >  2:556219/1‑148 (MQ=255)
                                                gtgtAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGt                                                                                                  <  2:615736/151‑1 (MQ=255)
                                                 tgtAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGGGGGAAAAATTAAAAAACGAAGACCCCCCATCAGAAAACGGGGGGCGCACACCCAgcggcg                                                                                                   >  2:1027381/1‑149 (MQ=255)
                                                  gtAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTa                                                                                                 <  2:612703/150‑1 (MQ=255)
                                                    aaCCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCg                                                                                              <  2:701107/151‑1 (MQ=255)
                                                       ccGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATTACACACGGTAGGCCCCCCATCAGAAAACGGGGTGCGCACACCGAGCTGCGTAGCGatc                                                                                           >  2:378368/1‑151 (MQ=255)
                                                       ccGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGatc                                                                                           <  2:1219998/151‑1 (MQ=255)
                                                           cGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGGGCGGGAAAAAATACAAACGGGAGCCCCCACCACCAAAAGCGGGGGGCGCACACCCCGCTGCGGAGGGGTCATcc                                                                                       >  1:392068/1‑151 (MQ=255)
                                                           cGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATCACACACCGGAGCCCCCCCACCAGAAAACGGGGGGGGCCCACCCAGTTGCGGAGCGatccccc                                                                                       >  1:444352/1‑147 (MQ=255)
                                                             gATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATTACAAACGGCAGCCCCCACAACAGAAAGCGGGGGGGGCACACCCAGCTGCGTAGCGGTCATccc                                                                                      >  2:976230/1‑150 (MQ=255)
                                                                cAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTg                                                                                   <  2:63232/150‑1 (MQ=255)
                                                                  aaaCTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATCACACACGGTAACCCCCCCAACACAAAACGGGGGGGGCACACCGAGCTGCGTGGCGGGCATCCCCTGCt                                                                                 >  1:742615/1‑150 (MQ=255)
                                                                   aaCTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGa                                                                               <  1:384089/151‑1 (MQ=255)
                                                                   aaCTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGa                                                                               <  2:616779/151‑1 (MQ=255)
                                                                      tCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGGGGGGCGCGGGAAAAATCACACACCGGAGACCCCACATCAGAAAGCGTGGGGCGCCCACCCAGCCGCGTAGCGATCATCCCTTTCTGGTa                                                                             >  2:487714/1‑150 (MQ=255)
                                                                        aTAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGAGGCGCGCGGGAAAATTTACAAACGGGAGCTCCCCCACCAGAAAAAGAGGGGCGCAAACCGAACTGCCGAGCGAACCTCCCCTGCTTATaaa                                                                           >  2:1240219/1‑150 (MQ=255)
                                                                          aGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGGGCGGGAAAAATTACACACGGTAGCCCCCACATCAAGAAACGGGGGGGGCACACCGAGACGCGTAGCGGGCATCCCCTTCTGATaaaa                                                                          >  2:1242128/1‑148 (MQ=255)
                                                                          aGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAAcc                                                                         <  1:545454/150‑1 (MQ=255)
                                                                              tgtCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATTAAAAACGGGAGGCCCCCCATCCCAAAACGGGGGGCGCACACCCAGCTGCGGAGCGGTCATCCCCTTCTGATAAAACAacac                                                                    >  1:1181790/1‑151 (MQ=255)
                                                                              tgtCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGAGGGCGGAAAAAAATACACACGGTAACCCCCCCATCAGAAAACGGGGGGCGCGCAACCAGCTGCGATGCGATCATCCCTTGCTGTTAAAacaaca                                                                     >  2:694414/1‑150 (MQ=255)
                                                                                   ccTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTg                                                                  <  2:1074788/148‑1 (MQ=255)
                                                                                    cTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGGGCGGGAAAAATCAAACACGGTAACCCCCCCATCAGAAAAGGGGGGGGGCACACCGAGCAGCGGAGCGATCATCCCTTCCTGATAAAACAACACTGGcac                                                              >  2:593348/1‑151 (MQ=255)
                                                                                       cAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATCACACACGGGAGCCCCCCCCTCAAAAAACGGGGGGCGCACACCCAGCTGCGGAGCGATCATCCCCTTCTGGTAAAACAACACCGGcacc                                                             >  2:582726/1‑149 (MQ=255)
                                                                                          aCCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTc                                                        <  1:918523/151‑1 (MQ=255)
                                                                                          aCCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTc                                                        <  1:227418/151‑1 (MQ=255)
                                                                                           ccGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGGGGGGGAAAAATTAAACACGGGGGCCCCCACAACAAAAAGAGGGGGGGGCACACCGAGCGGCCTTGCGGTCAACCCCTGCTGGTaaaac                                                                         >  2:613112/1‑131 (MQ=255)
                                                                                              gagggaggAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCgg                                                    <  1:258265/151‑1 (MQ=255)
                                                                                                   aggaAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGCGGGAAAAATCACAAACCGCAGCCCCCCCATCAGAAAGCGCGGGGGGCCCACCCCGCTGCGTAGCGATCATCCCCTTGTTTTAAACCACCACTCGCACCCCCTCCCGGGGAt                                                >  1:34653/1‑150 (MQ=255)
                                                                                                   aggaAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATCACACACCGGAGCCCCCACACCAAAAAGCGGGGTGCGGACACCGAGCTGCGGTGGGATCATTCCCTCCTGAGCAAACAACACTAGCACCACCTcgc                                                      >  1:1222893/1‑144 (MQ=255)
                                                                                                     gaaagaaaAGAGCGACCATTTTGTGGCGGTGGGGGGCGGGAAAAATTACAAACGGGAGCCCCCACATCAGAAAACGGGGGGGGCCCACCCAGCTGCGGAGCGGTCATCCCCTGCTGATAAAACAGCACCGGCACCCCTTCGCGGGGAtgg                                              >  1:534926/1‑150 (MQ=255)
                                                                                                      aaagaaaagaGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTAtggtg                                            <  1:982977/151‑1 (MQ=255)
                                                                                                       aagaaaagaGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTAtggtgg                                           <  2:1012575/151‑1 (MQ=255)
                                                                                                          aaaagaGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATTACACACAGGAGCCCCCCCAACAAAAAGCGGGGGGCGCACACCCAGCTGCGTAGCGATCATCCCCTGCTGATAAAACAACACCGGCACCACCTCCCGGGGAtggtgggt                                         >  2:547436/1‑148 (MQ=255)
                                                                                                          aaaagaGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCtg                                        <  1:108936/151‑1 (MQ=255)
                                                                                                          aaaagaGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCtg                                        <  1:740491/151‑1 (MQ=255)
                                                                                                            aagaGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCtgtg                                      <  2:642049/151‑1 (MQ=255)
                                                                                                            aagaGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCtgtg                                      <  1:424970/151‑1 (MQ=255)
                                                                                                              gagCGACCATTTTGTGGCGGTGGCGCGCGGAAAAACTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGAc                                    <  1:1208247/151‑1 (MQ=255)
                                                                                                                 cGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATCACACACCGTAGCCCCCACATCAGAAAGCGGGGGGCGCACACCGAGCTGCGTAGCGATCATCCCCTTCTGATAAAACAGCACTGGCACCACTTCGCGGGGATGGTGGCTGTGAcc                                   >  2:59997/1‑149 (MQ=255)
                                                                                                                 cGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCaa                                 <  1:53783/151‑1 (MQ=255)
                                                                                                                   aCCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATg                               <  2:534952/151‑1 (MQ=255)
                                                                                                                         ttGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCgg                          <  1:454282/150‑1 (MQ=255)
                                                                                                                          tgtgGCGGTGGCGCGCGGGAAAAATCACACACGGGAGCCCCCCCATCAGAAAGCGTGGGGCGCCCACCGCGCTGCGTAGCGATCATCCCTTGCTGATAAAACCACACTGGCACCACCTCCCGGGGATGGTGGCTGTGACCAATGGGCgg                          >  2:1113071/1‑149 (MQ=255)
                                                                                                                          tgtgGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGa                         <  1:855569/150‑1 (MQ=255)
                                                                                                                           gtggcggtggggCGCGGGAAAAATCACACACGGGTGCCCCCCCATCAGAAAGCGGGGGGCGCACACCGCGCTGCGCAGCGATCATCCCTTGCTGATAAACAAGCACTGGGACCCCCTCTCGGGGATGGTGGGTGTGACCAAAGGGCgg                          >  1:795227/1‑148 (MQ=255)
                                                                                                                           gtggcggtggcgCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATc                       <  1:208701/151‑1 (MQ=255)
                                                                                                                             ggcggtggcgCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCTCTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCg                      <  1:1008726/150‑1 (MQ=255)
                                                                                                                             ggcggtggcgCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGt                     <  2:404953/151‑1 (MQ=255)
                                                                                                                               cggtggcgCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTg                   >  1:380026/1‑151 (MQ=255)
                                                                                                                               cggtggcgCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTg                   >  2:283646/1‑151 (MQ=255)
                                                                                                                                    gcgcgcGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCg                <  1:820107/149‑1 (MQ=255)
                                                                                                                                         gcgGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATcc          >  1:1240898/1‑151 (MQ=255)
                                                                                                                                          cgGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCg         >  1:300289/1‑151 (MQ=255)
                                                                                                                                            ggAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCGc        >  2:1243790/1‑151 (MQ=255)
                                                                                                                                                  aTTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCGCCATCAc  <  2:654832/151‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
TTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCGCCATCAC  >  NC_000913/3508517‑3508811

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