Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,585,0580TC25.5% 111.5 / 11.6 47intergenic (+15/+23)yrhB/ggtputative heat shock chaperone/glutathione hydrolase proenzyme
Reads supporting (aligned to +/- strand):  ref base T (34/1);  new base C (0/12);  total (34/13)
Fisher's exact test for biased strand distribution p-value = 2.49e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GCAGATCTCCCGGTCGTCATTACTCTACATGGACGCTTTAGCCAGGGCTGGTATTTCTGTTTCGAAGCACGAGAATTTCTCGAAACTGGAGATGAGGCCGCGCGCTTAGCTGGTAACGCACCTTTTATTATTGATAAAGACAGT‑GGTGAAATTCATTCTCTGGGAAC‑GGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGG‑CT‑TGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTT‑A‑GTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTGAGCGTATCCGGGCTAAACCCTTTTTCGACACGCAGCTCGTCCGGCAACCACTGATGGTGGAAACGCGGCGCAT  >  NC_000913/3584810‑3585295
                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                               
gCAGATCTCCCGGTCGTCATTACTCTACATGGACGCTTTAGCCAGGGCTGGTATTTCTGTTTCGAAGCACGAGAATTTCTCGAAACTGGAGATGAGGCCGCGCGCTTAGCTGGTAACGCACCTTTTATTATTGATAAAGACAGT‑GGTGAAATTCATTCTCTGGGAAC‑GGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGG‑CT‑TGCCCTGACAAAATAGCCCTCTTcc                                                                                                                                                                                                                                               >  1:276394/1‑251 (MQ=255)
           ggTCGTCATTACTCTACATGGACGCTTTAGCCAGGGCTGGTATTTCTGTTTCGAAGCACGAGAATTTCTCGAAACTGGAGATGAGGCCGCGCGCTTAGCTGGTAACGCACCTTTTATTATTGATAAAGACAGT‑GGAGAAATTCATTCTCTGGGAAC‑GGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGG‑CT‑TGCCCTGACAAAATAGCCCTCTTCCCACGAAGAggg                                                                                                                                                                                                                                    >  2:68863/1‑251 (MQ=255)
                          aCATGGACGCTTTAGCCAGGGCTGGTATTTCTGTTTCGAAGCACGAGAATTTCTCGAAACTGGAGATGAGGCCGCGCGCTTAGCTGGTAACGCACCTTTTATTATTGATAAAGACAGT‑GGTGAAATTCATTCTCTGGGAAC‑GGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGG‑CT‑TGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTT‑A‑Gt                                                                                                                                                                                                                    <  2:24416/250‑1 (MQ=255)
                                      tAGCCAGGGCTGGTATTTCTGTTTCGAAGCACGAGAATTTCTCGAAACTGGAGATGAGGCCGCGCGCTTAGCTGGTAACGCACCTTTTATTATTGATAAAGACAGT‑GGTGAAATTCATTCTCTGGGAAC‑GGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGG‑CT‑TGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTT‑A‑GTACCCCGCCGTTa                                                                                                                                                                                                        >  1:30101/1‑250 (MQ=255)
                                                     tttCTGTTTCGAAGCACGAGAATTTCTCGAAACTGGAGATGAGGCCGCGCGCTTAGCTGGTAACGCACCTTTTATTATTGATAAAGACAGT‑GGTGAAATTCATTCTCTGGGAAC‑GGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGG‑CT‑TGCCCTGACAAAATAGCCCTCTTACCACGAAGAGGGCCGCTAACCTT‑A‑GTACCCCGCCGTTAAATCATCCACcgcgcg                                                                                                                                                                                        >  2:295308/1‑251 (MQ=255)
                                                                                               ggCCGCGCGCTTAGCTGGTAACGCACCTTTTATTATTGATAAAGACAGT‑GGTGAAATTCATTCTCTGGGAAC‑GGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGG‑CT‑TGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTT‑A‑GTACCCCGCCGTTAAat                                                                                                                                                                                                     >  1:231733/1‑196 (MQ=255)
                                                                                               gcccGCGCGTTTAGCTGGTAACGCCCTTTTTATTATTGAAAAAAACAGG‑GGTAAAATTCATTCCCGGGGAAC‑GAAAAAACGGCTGAAAAAATATCTACAG‑ATTAAAAAAAAAAAAAAGGTTCCTTCGGGTT‑TCCCCTGAAAAAATAGCCCTCTCCCCACGAAGAGGGCCGCTAACCTT‑A‑GTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCAt                                                                                                                                              <  1:292327/249‑1 (MQ=255)
                                                                                                       gcTTAGCTGGTAACGCACCTTTTATTATTGATAAAGACAGT‑GGTGAAATTCATTCTCTGGGAAC‑GGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGG‑CT‑TGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTT‑A‑GTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCttt                                                                                                                                      >  1:347107/1‑251 (MQ=255)
                                                                                                        ataagttGGTACCGCACATTTTATAATGTATAAAGAGAGT‑GTTGAAATTATTTCGC‑GGGAACCGGAAAACCCGCTGGAAAAATATCAATAGGATTGAAAAAAAAAAAGGGCT‑CCTTTGG‑CTGTGCCTTGAAAAAATCGCCCTCTCCCCAGGAAG‑GGGCCCTTTACCTTTA‑TTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGTATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTtg                                                                                                                                     <  2:126950/245‑1 (MQ=255)
                                                                                                          tAGCTGGTAACGCACCTTTTATTATTGATAAAGACAGT‑GGTGAAATTCATTCTCTGGGAAA‑GGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGG‑CT‑TGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTT‑A‑GTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTtgtg                                                                                                                                   >  1:549377/1‑251 (MQ=255)
                                                                                                          gaGCTGGTAACCCCCTTTTTATTATGGAAAAAG‑CAGTTGGTAAATTTCTTCCTCGGGGACC‑GGAAAACCCGCTGAAAAAATTTCTACAGGATTAAAAAAAAAAAAAGGCTACCTTCGG‑CT‑TCCCCTGAAAAAATAGCCCTCTCCCCAAGAAGAGGGCCGCTAACCTT‑A‑TTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTtgtg                                                                                                                                   <  1:396205/250‑1 (MQ=255)
                                                                                                                  aaCGCACCTTTTATTATTGATACAGAAAGT‑GGTGAAATTCATTCTCTGGGCAC‑GGCAAAACCGCTGGAAGAATATATACAGGATTACGAAATAAAAAAGGCTACCTTCGG‑CT‑TGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTT‑A‑GTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATCATGCTTTGCGTACTGcac                                                                                                                           >  2:425932/1‑249 (MQ=255)
                                                                                                                    cGCACCTTTTATTATTGATAAAGACAGT‑GGTGAAATTCATTCTCTGGGAAC‑GGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGG‑CT‑TGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTT‑A‑GTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCAt                                                                                                                         >  2:29667/1‑251 (MQ=255)
                                                                                                                             tattatTGATAAAGACAGT‑GGTGAAATTCATTCTCTGGGAAC‑GGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGG‑CT‑TGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTT‑A‑GTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCttt                                                                                                                >  2:117767/1‑251 (MQ=255)
                                                                                                                              attattGATAAAGACAGT‑GGTGAAATTCATTCTCTGGGAAC‑GGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGG‑CT‑TGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTT‑A‑GTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTc                                                                                                               >  1:310648/1‑251 (MQ=255)
                                                                                                                                    gATAAAGACAGT‑GGTGAAATTCATTCTCTGGGAAC‑GGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGG‑CT‑TGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTT‑A‑GTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGcc                                                                                                         >  1:335248/1‑251 (MQ=255)
                                                                                                                                       aaaGACAGT‑GGTGAAATTCATTCTCTGGGAAC‑GGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGG‑CT‑TGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTT‑A‑GTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACt                                                                                                      >  2:551185/1‑251 (MQ=255)
                                                                                                                                          gACAGT‑GGTGAAATTCATTCTCTGGGAAC‑GGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGG‑CT‑TGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTT‑A‑GTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTc                                                                                                   >  1:47693/1‑251 (MQ=255)
                                                                                                                                                                   ggAAC‑GGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGG‑CT‑TGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTT‑A‑GTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTg                                                                           >  2:42600/1‑251 (MQ=255)
                                                                                                                                                                     aaacgcaaacaCCGCGGGAAAAATATCTAAAGGATAAAAAAAAAAAAAAGCCTCCCTTCGG‑CT‑CCCCATAAAAAAATCCCCCTCTCCCCAGGAAGAGGGCCGCAACCCTT‑A‑GTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTGa                                                                          <  2:30101/241‑1 (MQ=255)
                                                                                                                                                                     aaC‑GGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGG‑CT‑TGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTT‑A‑GTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTg                                                                           >  1:495215/1‑249 (MQ=255)
                                                                                                                                                                            aaaaCCGCTGGAGGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGG‑CT‑TGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTT‑A‑GTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTGAGCGTATc                                                                   >  1:387503/1‑251 (MQ=255)
                                                                                                                                                                              aaCCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGG‑CT‑TGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTT‑A‑GTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTGAGCGTATCCg                                                                 >  2:59988/1‑251 (MQ=255)
                                                                                                                                                                                ccGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGG‑CT‑TGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTT‑A‑GTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTGAGCGTATCCggg                                                               >  2:327521/1‑251 (MQ=255)
                                                                                                                                                                                        agaataattacagaattaaaaaaaaaaaaaGGCTCCCTCCGG‑CT‑TCCCCGAAAAAAACACCCCTCTCCCCAGGAAGGGGCCCGC‑AACCTT‑ATGCCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTGAGCGTATCCGGGCTAAAccc                                                       <  2:310648/227‑1 (MQ=255)
                                                                                                                                                                                        agaaTATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGG‑CT‑TGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTT‑A‑GTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGCTGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCCGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTGAGCGTATCCGGGCTAAAccc                                                       >  2:150953/1‑251 (MQ=255)
                                                                                                                                                                                          aaTATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGG‑CT‑TGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTT‑A‑GTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTGAGCGTATCCGGGCTAAACCCtt                                                     >  2:174160/1‑251 (MQ=255)
                                                                                                                                                                                               caaagggataaagAAAAAAAAAAGCCTCCCTGCGG‑CT‑TCCCCTAAAAAAATCCCCCTCTCCCCACGAAGGGGCCCGCTACCCTT‑A‑GTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTGAGCGTATCCGGGCTAAACCCTTTTTc                                                 <  1:329983/238‑1 (MQ=255)
                                                                                                                                                                                                 ccagatttACAAAATAAAAAAGCCTCCCTTCCG‑CT‑TCCCCTGACAAATTACCCCCCTCCCCCCAAAGAGGCCCCCTAACTTT‑A‑TTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTGAGCGTATCCGGGCTAAACCCTTTTTCGac                                              <  1:490415/245‑1 (MQ=255)
                                                                                                                                                                                                 aCAGGATTACGAAATAAAAAAGGCTACCTTCGG‑CT‑TGCCCTGACAAAATAGCACTCTTCCCACGAAGAGGGCCGCTAACCTT‑A‑GTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTGAGCGTATCCGGGCTAAACCCTTTTTCGac                                              >  2:335727/1‑251 (MQ=255)
                                                                                                                                                                                                    ggATTACGAAATAAAAAAGGCTACCTTCGG‑CT‑TGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTT‑A‑GTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTGAGCGTATCCGGGCTAAACCCTTTTTCGACACg                                           >  1:551558/1‑251 (MQ=255)
                                                                                                                                                                                                             aaaaaaaaaGTCTCCCTGCGG‑CT‑TGCCCTAAAAAAAGCCCCCTCTCCCCACGAAGGGGCCCGCTACCCTT‑A‑GCCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTGAGCGTATCCGGGCTAAACCCTTTTTCGACACGCAGCTCGTc                                  <  1:31156/248‑1 (MQ=255)
                                                                                                                                                                                                              ntAAAAAAGGCTACCTTCGG‑CT‑TGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTT‑A‑GTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTGAGCGTATCCGGGCTAAACCCTTTTTCGACACGCAGCTCGTcc                                 >  2:245186/2‑251 (MQ=255)
                                                                                                                                                                                                              aaaaaaaaGGCTCCCTCCGG‑CT‑TCCCCTAAAAAAATCCCCCTCTCCCCACGAAGGGGGCCGCTACCCTT‑A‑GTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTGAGCGTATCCGGGCTAAACCCTTTTTCGACACGCAGCTCGTcc                                 <  1:245186/249‑1 (MQ=255)
                                                                                                                                                                                                              aTAAAAAAGGCTACCTTCGG‑CT‑TGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTT‑A‑GTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTGAGCGTAt                                                                    >  1:403507/1‑216 (MQ=255)
                                                                                                                                                                                                                  aaaaGGCTACCTTCGG‑CT‑TGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTT‑A‑GTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTGAGCGTATCCGGGCTAAACCCTTTTTCGACACGCAGCTCGTCCGGc                              >  2:282117/1‑250 (MQ=255)
                                                                                                                                                                                                                      ggCTACCTTCGG‑CT‑TGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTT‑A‑GTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTGAGCGTATCCGGGCTAAACCCTTTTTCGACACGCAGCTCGTCCGGCAAcc                          >  2:463454/1‑250 (MQ=255)
                                                                                                                                                                                                                      ggCTACCTTCGG‑CT‑TGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTT‑A‑GTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTGAGCGTATCCGGGCTAAACCCTTTTTCGACACGCAGCTCGTCCGGCAAcc                          >  1:553440/1‑250 (MQ=255)
                                                                                                                                                                                                                               tgcggcctccccgaaaaaaaTACCCCCCTCCCCACGAAGGGGCCCGCTACCCTT‑A‑GTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTGAGCGTATCCGGGCTAAACCCTTTTTCGACACGCAGCTCGTCCGGCAACCACTGAt                    <  1:59988/233‑1 (MQ=255)
                                                                                                                                                                                                                                 ggcctcccccTAACAAAAGACCCCTCCCCCCACGAAGGGGCCCGCTACCCTT‑A‑TTCCGCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTGAGCGTATCCGGGCTAAACCCTTTTTCGACACGCAGCTCGTCCGGCAACCACTGAtggtg                <  2:152500/244‑1 (MQ=255)
                                                                                                                                                                                                                                   ccttccccTGAAAAAATACCCCTCTCCCCACGAAGAGGCCCCCTACCCTT‑A‑GTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTGAGCGTATCCGGGCTAAACCCTTTTTCGACACGCAGCTCGTCCGGCAACCACTGATGGTGGa              <  1:174160/246‑1 (MQ=255)
                                                                                                                                                                                                                                         ccTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTT‑A‑GTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTGAGCGTATCCGGGCTAAACCCTTTTTCGACACGCAGCTCGTCCGGCAACCACTGATGGTGGAAAcgcg        >  1:131296/1‑251 (MQ=255)
                                                                                                                                                                                                                                          cTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTT‑A‑GTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTGAGCGTATCCGGGCTAAACCCTTTTTCGACACGCAGCTCGTCCGGCAACCACTGATGGTGGAAACgcgg       >  2:413322/1‑251 (MQ=255)
                                                                                                                                                                                                                                            gACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTT‑A‑GTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTGAGCGTATCCGGGCTAAACCCTTTTTCGACACGCAGCTCGTCCGGCAACCatcatcgaaatgtcgctcaac     >  1:431164/1‑231 (MQ=255)
                                                                                                                                                                                                                                             aCAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTT‑A‑GTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTGAGCGTATCCGGGCTAAACCCTTTTTCGACACGCAGCTCGTCCGGCAACCACTGATGGTGGAAACGCGgcgc    >  1:5024/1‑251 (MQ=255)
                                                                                                                                                                                                                                               aaaaTAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTT‑A‑GTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTGAGCGTATCCGGGCTAAACCCTTTTTCGACACGCAGCTCGTCCg                                >  2:56077/1‑221 (MQ=255)
                                                                                                                                                                                                                                               aaaaTAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTT‑A‑GTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTGAGCGTATCCGGGCTAAACCCTTTTTCGACACGCAGCTCGTCCGGCAACCACTGATGGTGGAAACGCGGCGCAt  >  2:118733/1‑251 (MQ=255)
                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                               
GCAGATCTCCCGGTCGTCATTACTCTACATGGACGCTTTAGCCAGGGCTGGTATTTCTGTTTCGAAGCACGAGAATTTCTCGAAACTGGAGATGAGGCCGCGCGCTTAGCTGGTAACGCACCTTTTATTATTGATAAAGACAGT‑GGTGAAATTCATTCTCTGGGAAC‑GGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGG‑CT‑TGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTT‑A‑GTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTGAGCGTATCCGGGCTAAACCCTTTTTCGACACGCAGCTCGTCCGGCAACCACTGATGGTGGAAACGCGGCGCAT  >  NC_000913/3584810‑3585295

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