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.

CGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGTCTGTTTGTTTTGAACAGTAAATGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACAATGAGTTGAAATTTCATATTGTTAATATTTATTAATGTATGCCAGGTGCGA‑TGAATCGTCATTGTATTCCCGGATTAACTATGTCCACAGCCCTGACGGGGAACTTCTCTGCGGGAGTGTCCGGGAATAATTAAAAACGATGCACACAGGGTTTAGCGCGTACATGTATTGTATTATGCCAACACCCCGGTGCTGACACGGAAGAAACCGGACGTTATGATTTAGCGTGGAAAGATTTGTGTAGTGTTCTGAATGCTCTCAGTAAATAGTAATGAATTATCAAAGGTATAGT  >  NC_000913/578767‑579254
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                           
cGAAGAAATGATGGGTGATGGTTTCCTTTTGTGCTACTGCTGCCGGTTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACAATGAGTTGAAATTTCATATTGTTAATATTTATTAATGTATGTCAGGTGCGA‑TGaa                                                                                                                                                                                                                                               <  2:90955/251‑1 (MQ=255)
cGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGTCTGTTTGTTTTGAACAGTAAATGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACAATGAGTTGAAATTTCATATTGTTAATATTTATTAATGTATGCCAGGTGCGA‑TGaa                                                                                                                                                                                                                                               <  1:374001/251‑1 (MQ=255)
     aaaTGATGGGTGATTGTTTCCTTTGGTGCTACTGCTGTCTGTTTGTTTTGAACAGTAAATGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACAATGAGTTGAAATTTCATATTGTTAATATTTATTAATGTATGCCAGGTGCGA‑TGAATCGTc                                                                                                                                                                                                                                          <  2:411173/251‑1 (MQ=255)
           tgGGTGATGGTTTCCTTTGGTGCTACTGCTGTCTGTTTGTTTTGAACAGTAAATGTCTGTTGCGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACAATGAGTTGAAATTTCATATTGTTAATATTTATTAATGTATGCCAGGTGCGA‑TGAATCGTCATTGt                                                                                                                                                                                                                                     >  2:276833/1‑250 (MQ=255)
              gTGATGGTTTCCTTTGGTGCTACTGCTGCCGGTTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACAATGAGTTGAAATTTCATATTGTTAATATTTATTAATGTATGTCAGGGGCGA‑TGAATCGTCATAGTATTc                                                                                                                                                                                                                                 <  1:124928/251‑1 (MQ=255)
                 aTGGTTTCCTTTGGTGCTACTGCTGCCGGTTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACAATGAGTTGAAATTTCATATTGTTAATATTTATTAATGTATGTCAGGTGCGA‑TGAATCGTCATTGTATTCCCg                                                                                                                                                                                                                              <  1:365805/251‑1 (MQ=255)
                               tGCTACTGCTGCCGGTTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACAATGAGTTGAAATTTCATATTGTTAATATTTATTAATGTATGTCAGGTGCGA‑TGAATCGTCATTGTATTCCCGGATTAACTATGTcc                                                                                                                                                                                                                >  2:118739/1‑251 (MQ=255)
                                 cTACTGCTGCCGGTTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACAATGAGTTGAAATTTCATATTGTTAATATTTATTAATGTATGTCAGGTGCGA‑TGAATCGTCATTGTATTCCCGGATTAACTATGTCcac                                                                                                                                                                                                              >  1:239161/1‑251 (MQ=255)
                                  tACTGCTGCCGGTTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACAATGAGTTGAAATTTCATATTGTTAATATTTATTAATGTATGTCAGGTGCGA‑TGAATCGTCATTGTATTCCCGGATTAACTATGTCcac                                                                                                                                                                                                              >  1:387102/1‑250 (MQ=255)
                                       ctgccggtttgtttTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACAATGAGTTGAAATTTCATATTGTTAATATTTATTAATGTATGTCAGGTGCGA‑TGAATCGTCATTGTATTCCCGGATTAACTATGTCCACAGCCCt                                                                                                                                                                                                        >  1:161530/7‑251 (MQ=255)
                                                     tGAACAGTAAACGTCTGTTGAGCACAGCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCAGTTTTTTCATGGTGTTATTCCCGATGCTGTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACAAGGAGTTTAAATTTCATAGTGTTAAGATTTATTAATGTATGTCAGGTGCGA‑TGAATCGTCATTGTATTCACGGATTAACTATGTCCACAGCCCTGACGGGGAACTtc                                                                                                                                                                                           <  2:253593/250‑1 (MQ=255)
                                                          aGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGTGGTTATTCCCGATTCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACAATGAGTTGAAATTTCATATTGTTAATATTTATTAATGTATGTCAGGTGCGA‑TGAATCGTCATTGTATTCCCGGATTAACTATGTCCACAGCCCTGCCGGGGACCTTCTCTGCg                                                                                                                                                                                     >  1:64516/1‑251 (MQ=255)
                                                           ttAAATGTCTGTTGAGCACATCCTGTAATAAGAAGGGCAAGCGCAGTAGCGCGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTAGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACAATGAGTTGAAATTTCATATTGTTAATATTTATTAATGTATGCCAGGTGCGA‑TGAATCGTCATTGTATTCCCGGATTAACTATGTCCACAGCCATGACGGGGAACTTCTCTGCgg                                                                                                                                                                                    >  2:402437/2‑251 (MQ=255)
                                                                 gTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGAAGCGAGTAGCATTTTTTTCATGGAGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACAATGCGTTGAAATTTCATATTGTTAATATTTATTAATGTATGTCAGGTGCGA‑TGCATCGTCATTGTATTCCCGGATTAACTATGTCCACAGCACTGACGGGGAACTTCTCTGCGGGAgtgt                                                                                                                                                                              >  2:199568/1‑251 (MQ=255)
                                                                     gTTGAGCACATCCTGGAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTGTCATGGTGGTAGTCCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACAATGAGTTGAAATTTCATATTGTTAATATGTATGAATGTTTGTCAGGGGCGG‑TGAAGCGTCATTGTATTCCTGTATGAACTGGGTCCTCAGCCCTGACGGGGAACTTCGCTTCGGGAGTGTCCgg                                                                                                                                                                          <  2:452261/251‑1 (MQ=255)
                                                                       tGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACAATGAGTTGAAATTTCATATTGTTAATATTTATTAATGTATGCCAGGTGCGA‑TGAATCGTCATTGTATTCCCGGATTAACTATGTCCACAGCCCTGACGGGGAACTTCTCTGCGGGAGTGTCCggg                                                                                                                                                                         <  1:10620/250‑1 (MQ=255)
                                                                       tGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACAATGAGTTGAAATTTCATATTGTTAATATTGAGTAATGTATGTCAGGTGCGA‑TGAATCGTCATTGTATTCCCGGATTAACTATGTCCACAGCCCTGACGGGGAACTTCTCTGCGGGAGTGTCCggg                                                                                                                                                                         <  1:255415/250‑1 (MQ=255)
                                                                        gAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGAATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTCAACAAACCCTAAACAATGAGTTGAAATTTCATATTGTTAATATTTATTAATGTATGTCAGGAGCGA‑TGAATCGTCATTGTATTCCCTGCTTAACTATGTCCACAGCCCTGACGGGGAACTTCTCTGCGGGAGTGTCCGGGa                                                                                                                                                                        >  1:451134/1‑250 (MQ=255)
                                                                         aGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTGTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACAATGAGTTGAAATTTCATATTGTTAAGATTTATTAATGTATGTCAGGTGCGA‑TGAATCGTCATTGTATTCCCGGATTAACTATGTCCACAGCCCTGACGGGGAACTTCTCTGCGGGAGTGTCCTGGaat                                                                                                                                                                      <  1:131707/251‑1 (MQ=255)
                                                                            acaTCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACAATGAGTTGAAATTTCATATTGTTAATATTTATTAATGTATGTCAGGTGCGA‑TGAATCGTCATTGTATTCCCGGATTAACTATGTCCACAGCCCTGACGGGGAACTTCTCTGCGGGAGTGTCCGGGaataat                                                                                                                                                                   >  1:52275/1‑251 (MQ=255)
                                                                                ccTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACAATGAGTTGAAATTTCATATTGTTAATATTTATTAATGTATGTCAGGTGCGA‑TGAATCGTCATTGTATTCCCGGATTAACTATGTCCACAGCCCTGACGGGGAACTTCTCTGCGGGAGTGTCCGGGAATATTTaaa                                                                                                                                                               <  1:343349/251‑1 (MQ=255)
                                                                                      ataaGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACAATGAGTTGAAATTTCATATTGTTAATATTTATTAATGTATGTCAGGTGCGA‑TGAATCGTCATTGTATTCCCGGATTAACTATGTCCACAGCCCTGACGGGGAACTTCTCTGCGGGAGTGTCCGGGAATAATTaagacgatg                                                                                                                                                         <  1:118739/251‑8 (MQ=255)
                                                                                              ggCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACAATGAGTTGAAATTTCATATTGTTAATATTTATTAATGTATGTCAGGTGCGA‑TGAATCGTCATTGTATTCCCGGATTAACTATGTCCACAGCCCTGACGGGGAACTTCTCTGCGGGAGTGTCCGGGAATAATTAAAA‑CGATGCACACAgg                                                                                                                                                <  1:342127/251‑1 (MQ=255)
                                                                                                ccAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACAATGAGTTGAAATTTCATATTGTTAATATTTATTAATTTATGTCAGGTGCGA‑TGAAACGGCATTGTATTCCCGGATTAACTATGTCCACAGCCCTGACGGGGAACTTCTCTGCGGGAGTGTCCGGGAAGAATTAAAA‑CGATGCACACAGGGt                                                                                                                                              <  1:344817/251‑1 (MQ=255)
                                                                                                     gcAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACAATGAGTTGAAATTTCATATTGTTAATATTTATTAATGTATGTCAGGTGCGA‑TGAATCGTCATTGTATTCCCGGATTAACTATGTCCACAGCCCTGACGGGGAACTTCTCTGCGGGAGTGTCCGGGAATAATTAAAA‑CGATGCACACAGGGTTTAgc                                                                                                                                         <  1:260912/251‑1 (MQ=255)
                                                                                                          aGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACAATGAGTTGAAATTTCATATTGTTAATATTTATTAATGTATGCCAGGTGCGA‑TGAATCGTCATTGTATTCCCGGATTAACTATGTCCACAGCCCTGACGGGGAACTTCTCTGCGGGAGTGTCCGGGAATAATTAAAAACGATGCACACAGGGTTTAGcgcg                                                                                                                                      <  1:140465/250‑1 (MQ=255)
                                                                                                                tAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACAATGAGTTGAAATTTCATATTGTTAATATTTATTAATGTATGTCAGGTGCGA‑TGAATCGTCATTGTATTCCCGGATTAACTATGTCCACAGCCCTGACGGGGAACTTCTCTGCGGGAGTGTCAGGGAATAATTAAAA‑CGATGCACACAGGGTTTAGCGCGTACacgt                                                                                                                               >  1:355800/1‑247 (MQ=255)
                                                                                                                   cATTTTTTTCATGGTGTTATTCCCGATGCTTTGTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACAATGAGTTGAAATTTCATATTGTTAATATTTATTAATGTATGTCAGGTGCGA‑TGAATCGTCATTGTATTCCCGGATTAACTATGTCCACAGCCCTGACGGGGAACTTCTCTGCGGGAGTGTCCGTGAATAATTAAAA‑CGATGCACACAGGGTTTAGCGCGTACTcgtattg                                                                                                                           <  1:344377/251‑1 (MQ=255)
                                                                                                                   cATTTTTTTCAAGGTGTTATTCCCGATGATTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTCAACAAACACTAAACACTGAGTTGAAATTTCATATTGTTAATATTTATTAATGTATGTCAGGTGCGA‑TGAATCGTCATTGTATTCCCGGATTAACTATGTCCACAGCCCTGACGGGGAACTTCTCTGCGGGAGTGTCCGGGAATAATTAAAA‑CGATGCAAACAGGGTTTAGCGCGTACAcgtattg                                                                                                                           >  2:211615/1‑251 (MQ=255)
                                                                                                                    aTTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACAATGAGTTGAAATTTCATATTGTTAATATTTATTAATGTATGTCAGGTGCGA‑TGAATCGTCATTGTATTCCCGGATTAACTATGTCCACAGCCCTGACGGGGAACTTCTCTGCGGGAGTGTCCGGGAATAATTAAAA‑CGATGCACACAGGGTTTAGCGCGTACAcgtattgc                                                                                                                          >  1:352749/1‑250 (MQ=255)
                                                                                                                       tttttCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACAATGAGTTGAAATTTCATATTGTTAATATTTATTAATGTATGCCAGGTGCGA‑TGAATCGTCATTGTATTCCCGGATTAACTATGTCCACAGCCCTGACGGGGAACTTCTCTGCGGGAGTGTCCGGGAATAATTAAAAACGATGCACACAGGGTTTAGCGCGTACAtgtattgta                                                                                                                         <  1:141710/250‑1 (MQ=255)
                                                                                                                                    tATTCCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAACATTAAACAAACCCTAAACAATGAGTTGCAATTTCATATTGTTAATATTTATTAATGTATGTCAGGTGCGA‑TGAATCGTCATTGTATTCCCGGATTAACTATGTCCACAGCCCTGACGGGGAACTTCTCTGCGGGAGTGTCCGGGAATAATTAAAA‑CGATGCACACAGGGTTTAGCGCGTACACGTATTGCATTATGCCAACGcccc                                                                                                          >  2:446367/1‑251 (MQ=255)
                                                                                                                                       tCCCGATGCTTTTTGAAGTTAGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACAATGAGTTGAAATTTCATATTGTTAATATTTATTAATGTATGTCAGGTGCGA‑TGAATCGTCATTGTATTCCCGGATTAACTATGTCCACAGCCCTGACGGGGAACTTCTCTGCGGGAGTGTCCGGGAATAATTAAAA‑CGATGCACACAGTGTTTAGCGCGTACACGTATTGCATTATGCCAACGCCCCgg                                                                                                        <  1:259991/250‑1 (MQ=255)
                                                                                                                                          cGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACAATGAGTTGAAATTTCATATTGTTAATATTTATTAATGTATGTCAGGTGCGA‑TGAATCGTCATTGTATTCCCGGATTAACTATGTCCACAGCCCTGACGGGGAACTTCTCTGCGGGAGTGTCCGGGAATAATTAAAA‑CGATGCACACATGGGTTAGCGCGTACACGTATTGCATTATGCCAACGCCCCGGTGCt                                                                                                    <  1:135359/251‑1 (MQ=255)
                                                                                                                                            aTGCTTGTTGAAGTTCTCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACAATGAGTTGAAATTTCATATTGTTAATATTTATTAATGTATGTCAGGTGCGA‑TGAATCGTCATTGTATTCCCGGATTAACTATGTCCAGATCCCTGACGGGGAACTTCTCTGCTGGAGTGTCCGGGAATAATTAAAA‑CGATGCACACAGGGTTTAGCGAGTACACGTATTGCATTATGCCAGCGCCCCGGTGCTGa                                                                                                  <  2:125565/251‑1 (MQ=255)
                                                                                                                                                       aGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACAATGAGTTGAAATTTCATATTGTTAATATTTATTAATGTATGCCAGGTGCGA‑TGAATCGTCATTGTATTCCCGGATTAACTATGTCCACAGCCCTGACGGGGAACTTCTCTGCGGGAGTGTCCGGGAATAATTAAAAACGATGCACACAGGGTTTAGCGCGTACATGTATTGTATTATGCCAACACCCCGGTGCTGACACGgaaga                                                                                         >  1:180743/1‑250 (MQ=255)
                                                                                                                                                            gCAGAATCGTATGTGAAGAACATTAAAGAAACCCTAAACAATGAGTTGCACTTTCATATTGTTAATATTTATTAATGTATGTCAGGTGCGA‑TGAATCGTCATTGTATTCCCGGATTAACTATGTCCACAGACCTGACGGGGAACTTCTCTGCGGGAGTGTCCGGGAATAATTAACA‑CGATGCACACAGGGTTTAGCGAGTACACGTATTGAATTATGCCAACCCCCCGGCGCTGACACGGAAGAAACCgg                                                                                   >  2:276751/1‑250 (MQ=255)
                                                                                                                                                                     tATGTGTAGAAAATTAAACAAACCCTAAACAATGAGTTGAAATTTCATATTGTTAATATTTATTAATGTATGTCAGGTGCGA‑TGAATCGTCATTGTATTCCCGGATTAACTATGTCCACAGCCCTGACGGGGAACTTCTCTGCGGGAGTGTCCGGGAATAATTAAAA‑CGATGCACACAGGGTTTAGCGCGTACACGTATTGCATTATGCCAACGCCCCGGTGCTGACACTGAAGAAACCGGACGTTATGAt                                                                         >  1:203726/1‑251 (MQ=255)
                                                                                                                                                                        gtgtAGAAAATTAAACAAGCCCTAAACAAGGAGTTGAAATTTCATATTGTTAATATTTAGTAATGTATGTCGGGT‑CGAAGGAATCTGCATTGTATTCCCGTATGAACTATGTCCACATCTCTGACGGGGATGTTTTCTGCGGGAGTGTCCTTGAATTATTAAGTACGATGCCAACTGGGATTAGCGCGTACAGGTATTGTATTGTGGGGACACCCAGGTGGTGACACTGAAGAAACCgggcatggtgattt                                                                       <  2:274361/251‑13 (MQ=255)
                                                                                                                                                                            aGAAAATTAAACAAACCCTAAACAATGAGTTGAAATTTCATATTGTTAATATTTATTAATGTATGCCCGGTGCGA‑TGAATCGTCATTGTATTCCCGGATTAACTATGTCCACAGCCCTGACGGGGAACTTCTCTGCGGGAGTGTCCGGGAATAATTAAAAACGATGCACACAGGGATTAGCGCGTACATGTATTGTATTATGCCAACACCCCGGTGCTGACACGGAAGAAACCGGACGTTATGATTTAGCg                                                                   >  2:439596/1‑251 (MQ=255)
                                                                                                                                                                                               aaaCAATGAGTTGAAATTTCATATTGTTAAGATGTATTAATGTATGCCAGGTGCGA‑TGAATCGTCATTGTATTCCCGGATTAACTATGTCCACAGCCCTGACGGGGAACTTCTCTGCGGGAGTGTCCGGGAGTAATTAAAAACGATGCACACATGGTTTAGCGCGTACATGTATTTTATTATGCCAACACCCCGGTGCTGACGCGGAAGAAGCCGGACGTTATTATTAAGCGTGGAAAGATTTGTGTAgt                                                 <  2:389795/250‑1 (MQ=255)
                                                                                                                                                                                                   aaTGAGTTGAAATTTCATATTGTTAATATTTATTAATGTATGCCAGGTGCGA‑TGAATCGTCATTGTATTCCCGGATTAACTATGTCCACAGCCCTGACGGGGAACTTCTCTGCGGGAGTGTCCGGGAATAATTAAAAACGATGCACACAGGGTTTAGCGCGTACATGTATTGTATTATGCCAACACCCCGGTGCTGACACGGAAGAAACCGGACGTTATGATTTAGCGTGGAAAGATTTGTGTAGTGTTCt                                            <  1:89786/251‑1 (MQ=255)
                                                                                                                                                                                                            aaaTTTCATATTGTTAATATTTATTAATGTATGTCAGGTGCGA‑TGAAGCGTCATTGTATTCCCGGATTAACTATGTCCACAGCCCTGACGGGGAACGTCTCTGCGGGAGTGTCCGGGAATAATTAAAA‑CGATGCACACAGGGTTTGGCGCGTACTCGTATTGCATTATGCCAACGCCCCGGGTCTGACACGGAAGAAACCGGACGGTAGGATATAGCGGGGAAAGATTTGTGTAGTGTTCTGAATGCtctc                                  <  1:433113/251‑1 (MQ=255)
                                                                                                                                                                                                                     atTGTTAATATTTATTAATGTATGCCAGGTGCGA‑TGAATCGTCATTGTATTCCCGGATTAACTATGTCCACAGCCCTGACGGGGAACTTCTCTGCGGGAGTGTCCGGGAATAATTAAAAACGATGCACACAGGGTTTAGCGCGTACATGTATTGTATTATGCCAACACCCCGGTGCTGACACGGAAGAAACCGGACGTTATGATTTAGCGTGGAAAGATTTGTGTAGTGTTCTGAATGCTCTCAGTAAATa                          <  1:276833/251‑1 (MQ=255)
                                                                                                                                                                                                                      ttGTTAATATTTATTAATGTATGACAGGTGCGA‑TGAATCGTCGTTGCATTCCCGACTTAACTATGTCCACAGCCCTGACGGG‑AACTTCTCTGCGGGAGCGTCCGGGAATACTTAAAAA‑GATGCACCCAGGGTTTAGCGCGTACACGGATTGCATTATGCCAACGCCCCGGTGCTGCCACGGAAGAAACCGGGCGTTATGATTTAGCGTGGAAAGATTTGTGTAGTGTTCTGAATGCTCTCAGTAAATAGt                        >  1:409097/1‑250 (MQ=255)
                                                                                                                                                                                                                            atatTTATTAATGTATGCCAGGTGCGA‑TGAATCGTCATTGTATTCCCGGATTAACTATGTCCACAGCCCTGACGGGGAACTTCTCTGCGGGAGTGTCCGGGAATAATTAAAAACGATGCACACAGGGTTTAGCGCGTACATGTATTGTATTATGCCAACACCCCGGTGCTGACACGGAAGAAACCGGACGTTATGATTTAGCGTGGAAAGATTTGTGTAGTGTTCTGAATGCTCTCAGTAAATAGTAATg                    >  1:76041/1‑250 (MQ=255)
                                                                                                                                                                                                                              atttattAATGTATGTCAGGTGCGA‑TGAATCGTCATTGTATTCCCGGATTAACTATGTCCACAGCCCTGACGGGGAACCTCTCTGCGGGAGTGTCCGGGAATAATTAAAA‑CGATGCACACAGGGTTTAGCGCGTACACGTATTGCATTATGCCAACGCCCCGGTGCTGACACGGAAGAAACCGGACGTTATGATTTAGCGTGGAAAGATTTGTGTAGTGTTCTGAATGCTCTCAGTAAATAGTAATGAAtt                >  1:62750/1‑251 (MQ=255)
                                                                                                                                                                                                                                   ttaATGTATGTCAGGTGCGA‑TGAATCGTCATTGTCTTCCCGGATTAACTATGTCCACAGCCCTGACGGGGAACTTCTCTGCGGGAGTGTCCGGGAATAATTAAAA‑CGATGCACACAGGGTTTAGCGCGTACACGTATTGCATTATGCCAACGCCCCGGTGCTGACACGGAAGAAACCGGACGTTATGATTTAGCGTGGAAAGATTTGTGTAGTGTTCTGAATGCTCTCAGTAAATAGTa                       >  1:349836/1‑239 (MQ=255)
                                                                                                                                                                                                                                   ttaATGTATGTCAGGTGCGA‑TGAATCGTCATTGTATTCCCGGATTAACTATGTCCACAGCCCTGACGGGGAACTTCTCTGCGGGAGTGTCCGGGAATAATTAAAA‑CGATGCACACAGGGTTTAGCGCGTACACGTATTGCATTATGCCAACGCCCCGGTGCTGACACGGAAGAAACCGGACGTTATGATTTAGCGTGGAAAGATTTGTGTAGTGTTCTGAATGCTCTCAGTAAATAGTa                       <  2:349836/239‑1 (MQ=255)
                                                                                                                                                                                                                                     aaTGTATGTCAGGTGCGA‑TGAATCGTCATTGTATTCCCGGATTAACTATGTCCACAGCCCTGACGGGGAACTTCTCTGCGGGAGTGTCCGGGAATAATTAAAT‑CGATGCACACAGGGTTTAGCTCGTACACGTATTGCATTATGCCAACGCCCCGGTGCTGACACGGAAGAAACCGGACGTTATGATTTAGCGTGGAAAGATTTGTGTAGTGTTCTGAATGCTCTCAGTAAATAGTAATGAATTATCAAAg         <  2:387102/251‑1 (MQ=255)
                                                                                                                                                                                                                                      atgtatgtCAGGTGCGA‑TGAATCGTCATTGTATTCCCGGATTAACTATGTCCACAGCCCTGACGGGGAACTTCTCTGCGGGAGTGTCCGGGAATAATTAAAA‑CGATGCACACAGGGTTTAGCGCGTACACGTATTGCATTATGCCAACGCCCCGGTGCTGACACGGAAGAAACCGGACGTTATGATTTAGCGTGGAAAGATTTGTGTAGTGTTCTGAATGCTCTCAGTAAATAGTAATGAATTATCAAAgg        <  1:199568/251‑1 (MQ=255)
                                                                                                                                                                                                                                       tGTATGCCAGGTGCGA‑TGAATCGTCATTGTATTCCCGGATTAACTATGTCCACAGCCCTGACGGGGAACTTCTCTGCGGGAGTGTCCGGGAATAATTAAAAACGATGCACACAGGGTTTAGCGCGTACATGTATTGTATTATGCCAACACCCCGGTGCTGACACGGAAGAAACCGGACGTTATGATTTAGCGTGGAAAGATTTGTGTAGTGTTCTGAATGCTCTCAGTAAATAGTAATGAATTATCAAAgg        >  1:309549/1‑251 (MQ=255)
                                                                                                                                                                                                                                             ccAGGTGCGA‑TGAATCGTCATTGTATTCCCGGATTAACTATGTCCACAGCCCTGACGGGGAACTTCTCTGCGGGAGTGTCCGGGAATAATTAAAAACGATGCACACAGGGTTTAGCGCGTACATGTATTGTATTATGCCAACACCCCGGTGCTGACACGGAAGAAACCGGACGTTATGATTTAGCGTGGAAAGCTTTGTGTAGTGTTCTGAATGCTCTCAGTAAATAGTAATGAATTATCAAAGGTATAGt  >  2:91939/1‑251 (MQ=255)
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                           
CGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGTCTGTTTGTTTTGAACAGTAAATGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACAATGAGTTGAAATTTCATATTGTTAATATTTATTAATGTATGCCAGGTGCGA‑TGAATCGTCATTGTATTCCCGGATTAACTATGTCCACAGCCCTGACGGGGAACTTCTCTGCGGGAGTGTCCGGGAATAATTAAAAACGATGCACACAGGGTTTAGCGCGTACATGTATTGTATTATGCCAACACCCCGGTGCTGACACGGAAGAAACCGGACGTTATGATTTAGCGTGGAAAGATTTGTGTAGTGTTCTGAATGCTCTCAGTAAATAGTAATGAATTATCAAAGGTATAGT  >  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: ‑