Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,649,2950TG47.9% 109.9 / 21.7 71intergenic (+254/‑314)ydfC/dicBuncharacterized protein, Qin prophage/Qin prophage; cell division inhibition protein
Reads supporting (aligned to +/- strand):  ref base T (4/33);  new base G (34/0);  total (38/33)
Fisher's exact test for biased strand distribution p-value = 4.64e-16
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGGATATTTCCCTGAAAAGTCAGACCATAACGCGAAAGCGCATGGCGAGGTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGC‑GGGGTTATTCCTT‑CCCCGTTGAGGACACC‑GGGTTGTCAGGTTGACCATACGCTTAAGTGACAACCCCGCTGCAACGCCCTCTGTTA‑TCAATTTTCTGGTGAC‑GTTTGGCGGTATCAGTTTTACTCCGTGACTGCTCTGCCGCCCTTTTT  >  NC_000913/1649148‑1649439
                                                                                                                                                    |                                                                                                                                                    
aGGATATTTCCCTGAAAAGTCAGACCATAACGCGAAAGCGCATGGCGAGGTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGC‑GGGGTTAt                                                                                                                                                   <  2:525483/151‑1 (MQ=255)
  gATATTTCCCTGAAAAGTCAGACCATAACGCGAAAGCGCATGGCGAGGTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGC‑GGGGTTAttc                                                                                                                                                 <  2:1259898/151‑1 (MQ=255)
  gATATTTCCCTGAAAAGTCAGACCATAACGCGAAAGCGCATGGCGAGGTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGC‑GGGGTTAttc                                                                                                                                                 <  1:775874/151‑1 (MQ=255)
     atTTCCCTGAAAAGTCAGACCATAACGCGAAAGCGCATGGCGAGGTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCNGNNNNGGCAACCCGCGAAANGNCGCGGCGGTAAGTATGGC‑GGGGGTAtttctt                                                                                                                                              >  2:233376/1‑151 (MQ=255)
     atTTCCCTGAAAAGTCAGACCATAACGCGAAAGCGCATGGCGAGGTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGC‑GGGGGTGtttctt                                                                                                                                              >  1:873783/1‑151 (MQ=255)
          ccTGAAAAGTCAGACCATAACGCGAAAGCGCATGGCGAGGTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGC‑GGGGTTATTCCTT‑CCCCg                                                                                                                                        <  1:853905/151‑1 (MQ=255)
           cTGAAAAGTCAGACCATAACGCGAAAGCGCATGGCGAGGTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGC‑GGGGGTTTTCCTT‑tccccc                                                                                                                                       >  2:518087/1‑149 (MQ=255)
             gAAAAGTCAGACCATAACGCGAAAGCGCATGGCGAGGTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGC‑GGGGGTTTTTCTT‑cccccgtg                                                                                                                                     >  1:636593/1‑147 (MQ=255)
             gAAAAGTCAGACCATAACGCGAAAGCGCATGGCGAGGTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGC‑GGGGGTTTTCCTT‑cccccctg                                                                                                                                     >  1:1085150/1‑147 (MQ=255)
             gAAAAGTCAGACCATAACGCGAAAGCGCATGGCGAGGTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGC‑GGGGGTATTTCTT‑CCCCCTTg                                                                                                                                     >  2:1193359/1‑151 (MQ=255)
                    cAGACCATAACGCGAAAGCGCATGGCGAGGTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGC‑GGGGTTATTCCTT‑CCCCGTTGAGGacac                                                                                                                              <  2:988103/151‑1 (MQ=255)
                    cAGACCATAACGCGAAAGCGCATGGCGAGGTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGC‑GGGGGTTTTTCTT‑CCCCCTTGAGGacac                                                                                                                              >  2:1238362/1‑151 (MQ=255)
                    cAGACCATAACGCGAAAGCGCATGGCGAGGTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGC‑GGGGGTATTTCTT‑TCCCGTTGAGGgcac                                                                                                                              >  2:1182106/1‑151 (MQ=255)
                     aGACCATAACGCGAAAGCGCATGGCGAGGTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGC‑GGGGGTTTTTCTT‑CCCCCTTGAggggacc                                                                                                                             >  1:284828/1‑146 (MQ=255)
                        ccATAACGCGAAAGCGCATGGCGAGGTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGC‑GGGGTTATTCCTT‑CCCCGTTGAGGACACC‑ggg                                                                                                                         <  1:1117958/151‑1 (MQ=255)
                        ccATAACGCGAAAGCGCATGGCGAGGTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGC‑GGGGGTTTTTCTT‑TCCCGTTGAGGacaaccgg                                                                                                                          >  2:293252/1‑146 (MQ=255)
                         cATAACGCGAAAGCGCATGGCGAGGTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGC‑GGGGTTATTCCTT‑CCCCGTTGAGGACACC‑GGGt                                                                                                                        <  1:337785/151‑1 (MQ=255)
                          aTAACGCGAAAGCGCATGGCGAGGTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGC‑GGGGTTATTCCTT‑CCCCGTTGAGGACACC‑GGGtt                                                                                                                       <  1:1263784/151‑1 (MQ=255)
                           tAACGCGAAAGCGCATGGCGAGGTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGC‑GGGGGTATTTCTT‑CCCCCGTGAGGGCACC‑CGGGTg                                                                                                                      >  2:607583/1‑151 (MQ=255)
                             aCGCGAAAGCGCATGGCGAGGTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGC‑GGGGTTATTCCTT‑CCCCGTTGAGGACACC‑GGGTTGTc                                                                                                                    <  2:185918/151‑1 (MQ=255)
                                cgAAAGCGCATGGCGAGGTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGC‑GGGGTTATTCCTT‑CCCCGTTGAGGACACC‑GGGTTGTCAgg                                                                                                                 <  2:592275/151‑1 (MQ=255)
                                 gAAAGCGCATGGCGAGGTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGC‑GGGGGTATTTCTT‑TCCCGTTTAGGGGACC‑GGGgtgttcggg                                                                                                                >  1:917253/1‑146 (MQ=255)
                                     gcgcATGGCGAGGTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGC‑GGGGTTATTCCTT‑CCCCGTTGAGGACACC‑GGGTTGTCAGGTTGAc                                                                                                            <  2:198266/151‑1 (MQ=255)
                                      cgcATGGCGAGGTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGC‑GGGGGTTTTCCTT‑CCCCCCTTTGGGGACC‑GGGGTGTCAGGGTGGcc                                                                                                           >  1:936392/1‑151 (MQ=255)
                                           ggCGAGGTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGC‑GGGGTTATTCCTT‑CCCCGTTGAGGACACC‑GGGTTGTCAGGTTGACCATAc                                                                                                       <  1:417245/150‑1 (MQ=255)
                                           ggCGAGGTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGC‑GGGGTTATTCCTT‑CCCCGTTGAGGACACC‑GGGTTGTCAGGTTGACCATAc                                                                                                       <  1:237650/150‑1 (MQ=255)
                                                 gTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGC‑GGGGTTATTCCTT‑CCCCGTTGAGGACACC‑GGGTTGTCAGGTTGACCATACGCTTAAg                                                                                                <  1:200340/151‑1 (MQ=255)
                                                          tCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGC‑GGGGGTATTTCCT‑CCCCCTTGAGGGCACC‑CGGGTGTCAGGGTTACCCTACGCTTTAGTGGCAAccc                                                                                       >  2:667803/1‑151 (MQ=255)
                                                              agatagCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGC‑GGGGTTATTCCTT‑CCCCGTTGAGGACACC‑GGGTTGTCAGGTTGACCATACGCTTAAGTGACAACCCCGc                                                                                    <  1:928652/150‑1 (MQ=255)
                                                               gatagCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGC‑GGGGTTATTCCTT‑CCCCGTTGAGGACACC‑ggggaaggaataac                                                                                                              <  1:536844/123‑12 (MQ=255)
                                                                  agCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGC‑GGGGGTATTTCTT‑CCCCCTTGAGGGCACC‑CGGGTGTCAGGGTGACCATACGCTTAAATGACAAACCCCCTGcc                                                                                >  2:1088168/1‑149 (MQ=255)
                                                                          gTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGC‑GGGGGTATTTCTT‑TCCCGTTGAGGAAACC‑CGGGTGTTAGGGTGACCATACGCTTAAATGACAACCCCCCTGCAAAGcccccc                                                                       >  2:722164/1‑148 (MQ=255)
                                                                            tAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGC‑GGGGTTATTCCTT‑CCCCGTTGAGGACACC‑GGGTTGTCAGGTTGACCATACGCTTAAGTGACAACCCCGCTGCAACGCCCTCTg                                                                      <  2:239749/150‑1 (MQ=255)
                                                                             aaaTTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGC‑GGGGTTATTCCTT‑CCCCGTTGAGGACACC‑GGGTTGTCAGGTTGACCATACGCTTAAGTGACAACCCCGCTGCAACGCCCTCTGt                                                                     <  1:91848/150‑1 (MQ=255)
                                                                             aaaTTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGC‑GGGGGTATTTCTT‑TCCCCTTGAGGGCACC‑GGGGTGTTAGGGTGAACATAAGCTTAAGTGACAACCCCGCCGCAACGCCCCCTGtt                                                                    >  2:677909/1‑151 (MQ=255)
                                                                                 tttCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGC‑GGGGTTATTCCTT‑CCCCGTTGAGGACACC‑GGGTTGTCAGGTTGACCATACGCTTAAGTGACAACCCCGCTGCAACGCCCTCTGTTA‑Tc                                                                <  2:1067115/150‑1 (MQ=255)
                                                                                 tttCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGC‑GGGGGTATTTCTT‑TCCCGCTGGGGGCACC‑GGGGTGTTCGGGTGGCCATAACCTTAAATGAAAACCCCGCTGCCACGCCcccttgtttt                                                                 >  2:176634/1‑144 (MQ=255)
                                                                                  ttCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGC‑GGGGTTATTCCTT‑CCCCGTTGAGGACACC‑GGGTTGTCAGGTTGACCATACGCTTAAGTGACAACCCCGCTGCAACGCCCTCTGTTA‑TCa                                                               <  1:934659/150‑1 (MQ=255)
                                                                                  ttCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGC‑GGGGGTTTTTCTT‑TCCCCCTGTGGGCACC‑CCGGTGTTCCGGTGACCACACCCTTTAGTGAAAAACCCCCCGCCACACCCCCTCTTT‑TCa                                                               >  1:1431749/1‑150 (MQ=255)
                                                                                    cgtcgACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGC‑GGGGTTATTTCTT‑TCCCCGTGTGGGCACC‑GGGGTGTCAGGGTGAGCATACGCCTAAAAGACAAACCCCCTGCCACGCCCCCTGTTA‑TCAAtt                                                            >  2:48020/1‑151 (MQ=255)
                                                                                    cgtcgACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGC‑GGGGGTATTCCTT‑TCCCGTTGAGGACACC‑GGGGTGTCAGGGTGACCATACCCTTAAGTGACAACCCCCCCGCAACGCCCCCTGTTT‑TCAAtt                                                            >  2:144747/1‑151 (MQ=255)
                                                                                       cgACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGC‑GGGGTTATTCCTT‑CCCCGTTGAGGACACC‑GGGTTGTCAGGTTGACCATACGCTTAAGTGACAACCCCGCTGCAACGCCCTCTGTTA‑TCAAtttt                                                          <  1:252280/150‑1 (MQ=255)
                                                                                        gACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGC‑GGGGGTATTTCTT‑CCCCCCTGTGGACACC‑CGGTTGTCAGGGTTTCCCTACGCTTTAGGGACAAACCCCCTGCAACCCCCCCTGTTT‑TTAAtttttt                                                        >  1:487898/1‑149 (MQ=255)
                                                                                        gACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGC‑GGGGGTATTTCTT‑CCCCCCTGGGGAGACC‑CGGGTGTTAGGGTGGCCATACGCTTAAGGGGCAACCCCCCCGCCACGCCCCCTTTTT‑TTCAtttttt                                                        >  1:170448/1‑149 (MQ=255)
                                                                                           cGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGC‑GGGGGTTTTTCTT‑CCCCCCTGTGGACACC‑CGGGTGTCAGGGTGACCATACGCTTAAGTGACAACCCCCCCGCCACCCCCCCTTTTT‑TTCAttttttttg                                                     >  1:1406688/1‑146 (MQ=255)
                                                                                            gTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGC‑GGGGTTATTCCTT‑CCCCGTTGAGGACACC‑GGGTTGTCAGGTTGACCATACGCTTAAGTGACAACCCCGCTGCAACGCCCTCTGTTA‑TCAATTTTCTGGTg                                                    <  2:1095700/151‑1 (MQ=255)
                                                                                                    ccGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGC‑GGGGTTATTCCTT‑CCCCGTTGAGGACACC‑GGGTTGTCAGGTTGACCATACGCTTAAGTGACAACCCCGCTGCAACGCCCTCTGTTA‑TCAATTTTCTGGTGAC‑GTTTg                                            <  2:806662/150‑1 (MQ=255)
                                                                                                     cGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGCGGGGGTTTTTCTTT‑CCCCCTTAGGGCCCCC‑GGGTT‑TTACGTTGACCCTACGCTTTAGTGGCAACCCCCCCGCAACGCCCTCTGTTT‑TTCATTTTTTTGTGACGGTTTgg                                           >  1:291814/1‑151 (MQ=255)
                                                                                                            ggCAACCCGCGAAATGGCGCGGCGGTAAGTATGGC‑GGGGTTATTCCTT‑CCCCGTTGAGGACACC‑GGGTTGTCAGGTTGACCATACGCTTAAGTGACAACCCCGCTGCAACGCCCTCTGTTA‑TCAATTTTCTGGTGAC‑GTTTGGCGGTATc                                    <  2:690577/150‑1 (MQ=255)
                                                                                                            ggCAACCCGCGAAATGGCGCGGCGGTAAGTATGGC‑GGGGTTATTCCTT‑CCCCGTTGAGGACACC‑GGGTTGTCAGGTTGACCATACGCTTAAGTGACAACCCCGCTGCAACGCCCTCTGTTA‑TCAATTTTCTGGTGAC‑GTTTGGCGGTATc                                    <  2:870341/150‑1 (MQ=255)
                                                                                                            ggCAACCCGCGAAATGGCGCGGCGGTAAGTATGGC‑GGGGTTATTCCTT‑CCCCGTTGAGGACACC‑GGGTTGTCAGGTTGACCATACGCTTAAGTGACAACCCCGCTGCAACGCCCTCTGTTA‑TCAATTTTCTGGTGAC‑GTTTGGCGGTATc                                    <  2:873784/150‑1 (MQ=255)
                                                                                                            ggCAACCCGCGAAATGGCGCGGCGGTAAGTATGGC‑GGGGGTATTCCTT‑CCCCCGTTAGGCGACC‑CGGGTGTGGGGTTGCCCCTATGCCTTAGAGACAAACCCGGTGCACCCCCCTCTGTTA‑TCTATTTTCTTGTGTC‑GATTTGCGGGATc                                    >  1:27768/1‑150 (MQ=255)
                                                                                                                 cccGCGAAATGGCGCGGCGGTAAGTATGGC‑GGGGTTATTCCTT‑CCCCGTTGAGGACACC‑GGGTTGTCAGGTTGACCATACGCTTAAGTGACAACCCCGCTGCAACGCCCTCTGTTA‑TCAATTTTCTGGTGAC‑GTTTGGCGGTATCAGtttt                              <  2:494625/151‑1 (MQ=255)
                                                                                                                  ccGCGAAATGGCGCGGCGGTAAGTATGGC‑GGGGGTTTTTCCT‑CCCCCTTGAGGGCACC‑CGGGTGTTAGGGTGACCATACGCCTAAATGACAAACCCCCCGCCACGCCCCCTGTTT‑TTAATTTTTTTGGGGC‑GGTTGGGGGGATCAGttttt                             >  1:1044569/1‑150 (MQ=255)
                                                                                                                  ccGCGAAATGGCGCGGCGGTAAGTATGGC‑GGGGGTATTCCTT‑CCCCGTTGAGGAAACC‑GGGGTGTCAGGTTGACCATACGCTTAAGTGACAACCCCGCTGCAACGCCCTCTGTTA‑TCAATTTTTTTGTGAC‑GGTTTGCGGGATTAGttttt                             >  2:722196/1‑150 (MQ=255)
                                                                                                                    gcgAAATGGCGCGGCGGTAAGTATGGC‑GGGGGTTTTTCTT‑CCCCCTTGAGGGCACC‑CGGGTGTTAGGGTGACCATACGCTTTAGTGACAAACCCCCCGCCAAGCCCCCTGTTT‑TTAATTTTTTTGGGGC‑GGTTGGGGGGATCAGgttttat                           >  2:960857/1‑148 (MQ=255)
                                                                                                                    gcgAAATGGCGCGGCGGTAAGTATGGC‑GGGGGTATTTCTT‑CCCCCCTTAGGGCACC‑CGGGTGTTAGGGTGGACATACGCTTTAGTGACAAACCCGCCGCCACGCCCCCCGTTT‑TTAAATTTTTTGTGGC‑GTTTTGCGGGATCAGgttttat                           >  1:923256/1‑148 (MQ=255)
                                                                                                                    gcgAAATGGCGCGGCGGTAAGTATGGC‑GGGGGTATTTCTT‑CCCCCCTGAGGACACC‑CGGGTGTCAGGGTGACCCTACGCTTTAGTGACAAACCCCCTGCAACGCCCCCTGTTA‑TTAATTTTTTGGTGAC‑GGTTGGCGGGATCAGGTTTACt                           >  2:955881/1‑151 (MQ=255)
                                                                                                                     cgAAATGGCGCGGCGGTAAGTATGGC‑GGGGGTATTTCTT‑CCCCCCTGAGGGAACC‑CGGGTGTCCGGGTGACCATACGCCTTAATGACAAACCCGCCGCCACCCCCCCTGTTA‑TTAATTTTTTGGTGAC‑GTTTGGCGGGATCAGTTTTtctc                          >  1:1291078/1‑151 (MQ=255)
                                                                                                                      gAAATGGCGCGGCGGTAAGTATGGC‑GGGGGTATTTCTT‑TCCCGTTGAGGACACC‑GGGGTGTCTGGGTGGGCATATGCTTTAATGAGAAAACCCCTGGCACGCCCTCTGTTT‑TCTATTTTTTTGGGAC‑GGTTGggggggttcatttttctcc                         >  1:933075/1‑136 (MQ=255)
                                                                                                                        aaTGGCGCGGCGGTAAGTATGGC‑GGGGGTATTTCTTCCCCCGTTG‑GGGCACC‑GGGGTGTCAGGGTTACCATACGCTTAAGTGACAACCCCGCTGCAACCCCCCCTTTTA‑TTAAATTTCTTGTGAC‑GGTTGGGGGGATCAGTTTTACTCCGt                       >  2:40412/1‑151 (MQ=255)
                                                                                                                         aTGGCGCGGCGGTAAGTATGGC‑GGGGGTATTTCTT‑CCCCCCTTTGGACACCGGGGTTGTTGGGTGAACCCAAC‑CCTAAGTGACAACCCCCCGGCCAACCCCCCCGTTATTCAATTTT‑TTGCGAC‑GTTTTGCGGGATCCGGTTTTCttcctt                      >  2:394930/1‑146 (MQ=255)
                                                                                                                          tGGCGCGGCGGTAAGTATGGC‑GGGGGTTTTTCTT‑TCCCGCTTGGGACACC‑CGGTTGTCCGGGTGGCCAAAAGCTTAAATGGCAAACCCCCCGCCACCCCCCCTGTTT‑TTAAATTTTTTGGGGC‑GGTTTGCGGGATCCAtttttatccctg                      >  2:132841/1‑142 (MQ=255)
                                                                                                                           ggCGCGGCGGTAAGTATGGC‑GGGGTTATTCCTT‑CCCCGTTGAGGACACC‑GGGTTGTCAGGTTGACCATACGCTTAAGTGACAACCCCGCTGCAACGCCCTCTGTTA‑TCAATTTTCTGGTGAC‑GTTTGGCGGTATCAGTTTTACTCCGTGAc                    <  2:1385278/151‑1 (MQ=255)
                                                                                                                            gcgcGGCGGTAAGTATGGC‑GGGGTTATTCCTT‑CCCCGTTGAGGACACC‑GGGTTGTCAGGTTGACCATACGCTTAAGTGACAACCCCGCTGCAACGCCCTCTGTTA‑TCAATTTTCTGGTGAC‑GTTTGGCGGTATCAGTTTTACTCCGTGACt                   <  2:1402770/151‑1 (MQ=255)
                                                                                                                                       aGTATGGC‑GGGGTTATTCCTT‑CCCCGTTGAGGACACC‑GGGTTGTCAGGTTGACCATACGCTTAAGTGACAACCCCGCTGCAACGCCCTCTGTTA‑TCAATTTTCTGGTGAC‑GTTTGGCGGTATCAGTTTTACTCCGTGACTGCTCTgccgcc        <  2:692734/151‑1 (MQ=255)
                                                                                                                                        gTATGGC‑GGGGTTATTCCTT‑CCCCGTTGAGGACACC‑GGGTTGTCAGGTTGACCATACGCTTAAGTGACAACCCCGCTGCAACGCCCTCTGTTA‑TCAATTTTCTGGTGAC‑GTTTGGCGGTATCAGTTTTACTCCGTGACTGCTCTGCCGccc       >  2:335275/1‑151 (MQ=255)
                                                                                                                                          aTGGC‑GGGGTTATTCCTT‑CCCCGTTGAGGACACC‑GGGTTGTCAGGTTGACCATACGCTTAAGTGACAACCCCGCTGCAACGCCCTCTGTTA‑TCAATTTTCTGGTGAC‑GTTTGGCGGTATCAGTTTTACTCCGTGACTGCTCTGCCGCCCtt     <  1:293258/151‑1 (MQ=255)
                                                                                                                                            ggC‑GGGGTTATTCCTT‑CCCCGTTGAGGACACC‑GGGTTGTCAGGTTGACCATACGCTTAAGTGACAACCCCGCTGCAACGCCCTCTGTTA‑TCAATTTTCTGGTGAC‑GTTTGGCGGTATCAGTTTTACTCCGTGACTGCTCTGCCGCCCtttt   >  1:378741/1‑151 (MQ=255)
                                                                                                                                            ggC‑GGGGTTATTCCTT‑CCCCGTTGAGGACACC‑GGGTTGTCAGGTTGACCATACGCTTAAGTGACAACCCCGCTGCAACGCCCTCTGTTA‑TCAATTTTCTGGTGAC‑GTTTGGCGGTATCAGTTTTACTCCGTGACTGCTCTGCCGCCCttt    <  2:885202/150‑1 (MQ=255)
                                                                                                                                             gC‑GGGGTTATTCCTT‑CCCCGTTGAGGACACC‑GGGTTGTCAGGTTGACCATACGCTTAAGTGACAACCCCGCTGCAACGCCCTCTGTTA‑TCAATTTTCTGGTGAC‑GTTTGGCGGTATCAGTTTTACTCCGTGACTGCTCTGCCGCCCttttt  <  2:225740/151‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
AGGATATTTCCCTGAAAAGTCAGACCATAACGCGAAAGCGCATGGCGAGGTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGC‑GGGGTTATTCCTT‑CCCCGTTGAGGACACC‑GGGTTGTCAGGTTGACCATACGCTTAAGTGACAACCCCGCTGCAACGCCCTCTGTTA‑TCAATTTTCTGGTGAC‑GTTTGGCGGTATCAGTTTTACTCCGTGACTGCTCTGCCGCCCTTTTT  >  NC_000913/1649148‑1649439

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