Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,585,0110GA23.8% 91.0 / 10.1 43E85K (GAA→AAA) yrhBputative heat shock chaperone
Reads supporting (aligned to +/- strand):  ref base G (27/5);  new base A (0/10);  total (27/16)
Fisher's exact test for biased strand distribution p-value = 2.04e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GATTACTTATCACGACGCATTCGCGAAAGCGAACCATTACCTTGATGATGCAGATCTCCCGGTCGTCATTACTCTACATGGACGCTTTAGCCAGGGCTGGTATTTCTGTTTCGAAGCACGAGAATTTCTCGAAACTGGAGATGAGGCCGCGCGCTTAGCTGGTAACGCACCTTTTATTATTGATAAAGACAGTGGTGAAATTCATTCTCTGGGAAC‑GGCAAAACC‑G‑CTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGC‑TACCTTCGGCTTGCCCT‑GACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTT‑AGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTGAGCGTATCCGGGCTAAACCCTTTTTCGACACGCAGCTCG  >  NC_000913/3584761‑3585261
                                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                             
gATTACTTCTCACGACGCATTCGCGAAAGCGAACCATTACCTTGATGATGCAGATCTCCCGGTCGTCATTACTCTACATGGACGCTTTAGCCAGGGCTGGTATTTCTGTTTCGAAGCACGAGAATTTCTCGAAACTGGAGATGAGGCCGCGCGCTTAGCTGGTAACGCACCTTTTATTATTGATAAAGACAGTGGTGAAATTCATTCTCTGGGAAC‑GGCAAAACC‑G‑CTGGAAGAATATCTACAGGATTACg                                                                                                                                                                                                                                                               >  1:562592/1‑251 (MQ=255)
       tATCACGACGCATTCGCGAAAGCGAACCATTACCTTGATGATGCAGATCTCCCGGTCGTCATTACTCTACATGGACGCTTTAGCCAGGGCTGGTATTTCTGTTTCGAAGCACGAGAATTTCTCGAAACTGGAGATGAGGCCGCGCGCTTAGCTGGTAACGCACCTTTTATTATTGATAAAGACAGTGGTGAAATTCATTCTCTGGGAAC‑GGCAAAACC‑G‑CTGGAAGAATATCTACAGGATTACGAAATaa                                                                                                                                                                                                                                                         >  1:370078/1‑250 (MQ=255)
        aTCACGACGCATTCGCGAAAGCGAACCATTACCTTGATGATGCAGATCTCCCGGTCGTCATTACTCTACATGGACGCTTTAGCCAGGGCTGGTATTTCTGTTTCGAAGCACGAGAATTTCTCGAAACTGGAGATGAGGCCGCGCGCTTAGCTGGTAACGCACCTTTTATTATTGATAAAGACAGTGGTGAAATTCATTCTCTGGGAAC‑GGCAAAACC‑G‑CTGGAAGAATATCTACAGGATTACGAAATaaa                                                                                                                                                                                                                                                        <  2:371800/250‑1 (MQ=255)
          cACGACGAATTCGCGAAAGCGAACCATTCCCTTGATGATGCAGATCTACCGGTCGTCATTACTCTACATGGACGCTTTAGCCAGGGCTGGTATTTCTGTTTCGAAGCACGAGAATTTCTCGAAACTGGAGATGAGGCCGCGCGCTTAGCTGGTAACGCACCTTTTATTATTGATAAAGACAGTGGTGAAATTCATTCTCTGGGAAC‑GGCAAAACA‑G‑CTGGAAGAATATCTACAGGATTACGAAATaaaaa                                                                                                                                                                                                                                                      >  2:296933/1‑250 (MQ=255)
            cgacgCATTCGCGAAAGCGAACCATTACCTTGATGATGCAGATCTCCCGGTCGTCATTACTCTACATGGACGCTTTAGCCAGGGCTGGTATTTCTGTTTCGAAGCACGAGAATTTCTCGAAACTGGAGATGAGGCCGCGCGCTTAGCTGGTAACGCACCTTTTATTATTGATAAAGACAGTGGTGAAATTCATTCTCTGGGAAC‑GGCAAAACC‑G‑CTGGAAGAATATCTACAGGATTACGAAATAAAAAAgg                                                                                                                                                                                                                                                   >  1:54124/1‑251 (MQ=255)
               cgCATTCGCGAAAGCGAACCATTACCTTGATGATGCAGATCTCCCGGTCGTCATTACTCTACATGGACGCTTTAGCCAGGGCTGGTATTTCTGTTTCGAAGCACGAGAATTTCTCGAAACTGGAGATGAGGCCGCGCGCTTAGCTGGTAACGCACCTTTTATTATTGATAAAGACAGTGGTGAAATTCATTCTCTGGGAAC‑GGCAAAACC‑G‑CTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGC‑t                                                                                                                                                                                                                                                >  2:494134/1‑250 (MQ=255)
                gCATTCGCGAAAGCGAACCATTACCTTGATGATGCAGATCTCCCGGTCGTCATTACTCTACATGGACGCTTTAGCCAGGGCTGGTATTTCTGTTTCGAAGCACGAGAATTTCTCGAAACTGGAGATGAGGCCGCGCGCTTAGCTGGTAACGCACCTTTTATTATTGATAAAGACAGTGGTGAAATTCATTCTCTGGGAAC‑GGCAAAACC‑G‑CTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGC‑TAc                                                                                                                                                                                                                                              >  2:240672/1‑251 (MQ=255)
                   ttCGCGAAAGCGAACCATTACCTTGATGATGCAGAGCTCCCGGTCGTCATTACTCTACATGGACGCTTTAGCCAGGGCTGGTATTTCTGTTTCGAAGCACGAGAATTTCTCGAAACTGGAGATGAGGCCGCGCGCTTAGCTGGTAACGCACCTTTTATTATTGATAAAGACAGTGGTGAAATTCATTCTCTGGGAAC‑GGCAAAACC‑G‑CTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGC‑TACCtt                                                                                                                                                                                                                                           >  2:227289/1‑251 (MQ=255)
                         aaaGCGAACCATTACCTTGATGATGCAGATCTCCCGGTCGTCATTACTCTACATGGACGCTTTAGCCAGGGCTGGTATTTCTGTTTCGAAGCACGAGAATTTCTCGAAACTGGAGATGAGGCCGCGCGCTTAGCTGGTAACGCACCTTTTATTATTGATAAAGACAGTGGTGAAATTCATTCTCTGGGAAC‑GGCAAAACC‑G‑CTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGC‑TACCTTCGGCtt                                                                                                                                                                                                                                     >  1:544521/1‑251 (MQ=255)
                         aaaGCGAACCATTACCTTGATGATGCAGATCTCCCGGTCGTCATTACTCTACATGGACGCTTTAGCCAGGGCTGGTATTTCTGTTTCGAAGCACGAGAATTTCTCGAAACTGGAGATGAGGCCGCGCGCTTAGCTGGTAACGCACCTTTTATTATTGATAAAGACAGTGGTGAAATTCATTCTCTGGGAAC‑GGCAAAACC‑G‑CTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGC‑TACCTTCGGCtt                                                                                                                                                                                                                                     <  2:33893/251‑1 (MQ=255)
                               aaCCATTACCTTGATGATGCAGATCTCCCGGTCGTCATTACTCTACATGGACGCTTTAGCCAGGGCTGGTATTTCTGTTTCGAAGCACGAGAATTTCTCGAAACTGGAGATGAGGCCGCGCGCTTAGCTGGTAACGCACCTTTTATTATTGATAAAGACAGTGGTGAAATTCATTCTCTGGGAAC‑GGCAAAACC‑G‑CTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGC‑TACCTTCGGCTTGCCCT‑g                                                                                                                                                                                                                              >  2:51614/1‑251 (MQ=255)
                                aCCATTACCTTGATGATGCAGATCTCCCGGTCGTCATTACTCTACATGGACGCTTTAGCCAGGGCTGGTATTTCTGTTTCGAAGCACGAGAATTTCTCGAAACTGGAGATGAGGCCGCGCGCTTAGCTGGTAACGCACCTTTTATTATTGATAAAGACAGTGGTGAAATTCATTCTCTGGGAAC‑GGCAAAACC‑G‑CTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGC‑TACCTTCGGCTTGCCCT‑g                                                                                                                                                                                                                              >  1:369330/1‑250 (MQ=255)
                                            atgatgCAGATCTCCCGGTCGTCATTACTCTACATGGACGCTTTAGCCAGGGCTGGTATTTCTGTTTCGACGCACGAGAATTTCTCGAAACTGGAGATGAGGCCGCGCGCTTATCTGGTAACGCACCTTTTATTATTGATAAAGACAGTGGTGAAATTCATTCTCTGGGAAC‑GGCAAAACC‑G‑CTGGAAGAATATCTACAGGATTACGAAATAAAAAGAGC‑TACCGTCAGCTTGCCCC‑GACAAAATAGCCct                                                                                                                                                                                                                 >  2:357605/1‑251 (MQ=255)
                                                   aGATCTCCCGGTCGTCATTACTCTACATGGACGCTTTAGCCAGGGCTGGTATTTCTGTTTCGAAGCACGAGAATTTCTCGAAACTGGAGATGAGGCCGCGCGCTTAGCTGGTAACGCACCTTTTATTATTGATAAAGACAGTGGTGAAATTCATTCTCTGGGAAC‑GGCAAAACC‑G‑CTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGC‑TACCTTCGGCTTGCCCT‑GACAAAATAGCCCTCTTccc                                                                                                                                                                                                           >  1:401292/1‑250 (MQ=255)
                                                   aGATCTCCCGGTCGTCATTACTCTACATGGACGCTTTAGCCAGGGCTGGTATTTCTGTTTCGAAGCACGAGAATTTCTCGAAACTGGAGATGAGGCCGCGCGCTTAGCTGGTAACGCACCTTTTATTATTGATAAAGACAGTGGTGAAATTCATTCTCTGGGAAC‑GGCAAAACC‑G‑CTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGC‑TACCTTCGGCTTGCCCT‑GACAAAATAGCCCTCTTCCCa                                                                                                                                                                                                          <  1:35402/251‑1 (MQ=255)
                                                    gATCTCCCGGTCGTCATTACTCTACATGGACGCTTTAGCCAGGGCTGGTATTTCTGTTTCGAAGCACGAGAATTTCTCGAAACTGGAGATGAGGCCGCGCGCTTAGCTGGTAACGCACCTTTTATTATTGATAAAGACAGTGGTGAAATTCATTCTCTGGGAAC‑GGCAAAACC‑G‑CTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGC‑TACCTTCGGCTTGCCCT‑GACAAAATAGCCCTCTTCCCAc                                                                                                                                                                                                         >  1:286039/1‑251 (MQ=255)
                                                                gtcATTACTCTACATGGACGCTTTAGCCAGGGCTGGTATTTCTGTTTCGAAGCACGAGAATTTCTCGAAACTGGAGATGAGGCCGCGCGCTTAGCTGGTAACGCACCTTTTATTATTGATAAAGACAGTGGTGAAATTCATTCTCTGGGAAC‑GGCAAAACC‑G‑CTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGC‑TACCTTCGGCTTGCCCT‑GACAAAATAGCCCTCTTCCCACGAAGAGGGCCGc                                                                                                                                                                                             >  1:111448/1‑251 (MQ=255)
                                                                      aCTCTACATGGACGCTTTAGCCAGGGCTGGTATTTCTGTTTCGAAGCACGAGAATTTCTCGAAACTGGAGATGAGGCCGCGCGCTTAGCTGGTAACGCACCTTTTATTATTGATAAAGACAGTGGTGAAATTCATTCTCTGGGAAC‑GGCAAAACC‑G‑CTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGC‑TACCTTCGGCTTGCCCT‑GACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAAcc                                                                                                                                                                                        <  2:370078/250‑1 (MQ=255)
                                                                                     tttAGCCAGGGCTGGTATTTCTGTTTCGAAGCACGAGAATTTCTCGAAACTGGAGATGAGGCCGCGAGCTTAGCTGGTAACGCACCTTTTATTATTGATAAAGACAGTGGTGAAATTCATTCTCTGGGAAC‑GGCAAAACC‑G‑CTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGC‑TACCTTCGGCTTGCCCT‑GACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTT‑AGTACCCCGCCGtt                                                                                                                                                                       >  1:79230/1‑251 (MQ=255)
                                                                                                         cTGTTTCGAAGCACGAGAATTTCTCGAAACTGGAGATGAGGCCGCGCGCTTAGCTGGTAACGCACCTTTTATTATTGATAAAGACAGTGGTGAAATTCATTCTCTGGGAAC‑GGCAAAACC‑G‑CTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGC‑TACCTTCGGCTTGCCCT‑GACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTT‑AGTACCCCGCCGTTAAATCATCCACCGAGCGCgg                                                                                                                                                   >  2:98824/1‑251 (MQ=255)
                                                                                                             ttCGAAGCAGGAGATTTTCTCGAAACTGGAGTGGGGGCCGAGCGTTTAGTGGGTACCGCACTTTTTATTATTGAAAAAGAAAGTGGGAAAATTCTTTCTC‑GGGAAA‑GGGAAAACCCG‑CTGGAAAATTATCTACAGAATTACAAAAAAAAAAGGGC‑TCCTTTCGG‑TTGCCCTTGAAAAAATAGCCCTTTTCCCAGGAAGTGGGCCGTTACCTTT‑AGTTCCCCGCCGTTAAATCATCCACCGAGCGCGGGTc                                                                                                                                                <  1:202865/250‑1 (MQ=255)
                                                                                                                   gcaggagaTTTTCTCGAAACTGGAGAGGGGGCCGCGCGTTTGGCTGTAAACGCACTTTTTATTATTAATAAAGACAGTGGTAAAATTCATCCTCTGGGAAC‑GGCAAACCC‑G‑CGGGAAGAATTTCTACAGATTTACAAAAAAAAAAAGGC‑TCCCTTCGGCTTCCCCT‑GACAAAATACCCCCCTCCCCACGAAGGGGCCCGCAACCCTT‑AGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATgcg                                                                                                                                         <  1:86569/247‑1 (MQ=255)
                                                                                                                                  gAAACTGGAGATGAGGCCGCGCGCTTAGCTGGTAACGCACCTTTTATTATTGATAAAGACAGTGGTGAAATTCATTCTCTGGGAAC‑GGCAAAACC‑G‑CTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGC‑TACCTTCGGCTTGCCCT‑GACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTT‑AGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCg                                                                                                                          >  2:318568/1‑251 (MQ=255)
                                                                                                                                      cgggagaggggcccGCGCGCTTAGCTGGTCCCCCCCTTTTTATTATTGAAAAAAACGGTGGTAAAATTCATTCTC‑GGGAACCGGCAAACCC‑G‑CTGAAAAAATACCTACAGGATTACAAAAAAAAAAAGGC‑TCCCTCCGCCTTCCCCT‑ACCAAAATACCCCCCTCCCCCCGAAGGGGCCCGCTACCCTT‑AGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCg                                                                                                                      <  1:291576/239‑1 (MQ=255)
                                                                                                                                              ggggccggggtcTTAGCTGGTACCGCACCTTTTATTATGGAAAAAAACGGTGGAAAAATTATTTCTGTGGGAAC‑GACAAACCC‑G‑GTGAAAAAATATCTACAGAATGAC‑AAAAAAAAAAGGCCTCCTTTCGCCTTCCCCG‑GAAAAAATACCCCTCTCCCCACGAAGGGGCCCGCAACCCTT‑GTTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAAcc                                                                                                              <  1:404639/240‑1 (MQ=255)
                                                                                                                                                ggCCGCGCGCTTAGCTGGTAACGCACCTTTTATTATTGATAAAGACAGTGGTGAAATTCATTCTCTGGGAAC‑GGCAAAACC‑G‑CTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGC‑TACCTTCGGCTTGCCCT‑GACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTT‑AGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCAt                                                                                                            >  1:205488/1‑251 (MQ=255)
                                                                                                                                                        gcTTAGCTGGTAACGCACCTTTTATTATTGATAAAGACAGTGGTGAAATTCATTCTCTGGGAAC‑GGCAAAACC‑G‑CTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGC‑TACCTTCGGCTTGCCCT‑GACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTT‑AGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCttt                                                                                                    >  1:357795/1‑251 (MQ=255)
                                                                                                                                                            aggtGGTACCGCACTTTTTATTTTTAAAAAAGACGGTGGTAAATTTCATTCGATGG‑AAC‑CGAAAAACCCG‑TTGAAAAAAAATCTCCGGGATTAAGAAAAAAAAAGGGC‑TTCTTTCGCCTGCCCCT‑GACAAAATCGCCCTCCCCCCAGGAGGGGGCCCGCAAACCTT‑ATTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTtgtg                                                                                                 <  1:446682/247‑1 (MQ=255)
                                                                                                                                                                  aaaCGCACTTTTTATTATGAAAAAAAACGGTGTTGAATTTAATTCTC‑GGGAAC‑CGGAAAACC‑GGCGGAAAAAATTTCTACAGGTTTAAAAAAAAAAAAAGCC‑TCCCCCCGGCTTCCCCA‑AAAAAAACACCCCTCTCCCCACGAAGAGGGCCGCTACCCTT‑AGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGcc                                                                                          <  2:357795/250‑1 (MQ=255)
                                                                                                                                                                           ttttATTATTAAAAAAGAGAGTGGTAAAATTCATTCTCGGGGACC‑GACAAAACC‑G‑CGGGAAGAATATCTACAGATTTAAAAAAAAAAAAAGCC‑TACCTTCGGCTTCCCCG‑AACAAATTACCCCTCCCCCCCCGAGGAGGGCCCCTACCCTT‑AGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCct                                                                                  <  2:54124/250‑1 (MQ=255)
                                                                                                                                                                                   ttGATAAAGACAGTGGTGAAAGTCATTCTCTGGGAAC‑GGCAAAACC‑G‑CTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGC‑TACCTTCGGCTTGCCCT‑GACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTT‑AGTACACTGCCGTTAAATCATCCACCGAGCGCGGGTCGGCTGCGCCGTACAACTCACCGTCCGGCCCAACCATACAGCTTTGTGTACTGCACATCTCCTCTTTCAgcc                                                                         >  2:258754/1‑250 (MQ=255)
                                                                                                                                                                                           gACAGTGGTGAAATTCATTCTCGGGAAAC‑GGAAAACCC‑G‑CTGAAAGAATATCTACAGTATTACAAAAAAAAAAAGGC‑TCCCTTCGGCTTCCCCT‑GAAAAAATACCCCTCTTCCCACGAAGAGGCCCGCTACCCTT‑AGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTc                                                                 <  1:285624/251‑1 (MQ=255)
                                                                                                                                                                                            aaaGTGGTAAAATTCTTTCTCGGGGAAC‑GACAAACCC‑G‑CTGAAAGATTATCTACAGGATTACAAAAAAAAAAAGGC‑TCCCGCGGGCTCCCCCT‑GAAAAAACCGCCTTCTCCCCACGAAGAGG‑CCGGTTACCTTTAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCt                                                                <  1:430352/249‑1 (MQ=255)
                                                                                                                                                                                               gtggtgAAATTCATTCTCTGGGAAC‑GGCAAAACC‑G‑CTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGC‑TACCTTCGGCTTGCCCT‑GACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTT‑AGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCTCCTCTTTCAGCGCCACTTTCTGAc                                                             >  2:340750/1‑251 (MQ=255)
                                                                                                                                                                                                         tcattcTCGGGGAGC‑GGCAAACCC‑G‑CTGGAAGAATATCTACAGGATTACAAAAAAAAAAAGGC‑TACCTTCGGCTTCCCCT‑GACAAAATACCCCTCTCCCCACAAAGAGGCCCGCTACCCTT‑AGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCtt                                                    <  1:353823/250‑1 (MQ=255)
                                                                                                                                                                                                                     aaC‑GGAAAAACC‑G‑CGGAAAGAATATCTACAGGATTACAAAAAAAAAAAGGC‑TCCCTTCGGCTGGCCCT‑GACAAAATAGCCCTCTCCCCACGAAGGGGCCCGCTACCCTT‑AGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTGAg                                       <  1:548709/251‑1 (MQ=255)
                                                                                                                                                                                                                         ggCAAAACC‑G‑CTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGC‑TACCTTCGGCTTGCCCT‑GACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTT‑AGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTGAGCGt                                    >  1:266571/1‑251 (MQ=255)
                                                                                                                                                                                                                            aaaccc‑G‑CTGAAAAAATATCTACATGATTACAAAATAAAAAAGCC‑TCCCTTCGGCTTCCCCT‑GAAAAAATCCCCCCCTCCCCACGAAGGGGCCCGCCACCCTT‑AGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTGAGCGTATc                                 <  1:213136/247‑1 (MQ=255)
                                                                                                                                                                                                                                      tGGAAGAATATCTACAGGATTACGAAATAAAAAAGGC‑TACCTTCGGCTTGCCCT‑GACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTT‑AGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTGAGCGTATCAGGGCTa                          >  2:203285/1‑250 (MQ=255)
                                                                                                                                                                                                                                       ggAAGAATATCTACAGGATTACGAAATAAAAAAGGC‑TACCTTCGGCTTGCCCT‑GACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTT‑AGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTGAGCGTATCCGGGCTaa                         >  2:234114/1‑250 (MQ=255)
                                                                                                                                                                                                                                           gaaTATCTACAGGATTACGAAATAAAAAAGGC‑TACCTTCGGCTTGCCCT‑GACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTT‑AGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTGAGCGTATCCGGGCTAAACCCt                    >  1:99220/1‑251 (MQ=255)
                                                                                                                                                                                                                                             atatCTACAGGATTACGCAATAAAAAAGGC‑TACCTTCGGCTTGCCTT‑GACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTT‑AGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTGAGCGTATCCGGGCTAAACCCttt                  >  1:487077/1‑251 (MQ=255)
                                                                                                                                                                                                                                                   ccAGGATTACGAAATAAAAAAGGC‑TACCTTCGGCTTGCCCT‑GACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTT‑AGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTGAGCGTATCCGGGCTAAACCCTTTTTCGac            >  2:324718/2‑251 (MQ=255)
                                                                                                                                                                                                                                                    cAGGATTACGAAATAAAAAAGGC‑TACCTTCGGCTTGCCCT‑GACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTT‑AGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTGAGCGTATCCGGGCTAAACCCTTTTTCGac            >  2:413453/1‑250 (MQ=255)
                                                                                                                                                                                                                                                             gAAATAAAAAAGGC‑TACCTTCGGCTTGCCCT‑GACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTT‑AGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTGAGCGTATCCGGGCTAAACCCTTTTTCGACACGCAGCTCg  >  1:155715/1‑251 (MQ=255)
                                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                             
GATTACTTATCACGACGCATTCGCGAAAGCGAACCATTACCTTGATGATGCAGATCTCCCGGTCGTCATTACTCTACATGGACGCTTTAGCCAGGGCTGGTATTTCTGTTTCGAAGCACGAGAATTTCTCGAAACTGGAGATGAGGCCGCGCGCTTAGCTGGTAACGCACCTTTTATTATTGATAAAGACAGTGGTGAAATTCATTCTCTGGGAAC‑GGCAAAACC‑G‑CTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGC‑TACCTTCGGCTTGCCCT‑GACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTT‑AGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTGAGCGTATCCGGGCTAAACCCTTTTTCGACACGCAGCTCG  >  NC_000913/3584761‑3585261

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