Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,560,9850CA25.4% 109.9 / 11.2 63intergenic (+87/‑383)yffL/yffMCPZ‑55 prophage; uncharacterized protein/CPZ‑55 prophage; uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base C (23/24);  new base A (6/10);  total (29/34)
Fisher's exact test for biased strand distribution p-value = 5.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GATCGGCAGAAAGGGAGTGCAACGAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACA  >  NC_000913/2560837‑2561128
                                                                                                                                                    |                                                                                                                                               
gATCGGCAGAAAGGGAGTGCAACGAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTc                                                                                                                                                 <  3:27452/149‑1 (MQ=255)
      cAGAAAGGGAGTGCAACGAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTGAGTGGTGGGTCAAAAGt                                                                                                                                           >  7:347290/1‑149 (MQ=255)
       aGAAAGGGAGTGCAACGAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTGGCCGTGAGTGGTTTCTCAAAAGtt                                                                                                                                          >  3:331918/1‑149 (MQ=255)
          aaGGGAGTGCAACGAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGcc                                                                                                                                       <  7:332855/149‑1 (MQ=255)
            gggAGTGCAACGAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGt                                                                                                                                     <  1:218323/149‑1 (MQ=255)
            gggAGTGCAACGAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAATGGTGGGTAAAAATTTGcct                                                                                                                                      >  8:385595/1‑147 (MQ=255)
                gTGCAACGAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTATTGGTGGGTCAAAAGTTGCCGTTAAg                                                                                                                                 >  7:562947/1‑149 (MQ=255)
                   cAACGAATACCAGTTATTATTCAAGCCGTCCAACATGAAATAACGGAATCAGGCGTGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGt                                                                                                                               <  5:91707/148‑1 (MQ=255)
                   cAACGAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGt                                                                                                                               <  7:392862/148‑1 (MQ=255)
                    aaCGAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTgg                                                                                                                             >  4:358924/1‑149 (MQ=255)
                    aaCGAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGGCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTgg                                                                                                                             <  1:158206/149‑1 (MQ=255)
                      cGAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGAGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGt                                                                                                                           <  2:566412/149‑1 (MQ=255)
                            ccAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAATTTTGACTTTAAGGTGGGTCAAAAg                                                                                                                     >  7:334825/1‑149 (MQ=255)
                            ccAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCaaaa                                                                                                                      <  6:540138/148‑1 (MQ=255)
                             cAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGACGTGTAGGGGGGGCaaaatt                                                                                                                    >  2:433705/1‑147 (MQ=255)
                               gTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTATTGGTGGGTCAAATTTTGCCGTGTAGGAGTGTCAAttgttg                                                                                                                  >  3:96097/1‑149 (MQ=255)
                                  attattCAAGCCGTACAACATGAAATAACGGAATCAGTCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGGCAAAAGTTGcc                                                                                                                <  1:146569/148‑1 (MQ=255)
                                  attattCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAAAGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGAGGGTCAAAAGTTGCCGTTAGTGAGGGGGCAATAGTTGACATTTAGGGGTGTGAAAAGATGACg                                                                                                               >  3:57408/1‑149 (MQ=255)
                                     attCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTTAGGTGGGTCAAAAGTTgaagtg                                                                                                             >  1:97757/1‑143 (MQ=255)
                                        taaGCCGTACAACATGAAATAAAGGCATCCGGCTCGCAATTTTCACGCATTTTCTTTTCCTTTATTCCATGAATTCAGGTGGGGAATTTGTTGGCGTTATTTTTGGGTCAAATCTTTCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAg                                                                                                          <  6:269361/147‑1 (MQ=255)
                                         aaGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAGTTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAATGTTGCCGTTTAGGTGGGTCAAAAGTTGCCGTTAAGGt                                                                                                        >  1:271304/1‑149 (MQ=255)
                                         aaGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTgg                                                                                                                             <  1:86051/128‑1 (MQ=255)
                                         aaGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTgg                                                                                                                             <  1:77981/128‑1 (MQ=255)
                                         aaGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTgg                                                                                                                             >  2:86051/1‑128 (MQ=255)
                                         aaGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTgg                                                                                                                             >  2:77981/1‑128 (MQ=255)
                                         aaGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTGAGTGGTGGTTCAACAGTTGCCGTTAAGGTTGGTCAAAAATTCCCGtgtat                                                                                                          >  1:431630/1‑143 (MQ=255)
                                         aaGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGATAGTGTTGGATAAAAAGATGAATTTAAGGATTGTCAAATTTTGACGTCaaacg                                                                                                        >  8:70057/1‑146 (MQ=255)
                                           gCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTgg                                                                                                      <  3:358924/149‑1 (MQ=255)
                                                acaacaTGAAATAACGGAATCAGTCTGGCAATGTTCCCGCCTTTTGTTTGCCTTAATTCCATGAATTCAGGTGGGGCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGt                                                                                                    <  5:527835/146‑1 (MQ=255)
                                                acaacaTGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGGGGGGGTTCAAAAGTTGCCGAGAAGGTGGGTCAAAAGTTGGCGGTAAGGTTGGTCa                                                                                                  >  6:527835/1‑148 (MQ=255)
                                                 caacaTGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTAGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCaa                                                                                                 <  1:71070/148‑1 (MQ=255)
                                                  atcctGAAATTACGGCACCCGTAGGCCATTGTTCCCGCCTTTTCTTTGCCTTAATTCCATGCATTCAGGTGGGCCAAAAGATGGCCTTTGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCaaaa                                                                                               <  8:452758/145‑1 (MQ=21)
                                                  aacaTGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGGGGTGCGGCAAAAGTGGCCGGTAAGGTGGGGCAAAAGTTGCCGTGaa                                                                                                           >  5:84381/1‑137 (MQ=255)
                                                    caTGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGAGTGTCAAAAGTTGGCGTTAAGGTGGGTGAAATGt                                                                                             >  7:559028/1‑149 (MQ=39)
                                                      tGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTaa                                                                                                                                  >  5:362004/1‑110 (MQ=255)
                                                      tGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGACGTTAGTGGTGGGTCAAAAGTTGCCGTTaa                                                                                                                                  <  6:362004/110‑1 (MQ=255)
                                                       gAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGTTGGGTAAAATGTTGCCTTTAAGGGGTGTCAAAATATGACGTTAAGGGGGGTGaaatgatgg                                                                                          >  1:564802/1‑143 (MQ=255)
                                                         aaTAACGGAATCAGGCGGGCAATGTTCCCGCGTTTTCTTTGCCTTAATTCCATGAATTCAGGAGGGTCAAAAGTTGCCGTTTGTGGTGGGTAATATGTTGCCTTTAAGGTGTGTCAAAAGTTGCAGTTAAGGTGTGTGAATAGttt                                                                                           >  7:209068/1‑145 (MQ=255)
                                                                                     cGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGTTGGGTCAAAAGTTGCAGTTAAGGAGGGTCAAAAGTTTCAGTGAACGTGGGGGAAAATTTGCAGAGAAAGAATTCTACAATTAAGAATTCTa                                                            >  7:234146/1‑149 (MQ=25)
                                                                                     agCCATTTCTTTTCCTTAAATCCATGAATTCATTGTTCTCAAAATTAGCCCTTATGGGTGGGTCAAAAGTTGCCTTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTa                                                            <  8:100993/148‑1 (MQ=21)
                                                                                       ccTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTAAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCAGTTAAGGTGGGTGAAAAGTTGCCGAGTAAGAATGCGAGTATTAATAATTCTAgg                                                          >  7:458496/1‑149 (MQ=37)
                                                                                        cTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGAGGGTCAAAATTTGCAGTTAAGGTGGGGGAAAATTTGAAGAGAAAGAATTCTACTATTAATAATTCTAg                                                           >  8:562629/1‑146 (MQ=25)
                                                                                                ttccTTAATTTAATGAATTCAGGGGGGTTAAAAGAAGGCTTTATGGGGGGGTAAAAAGCTGCCTTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGa                                                  <  8:334825/146‑1 (MQ=25)
                                                                                                tGCCTTAATTCAATGAATTCAGGTGGGTTCAAATTTGCCCTTATTGGTGGGTAAAAAGATGCCGTTAAGGTGGGTCAAAAGTTGCCGTTCAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGa                                                  <  4:436548/148‑1 (MQ=37)
                                                                                                 gCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCAGTTAAGGTGGGTGAAATTTTGAAGAGAAAGAATTCTTCTATTAAGAATACTAGGCGACGGGAGt                                                >  7:481914/1‑149 (MQ=25)
                                                                                                       attgaaTGAATTCAGGTGGGTATTAATTAGGCGTTATTTGTGGCTAATTTGTTGCCTTTTAGGTGGGTAAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTa                                          <  4:871/144‑1 (MQ=25)
                                                                                                       aTTCAATGAATTCTGCGTGGAAAAAATATGGCGTTACTGGGGGCTCAAAAGAGGCCGTTAATGTGCGTCAAAAGTTTCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTa                                          <  6:6572/149‑1 (MQ=21)
                                                                                                        atcaatcaaTTCATGTCGGTCAAAAGACGCCCTTACTGGGGGGTAAAAAGGTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTa                                          <  1:433705/141‑1 (MQ=25)
                                                                                                          ccATGAATTCAGGGGGGTCAAATGCTGCCGTTAGTGGTGGGTAAAAAGTTGCCTTTAAGGTGGGTCAAAACTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGc                                        <  7:385595/148‑1 (MQ=38)
                                                                                                            aTGAATTCAGGTGGGTCAAAAGTTGCCGTTAATGTTGTGTAAAAATTTGCATTATAGGAGACTCAAATTTTGAAGTATAGGTGGCTCAAATTTTGACTAGTAAGAAGTCTAAAATTAACATTAATAGGCGACGGGAGtgttattgact                                      >  8:99922/1‑141 (MQ=255)
                                                                                                                         gctCAAAAGCAGCCCTTAGGGGTGGCTAAAACGGTGCCGTTAAGGTGGGTCAAACCTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATa                        <  1:484966/147‑1 (MQ=37)
                                                                                                                          gTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTGTAGGTGGGTGAAAATTTGCCGTTAAGGTGGCGCAAAGGTTGCGGAGAAAGAAGTCTAAAATTAAGAATAGTAGGCGACTGGAGTTTTGTAGACTAACGGTCGGTCa                          >  2:339201/1‑146 (MQ=19)
                                                                                                                             caaaGTAGCCCTTATTGTTGGTTATAAAGCTGCCGTTAAGTTGGGTCAAAATTTGCGTTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATagaga                    <  5:216942/148‑1 (MQ=255)
                                                                                                                                tgTTACCGTTATTGGTGGGTCAAATCTTGCCTTTAAGGTGGGTCAAAAGTTGCCGTTAAGGGGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGTCGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTc                 <  3:387562/148‑1 (MQ=38)
                                                                                                                                       tttATTGGTGGGTAAAACGTTGCAGTTATGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTaga          <  8:481914/148‑1 (MQ=38)
                                                                                                                                       gTTAGTGTTGGGTCAAAAGTTGCCGTTAAGGTGGGGCAAACGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTaga          <  1:529385/149‑1 (MQ=255)
                                                                                                                                       cttacgggtggGTAAAAACGTGCCGTTAAGGTGGGTAAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTaga          <  6:97155/143‑1 (MQ=38)
                                                                                                                                         tattggtggGTAAAAACTTGCCTTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTaga          <  3:165775/144‑1 (MQ=255)
                                                                                                                                          agaggtggGTCAAAACATGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGGCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGt        <  3:59874/145‑1 (MQ=255)
                                                                                                                                           gtggtgGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAAGTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCATCCAGGTCGGTCATAGAGAGTCGCCTagagat       >  1:522310/1‑146 (MQ=255)
                                                                                                                                            tggtggGTCAAAAGTTCCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTa      <  7:562629/148‑1 (MQ=255)
                                                                                                                                               tggGACAATAGTTGCCGTTTAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACa  <  8:562947/149‑1 (MQ=255)
                                                                                                                                                  gTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCa                          >  4:405091/1‑122 (MQ=255)
                                                                                                                                                  gTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCa                          <  3:405091/122‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                               
GATCGGCAGAAAGGGAGTGCAACGAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACA  >  NC_000913/2560837‑2561128

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