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

AGGATATTTCCCTGAAAAGTCAGACCATAACGCGAAAGCGCATGGCGAGGTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGCGGGGTTATTCCTTCCCCGTTGAGGACACCGGGTTGTCAGGTTGACCATACGCTTAAGTGACAACCCCGCTGCAACGCCCTCTGTTATCAATTTTCTGGTGACGTTTGGCGGTATCAGTTTTACTCCGTGACTGCTCTGCCGCCCTTTTT  >  NC_000913/1649148‑1649439
                                                                                                                                                    |                                                                                                                                                    
aGGATATTTCCCTGAAAAGTCAGACCATAACGCGAAAGCGCATGGCGAGGTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGCGGGGTTAt                                                                                                                                                   <  2:525483/151‑1 (MQ=255)
  gATATTTCCCTGAAAAGTCAGACCATAACGCGAAAGCGCATGGCGAGGTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGCGGGGTTAttc                                                                                                                                                 <  2:1259898/151‑1 (MQ=255)
  gATATTTCCCTGAAAAGTCAGACCATAACGCGAAAGCGCATGGCGAGGTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGCGGGGTTAttc                                                                                                                                                 <  1:775874/151‑1 (MQ=255)
     atTTCCCTGAAAAGTCAGACCATAACGCGAAAGCGCATGGCGAGGTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCNGNNNNGGCAACCCGCGAAANGNCGCGGCGGTAAGTATGGCGGGGGTAtttctt                                                                                                                                              >  2:233376/1‑151 (MQ=255)
     atTTCCCTGAAAAGTCAGACCATAACGCGAAAGCGCATGGCGAGGTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGCGGGGGTGtttctt                                                                                                                                              >  1:873783/1‑151 (MQ=255)
          ccTGAAAAGTCAGACCATAACGCGAAAGCGCATGGCGAGGTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGCGGGGTTATTCCTTCCCCg                                                                                                                                        <  1:853905/151‑1 (MQ=255)
           cTGAAAAGTCAGACCATAACGCGAAAGCGCATGGCGAGGTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGCGGGGGTTTTCCTTtccccc                                                                                                                                       >  2:518087/1‑149 (MQ=255)
             gAAAAGTCAGACCATAACGCGAAAGCGCATGGCGAGGTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGCGGGGGTTTTTCTTcccccgtg                                                                                                                                     >  1:636593/1‑147 (MQ=255)
             gAAAAGTCAGACCATAACGCGAAAGCGCATGGCGAGGTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGCGGGGGTTTTCCTTcccccctg                                                                                                                                     >  1:1085150/1‑147 (MQ=255)
             gAAAAGTCAGACCATAACGCGAAAGCGCATGGCGAGGTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGCGGGGGTATTTCTTCCCCCTTg                                                                                                                                     >  2:1193359/1‑151 (MQ=255)
                    cAGACCATAACGCGAAAGCGCATGGCGAGGTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGCGGGGTTATTCCTTCCCCGTTGAGGacac                                                                                                                              <  2:988103/151‑1 (MQ=255)
                    cAGACCATAACGCGAAAGCGCATGGCGAGGTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGCGGGGGTTTTTCTTCCCCCTTGAGGacac                                                                                                                              >  2:1238362/1‑151 (MQ=255)
                    cAGACCATAACGCGAAAGCGCATGGCGAGGTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGCGGGGGTATTTCTTTCCCGTTGAGGgcac                                                                                                                              >  2:1182106/1‑151 (MQ=255)
                     aGACCATAACGCGAAAGCGCATGGCGAGGTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGCGGGGGTTTTTCTTCCCCCTTGAggggacc                                                                                                                             >  1:284828/1‑146 (MQ=255)
                        ccATAACGCGAAAGCGCATGGCGAGGTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGCGGGGTTATTCCTTCCCCGTTGAGGACACCggg                                                                                                                         <  1:1117958/151‑1 (MQ=255)
                        ccATAACGCGAAAGCGCATGGCGAGGTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGCGGGGGTTTTTCTTTCCCGTTGAGGacaaccgg                                                                                                                          >  2:293252/1‑146 (MQ=255)
                         cATAACGCGAAAGCGCATGGCGAGGTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGCGGGGTTATTCCTTCCCCGTTGAGGACACCGGGt                                                                                                                        <  1:337785/151‑1 (MQ=255)
                          aTAACGCGAAAGCGCATGGCGAGGTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGCGGGGTTATTCCTTCCCCGTTGAGGACACCGGGtt                                                                                                                       <  1:1263784/151‑1 (MQ=255)
                           tAACGCGAAAGCGCATGGCGAGGTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGCGGGGGTATTTCTTCCCCCGTGAGGGCACCCGGGTg                                                                                                                      >  2:607583/1‑151 (MQ=255)
                             aCGCGAAAGCGCATGGCGAGGTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGCGGGGTTATTCCTTCCCCGTTGAGGACACCGGGTTGTc                                                                                                                    <  2:185918/151‑1 (MQ=255)
                                cgAAAGCGCATGGCGAGGTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGCGGGGTTATTCCTTCCCCGTTGAGGACACCGGGTTGTCAgg                                                                                                                 <  2:592275/151‑1 (MQ=255)
                                 gAAAGCGCATGGCGAGGTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGCGGGGGTATTTCTTTCCCGTTTAGGGGACCGGGgtgttcggg                                                                                                                >  1:917253/1‑146 (MQ=255)
                                     gcgcATGGCGAGGTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGCGGGGTTATTCCTTCCCCGTTGAGGACACCGGGTTGTCAGGTTGAc                                                                                                            <  2:198266/151‑1 (MQ=255)
                                      cgcATGGCGAGGTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGCGGGGGTTTTCCTTCCCCCCTTTGGGGACCGGGGTGTCAGGGTGGcc                                                                                                           >  1:936392/1‑151 (MQ=255)
                                           ggCGAGGTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGCGGGGTTATTCCTTCCCCGTTGAGGACACCGGGTTGTCAGGTTGACCATAc                                                                                                       <  1:417245/150‑1 (MQ=255)
                                           ggCGAGGTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGCGGGGTTATTCCTTCCCCGTTGAGGACACCGGGTTGTCAGGTTGACCATAc                                                                                                       <  1:237650/150‑1 (MQ=255)
                                                 gTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGCGGGGTTATTCCTTCCCCGTTGAGGACACCGGGTTGTCAGGTTGACCATACGCTTAAg                                                                                                <  1:200340/151‑1 (MQ=255)
                                                          tCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGCGGGGGTATTTCCTCCCCCTTGAGGGCACCCGGGTGTCAGGGTTACCCTACGCTTTAGTGGCAAccc                                                                                       >  2:667803/1‑151 (MQ=255)
                                                              agatagCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGCGGGGTTATTCCTTCCCCGTTGAGGACACCGGGTTGTCAGGTTGACCATACGCTTAAGTGACAACCCCGc                                                                                    <  1:928652/150‑1 (MQ=255)
                                                               gatagCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGCGGGGTTATTCCTTCCCCGTTGAGGACACCggggaaggaataac                                                                                                              <  1:536844/123‑12 (MQ=255)
                                                                  agCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGCGGGGGTATTTCTTCCCCCTTGAGGGCACCCGGGTGTCAGGGTGACCATACGCTTAAATGACAAACCCCCTGcc                                                                                >  2:1088168/1‑149 (MQ=255)
                                                                          gTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGCGGGGGTATTTCTTTCCCGTTGAGGAAACCCGGGTGTTAGGGTGACCATACGCTTAAATGACAACCCCCCTGCAAAGcccccc                                                                       >  2:722164/1‑148 (MQ=255)
                                                                            tAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGCGGGGTTATTCCTTCCCCGTTGAGGACACCGGGTTGTCAGGTTGACCATACGCTTAAGTGACAACCCCGCTGCAACGCCCTCTg                                                                      <  2:239749/150‑1 (MQ=255)
                                                                             aaaTTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGCGGGGTTATTCCTTCCCCGTTGAGGACACCGGGTTGTCAGGTTGACCATACGCTTAAGTGACAACCCCGCTGCAACGCCCTCTGt                                                                     <  1:91848/150‑1 (MQ=255)
                                                                             aaaTTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGCGGGGGTATTTCTTTCCCCTTGAGGGCACCGGGGTGTTAGGGTGAACATAAGCTTAAGTGACAACCCCGCCGCAACGCCCCCTGtt                                                                    >  2:677909/1‑151 (MQ=255)
                                                                                 tttCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGCGGGGTTATTCCTTCCCCGTTGAGGACACCGGGTTGTCAGGTTGACCATACGCTTAAGTGACAACCCCGCTGCAACGCCCTCTGTTATc                                                                <  2:1067115/150‑1 (MQ=255)
                                                                                 tttCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGCGGGGGTATTTCTTTCCCGCTGGGGGCACCGGGGTGTTCGGGTGGCCATAACCTTAAATGAAAACCCCGCTGCCACGCCcccttgtttt                                                                 >  2:176634/1‑144 (MQ=255)
                                                                                  ttCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGCGGGGTTATTCCTTCCCCGTTGAGGACACCGGGTTGTCAGGTTGACCATACGCTTAAGTGACAACCCCGCTGCAACGCCCTCTGTTATCa                                                               <  1:934659/150‑1 (MQ=255)
                                                                                  ttCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGCGGGGGTTTTTCTTTCCCCCTGTGGGCACCCCGGTGTTCCGGTGACCACACCCTTTAGTGAAAAACCCCCCGCCACACCCCCTCTTTTCa                                                               >  1:1431749/1‑150 (MQ=255)
                                                                                    cgtcgACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGCGGGGTTATTTCTTTCCCCGTGTGGGCACCGGGGTGTCAGGGTGAGCATACGCCTAAAAGACAAACCCCCTGCCACGCCCCCTGTTATCAAtt                                                            >  2:48020/1‑151 (MQ=255)
                                                                                    cgtcgACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGCGGGGGTATTCCTTTCCCGTTGAGGACACCGGGGTGTCAGGGTGACCATACCCTTAAGTGACAACCCCCCCGCAACGCCCCCTGTTTTCAAtt                                                            >  2:144747/1‑151 (MQ=255)
                                                                                       cgACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGCGGGGTTATTCCTTCCCCGTTGAGGACACCGGGTTGTCAGGTTGACCATACGCTTAAGTGACAACCCCGCTGCAACGCCCTCTGTTATCAAtttt                                                          <  1:252280/150‑1 (MQ=255)
                                                                                        gACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGCGGGGGTATTTCTTCCCCCCTGTGGACACCCGGTTGTCAGGGTTTCCCTACGCTTTAGGGACAAACCCCCTGCAACCCCCCCTGTTTTTAAtttttt                                                        >  1:487898/1‑149 (MQ=255)
                                                                                        gACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGCGGGGGTATTTCTTCCCCCCTGGGGAGACCCGGGTGTTAGGGTGGCCATACGCTTAAGGGGCAACCCCCCCGCCACGCCCCCTTTTTTTCAtttttt                                                        >  1:170448/1‑149 (MQ=255)
                                                                                           cGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGCGGGGGTTTTTCTTCCCCCCTGTGGACACCCGGGTGTCAGGGTGACCATACGCTTAAGTGACAACCCCCCCGCCACCCCCCCTTTTTTTCAttttttttg                                                     >  1:1406688/1‑146 (MQ=255)
                                                                                            gTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGCGGGGTTATTCCTTCCCCGTTGAGGACACCGGGTTGTCAGGTTGACCATACGCTTAAGTGACAACCCCGCTGCAACGCCCTCTGTTATCAATTTTCTGGTg                                                    <  2:1095700/151‑1 (MQ=255)
                                                                                                    ccGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGCGGGGTTATTCCTTCCCCGTTGAGGACACCGGGTTGTCAGGTTGACCATACGCTTAAGTGACAACCCCGCTGCAACGCCCTCTGTTATCAATTTTCTGGTGACGTTTg                                            <  2:806662/150‑1 (MQ=255)
                                                                                                     cGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGCGGGGGTTTTTCTTTCCCCCTTAGGGCCCCCGGGTTTTACGTTGACCCTACGCTTTAGTGGCAACCCCCCCGCAACGCCCTCTGTTTTTCATTTTTTTGTGACGGTTTgg                                           >  1:291814/1‑151 (MQ=255)
                                                                                                            ggCAACCCGCGAAATGGCGCGGCGGTAAGTATGGCGGGGTTATTCCTTCCCCGTTGAGGACACCGGGTTGTCAGGTTGACCATACGCTTAAGTGACAACCCCGCTGCAACGCCCTCTGTTATCAATTTTCTGGTGACGTTTGGCGGTATc                                    <  2:690577/150‑1 (MQ=255)
                                                                                                            ggCAACCCGCGAAATGGCGCGGCGGTAAGTATGGCGGGGTTATTCCTTCCCCGTTGAGGACACCGGGTTGTCAGGTTGACCATACGCTTAAGTGACAACCCCGCTGCAACGCCCTCTGTTATCAATTTTCTGGTGACGTTTGGCGGTATc                                    <  2:870341/150‑1 (MQ=255)
                                                                                                            ggCAACCCGCGAAATGGCGCGGCGGTAAGTATGGCGGGGTTATTCCTTCCCCGTTGAGGACACCGGGTTGTCAGGTTGACCATACGCTTAAGTGACAACCCCGCTGCAACGCCCTCTGTTATCAATTTTCTGGTGACGTTTGGCGGTATc                                    <  2:873784/150‑1 (MQ=255)
                                                                                                            ggCAACCCGCGAAATGGCGCGGCGGTAAGTATGGCGGGGGTATTCCTTCCCCCGTTAGGCGACCCGGGTGTGGGGTTGCCCCTATGCCTTAGAGACAAACCCGGTGCACCCCCCTCTGTTATCTATTTTCTTGTGTCGATTTGCGGGATc                                    >  1:27768/1‑150 (MQ=255)
                                                                                                                 cccGCGAAATGGCGCGGCGGTAAGTATGGCGGGGTTATTCCTTCCCCGTTGAGGACACCGGGTTGTCAGGTTGACCATACGCTTAAGTGACAACCCCGCTGCAACGCCCTCTGTTATCAATTTTCTGGTGACGTTTGGCGGTATCAGtttt                              <  2:494625/151‑1 (MQ=255)
                                                                                                                  ccGCGAAATGGCGCGGCGGTAAGTATGGCGGGGGTTTTTCCTCCCCCTTGAGGGCACCCGGGTGTTAGGGTGACCATACGCCTAAATGACAAACCCCCCGCCACGCCCCCTGTTTTTAATTTTTTTGGGGCGGTTGGGGGGATCAGttttt                             >  1:1044569/1‑150 (MQ=255)
                                                                                                                  ccGCGAAATGGCGCGGCGGTAAGTATGGCGGGGGTATTCCTTCCCCGTTGAGGAAACCGGGGTGTCAGGTTGACCATACGCTTAAGTGACAACCCCGCTGCAACGCCCTCTGTTATCAATTTTTTTGTGACGGTTTGCGGGATTAGttttt                             >  2:722196/1‑150 (MQ=255)
                                                                                                                    gcgAAATGGCGCGGCGGTAAGTATGGCGGGGGTTTTTCTTCCCCCTTGAGGGCACCCGGGTGTTAGGGTGACCATACGCTTTAGTGACAAACCCCCCGCCAAGCCCCCTGTTTTTAATTTTTTTGGGGCGGTTGGGGGGATCAGgttttat                           >  2:960857/1‑148 (MQ=255)
                                                                                                                    gcgAAATGGCGCGGCGGTAAGTATGGCGGGGGTATTTCTTCCCCCCTTAGGGCACCCGGGTGTTAGGGTGGACATACGCTTTAGTGACAAACCCGCCGCCACGCCCCCCGTTTTTAAATTTTTTGTGGCGTTTTGCGGGATCAGgttttat                           >  1:923256/1‑148 (MQ=255)
                                                                                                                    gcgAAATGGCGCGGCGGTAAGTATGGCGGGGGTATTTCTTCCCCCCTGAGGACACCCGGGTGTCAGGGTGACCCTACGCTTTAGTGACAAACCCCCTGCAACGCCCCCTGTTATTAATTTTTTGGTGACGGTTGGCGGGATCAGGTTTACt                           >  2:955881/1‑151 (MQ=255)
                                                                                                                     cgAAATGGCGCGGCGGTAAGTATGGCGGGGGTATTTCTTCCCCCCTGAGGGAACCCGGGTGTCCGGGTGACCATACGCCTTAATGACAAACCCGCCGCCACCCCCCCTGTTATTAATTTTTTGGTGACGTTTGGCGGGATCAGTTTTtctc                          >  1:1291078/1‑151 (MQ=255)
                                                                                                                      gAAATGGCGCGGCGGTAAGTATGGCGGGGGTATTTCTTTCCCGTTGAGGACACCGGGGTGTCTGGGTGGGCATATGCTTTAATGAGAAAACCCCTGGCACGCCCTCTGTTTTCTATTTTTTTGGGACGGTTGggggggttcatttttctcc                         >  1:933075/1‑136 (MQ=255)
                                                                                                                        aaTGGCGCGGCGGTAAGTATGGCGGGGGTATTTCTTCCCCCGTTGGGGCACCGGGGTGTCAGGGTTACCATACGCTTAAGTGACAACCCCGCTGCAACCCCCCCTTTTATTAAATTTCTTGTGACGGTTGGGGGGATCAGTTTTACTCCGt                       >  2:40412/1‑151 (MQ=255)
                                                                                                                         aTGGCGCGGCGGTAAGTATGGCGGGGGTATTTCTTCCCCCCTTTGGACACCGGGGTTGTTGGGTGAACCCAACCCTAAGTGACAACCCCCCGGCCAACCCCCCCGTTATTCAATTTTTTGCGACGTTTTGCGGGATCCGGTTTTCttcctt                      >  2:394930/1‑146 (MQ=255)
                                                                                                                          tGGCGCGGCGGTAAGTATGGCGGGGGTTTTTCTTTCCCGCTTGGGACACCCGGTTGTCCGGGTGGCCAAAAGCTTAAATGGCAAACCCCCCGCCACCCCCCCTGTTTTTAAATTTTTTGGGGCGGTTTGCGGGATCCAtttttatccctg                      >  2:132841/1‑142 (MQ=255)
                                                                                                                           ggCGCGGCGGTAAGTATGGCGGGGTTATTCCTTCCCCGTTGAGGACACCGGGTTGTCAGGTTGACCATACGCTTAAGTGACAACCCCGCTGCAACGCCCTCTGTTATCAATTTTCTGGTGACGTTTGGCGGTATCAGTTTTACTCCGTGAc                    <  2:1385278/151‑1 (MQ=255)
                                                                                                                            gcgcGGCGGTAAGTATGGCGGGGTTATTCCTTCCCCGTTGAGGACACCGGGTTGTCAGGTTGACCATACGCTTAAGTGACAACCCCGCTGCAACGCCCTCTGTTATCAATTTTCTGGTGACGTTTGGCGGTATCAGTTTTACTCCGTGACt                   <  2:1402770/151‑1 (MQ=255)
                                                                                                                                       aGTATGGCGGGGTTATTCCTTCCCCGTTGAGGACACCGGGTTGTCAGGTTGACCATACGCTTAAGTGACAACCCCGCTGCAACGCCCTCTGTTATCAATTTTCTGGTGACGTTTGGCGGTATCAGTTTTACTCCGTGACTGCTCTgccgcc        <  2:692734/151‑1 (MQ=255)
                                                                                                                                        gTATGGCGGGGTTATTCCTTCCCCGTTGAGGACACCGGGTTGTCAGGTTGACCATACGCTTAAGTGACAACCCCGCTGCAACGCCCTCTGTTATCAATTTTCTGGTGACGTTTGGCGGTATCAGTTTTACTCCGTGACTGCTCTGCCGccc       >  2:335275/1‑151 (MQ=255)
                                                                                                                                          aTGGCGGGGTTATTCCTTCCCCGTTGAGGACACCGGGTTGTCAGGTTGACCATACGCTTAAGTGACAACCCCGCTGCAACGCCCTCTGTTATCAATTTTCTGGTGACGTTTGGCGGTATCAGTTTTACTCCGTGACTGCTCTGCCGCCCtt     <  1:293258/151‑1 (MQ=255)
                                                                                                                                            ggCGGGGTTATTCCTTCCCCGTTGAGGACACCGGGTTGTCAGGTTGACCATACGCTTAAGTGACAACCCCGCTGCAACGCCCTCTGTTATCAATTTTCTGGTGACGTTTGGCGGTATCAGTTTTACTCCGTGACTGCTCTGCCGCCCtttt   >  1:378741/1‑151 (MQ=255)
                                                                                                                                            ggCGGGGTTATTCCTTCCCCGTTGAGGACACCGGGTTGTCAGGTTGACCATACGCTTAAGTGACAACCCCGCTGCAACGCCCTCTGTTATCAATTTTCTGGTGACGTTTGGCGGTATCAGTTTTACTCCGTGACTGCTCTGCCGCCCttt    <  2:885202/150‑1 (MQ=255)
                                                                                                                                             gCGGGGTTATTCCTTCCCCGTTGAGGACACCGGGTTGTCAGGTTGACCATACGCTTAAGTGACAACCCCGCTGCAACGCCCTCTGTTATCAATTTTCTGGTGACGTTTGGCGGTATCAGTTTTACTCCGTGACTGCTCTGCCGCCCttttt  <  2:225740/151‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
AGGATATTTCCCTGAAAAGTCAGACCATAACGCGAAAGCGCATGGCGAGGTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGCGGGGTTATTCCTTCCCCGTTGAGGACACCGGGTTGTCAGGTTGACCATACGCTTAAGTGACAACCCCGCTGCAACGCCCTCTGTTATCAATTTTCTGGTGACGTTTGGCGGTATCAGTTTTACTCCGTGACTGCTCTGCCGCCCTTTTT  >  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: