Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,383,7220GA26.1% 87.4 / 15.9 46F80F (TTC→TTT) ycjWDNA‑binding transcriptional repressor YcjW
Reads supporting (aligned to +/- strand):  ref base G (30/4);  new base A (0/12);  total (30/16)
Fisher's exact test for biased strand distribution p-value = 4.68e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CGACAGAATAAACATGGGCATATTGACCTTCAACTTTGCCAATCACCACCACCGGAATATCATATTTATCGAGTTGGGCAAAAAATGACTCATCCGCTGGCGAACTGAGCATAATAATGCCTTTAAT‑CATTTTCTGCTTAATTTTGCTTTCGCATTTTTGCAGATCTTCTGCCGGATTGTG‑CGACGTTTGTAA‑TATCACGTCGAAACCTTCTTCTTCAGCTTTGGCGGTG‑ATGGCATGCAAAACTTC‑AGAGAAAAAC‑GGATTACCCGCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCATCAAAACCCGAAGAGGTCAGCGCGCGGGCCAGCTTGTTTGGCTGATACTGTAATTCTTCAATGGCCCGTAACACTTTTTCGCGCGCTTCCGGGGAGATATTAGTTTGCTTATTCAGCACGCGTGATACGGTGGATTTTGATACGCCTGCAACCCTGGCAATATCATAAATAGTAGGGGACATAGGTCAGGA  >  NC_000913/1383472‑1383970
                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                         
cGACAGAATAAACATGGGCATATTGACCTTCAACTTTGCCAATCACCACCACCGGAATATCATATTTATCGAGTTGGGCAAAAAATGACTCATCCGCTGGCGAACTGAGCATAATAATGCCTTTAAT‑CATTTTCTGCTTAATTTTGCTTTCGCATTTTTGCAGATCTTCTGCCGGATTGTG‑CGACGTTTGTAA‑TATCACGTCGAAACCGTCTTCTTCAGCTTTGGCGGTG‑ATGGCATGCAAAACTTC‑agag                                                                                                                                                                                                                                                           <  2:137493/251‑1 (MQ=255)
     gAATAAACATGGGCATATTGACCTTCAACTTTGCCAATCACCACCACCGGAATATCATATTTATCGAGTTGGGCAAAAAATGACTCATCCGCTGGCGAACTGAGCATAATAATGCCTTTAAT‑CATTTTCTGCTTAATTTTGCTTTCGCATTTTTGCAGATCTTCTGCCGGATTGTG‑CGACGTTTGTAA‑TATCACGTCGAAACCTTCTTCTTCAGCTTTGGCGGTG‑ATGGCATGCAAAACTTC‑AGAGaaaaa                                                                                                                                                                                                                                                      <  1:213796/251‑1 (MQ=255)
          aaCATGGGCATATTGAACTTCAACTTTGCCAATCACCCCCACCGGCATCTCATATTTATCGAGTTGGGCAAAAAATGACTCATCCGCTGGCGAACTGAGCATAATAATGCCTTTAAT‑CATTTTCTGCTTAATTTTGCTTTCGCATTTTTGCAGATCTTCTGCCGGATCGTG‑CGACGTTTGTCA‑TATCACGTCGAAACCTTCTTCTTCAGCTTTGGCGGTG‑ATGGCATGCAAAACTTC‑AGAGAAAAAC‑GGAt                                                                                                                                                                                                                                                >  1:341841/1‑251 (MQ=255)
                       tGACCTTCAACTTTGCCAATCACCACCACCGGAATATCATATTTATCGAGTTGGGCAAAAAATGACTCATCCGCTGGCGAACTGAGCATAATAATGCCTTTAAT‑CATTTTCTGCTTAATTTTGCTTTCGCATTTTTGCAGATCTTCTGCCGGATTGTG‑CGACGTTTGTAA‑TATCACGTCGAAACCTTCTTCTTCAGCTTTGGCGGTG‑ATGGCATGCAAAACTTC‑AGAGAAAAAC‑GGATTACCCGCCGTAGt                                                                                                                                                                                                                                   >  2:116185/1‑251 (MQ=255)
                                 cTTTGCCAATCACCAACCCCGGAATATCATATTTATCGCGTTGGGCAAAAAATGACTCATCCGCTGGCGAACTGCGCATAATAATGCCTTTAAT‑CATTTTCTGCTTAATTTTGCTTTCGCATTTTTGCAGATCTTCTGCCGGATTGTG‑CGACGTTTGTAC‑TATCACGTCGAAACCTTCTTCTTCAGCTTTGGCGGTG‑ATGGCATGCCAAACTTC‑AGAGAAACAC‑GGATTACCCGCCGTAGTTTTGGTCGa                                                                                                                                                                                                                          >  2:454421/1‑250 (MQ=255)
                                                          aaaaatttttgcgattgggaaaaaaaaTGACTCCACCGCTGGAGAATAGAAAAAAATAATGCTTTTAAT‑CATTTTCTGTTTATTTTTTTTTTCTCATTTTTGCAGATCTTCTGC‑GGATTGTGGCGCAGTTTGTAA‑AACCACGACAAACCCTTTTTCGCCGGTTTTGGCGGGG‑AGGGCATGAAAAATTGA‑AAAAAAAAAG‑GGATACCCCGCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCAt                                                                                                                                                                                                <  2:349680/230‑1 (MQ=255)
                                                            cATATTTATCGAGTTGGGCAAAAAATGACTCATCCGCTGGCGAACTGAGCCTAATAATGCCTTTAAT‑CATTTTCTGCTTAATTTTGCTTTCGCATTTTTGCAGATCTTCTGCCGGATTGTG‑CGACGTTTGTAA‑TATCACGTCGAAACCTTCTTCTTCAGCTTTGGCGGTG‑ATGGCATGCAAAACTTC‑AGAGAAAAAC‑GGATTACCCGCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCATc                                                                                                                                                                                               >  2:617095/1‑250 (MQ=255)
                                                               atTTATCGAGTTGGGCAAAAAATGACTCATCCGCTGGCGAACTGAGCCTAATAATGCCTTTAAT‑CATTTTCTGCTTAATTTTGCTTTCGCATTTTTGCAGATCTTCTGCCGGATTGTG‑CGACGTTTGTAA‑TATCACGTCGAAACCTTCTTCTTCAGCTTTGGCGGTG‑ATGTCATGCAAAACTTC‑AGAGAAAAAC‑GGATTACCCGCCGTCGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCATCaaa                                                                                                                                                                                            >  1:665128/1‑250 (MQ=255)
                                                                          tGGGCAAAAAATGACTCATCCGCTGGCGAACTGAGCATAATAATGCCTTTAAT‑CATTTTCTGCTTAATTTTGCTTTCGCATTTTTGCAGATCTTCTGCCGGATTGTG‑CGACGTTTGTAA‑TATCACGTCGAAACCTTCTTCTTCAGCTTTGGCGGTG‑ATGGCATGCAAAACTTC‑AGAGAAAAAC‑GGATTACCCGCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCATCAAAACCCGAAGAGGt                                                                                                                                                                                >  1:166121/1‑251 (MQ=255)
                                                                               aaaaaaTGACTCATCCGCTGGCGCAATGCGCATCATAATGCCTTTAAT‑CATTTTCTGCTTAATTTTGCTTTCGCATTTTTGCAGATCTTCTGCCGGATTGTG‑CGACGTTTGTAA‑TATCACGTCGAAACATTCTTCTTCAGCTTTGGCGGTG‑ATGGCATGCAAAACTTC‑AGAGAAAAAC‑GGCTTACCCGACGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCATCAAAACCCGAAGAGGTCAgcg                                                                                                                                                                           >  2:356942/1‑251 (MQ=255)
                                                                               aaaaaaTGACTCATCCGCTGGCGAACTGAGCATAATAATGCCTTTAAT‑CATTTTCTGCTTAATTTTGCTTTCGCATTTTTGCAGATCTTCTGCCGGATTGTG‑CGACGTTTGTAA‑TATCACGTCGAAACCTTCTTCTTCAGCTTTGGCGGTG‑ATGGCATGCAAAACTTC‑AGAGAAAAAC‑GGATTACCCGCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCATCAAAACCCGAAGAGGTCAgcg                                                                                                                                                                           >  2:383900/1‑251 (MQ=255)
                                                                                 aaaaTGACTCACCCGCGGGCGAAGTGAGCATAATAAGTCATTAAAA‑TATTTTCTGCTAATTTTTGCTTTCGCTTTTTTGGAGATCTTCTGCGGGATGGTG‑CGACGTTTGTAA‑TACCACGTCGAACCCTTCTTCGTCAGCTTTGGGGGTG‑AGGGCAGGAAAAACTTG‑AGAAAAAACC‑GGTTTACCCGCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCATCAAAACCCGAAGAGGTCAgcgcg                                                                                                                                                                         <  1:341673/251‑1 (MQ=255)
                                                                                    aTGACTCATCCGCTGGCGAACTGAGCATAATAATGCCTTTAAT‑CATTTTCTGCTTAATTTTGCTTTCGCATTTTTGCAGATCTTCTGCCGGATTGTG‑CGACGTTTGTAA‑TATCACGTCGAAACCTTCTTCTTCAGCTTTGGCGGTG‑ATGGCATGCAAAACTTC‑AGAGAAAAAC‑GGATTACCCGCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCATCAAAACCCGAAGAGGTCAGCGCGCgg                                                                                                                                                                      >  2:215526/1‑251 (MQ=255)
                                                                                      gACTCATCCGCTGGCGACCTGAGAATAATATTGC‑TTTAATTCATTTTCTGTTTAATTTTGCTTCCGCATTTTTGCAGATCTTCTGCCGGTTTGGG‑CGACTTTTGTAA‑TATCACGTCAAAACCTTCTTCTTCACCTTTGGCGGTG‑ATGGCATGCAAAATTTA‑AGAAAAAAAC‑GGATTCCCCGCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCATCAAAACCCGAAGAGGTCAGCGCGCGGGc                                                                                                                                                                    <  1:603079/251‑1 (MQ=255)
                                                                                       cacacccctcttccgatcggagaaaaatattGCTTTTAAC‑CTTTTTCTGCTTAATTTTGCTTTCGCTTTTGTGCAGATCTTCTGCGGGTTTGGG‑CGCCGTTTGAAA‑AACCCCGACAAACCCTTCTTCTCCAGCTTGGGGGGGA‑AGGGCA‑GAAAAACCTA‑AAAAAAAAACCGTATTCCCCGCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCATCAAAACCCGAAGAGGTCAGCGCGCGGGcc                                                                                                                                                                   <  2:339029/226‑1 (MQ=255)
                                                                                          cccccGCTGGCGAACTGAGTATAATAATGCCTTTAAT‑CATTTTGTGTTTAATTTTGTTTTCTCTTTTTTGCAGATCTTCTGCGGGATGGTG‑CGATTTTGGTAA‑TATCACGTCGAAACTTCCTCCTTCAG‑TTGTGCGGTGGTTGGAAGACAAAACTAC‑AGAAAAAATC‑AGATCCCCCGCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCATCAAAACCCGAAGAGGTCAGCGCGCGGGCCa                                                                                                                                                                  <  1:278264/246‑1 (MQ=255)
                                                                                           ccccGCTGGCGAACTGAGAATAATAATGCCTTTAAT‑TTTTTTCTTTTTATTTTTGTTTTCGAATTTTTGCAGATCTTCTGCGGGATTTTG‑CGACTTTTGAAA‑TATCACGTAAAACCCTTCTTCTTCAGTTTGGGGGGTG‑AGGGCATGAAAAACTAC‑AGAAAAAAAC‑GTATTACCCGCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCATCAAAACCCGAAGAGGTCAGCGCGCGGGCCAGCt                                                                                                                                                               <  1:84409/249‑1 (MQ=255)
                                                                                                  ggcaaaCTGAGAATAATAATGCCTTTAAT‑CATTTTCTGCTAATTTTTGCTTTCGCATTTTTGCAGATTTTCGGCGGGATTGTG‑CGATGTTTGAAA‑TATCCCTTCGAACCCTTCTTCTCCAGCTTTGGCGGGG‑ATGGCATGAAAAACTAA‑AGAGAAAAAC‑GGATTCCCCGCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCATCAAAACCCGAAGAGGTCAGCGCGCGGGCCAGCTTGTTTgg                                                                                                                                                        <  1:215526/247‑1 (MQ=255)
                                                                                                       acggatcATAATAATGCCTTTAAT‑CATTTTCTGCTTAATTTTGCTTTCGCATTTTTGCAGATCTTCGGCCGGATTGTG‑CGACGTTTGTAA‑TACCACGCCGAAACCTTCTTCTTCAGCTTTGGCGGTG‑ATGGCATGAAAAACTTC‑AGAGAAAAAC‑GAATCACCCGCCGTAGTTTTg                                                                                                                                                                                                                               <  1:171744/169‑1 (MQ=255)
                                                                                                       aCTGAGCATAATAATGCCTTTAAT‑CATTTTCTGCTTAATTTTGCTTTCGCATTTTTGCAGCTCTTCTGCCGGATTGTG‑CGACGTTTGTAA‑TATCACGTCGAAACCTTCTTCTTCAGCTTTGGCGGTG‑ATGGCATGCAAAACTTC‑AGAGAAAAAC‑GGATTACCCGCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCATCAAAACCCGAAGAGGTCAGCGCGCGGGCCAGCTTGTTTGGCTGAt                                                                                                                                                   >  2:677764/1‑251 (MQ=255)
                                                                                                       aCTGAGCATAATAATGCCTTTAAT‑CATTTTCTGCTTAATTTTGCTTTCGCATTTTTGCAGATCTTCTGCCGGATTGTG‑CGACGTTTGTAA‑TATCACGTCGAAACCTTCTTCTTCAGCTTTGGCGGTG‑ATGGCATGCAAAACTTC‑AGAGAAAAAC‑GGATTACCCGCCGTAGTTTTg                                                                                                                                                                                                                               >  2:171744/1‑175 (MQ=255)
                                                                                                       aCTGAGCATAATAATGCCTTTAAT‑CATTTTCTGCTTAATTTTGCTTTCGCATTTTTGCAGATCTTCTGCCGGATTGTG‑CGACGTTTGTAA‑TATCACGTCGAAACCTTCTTCTTCAGCTTTGGCGGTG‑ATGGCATGCAAAACTTC‑AGAGAAAAAC‑GGATTACCCGCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCATCAAAACCCGAAGAGGTCAGCGCGCGGGCCAGCTTGTTTGGCTGAt                                                                                                                                                   >  1:69129/1‑251 (MQ=255)
                                                                                                          gAGCATAATAATGCCTTTAAT‑CATTTTCTGCTTAATTTTGCTTTCGCATTTTTGCAGATCTTCTGCCGGATTGTG‑CGACGTTTGTAA‑TATCACGTCGAAACCTTCTTCTTCAGCTTTGGCGGTG‑ATGGCATGCAAAACTTC‑AGAGAAAAAC‑GGATTACCCGCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCATCAAAACCCGCAGAGGTCAGCGCGCGGGCCAGCTTGTTTGGCTGATACt                                                                                                                                                >  2:659692/1‑251 (MQ=255)
                                                                                                          gAGCATAATAATGCCTTTAAT‑CATTTTCTGCTTAATTTTGCTTTCGCATTTTTGCAGATCTTCTGCCGGATTGTG‑CGACGTTTGTAA‑TATCACGTCGAAACCTTCTTCTTCAGCTTTGGCGGTG‑ATGGCATGCAAAACTTC‑AGAGAAAAAC‑GGATTACCCGCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCATCAAAACCCGAAGAGGTCAGCGCGCGGGcc                                                                                                                                                                   >  1:339029/1‑232 (MQ=255)
                                                                                                                aataatGCCTTTAAT‑CATTTTCTGCTTAATTTTGCTTTCGCATTTTTGCAGATCTTCTGCCGGATTGTG‑CGACGTTTGTAA‑TATCACGTCGAAACCTTCTTCTTCAGCTTTGGCGGTG‑ATGGCATGCAAAACTTC‑AGAGAAAAAC‑GGATTACCCGCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCATCAAAACCCGAAGAGGTCAGCGCGCGGGCCAGCTTGTTTGGCTGATACTGTAAtt                                                                                                                                          >  2:475779/1‑251 (MQ=255)
                                                                                                                            aaT‑CATTTTCTGCTTAATTTTGCTTTCGCATTTTTGCAGATCTTCTGCCGGATTGTG‑CGACGTTTGTAA‑TATCACGTCGAAACCTTCTTCTTCAGCTTTGGCGGTG‑ATGGCATGCAAAACTTC‑AGAGAAAAAC‑GGATTACCCGCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCATCAAAACCCGAAGAGGTCAGCGCGCGGGCCAGCTTGTTTGGCTGATACTGTAATTCTTCAATGGccc                                                                                                                              >  2:574818/1‑251 (MQ=255)
                                                                                                                             aT‑CATTTTCTGCTTAATTTTGCTTTAGCATTTTTGCAGATCTTCTGCCGGATTGTG‑CGACGTTTGTAA‑TATCACGTCGAAACCTTCTTCTTCAGCTTTGGCGGTG‑ATGGCATGCAAAACTTC‑AGAGAAAAAC‑GGATTACCCGCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCATCAAAACCCGCAGAGGGCCGCGCGCGGGCCAGCTTGTTTGGCTGATACTGTAATTTTTCAATGGCCCg                                                                                                                             >  2:397561/1‑251 (MQ=255)
                                                                                                                                        gCTTAATTTTGCTTTCGCATTTTTGCAGATCTTCTGCCGGATTGTG‑CGACGTTTGTAA‑TATCACGTCGAAACCTTCTTCTTCAGCTTTGGCGGTG‑ATGGCATGCAAAACTTC‑AGAGAAAAAC‑GGATTACCCGCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCATCAAAACCCGAAGAGGTCAGCGCGCGGGCCAGCTTGTTTGGCTGATACTGTAATTCTTCAATGGCCCGTAACACtttt                                                                                                                   >  1:391036/1‑251 (MQ=255)
                                                                                                                                          ttAATTTTGCTTTCGCATTTTTGCAGATCTTCTGCCGGATTGTG‑CGACGTTTGTAA‑TATCACGTCGAAACCTTCTTCTTCAGCTTTGGCGGTG‑ATGGCATGCAAAACTTC‑AGAGAAAAAC‑GGAATACCCGCCGTAGTTTTGGTCGGACGGGTAGAAATCACCATAATGGCATCAAAACCCGAAGAGGTCAGCGCGCGGGCCAGCTTGTTTGGCTGATACTGGAATTCTTCAATGGCCCGTAACACtttttt                                                                                                                 >  2:472617/1‑250 (MQ=255)
                                                                                                                                            atttttGTTTTCGCTTTTTGGCAGATTTGCTGCGGGTTTGGG‑CGCCGTTTGAAA‑TACCACGTCAAACCCTTCTTCTTCAGTTTGGGCGTTG‑ATGGCATGAAAAACTAA‑AGAAAAAAAC‑GGATTCCCCGCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCATCAAAACCCGAAGAGGTCAGCGCGCGGGCCAGCTTGTTTGGCTGATACTGTAATTCTTCAATGGCCCGTAACACTTTTTcg                                                                                                                <  1:383900/248‑1 (MQ=255)
                                                                                                                                            aaTTTTGCTTTCGCATTTTTGCAGATCTTCTGCCGGATTGTG‑CGACGTTTGTAA‑TATCACGTCGAAACCTTCTTCTTCAGCTTTGGCGGTG‑ATGGCATGCAAAACTTC‑AGAGAAAAAC‑GGATTACCCGCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCATCAAAACCCGAAGAGGTCAGCGCGCGGGCCAGCTTGTTTGGCTGATACTGTAATTCTTCAATGGCCCGTAACACTTTTTcgc                                                                                                               >  1:684233/1‑251 (MQ=255)
                                                                                                                                               tttGCTTTCGCATTTTTGCAGATCTTCTGCCGGATTGTG‑CGACGTTTGTAA‑TATCACGTCGAAACCTTCTTCTTCAGCTTTGGCGGTG‑ATGGCATGCAAAACTTC‑AGAGAAAAAC‑GGATTACCCGCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCATCAAAACCCGAAGAGGTTAGCGCGCGGGCCAGCTTGTTTGGCTGATACTGTAATTCTTCAATGGCCCGTAACACTTTTTcgcgcg                                                                                                            >  2:472663/1‑251 (MQ=255)
                                                                                                                                                     ttCGCATTTTTGCAGATCTTCTGCCGGATTGTG‑CGACGTTTGTAA‑TATCACGTCGAAACCTTCTTCTTCAGCTTTGGCGGTG‑ATGGCATGCAAAACTTC‑AGAGAAAAAC‑GGATTACCCGCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCATCAAAACCCGAag                                                                                                                                                                                    >  1:216948/1‑173 (MQ=255)
                                                                                                                                                            ttttGCAGATCTTCTGCCGGATTGTG‑CGACGTTTGTAA‑TATCACGTCGAAACCTTCTTCTTCAGCTTTGGCGGTG‑ATGGCATGCAAAACTTC‑AGAGAAAAAC‑GGATTACCCGCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCATCAAAACCCGAAGAGGTCAGCGCGCGGGCCAGCTTGTTTGGCTGATACTGTAATTCTTCAATGGCCCGTAACACTTTTTCGCGCGCTTCCgggg                                                                                                   >  1:579419/1‑247 (MQ=255)
                                                                                                                                                                  aGATCTTCTGCCGGATTGTG‑CGACGTTTGTAA‑TATCACGTCGAAACCTTCTTCTTCAGCTTTGGCGGTG‑ATGGCATGCAAAACTTC‑AGAGAAAAAC‑GGATTACCCGCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCATCAAAACCCGAAGAGGTCAGCGCGCGGGCCAGCTTGTTTGGCTGATACTGTAATTCTTCAATGGCCCGTAACACTTTTTCGCGCGCTTCCGGGGAGATATTAGt                                                                                         >  1:164085/1‑251 (MQ=255)
                                                                                                                                                                      cttctGCCGGATTGTG‑CGACGTTTGTAA‑TATCACGTCGAAACCTTCTTCTTCAGCTTTGGCGGTG‑ATGGCATGCAAAACTTC‑AGAGAAAAAC‑GGATTACCCGCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCATCAAAACCCGAAGAGGGCAGCGCGCGGGCCAGCTTGTTTGGCTGATACTGTAATTCTTCAATGGCCCGTAACACTTTTTCGCGCGCTTCCGGGGTGATATTAGTTTGc                                                                                     >  2:509312/1‑251 (MQ=255)
                                                                                                                                                                       ttctGCCGGATTGTG‑CGACGTTTGTAA‑TATCACGTCGAAACCTTCTTCTTCAGCTTTGGCGGTG‑ATGGCATGCAAAACTTC‑AGAGAAAAAC‑GGATTACCCGCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCATCAAAACCCGAAGAGGTCAGCGCGCGGGCCAGCTTGTTTGGCTGATACTGTAATTCTTCAATGGCCCGTAACACTTTTTCGCGCGCTTCCGGGGAGATATTAGTTTGCt                                                                                    >  1:675350/1‑251 (MQ=255)
                                                                                                                                                                        tcggCCGGATTGTG‑CGTCGTT‑GTAAATACCAGCTCGAACCCTTCTACTTCGGTTTTGGCGGTG‑ATGGCATGAAAAACTTC‑AGAAAAAAGC‑GGATTGCCCGCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCATCAAAACCCGAAGAGGTCAGCGCGCGGGCCAGCTTGTTTGGCTGATACTGTAATTCTTCAATGGCCCGTAACACTTTTTCGCGCGCTTCCGGGGAGATATTAGttt                                                                                       <  2:382705/244‑1 (MQ=255)
                                                                                                                                                                            ccGGATTGTG‑CGACGTTTGTAA‑TATCACGTCGAAACCTTATTCTTCAGCTTTTGCGGTG‑ATGGCATGCAAACCTTC‑AGAGAAAAAC‑GGCTTACCCGCCGTAGTTTTGGTCGAAGGGGTAGAAATCACCATAATGGCATCAAAACCCGAAGAGGTCAGCGCGCGGGCCAGCTTGTTTGGCTGATACTGTAATTATTCAATGGCCCGTAACACTTTTTCGCGCGCTTCCGGGGAGATATTAGTTTGATTATTc                                                                               >  1:364815/1‑251 (MQ=255)
                                                                                                                                                                             cGGATTGTG‑CGACGTTTTTAA‑TACCACGTCAAACCCTTCTTCTTCAGCTTGGGGGGTG‑ATGGCATCCAAAACTTC‑AGAAAAAAAG‑GGATTACCCGCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCATCAAAACCCGAAGAGGTCAGCGCGCGGGCCAGCTTGTTTGGCTGATACTGTAAtt                                                                                                                                          <  1:491706/191‑1 (MQ=255)
                                                                                                                                                                             cGGATTGTG‑CGACGTTTGTAA‑TATCACGTCGAAACCTTCTTCTTCAGCTTTGGCGGTG‑ATGGCATGCAAAACTTC‑AGAGAAAAAC‑GGATTACCCGCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCATCAAAACCCGAAGAGGTCAGCGCGCGGGCCAGCTTGTTTGGCTGATACTGTAAtt                                                                                                                                          >  2:491706/1‑191 (MQ=255)
                                                                                                                                                                              ggtttttgccGACTTTTGTAA‑TACCACTACGAACCCTTTTTCTCCATCTTTGGCGGGT‑ATGGCAT‑AAAAACTTCCAAAAAAAAAC‑GTATCCCCCGCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCATCAAAACCCGAAGAGGTCAGCGCGCGGGCCAGCTTGTTTGGCTGATACTGTAATTCTTCAATGGCCCGTAACACTTTTTCGCGCGCTTCCGGTGAGATATTAGTTTGCTTATTCa                                                                              <  1:599202/242‑1 (MQ=255)
                                                                                                                                                                                         gcgttttgtaatacccctccgataCCTTCTTTTCCAGCTTTGGCGGTG‑ATGGCTTGAAAAACTTA‑AAAGAAAACC‑GAATTACCCGCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCATCAAAACCCGAAGAGGTCAGCGCTCGGGCCAGCTTGTTTGGCTGATACTGTAATTCTTCAATGGCCCGTAACACTTTTTCGCGCTCTTCCGGGGAGATATTAGTTTGCTTATTCAGCACGCGTGa                                                                    <  1:356942/228‑1 (MQ=255)
                                                                                                                                                                                                     atCACGTCGAAACCTTCTTCTTCAGCTTTGGCGGTG‑ATGGCATGCAAAACTTC‑AGAGAAAAAC‑GGATTACCCGCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCATCAAAACCCGAAGAGGTAAGCGCGCGGGCCAGCTTGTTTGGCTGATACTGTAATTCTTCAATGGCCCGTAACACTTTTTCGCGCGCTTCCGGGGAGATATTAGTTTGCTTATTCAGCACGCGGGATACGGGGGAttt                                                        >  2:32117/1‑251 (MQ=255)
                                                                                                                                                                                                      ccACGTAGAACCCTTCTTCTTCAGCTTGGGCGGGG‑ATGGCAGAAAAAACTAC‑AAAAAAAAAC‑GGTTTCCCCGCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCATCAAAACCCGAAGAGGTCAGCGCGCGGGCCAGCTTGTTTGGCTGATACTGTAATTCTTCAATGGCCCGTAACACTTTTTCGCGCGCTTCCGGGGAGATATTAGTTTGCTTATTCAGCACGCGTGATACGGTGGAtttt                                                       <  1:570914/250‑1 (MQ=255)
                                                                                                                                                                                                              aacccTTCTTCTTCAGCTTTGGCGGTG‑ATGGATTGCAAAACTTC‑AGAAAAAACC‑GGATTCCCCGCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCATCAAAACCCGAAGAGGTCAGCGCGCGGGCCAGCTTGTTTGGCTGATACTGTAATTCTTCAATGGCCCGTAACACTTTTTCGCGCGCTTCCGGGGAGATATTAGTTTGCTTATTCAGCACGCGTGATACGGTGGATTTTGATACGcc                                               <  1:20115/248‑1 (MQ=255)
                                                                                                                                                                                                                                     ggTG‑ATGGCATGCAAAACTTC‑AGAGAAAACC‑GGATTACCCGCCGCAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCATCAAAACCCGAAGAGGTCAGCGCGCGGGCCAGCTTGTTTGGCTGATACTGTAATTCTTCAATGGCCCGTAACACTTTTTCGCGCGCTTCCGGGGAGATATTAGTTTGCTTATTCAGCACGCGTGATACGGTGGATTTTGATACGCCTGCAACCCTGGCAATATCATaa                         >  2:642150/1‑250 (MQ=255)
                                                                                                                                                                                                                                                              agaAAAAC‑GGATTACCCGCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCATCAAAACCCGTAGAGGTCAGCGCGCGGGCCAGCTTGTTTGGCTGATACTGTAATTCTTCAATGGCCCGTAACACTTTTTCGCGCGCTTCCGGGGAGATATTAGTTTGCTTATTCAGCACGCGTGATACGGTGGATTTTGATACGCCTGCAACCCTGGCAATATCATAAATAGTAGGGGACATAGGTCAGGa  >  2:47670/1‑250 (MQ=255)
                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                         
CGACAGAATAAACATGGGCATATTGACCTTCAACTTTGCCAATCACCACCACCGGAATATCATATTTATCGAGTTGGGCAAAAAATGACTCATCCGCTGGCGAACTGAGCATAATAATGCCTTTAAT‑CATTTTCTGCTTAATTTTGCTTTCGCATTTTTGCAGATCTTCTGCCGGATTGTG‑CGACGTTTGTAA‑TATCACGTCGAAACCTTCTTCTTCAGCTTTGGCGGTG‑ATGGCATGCAAAACTTC‑AGAGAAAAAC‑GGATTACCCGCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCATCAAAACCCGAAGAGGTCAGCGCGCGGGCCAGCTTGTTTGGCTGATACTGTAATTCTTCAATGGCCCGTAACACTTTTTCGCGCGCTTCCGGGGAGATATTAGTTTGCTTATTCAGCACGCGTGATACGGTGGATTTTGATACGCCTGCAACCCTGGCAATATCATAAATAGTAGGGGACATAGGTCAGGA  >  NC_000913/1383472‑1383970

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