Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913579,0040CT72.0% 67.7 / 38.4 51intergenic (‑111/+180)borD/ybcVprophage lipoprotein BorD/DUF1398 domain‑containing protein YbcV
Reads supporting (aligned to +/- strand):  ref base C (6/8);  new base T (16/20);  total (23/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.83e-01
Rejected as consensus: Frequency below/above cutoff threshold.

CGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGTCTGTTTGTTTTGAACAGTAAATGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACAATGAGTTGAAATTTCATATTGTTAATATTTATTAATGTATGCCAGGTGCGATGAATCGTCATTGTATTCCCGGATTAACTATGTCCACAGCCCTGACGGGGAACTTCTCTGCGGGAGTGTCCGGGAATAATTAAAAACGATGCACACAGGGTTTAGCGCGTACATGTATTGTATTATGCCAACACCCCGGTGCTGACACGGAAGAAACCGGACGTTATGATTTAGCGTGGAAAGATTTGTGTAGTGTTCTGAATGCTCTCAGTAAATAGTAATGAATTATCAAAGGTATAGT  >  NC_000913/578767‑579254
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                           
cGAAGAAATGATGGGTGATGGTTTCCTTTTGTGCTACTGCTGCCGGTTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACAATGAGTTGAAATTTCATATTGTTAATATTTATTAATGTATGTCAGGTGCGATGaa                                                                                                                                                                                                                                               <  2:90955/251‑1 (MQ=255)
cGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGTCTGTTTGTTTTGAACAGTAAATGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACAATGAGTTGAAATTTCATATTGTTAATATTTATTAATGTATGCCAGGTGCGATGaa                                                                                                                                                                                                                                               <  1:374001/251‑1 (MQ=255)
     aaaTGATGGGTGATTGTTTCCTTTGGTGCTACTGCTGTCTGTTTGTTTTGAACAGTAAATGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACAATGAGTTGAAATTTCATATTGTTAATATTTATTAATGTATGCCAGGTGCGATGAATCGTc                                                                                                                                                                                                                                          <  2:411173/251‑1 (MQ=255)
           tgGGTGATGGTTTCCTTTGGTGCTACTGCTGTCTGTTTGTTTTGAACAGTAAATGTCTGTTGCGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACAATGAGTTGAAATTTCATATTGTTAATATTTATTAATGTATGCCAGGTGCGATGAATCGTCATTGt                                                                                                                                                                                                                                     >  2:276833/1‑250 (MQ=255)
              gTGATGGTTTCCTTTGGTGCTACTGCTGCCGGTTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACAATGAGTTGAAATTTCATATTGTTAATATTTATTAATGTATGTCAGGGGCGATGAATCGTCATAGTATTc                                                                                                                                                                                                                                 <  1:124928/251‑1 (MQ=255)
                 aTGGTTTCCTTTGGTGCTACTGCTGCCGGTTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACAATGAGTTGAAATTTCATATTGTTAATATTTATTAATGTATGTCAGGTGCGATGAATCGTCATTGTATTCCCg                                                                                                                                                                                                                              <  1:365805/251‑1 (MQ=255)
                               tGCTACTGCTGCCGGTTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACAATGAGTTGAAATTTCATATTGTTAATATTTATTAATGTATGTCAGGTGCGATGAATCGTCATTGTATTCCCGGATTAACTATGTcc                                                                                                                                                                                                                >  2:118739/1‑251 (MQ=255)
                                 cTACTGCTGCCGGTTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACAATGAGTTGAAATTTCATATTGTTAATATTTATTAATGTATGTCAGGTGCGATGAATCGTCATTGTATTCCCGGATTAACTATGTCcac                                                                                                                                                                                                              >  1:239161/1‑251 (MQ=255)
                                  tACTGCTGCCGGTTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACAATGAGTTGAAATTTCATATTGTTAATATTTATTAATGTATGTCAGGTGCGATGAATCGTCATTGTATTCCCGGATTAACTATGTCcac                                                                                                                                                                                                              >  1:387102/1‑250 (MQ=255)
                                       ctgccggtttgtttTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACAATGAGTTGAAATTTCATATTGTTAATATTTATTAATGTATGTCAGGTGCGATGAATCGTCATTGTATTCCCGGATTAACTATGTCCACAGCCCt                                                                                                                                                                                                        >  1:161530/7‑251 (MQ=255)
                                                     tGAACAGTAAACGTCTGTTGAGCACAGCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCAGTTTTTTCATGGTGTTATTCCCGATGCTGTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACAAGGAGTTTAAATTTCATAGTGTTAAGATTTATTAATGTATGTCAGGTGCGATGAATCGTCATTGTATTCACGGATTAACTATGTCCACAGCCCTGACGGGGAACTtc                                                                                                                                                                                           <  2:253593/250‑1 (MQ=255)
                                                          aGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGTGGTTATTCCCGATTCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACAATGAGTTGAAATTTCATATTGTTAATATTTATTAATGTATGTCAGGTGCGATGAATCGTCATTGTATTCCCGGATTAACTATGTCCACAGCCCTGCCGGGGACCTTCTCTGCg                                                                                                                                                                                     >  1:64516/1‑251 (MQ=255)
                                                           ttAAATGTCTGTTGAGCACATCCTGTAATAAGAAGGGCAAGCGCAGTAGCGCGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTAGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACAATGAGTTGAAATTTCATATTGTTAATATTTATTAATGTATGCCAGGTGCGATGAATCGTCATTGTATTCCCGGATTAACTATGTCCACAGCCATGACGGGGAACTTCTCTGCgg                                                                                                                                                                                    >  2:402437/2‑251 (MQ=255)
                                                                 gTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGAAGCGAGTAGCATTTTTTTCATGGAGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACAATGCGTTGAAATTTCATATTGTTAATATTTATTAATGTATGTCAGGTGCGATGCATCGTCATTGTATTCCCGGATTAACTATGTCCACAGCACTGACGGGGAACTTCTCTGCGGGAgtgt                                                                                                                                                                              >  2:199568/1‑251 (MQ=255)
                                                                     gTTGAGCACATCCTGGAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTGTCATGGTGGTAGTCCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACAATGAGTTGAAATTTCATATTGTTAATATGTATGAATGTTTGTCAGGGGCGGTGAAGCGTCATTGTATTCCTGTATGAACTGGGTCCTCAGCCCTGACGGGGAACTTCGCTTCGGGAGTGTCCgg                                                                                                                                                                          <  2:452261/251‑1 (MQ=255)
                                                                       tGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACAATGAGTTGAAATTTCATATTGTTAATATTTATTAATGTATGCCAGGTGCGATGAATCGTCATTGTATTCCCGGATTAACTATGTCCACAGCCCTGACGGGGAACTTCTCTGCGGGAGTGTCCggg                                                                                                                                                                         <  1:10620/250‑1 (MQ=255)
                                                                       tGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACAATGAGTTGAAATTTCATATTGTTAATATTGAGTAATGTATGTCAGGTGCGATGAATCGTCATTGTATTCCCGGATTAACTATGTCCACAGCCCTGACGGGGAACTTCTCTGCGGGAGTGTCCggg                                                                                                                                                                         <  1:255415/250‑1 (MQ=255)
                                                                        gAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGAATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTCAACAAACCCTAAACAATGAGTTGAAATTTCATATTGTTAATATTTATTAATGTATGTCAGGAGCGATGAATCGTCATTGTATTCCCTGCTTAACTATGTCCACAGCCCTGACGGGGAACTTCTCTGCGGGAGTGTCCGGGa                                                                                                                                                                        >  1:451134/1‑250 (MQ=255)
                                                                         aGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTGTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACAATGAGTTGAAATTTCATATTGTTAAGATTTATTAATGTATGTCAGGTGCGATGAATCGTCATTGTATTCCCGGATTAACTATGTCCACAGCCCTGACGGGGAACTTCTCTGCGGGAGTGTCCTGGaat                                                                                                                                                                      <  1:131707/251‑1 (MQ=255)
                                                                            acaTCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACAATGAGTTGAAATTTCATATTGTTAATATTTATTAATGTATGTCAGGTGCGATGAATCGTCATTGTATTCCCGGATTAACTATGTCCACAGCCCTGACGGGGAACTTCTCTGCGGGAGTGTCCGGGaataat                                                                                                                                                                   >  1:52275/1‑251 (MQ=255)
                                                                                ccTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACAATGAGTTGAAATTTCATATTGTTAATATTTATTAATGTATGTCAGGTGCGATGAATCGTCATTGTATTCCCGGATTAACTATGTCCACAGCCCTGACGGGGAACTTCTCTGCGGGAGTGTCCGGGAATATTTaaa                                                                                                                                                               <  1:343349/251‑1 (MQ=255)
                                                                                      ataaGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACAATGAGTTGAAATTTCATATTGTTAATATTTATTAATGTATGTCAGGTGCGATGAATCGTCATTGTATTCCCGGATTAACTATGTCCACAGCCCTGACGGGGAACTTCTCTGCGGGAGTGTCCGGGAATAATTaagacgatg                                                                                                                                                         <  1:118739/251‑8 (MQ=255)
                                                                                              ggCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACAATGAGTTGAAATTTCATATTGTTAATATTTATTAATGTATGTCAGGTGCGATGAATCGTCATTGTATTCCCGGATTAACTATGTCCACAGCCCTGACGGGGAACTTCTCTGCGGGAGTGTCCGGGAATAATTAAAACGATGCACACAgg                                                                                                                                                <  1:342127/251‑1 (MQ=255)
                                                                                                ccAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACAATGAGTTGAAATTTCATATTGTTAATATTTATTAATTTATGTCAGGTGCGATGAAACGGCATTGTATTCCCGGATTAACTATGTCCACAGCCCTGACGGGGAACTTCTCTGCGGGAGTGTCCGGGAAGAATTAAAACGATGCACACAGGGt                                                                                                                                              <  1:344817/251‑1 (MQ=255)
                                                                                                     gcAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACAATGAGTTGAAATTTCATATTGTTAATATTTATTAATGTATGTCAGGTGCGATGAATCGTCATTGTATTCCCGGATTAACTATGTCCACAGCCCTGACGGGGAACTTCTCTGCGGGAGTGTCCGGGAATAATTAAAACGATGCACACAGGGTTTAgc                                                                                                                                         <  1:260912/251‑1 (MQ=255)
                                                                                                          aGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACAATGAGTTGAAATTTCATATTGTTAATATTTATTAATGTATGCCAGGTGCGATGAATCGTCATTGTATTCCCGGATTAACTATGTCCACAGCCCTGACGGGGAACTTCTCTGCGGGAGTGTCCGGGAATAATTAAAAACGATGCACACAGGGTTTAGcgcg                                                                                                                                      <  1:140465/250‑1 (MQ=255)
                                                                                                                tAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACAATGAGTTGAAATTTCATATTGTTAATATTTATTAATGTATGTCAGGTGCGATGAATCGTCATTGTATTCCCGGATTAACTATGTCCACAGCCCTGACGGGGAACTTCTCTGCGGGAGTGTCAGGGAATAATTAAAACGATGCACACAGGGTTTAGCGCGTACacgt                                                                                                                               >  1:355800/1‑247 (MQ=255)
                                                                                                                   cATTTTTTTCATGGTGTTATTCCCGATGCTTTGTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACAATGAGTTGAAATTTCATATTGTTAATATTTATTAATGTATGTCAGGTGCGATGAATCGTCATTGTATTCCCGGATTAACTATGTCCACAGCCCTGACGGGGAACTTCTCTGCGGGAGTGTCCGTGAATAATTAAAACGATGCACACAGGGTTTAGCGCGTACTcgtattg                                                                                                                           <  1:344377/251‑1 (MQ=255)
                                                                                                                   cATTTTTTTCAAGGTGTTATTCCCGATGATTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTCAACAAACACTAAACACTGAGTTGAAATTTCATATTGTTAATATTTATTAATGTATGTCAGGTGCGATGAATCGTCATTGTATTCCCGGATTAACTATGTCCACAGCCCTGACGGGGAACTTCTCTGCGGGAGTGTCCGGGAATAATTAAAACGATGCAAACAGGGTTTAGCGCGTACAcgtattg                                                                                                                           >  2:211615/1‑251 (MQ=255)
                                                                                                                    aTTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACAATGAGTTGAAATTTCATATTGTTAATATTTATTAATGTATGTCAGGTGCGATGAATCGTCATTGTATTCCCGGATTAACTATGTCCACAGCCCTGACGGGGAACTTCTCTGCGGGAGTGTCCGGGAATAATTAAAACGATGCACACAGGGTTTAGCGCGTACAcgtattgc                                                                                                                          >  1:352749/1‑250 (MQ=255)
                                                                                                                       tttttCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACAATGAGTTGAAATTTCATATTGTTAATATTTATTAATGTATGCCAGGTGCGATGAATCGTCATTGTATTCCCGGATTAACTATGTCCACAGCCCTGACGGGGAACTTCTCTGCGGGAGTGTCCGGGAATAATTAAAAACGATGCACACAGGGTTTAGCGCGTACAtgtattgta                                                                                                                         <  1:141710/250‑1 (MQ=255)
                                                                                                                                    tATTCCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAACATTAAACAAACCCTAAACAATGAGTTGCAATTTCATATTGTTAATATTTATTAATGTATGTCAGGTGCGATGAATCGTCATTGTATTCCCGGATTAACTATGTCCACAGCCCTGACGGGGAACTTCTCTGCGGGAGTGTCCGGGAATAATTAAAACGATGCACACAGGGTTTAGCGCGTACACGTATTGCATTATGCCAACGcccc                                                                                                          >  2:446367/1‑251 (MQ=255)
                                                                                                                                       tCCCGATGCTTTTTGAAGTTAGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACAATGAGTTGAAATTTCATATTGTTAATATTTATTAATGTATGTCAGGTGCGATGAATCGTCATTGTATTCCCGGATTAACTATGTCCACAGCCCTGACGGGGAACTTCTCTGCGGGAGTGTCCGGGAATAATTAAAACGATGCACACAGTGTTTAGCGCGTACACGTATTGCATTATGCCAACGCCCCgg                                                                                                        <  1:259991/250‑1 (MQ=255)
                                                                                                                                          cGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACAATGAGTTGAAATTTCATATTGTTAATATTTATTAATGTATGTCAGGTGCGATGAATCGTCATTGTATTCCCGGATTAACTATGTCCACAGCCCTGACGGGGAACTTCTCTGCGGGAGTGTCCGGGAATAATTAAAACGATGCACACATGGGTTAGCGCGTACACGTATTGCATTATGCCAACGCCCCGGTGCt                                                                                                    <  1:135359/251‑1 (MQ=255)
                                                                                                                                            aTGCTTGTTGAAGTTCTCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACAATGAGTTGAAATTTCATATTGTTAATATTTATTAATGTATGTCAGGTGCGATGAATCGTCATTGTATTCCCGGATTAACTATGTCCAGATCCCTGACGGGGAACTTCTCTGCTGGAGTGTCCGGGAATAATTAAAACGATGCACACAGGGTTTAGCGAGTACACGTATTGCATTATGCCAGCGCCCCGGTGCTGa                                                                                                  <  2:125565/251‑1 (MQ=255)
                                                                                                                                                       aGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACAATGAGTTGAAATTTCATATTGTTAATATTTATTAATGTATGCCAGGTGCGATGAATCGTCATTGTATTCCCGGATTAACTATGTCCACAGCCCTGACGGGGAACTTCTCTGCGGGAGTGTCCGGGAATAATTAAAAACGATGCACACAGGGTTTAGCGCGTACATGTATTGTATTATGCCAACACCCCGGTGCTGACACGgaaga                                                                                         >  1:180743/1‑250 (MQ=255)
                                                                                                                                                            gCAGAATCGTATGTGAAGAACATTAAAGAAACCCTAAACAATGAGTTGCACTTTCATATTGTTAATATTTATTAATGTATGTCAGGTGCGATGAATCGTCATTGTATTCCCGGATTAACTATGTCCACAGACCTGACGGGGAACTTCTCTGCGGGAGTGTCCGGGAATAATTAACACGATGCACACAGGGTTTAGCGAGTACACGTATTGAATTATGCCAACCCCCCGGCGCTGACACGGAAGAAACCgg                                                                                   >  2:276751/1‑250 (MQ=255)
                                                                                                                                                                     tATGTGTAGAAAATTAAACAAACCCTAAACAATGAGTTGAAATTTCATATTGTTAATATTTATTAATGTATGTCAGGTGCGATGAATCGTCATTGTATTCCCGGATTAACTATGTCCACAGCCCTGACGGGGAACTTCTCTGCGGGAGTGTCCGGGAATAATTAAAACGATGCACACAGGGTTTAGCGCGTACACGTATTGCATTATGCCAACGCCCCGGTGCTGACACTGAAGAAACCGGACGTTATGAt                                                                         >  1:203726/1‑251 (MQ=255)
                                                                                                                                                                        gtgtAGAAAATTAAACAAGCCCTAAACAAGGAGTTGAAATTTCATATTGTTAATATTTAGTAATGTATGTCGGGTCGAAGGAATCTGCATTGTATTCCCGTATGAACTATGTCCACATCTCTGACGGGGATGTTTTCTGCGGGAGTGTCCTTGAATTATTAAGTACGATGCCAACTGGGATTAGCGCGTACAGGTATTGTATTGTGGGGACACCCAGGTGGTGACACTGAAGAAACCgggcatggtgattt                                                                       <  2:274361/251‑13 (MQ=255)
                                                                                                                                                                            aGAAAATTAAACAAACCCTAAACAATGAGTTGAAATTTCATATTGTTAATATTTATTAATGTATGCCCGGTGCGATGAATCGTCATTGTATTCCCGGATTAACTATGTCCACAGCCCTGACGGGGAACTTCTCTGCGGGAGTGTCCGGGAATAATTAAAAACGATGCACACAGGGATTAGCGCGTACATGTATTGTATTATGCCAACACCCCGGTGCTGACACGGAAGAAACCGGACGTTATGATTTAGCg                                                                   >  2:439596/1‑251 (MQ=255)
                                                                                                                                                                                               aaaCAATGAGTTGAAATTTCATATTGTTAAGATGTATTAATGTATGCCAGGTGCGATGAATCGTCATTGTATTCCCGGATTAACTATGTCCACAGCCCTGACGGGGAACTTCTCTGCGGGAGTGTCCGGGAGTAATTAAAAACGATGCACACATGGTTTAGCGCGTACATGTATTTTATTATGCCAACACCCCGGTGCTGACGCGGAAGAAGCCGGACGTTATTATTAAGCGTGGAAAGATTTGTGTAgt                                                 <  2:389795/250‑1 (MQ=255)
                                                                                                                                                                                                   aaTGAGTTGAAATTTCATATTGTTAATATTTATTAATGTATGCCAGGTGCGATGAATCGTCATTGTATTCCCGGATTAACTATGTCCACAGCCCTGACGGGGAACTTCTCTGCGGGAGTGTCCGGGAATAATTAAAAACGATGCACACAGGGTTTAGCGCGTACATGTATTGTATTATGCCAACACCCCGGTGCTGACACGGAAGAAACCGGACGTTATGATTTAGCGTGGAAAGATTTGTGTAGTGTTCt                                            <  1:89786/251‑1 (MQ=255)
                                                                                                                                                                                                            aaaTTTCATATTGTTAATATTTATTAATGTATGTCAGGTGCGATGAAGCGTCATTGTATTCCCGGATTAACTATGTCCACAGCCCTGACGGGGAACGTCTCTGCGGGAGTGTCCGGGAATAATTAAAACGATGCACACAGGGTTTGGCGCGTACTCGTATTGCATTATGCCAACGCCCCGGGTCTGACACGGAAGAAACCGGACGGTAGGATATAGCGGGGAAAGATTTGTGTAGTGTTCTGAATGCtctc                                  <  1:433113/251‑1 (MQ=255)
                                                                                                                                                                                                                     atTGTTAATATTTATTAATGTATGCCAGGTGCGATGAATCGTCATTGTATTCCCGGATTAACTATGTCCACAGCCCTGACGGGGAACTTCTCTGCGGGAGTGTCCGGGAATAATTAAAAACGATGCACACAGGGTTTAGCGCGTACATGTATTGTATTATGCCAACACCCCGGTGCTGACACGGAAGAAACCGGACGTTATGATTTAGCGTGGAAAGATTTGTGTAGTGTTCTGAATGCTCTCAGTAAATa                          <  1:276833/251‑1 (MQ=255)
                                                                                                                                                                                                                      ttGTTAATATTTATTAATGTATGACAGGTGCGATGAATCGTCGTTGCATTCCCGACTTAACTATGTCCACAGCCCTGACGGGAACTTCTCTGCGGGAGCGTCCGGGAATACTTAAAAAGATGCACCCAGGGTTTAGCGCGTACACGGATTGCATTATGCCAACGCCCCGGTGCTGCCACGGAAGAAACCGGGCGTTATGATTTAGCGTGGAAAGATTTGTGTAGTGTTCTGAATGCTCTCAGTAAATAGt                        >  1:409097/1‑250 (MQ=255)
                                                                                                                                                                                                                            atatTTATTAATGTATGCCAGGTGCGATGAATCGTCATTGTATTCCCGGATTAACTATGTCCACAGCCCTGACGGGGAACTTCTCTGCGGGAGTGTCCGGGAATAATTAAAAACGATGCACACAGGGTTTAGCGCGTACATGTATTGTATTATGCCAACACCCCGGTGCTGACACGGAAGAAACCGGACGTTATGATTTAGCGTGGAAAGATTTGTGTAGTGTTCTGAATGCTCTCAGTAAATAGTAATg                    >  1:76041/1‑250 (MQ=255)
                                                                                                                                                                                                                              atttattAATGTATGTCAGGTGCGATGAATCGTCATTGTATTCCCGGATTAACTATGTCCACAGCCCTGACGGGGAACCTCTCTGCGGGAGTGTCCGGGAATAATTAAAACGATGCACACAGGGTTTAGCGCGTACACGTATTGCATTATGCCAACGCCCCGGTGCTGACACGGAAGAAACCGGACGTTATGATTTAGCGTGGAAAGATTTGTGTAGTGTTCTGAATGCTCTCAGTAAATAGTAATGAAtt                >  1:62750/1‑251 (MQ=255)
                                                                                                                                                                                                                                   ttaATGTATGTCAGGTGCGATGAATCGTCATTGTCTTCCCGGATTAACTATGTCCACAGCCCTGACGGGGAACTTCTCTGCGGGAGTGTCCGGGAATAATTAAAACGATGCACACAGGGTTTAGCGCGTACACGTATTGCATTATGCCAACGCCCCGGTGCTGACACGGAAGAAACCGGACGTTATGATTTAGCGTGGAAAGATTTGTGTAGTGTTCTGAATGCTCTCAGTAAATAGTa                       >  1:349836/1‑239 (MQ=255)
                                                                                                                                                                                                                                   ttaATGTATGTCAGGTGCGATGAATCGTCATTGTATTCCCGGATTAACTATGTCCACAGCCCTGACGGGGAACTTCTCTGCGGGAGTGTCCGGGAATAATTAAAACGATGCACACAGGGTTTAGCGCGTACACGTATTGCATTATGCCAACGCCCCGGTGCTGACACGGAAGAAACCGGACGTTATGATTTAGCGTGGAAAGATTTGTGTAGTGTTCTGAATGCTCTCAGTAAATAGTa                       <  2:349836/239‑1 (MQ=255)
                                                                                                                                                                                                                                     aaTGTATGTCAGGTGCGATGAATCGTCATTGTATTCCCGGATTAACTATGTCCACAGCCCTGACGGGGAACTTCTCTGCGGGAGTGTCCGGGAATAATTAAATCGATGCACACAGGGTTTAGCTCGTACACGTATTGCATTATGCCAACGCCCCGGTGCTGACACGGAAGAAACCGGACGTTATGATTTAGCGTGGAAAGATTTGTGTAGTGTTCTGAATGCTCTCAGTAAATAGTAATGAATTATCAAAg         <  2:387102/251‑1 (MQ=255)
                                                                                                                                                                                                                                      atgtatgtCAGGTGCGATGAATCGTCATTGTATTCCCGGATTAACTATGTCCACAGCCCTGACGGGGAACTTCTCTGCGGGAGTGTCCGGGAATAATTAAAACGATGCACACAGGGTTTAGCGCGTACACGTATTGCATTATGCCAACGCCCCGGTGCTGACACGGAAGAAACCGGACGTTATGATTTAGCGTGGAAAGATTTGTGTAGTGTTCTGAATGCTCTCAGTAAATAGTAATGAATTATCAAAgg        <  1:199568/251‑1 (MQ=255)
                                                                                                                                                                                                                                       tGTATGCCAGGTGCGATGAATCGTCATTGTATTCCCGGATTAACTATGTCCACAGCCCTGACGGGGAACTTCTCTGCGGGAGTGTCCGGGAATAATTAAAAACGATGCACACAGGGTTTAGCGCGTACATGTATTGTATTATGCCAACACCCCGGTGCTGACACGGAAGAAACCGGACGTTATGATTTAGCGTGGAAAGATTTGTGTAGTGTTCTGAATGCTCTCAGTAAATAGTAATGAATTATCAAAgg        >  1:309549/1‑251 (MQ=255)
                                                                                                                                                                                                                                             ccAGGTGCGATGAATCGTCATTGTATTCCCGGATTAACTATGTCCACAGCCCTGACGGGGAACTTCTCTGCGGGAGTGTCCGGGAATAATTAAAAACGATGCACACAGGGTTTAGCGCGTACATGTATTGTATTATGCCAACACCCCGGTGCTGACACGGAAGAAACCGGACGTTATGATTTAGCGTGGAAAGCTTTGTGTAGTGTTCTGAATGCTCTCAGTAAATAGTAATGAATTATCAAAGGTATAGt  >  2:91939/1‑251 (MQ=255)
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                           
CGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGTCTGTTTGTTTTGAACAGTAAATGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACAATGAGTTGAAATTTCATATTGTTAATATTTATTAATGTATGCCAGGTGCGATGAATCGTCATTGTATTCCCGGATTAACTATGTCCACAGCCCTGACGGGGAACTTCTCTGCGGGAGTGTCCGGGAATAATTAAAAACGATGCACACAGGGTTTAGCGCGTACATGTATTGTATTATGCCAACACCCCGGTGCTGACACGGAAGAAACCGGACGTTATGATTTAGCGTGGAAAGATTTGTGTAGTGTTCTGAATGCTCTCAGTAAATAGTAATGAATTATCAAAGGTATAGT  >  NC_000913/578767‑579254

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