Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,287,420 G→C 7.5% intergenic (‑92/+118) tyrV ← / ← tyrT tRNA‑Tyr/tRNA‑Tyr

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,287,4200GC 7.5% 159.0 / 2.8 51intergenic (‑92/+118)tyrV/tyrTtRNA‑Tyr/tRNA‑Tyr
Reads supporting (aligned to +/- strand):  ref base G (29/18);  new base C (2/2);  total (31/20)
Fisher's exact test for biased strand distribution p-value = 6.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.60e-01

GCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTACAGTCTGCTCCCTTTGGCCGCTCGGGAACCCCACCAGGGGTAATTCAAATTTTGAGGTAATGCTTGAGATGGTGGTGGGGGAAGGATTATTCGTCGCTTCGCTCCTCACCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTTCTCACCCTTCCCCGATGAGTGCAAACTTTCACAATCTCACCGAAGTTACCACATCGCTGTGGTGAATTATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTACAGTCTGCTCCCTTTGGCCGCTCGGGAACCCCACCACGGGGTAATGCTTTTTACTGGCCTGCTCCCTTATCGGGAAGCGGGGCGCATCATATCAAATGACGCGCCGCTGTAAAGTGTTACGTTGAGAAAAATG  >  NC_000913/1287131‑1287720
                                                                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                                                                            
gacgtgtgcttttccgattgtttgAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTACAGTCTGCTCCCTTTGGCCGCTCGGGAACCCCACCAGGGGTAATTCAAATTTTGAGGTAATGCTTGAGATGGTGGTGGGGGAAGGATTATTCGTCGCTTCGCTCCTCACCCTTCGGGCCGTTGCCACAGGCAACggcc                                                                                                                                                                                                                                                                                                   <  2:595204/278‑4 (MQ=38)
     ctcGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTACAGTCTGCTCCCTTTGGCCGCTCGGGAACCCCACCAGGGGTAATTCAAATTTTGAGGTAATGCTTGAGATGGTGGTGGGGGAAGGATTATTCGTCGCTTCGCTCCTCACCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCg                                                                                                                                                                                                                                                                                              >  2:176202/1‑301 (MQ=255)
                    gTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTACAGTCTGCCCCCTTTGGCCGCTCGGGAACCCCACCAGGGGTAATTCAAATTTTGAGGTAATGCTTGAGATGGTGGTGGGGGAAGGATTATTCGTCGCTTCGCTCCTCACCCTTCGGGCCGTAGCCACAGGCAACggcaagaacggaacgacgacgggg                                                                                                                                                                                                                                                                               >  1:594443/1‑278 (MQ=38)
                        aaCCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGGGATGGTGGTGGGGGAAGGAGTCGACCCTTCGAAGTCGATGACGGCAGATTTACAGTCTGCTCACTTTGGCCGCTCGGGAACCCCACCAGGGGTAATCCAAATGATGAGGTAATGCGTGAGATGGTGGTGGGGGCACGATTATTCGCAGCTTCGCTCCTCACCCCTCGGGCCGTAGCCAGAGGCAACGCACTCCCGCATCCCCCCCAACCGAAcc                                                                                                                                                                                                                                                                           >  1:285616/1‑301 (MQ=25)
                                     gCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTACAGTCTGCTCCCTTTGGCCGCTCGGGAACCCCACCAGGGGTAATTCAAATTTTGAGGTAATGCTTGAGATGGTGGTGGGGGAAGGATTATTCGTCGCTTCGCTCCTCACCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGtt                                                                                                                                                                                                                                                              >  2:621593/1‑301 (MQ=255)
                                       ttCTCCTCCTTCCCCGCTTGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTTAAAGTGTTGGTGGTGGTGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTACAGTCTGCTCCCTTTGGCCGCTCGGGAACCCCACCAGGGGTAATTCAAATTTTGAGGTAATGCTTGAGATGGTGGTGGGGGAAGGATTATTCGTCGCTTCGCTCCTCACCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTtct                                                                                                                                                                                                                                                            <  1:620735/301‑1 (MQ=255)
                                           aaTCCTTCCCGGCATGGTCAGAATATTTAATTGCGTGTTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGAGGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTACAGTCTGCTCCCTTTGGCCGCTCGGGAACCCCACCAGGGGTAATTCAAATTTTGAGGTAATGCTTGAGATGGTGGTGGGGGAAGGATTATTCGTCGCTTCGCTCCTCACCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTTCTCAcc                                                                                                                                                                                                                                                        <  1:156276/300‑1 (MQ=255)
                                                  ccccGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTACAGTCTGCTCCCTTTGGCCGCTCGGGAACCCCACCAGGGGTAATTCAAATTTTGAGGTAATGCTTGAGATGGTGGTGGGGGAAGGATTATTCGTCGCTTCGCTCCTCACCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCCAGCTTTCGCCCGAATCGAACCATAGTCGAAAGTTACCACCCTTcccc                                                                                                                                                                                                                                                 >  1:547223/1‑301 (MQ=255)
                                                                                        gAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTACAGTCTGCTCCCTTTGGCCGCTCGGGAACCCCACCAGGGGTAATTCAAATTTTGAGGTAATGCTTGAGATGGTGGTGGGGGAAGGATTATTCGTCGCTTCGCTCCTCACCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTTCTCACCCTTCCCCGATGAGTGCAAACTTTCACAATCTCACCGAAGTTACca                                                                                                                                                                                                           <  2:547896/301‑1 (MQ=255)
                                                                                           gTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTACAGTCTGCTCCCTTTGGCCGCTCGGGAACCCCACCAGGGGTAATTCAAATTTTGAGGTAATGCTTGAGATGGTGGTGGGGGAAGGATTATTCGTCGCTTCGCTCCTAACCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGCTCTCACCCTTCCCACAAGAGTGCAAACTTTCACAACCCAACCAAAGTTCCcacaa                                                                                                                                                                                                        >  1:282090/1‑300 (MQ=255)
                                                                                                   atttttGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTACAGTCTGCTCCCTTTGGCCGCTCGGGAACCCCACCAGGGGTAATTCAAATTTTGAGGTAATGCTTGAGATGGTGGTGGGGGAAGGATTATTCGTCGCTTCGCTCCTCACCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTTCTCACCCTTCCCCGATGAGTGCAAACTTTCACAATCTCACCGAAGTTACCACATCGCTgtgg                                                                                                                                                                                                <  2:285919/299‑1 (MQ=255)
                                                                                                              tGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTACAGTCTGCTCCCTTTGGCCGCTCGGGAACCCCACCAGGGGTAATTCAAATTTTGAGGTAATGCTTGAGATGGTGGTGGGGGAAGGATTATTCGTCGCTTCGCTCCTCCCCCTTCGGGCCGTTTCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCCAAGGTTCTCACCCTTCCCCGATGAGTGCAAACTTTCACAATCTCACCGAAGTTACCACATCGCTGTGGTAAATTAtggt                                                                                                                                                                                     >  2:229611/1‑301 (MQ=255)
                                                                                                                         ggggAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTACAGTCTGCTCCCTTTGGCCGCTCGGGAACCCCACCAGGGGTAATTCAAATTTTGAGGTAATGCTTGAGATGGTGGTGGGGGAAGGATTATTCGTCGCTTCGCTCCTCACCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTACTCACCCTTCCCCGATGAGTGCAAACTTTCCCAATCTCACCGAAGTTACCACATCGCTGTGGTGAATTATGGTGGTGGGGGaaa                                                                                                                                                                          >  2:178760/1‑300 (MQ=255)
                                                                                                                              aGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTACAGTCTGCTCCCTTTGGCCGCTCGGGAACCCCACCAGGGGTAATTCAAATTTTGAGGTAATGCTTGAGATGGTGGTGGGGGAAGGATTATTCGTCGCTTCGCTCCTCACCATTCGGTCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCTAACCTTAGACGAAGGTTCTCACCCTTCCCCGATGGATACAAGATTTCACAATCTAACCGAATATACCACATCGCTGTGGTCCATTAAGGTTGTtgggtacagcatt                                                                                                                                                                     >  1:258733/1‑292 (MQ=38)
                                                                                                                                 tttCGAACCTTCGAAGTCGATGACGGCAGATTTACAGTCTGCTCCCTTTGGCCGCTCGGGAACCCCACCAGGGGTAATTCAAATTTTGAGGTAATGCTTGAGATGGTGGTGGGGGAAGGATTATTCGTCGCTTCGCTCCTCACCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTTCTCACCCTTCCCCGATGAGTGCAAACTTTCACAATCTCACCGAAGTTACCACATCGCTGTGGTGAATTATGGTGGTGGGGGAAGGATTCGaa                                                                                                                                                                  <  2:282393/300‑1 (MQ=255)
                                                                                                                                 aTTCGAACCTTCGAAGTCGATGACGGCAGATTTACAGTCTGCTCCCTTTGGCCGCTCGGGAACCCCACCAGGGGTAATTCAAATTTTGAGGTAATGCTTGAGATGGTGGTGGGGGAAGGATTATTCGTCGCTTCGCTCCTCACCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTTCTCACCCTTCCCCGATGAGTGCAAACTTTCACAATCTCACCGAAGTTACCACATCGCTGTGGTGAATTATGGTGGTGGGGGAAGGATTCGaa                                                                                                                                                                  >  2:589059/1‑301 (MQ=255)
                                                                                                                                                      gtcGGTAGATTTACTGTCTGCTCTCTTTGGTCGATCGTGAACCCTACCTGGTGTAAGTCAAATTTTGAGTTAATGTTTGGGATGGTGGTGGCGGAAGGATTATTCGTCGCTTCGCTCCTCACCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTTCTCACCCTTCCCCGATGAGTGCAAACTTTCACAATCTCACCGAAGTTACCACATCGCTGTGGTGAATTATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCa                                                                                                                                             <  1:176005/299‑1 (MQ=38)
                                                                                                                                                       aCGGCAGATTTACAGTCTGCTCCCTTTGGCCGCTCGGGAACCCCACCAGGGGTAATTCAAATTTTGAGGTAATGCTTGAGATGGTGGTGGGGGAAGGATTATTCGTCGCTTCGCTCCTCACCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTTCTCACCCTTCCCCGATGAGTGCAAACTTTCACAATCTAACCGAAGTTACCACATCGCTGTGGTGAATTATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCaa                                                                                                                                            >  2:113475/1‑300 (MQ=255)
                                                                                                                                                          gCAGATTTACAGTCTGCTCCCTTTGGCCGCTCGGGAACCCCACCAGTGGTAATTCAACTTTTTAGGTAATGCTTGAGATGGCGGTGGGGGAAGGATTATTCGTCGCTTCGCTCCTCACCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGCTCTCACCCTTCCCCGATGAGTGCAAACTTTCACAATCTCACCGAAGTTACAACATCGCTGTGGTGACTAATGGTGGTGGGGGAAAGACTAGAACCTTCCAAGACAAAGACGGCAGaat                                                                                                                                         >  1:109361/1‑299 (MQ=255)
                                                                                                                                                               ttttCAGTCTGCTCCTGTGGGTCGCTCGGGAATCCCACCTGGGGGTATTCAAATTTTGGGTTAATGCTTGAGATGGTGGTGGGGGAATGATTATTCGTCGCTTCGCTCCTCACCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCGCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTTCTCACCCTTCCCTGATGAGTGCAAACTTGCACAATCTCACCGAAGTTACTACATCGCTGTGGTGAACTATGCTGGCGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTACAg                                                                                                                                    <  1:588335/301‑1 (MQ=38)
                                                                                                                                                                   cAGTCTGCTCCCTTTGGCCGCTCGGGAACCCCACCAGGGGTAATTCAAATTTTGAGGTCATGCTTGAGATGGTGGTGGGGGAAGGATTATTCGTCGCTTCGCTCCTCACCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAATGTTCACACCCTTCCCCGATGAGTGCAAACTTTCACAATCTCACCGACGTTACCACATCGCTGTGGTGAATTCTGGTGGTAGGGGAAGGATTCGAACCATCGAAATCGATACCGCCAGACTTACaacctt                                                                                                                                >  1:633321/1‑296 (MQ=255)
                                                                                                                                                                             ccTTTGGCCGCTCGGGAACCCCACCAGGGGTAATTCAAATTTTGAGGTAATGCTTGAGATGGTGGTGGGGGAAGGATTATTCGTCGCTTCGCTCCTCACCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTTCTCACCCTTCCCCGATGAGTGCAAACTTTCACAATCTAACCGAAGTTCCCACATAGCTGTGGTGAATTATGGTGGTGGGGAAAGGATTCAAACCTTCGAAATAGAAGACAGCAGATTTACAGTAAGCTCCCTTTgg                                                                                                                      >  2:287708/1‑301 (MQ=255)
                                                                                                                                                                                 caaaatttgaattacccctggtggggTAATTCAAATTTTGAGGTAATGCTTGAGATGGTGGTGGGGGAAGGATTATTCGTCGCTTCGCTCCTCACCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTTCTCACCCTTCCCCGATGAGTGCAAACTTTCACAATCTCACCGCAGTTACCCCATCTCTGTGGTGAATTATGGTGGTGGGGGAAGGATCCGAACCTTCGAAGTCGATACCAGCCAATTTACAATCTCATCCCTTTGGcccc                                                                                                                  >  1:574206/23‑299 (MQ=255)
                                                                                                                                                                                  ggCCGCTCGGGAACCCCACCAGGGGTAATTCAAATTTTGAGGTAATGCTTGAGATGGTGGTGGGGGAAGGATTATTCGTCGCTTCGCTCCTCCCCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGACCCTTAGTCGAAGGTTCTCACCCTTCCCCGATGATTGCCACCTTTAACCATCTCACCGCAGTTACCACATCGCCTCGGAGAATTATGGCGGCGGGGGAACGATTCCACCCTCCCAAGCAGATGACTGCAGATTCACAGTCTGCTCCCCTTAGccccc                                                                                                                 >  1:583503/1‑298 (MQ=255)
                                                                                                                                                                                   gCCGCTCGGGAACCCCACCAGGGGTAATTCAAATTTTGAGGTAATGCTTGAGATGGTGGTGGGGGAAGGATTATTCGTCGCTTCGCTCCTCACCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTTCTCACCCTTCCCCGATGAGTGCAAACTTTCACAATCTCACCGAAGTTACCACATCGCTGTGGTGAATTATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTACAGTCTGCTCCCTTTGGCCGCTc                                                                                                                <  2:109485/301‑1 (MQ=255)
                                                                                                                                                                                          ggtatccctACCTGTGGTAATTCAAATTTTGAGGTAATGCTTGAGATGGTGGTGGGGGAAGGATTATTCGTCGCTTCGCTCCTCCCCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTTCTCACCCTTCCCCGATGAGTGCAAACTTTCACAATCTCACCGAAGTTACCACATCGCTGTGGTGAATTATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTACAGTCTGCTCCCTTTGGCCGCTCGGGAAcc                                                                                                         <  2:99129/296‑1 (MQ=255)
                                                                                                                                                                                           atctccccACCAGGGGTAATTCAAATTTTGAGGTAATGCTTGGGATGGTGGTGGGGGAAGGATTATTCGTCGCTTCGCTCCTCACCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTTCTCACCCTTCCCCGATGAGTGCAAACTTTCACAATCTCACCGAAGTTACCACATCGCTGTGGTGAATTATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTACAGTCTGCTCCCTTTGGCCGCTCGGGAAccc                                                                                                        <  1:161664/297‑1 (MQ=255)
                                                                                                                                                                                             ttctcctccaGGGGTAATTCTAATTTGTAGGTAATGCTTGAGATGGTGGTGGGTGAATGATTATTCGTCGCTGCGCTCCTCACCCTTCGTGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCGTAGTCGAAGGTTCTCACCCTTCCCCGAGGAGTGAAAACTTTCTCAGTCTCACCGAAGTTACCACATCGCTGTGGTGAATTATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTACAGTCTGCTCCCTTTGGCCGCTCGGGAACCcca                                                                                                      <  2:5740/294‑1 (MQ=38)
                                                                                                                                                                                             aaCCCCACCAGGGGTAATTCAAATTTTGAGGTAATGCTTGAGATGGTGGTGGGGGAAGGATTATTCGTCGCTTCGCTCCTCACCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTTCTCACCCTTCCCCGATGAGTGCAAACTTTCACCATCTCACCGAAGTTACCACATCGCTGTGGTGAATTATGGTGGTGGGGGAATGATTCGAACCTTCGAAGTCGATGACGGCAGATTTACAGTCTGCTCCCTTTGGCCGCTCGGACAccccc                                                                                                      >  1:317466/1‑300 (MQ=255)
                                                                                                                                                                                               ccccACCAGGGGTAATTCAAATTTTGAGGTAATGCTTGAGATGGTGGTGGGGGAAGGATTATTCGTCGCTTCGCTCCTCACCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTTCTCACCCTTCCCCGATGAGTGCAAACTTTCACAATCTCACCGAAGTTACCACATCGCTGTGGTGAATTATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTACAGTCTGCTCCCTTTGGCCGCTCGGGAAcccagaa                                                                                                    >  2:161815/1‑297 (MQ=255)
                                                                                                                                                                                               ccccACCAGGGGTAATTCAAATTTTGAGGTAATGCTTGAGATGGTGGTGGGGGAAGGATTATTCGTCGCTTCGCTCCTCACCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTTCTCACCCTTCCCCGATGAGTGCAAACTTTCACAATCTCACCGAAGTTACCACATCGCTGTGGTGAATTATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTACAGTCTGCTCCCTTTGGCCGCTCGGGAACCccacc                                                                                                    >  2:375993/1‑301 (MQ=255)
                                                                                                                                                                                                                    ttttGAGGTAATGCTTGAGATGGTGGTGGGGGAAGGATTATTCGTCGATTCGCGCCTCACCCTTCGCGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAAACTTAGTCGAAGGTTCTCACCCTTCCCCGCTGAGTGCAAACTTTCACAATCTCACCGAAGTTACCACATCGCTGAGTTGAATTATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGCCGGAAGATTTACAGTCTGCTCCATTTAGCCACTCGTGAACCCCACCACGGGACAATGCTTTTTACtt                                                                               >  1:426133/1‑300 (MQ=38)
                                                                                                                                                                                                                       tGAGGTATTGCTTGAGATGGTGGTGGGGGAAGGATTATTCGTCGCTTCGCTCCTCCCCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGACCCTTAGTCGAAGGTTCTCACCCTTCCCCGATGAGTGCAAACTTTCACAATCTCACCGAAGTTACCACATCGCTGTGGTGAATTATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTACAGTCTGCTCCCTTTGGCCGCTCGGGAACCCCACCACGGGGTAATGCTTTTTACTGGcc                                                                            <  2:17256/301‑1 (MQ=255)
                                                                                                                                                                                                                               tGCTTGAGATGGTGGTGGGGGAAGGATTATTCGTCGCTTCGCTCCTCACCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTTCTCACCCTTCCCCGATGAGTGCAAACTTTCACAATCTCACCGAAGTTACCACATCGCTGTGGTGAATTATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTACAGTCTGCTCCCTTTGGCCGCTCGGGAACCCCACCACGGGGTAATGCTTTTTACTGGCCTGCTCCCt                                                                    >  2:301576/1‑301 (MQ=255)
                                                                                                                                                                                                                                gCTTGAGATGGTGGTGGGGGAAGGATTATTCGTCGCTTCGCTCCTCACCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTTCTCACCCTTCCCCGATGAGTGCAAACTTTCACAATCTCACCGAAGTTACCACATCGCTGTGGTGAATTATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTACAGTCTGCTCCCTTTGGCCGCCCGGGAACCCCACCACGGGGTAATGCTTTTAACTGGCCCGCTCCCtt                                                                   >  1:405904/1‑301 (MQ=255)
                                                                                                                                                                                                                                     tgaTGGTGGTGGGGGAAGGATTATTCGTCGCTTCGCTCCTCACCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTTCTCACCCTTCCCCGATGAGTGCAAACTTTCACAATCTCACCGAAGTTACCACATCGCTGTGGTGAATTATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTACAGTCTGCTCCCTTTGGCCGCTCGGGAACCCCACCACGGGGTAATGCTTTTTACTGGCCTGCTCCCTTATCgg                                                              <  2:574921/300‑1 (MQ=255)
                                                                                                                                                                                                                                      ggtggtggtggGGGGTGGTTTTTTCTGCGTTTCGCTCCTCACCATTTGGTCCTTTGCCGGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGATGGTTCGCACCCGTCCCCGATGAGATCTACCCTTCTAAACCCCACTGACCTTAACGCAGCGCTGCCGATCTTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTACAGTCTGCTCCCTTTGGCCGCTCGGGAACCCCACCACGGGGTAATGCTTTTTACTGGCCTGCTCCCTTATCggg                                                             <  1:229368/299‑1 (MQ=21)
                                                                                                                                                                                                                                            ggtggGGGAAGGATTATTCGTCGCTTCGCTCCTCACCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTTCTCACCCTTCCCCGATGAGTGCAAACTTTCACAATCTCACCGAAGTTACCACATCGCTGTGGTGAATTATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTACAGTCTGCTCCCTTTGGCCGCTCGGGAACCCCACCACGGGGTAATGCTTTTTACTGGCCTGCTCCCTTATCGGGAAGCgg                                                       <  2:406467/301‑1 (MQ=255)
                                                                                                                                                                                                                                               gggggAAGGATTATTCGTCGCTTCGCTCCTCACCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTTCTCACCCTTCCCCGATGAGTGCAAACTTTCACAATCTCACCGAAGTTACCACATCGCTGTGGTGAATTATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTACAGTCTGCTCCCTTTGGCCGCTCGGGAACCCCACCACGGGGTAATGCTTTTTACTGGCCTGCTCCCTTATCGGGAAGCggggg                                                    >  2:222800/1‑300 (MQ=255)
                                                                                                                                                                                                                                                       gATTATTCGTCGCTTCGCTCCTCACCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTTCTCACCCTTCCCCGATGAGTGCAAACTTTCACAATCTCACCGAAGTTACCACATCGCTGTGGTGAATTATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTACAGTCTGCTCCCTTTGGCCGCTCGGGAACCCCACCACGGGGTAATGCTTTTTACTGGCCTGCTCCCTTATCGGGAAGCGGGGCGCATCata                                            >  2:70527/1‑301 (MQ=255)
                                                                                                                                                                                                                                                       gATTATTCGTCGCTTCGCTCCTCACCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTTCTCACCCTTCCCCGATGAGTGCAAACTTTCACAATCTCACCGAAGTTACCACATCGCTGTGGTGAATTATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTACAGTCTGCTCCCTTTGGCCGCTCGGGAACCCCACCACGGGGTAATGCTTTTTACTGGCCTGCTCCCTTATCGGGAAACGGGGCGcatcaaa                                            >  2:482707/1‑299 (MQ=255)
                                                                                                                                                                                                                                                           tttCGTCGCTTCGCTCCTCTCCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTTCTCACCCTTCCCCGATGAGTGCAAACTTTCACAATCTCACCGAAGTTACCACATCGCTGTGGTGAATTATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTACAGTCTGCTCCCTTTGGCCGCTCGGGAACCCCACCACGGGGTAATGCTTTTTACTGGCCTGCTCCCTTATCGGGAAGCGGGGCGCATCATATCaa                                        <  1:222557/300‑1 (MQ=255)
                                                                                                                                                                                                                                                                       gctCCTCACCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTTCTCACCCTTCCCCGATGAGTGCAAACTTTCACAATCTCACCGAAGTTACCACATCGCTGTGGTGAATTATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTACAGTCTGCTCCCTTTGGCCGCTCGGGAACCCCACCACGGGGTAATGCTTTTTACTGGCCTGCTCCCTTATCGGGAAGCGGTGCGCATCATATCAAATGAAGcgccgc                            >  2:137218/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                         ctctgctctgttgtgtggTTGCCTGTGGCAATGTTCTCTCGCTTTCGCTCGCATCGAACCTTATTCTAATGTTCTCTGCTTTCCTCGATGGGGGCAAACTGTCACATACACACCGAAGCTACCTCATCGCTGACGTGAATTATGCGGTTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTACAGTCTGCTCCCTTTGGCCGCTCGGGAACCCCACCACGGGGTAATGCTTTTTACTGGCCTGCTCCCTTATCGGGAAGCGGGGCGCATCATATCAAATGACGCGCCGCTg                          <  1:178563/284‑1 (MQ=25)
                                                                                                                                                                                                                                                                               cccTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTCGTCGAAGGTTCTCACCCTTCCCCGATGAGTGCCAACTTTCACAATCTCACCGAAGTTACCACATCGCTGTGGTGAATTATGGTGGTGGGGGAAGGATTCGAACCTTCGCAGTCGCTGACGGCAGATTTACCGTCTGCTCCCTTTGGCCCCCCGGGCACCCCACCACGGGGTCATGCTTTTTACTGGCCTGCTCCCCTAACGGGAACCGGGGAGCATCatacacaaagccgccccccagcacaat                    >  1:503924/1‑277 (MQ=255)
                                                                                                                                                                                                                                                                               cccTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTTCTCACCCTTCCCCGATGAGTGCAAACTTTCACAATCTCACCGAAGTTACCACATCGCTGTGGTGAATTATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCCGATTTACAGTCTGCTCCCTTTGGCCGCTCGGGAACCCCACCACGGGGTAATGCCTTTTACTGCACTGCCACCTTATCAGAAAACCGGGCCCaacaaaacaacaaacgacccgcattaacat                    >  1:247300/1‑281 (MQ=37)
                                                                                                                                                                                                                                                                               cccTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTTCTCACCCTTCCCCGATGAGTGCAAACTTTCACAATCTCACCGAAGTTACCACATCGCTGTGGTGAATTATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTACAGTCTGCTCCCTTTTGCCGCTCGGAAACCCCACCACCGGGTAATGCTTTTTACTCGCCTGCTCCCTTATCGGAAAGCGGGGCGCATCATACCAAATGACTCACCACTGtccagt                    >  1:574209/1‑296 (MQ=38)
                                                                                                                                                                                                                                                                                   tCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTTCTCACCCTTCCCCGATGAGTGCAAACTTTCACAATCTCACCGAAGTTACCACATCGCTGTGGTGAATTATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTACAGTCTGCTCCCTTTGGCCGCTCGGGAACCCCACCACGGGGTAATGCTTTTTACTGGCCTGCTCCCTTATCGGGAAGCGGGGCGCATCATATCAAATGACGCGCCGCTGTAAAGTGTTa                <  2:317870/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                                     agaacGTTGCCACAGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTTCTCACCCTTCCCCGATGAGTGCAAACTTTCACAATCTCACCGAAGTTACCACATCGCTGTGGTGAATTATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTACAGTCTGCTCCCTTTGGCCGCTCGGGAACCCCACCACGGGGTAATGCTTTTTACTGGCCTGCTCCCTTATCGGGAAGCGGGGCGCATCATATCAAATGACGCGCCGCTGTAAAGTGTTACg              >  2:392497/5‑301 (MQ=255)
                                                                                                                                                                                                                                                                                      ggCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTTCTCACCCTTCCCCGATGAGTGCAAACTTTCACAATCTCACCGAAGTTACCACATCGCTGTGGTGAATTATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTACAGTCTGCTCCCTTTGGCCGCTCGGGAACCCCACCACGGGGTAATGCTTTTTACTGGCCTGCTCCCTTATCGGGAAGCGGGGCGCATCATATCAAATGACGCGCCGCTGTAAAGTGTTCCGt             >  1:179919/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                                        acGTTGCCACAGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTTCTCACCCTTCCCCGATGAGTGCAAACTTTCACAATCTCACCGAAGTTACCACATCGCTGTGGTGAATTATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTACAGTCTGCTCCCTTTGGCCGCTCGGGAACCCCACCACGGGGTAATGCTTTTTACTGGCCTGCTCCCTTATCGGGAAGCGGGGCGCATCATATCAAATGACGCGCCGCTGTAAAGTGTTACGTTg           <  1:391983/300‑1 (MQ=255)
                                                                                                                                                                                                                                                                                                 gtgGCTACGTTCTTTTGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTTCTCACCCTTCCCCGATGAGTGCAAACTTTCACAATCTCACCGAAGTTACCACATCGCTGTGGTGAATTATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTACAGTCTGCTCCCTTTGGCCGCTCGGGAACCCCACCACGGGGTAATGCTTTTTACTGGCCTGCTCCCTTATCGGGAAGCGGGGCGCATCATATCAAATGACGCGCCGCTGTAAAGTGTTACGTTGAGAAAAATg  <  1:375497/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                                                                            
GCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTACAGTCTGCTCCCTTTGGCCGCTCGGGAACCCCACCAGGGGTAATTCAAATTTTGAGGTAATGCTTGAGATGGTGGTGGGGGAAGGATTATTCGTCGCTTCGCTCCTCACCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTTCTCACCCTTCCCCGATGAGTGCAAACTTTCACAATCTCACCGAAGTTACCACATCGCTGTGGTGAATTATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTACAGTCTGCTCCCTTTGGCCGCTCGGGAACCCCACCACGGGGTAATGCTTTTTACTGGCCTGCTCCCTTATCGGGAAGCGGGGCGCATCATATCAAATGACGCGCCGCTGTAAAGTGTTACGTTGAGAAAAATG  >  NC_000913/1287131‑1287720

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