Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,816,036 2 bp→TT coding (513‑514/687 nt) rph ← truncated RNase PH

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,816,0360CT98.5% 234.5 / NA 67V172I (GTT→ATT) rphtruncated RNase PH
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (38/28);  total (38/29)
Fisher's exact test for biased strand distribution p-value = 4.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
*NC_0009133,816,0370CT98.5% 238.1 / NA 67A171A (GCG→GCArphtruncated RNase PH
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (37/28);  minor base G (0/1);  total (38/29)
Fisher's exact test for biased strand distribution p-value = 4.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

GACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTACCAGCGCCACGCAGGCACCCGTAATCGACGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCCGCGCGAAGAGCACGGGCGATCAGACGCTGGAT  >  NC_000913/3815798‑3816282
                                                                                                                                                                                                                                                ||                                                                                                                                                                                                                                                     
gACAAAAAAAAGGCGACTCATCAGTCGACTTAAAAATCAGTTTGCCCGCGCCGCCTTCTGCGTCGCTACAATGGATGCGATTCCCCCTCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTTAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGt                                                                                                                                                                                                                                             <  1:266748/251‑1 (MQ=255)
  cAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTTAGTCGCTACAATGGATTCGATTCCCCCCCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGGGAACGGCTCGCCTTCTGCCGTCCACTGCACTTCAATGATGCGCCCGTCTTCGTTCATCACTACGTTCCTGTCGGTCTCTGAGGCTGAGTCCTCAACGGATTCCAGATCGCAAATTGCTTCACCGttt                                                                                                                                                                                                                                           >  2:345905/1‑250 (MQ=255)
          aGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCCTCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAAtacc                                                                                                                                                                                                                                   >  1:42203/1‑248 (MQ=255)
           ggCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCCTCGGGCCAGCGCCAACAAGGGGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGCTCGCAAATTGCTTCACCGTTTACAAtaccc                                                                                                                                                                                                                                  >  2:324339/1‑247 (MQ=255)
               aCTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCCTCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGAc                                                                                                                                                                                                                                >  2:122259/1‑249 (MQ=255)
                  catcaGTCGCCTTAAAAATCCGTTTGCCAGCGCCGCCTTCTGCGTCTCTACAATGGATTCGATTCCCCCTCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATTGGTGACCGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTTCTTCACCGTTTACAATACCGACAGaa                                                                                                                                                                                                                            <  1:128244/250‑1 (MQ=255)
                        tCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCCTCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCgg                                                                                                                                                                                                                     >  1:48131/1‑251 (MQ=255)
                            cTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCCTCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTAc                                                                                                                                                                                                                 >  2:314612/1‑251 (MQ=255)
                              tAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCCTCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTAcc                                                                                                                                                                                                                <  1:356191/250‑1 (MQ=255)
                                          tGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCCTCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTACCATCCCTTTCATCg                                                                                                                                                                                                   >  2:106038/1‑251 (MQ=255)
                                            ccAGCGCCGCCTTCTGCGACGCTACAATGGATTCGCTTCCCCCTCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGCCTTCGGTCATCACTACGTTCCTGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAACACCGACAGAAACTGCGGCTACCATCCCTTTCATCggg                                                                                                                                                                                                 >  2:386509/1‑250 (MQ=255)
                                                 gccgccTTCTGCGTCGCTACAATGGATTCGATTCCCCCTCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTACCATCCCTTTCATCGGGTTGGt                                                                                                                                                                                            >  1:333844/1‑251 (MQ=255)
                                                 gccgccTTCTGCGTCGCTACAATGGATTCGATTCCCCCTCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTACCATCCCTTTCATCGGGTTGGt                                                                                                                                                                                            >  1:293302/1‑251 (MQ=255)
                                                    gccTTCTGCGTCGCTACAATGGATTCGATTCCCCCTCGGGCCAGCGCCAACAAGGTGAGTCGCTATTCATGGGTGAACGGCTCGCCTTCAGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCCTCACTTCGTTCCTGTCGGTCTCTGCGGCTGGGTCTTCAACGTATTCCAGAGCGCCAATAGCTTCACCGGTTACAATACCGACAGAAACTGCGGACACCATCCCTTTCaccaggttgggttt                                                                                                                                                                                         >  2:337481/1‑238 (MQ=255)
                                                                  tACAATGGATTCGATTCCCCCTCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTACCATCCCTTTCATCGGGTTGGTTTTCAGCTTGCCGtttt                                                                                                                                                                           >  2:128501/1‑251 (MQ=255)
                                                                       tGGATTCGATTCCCCCTCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTACCATCCCTTTCATCGGGTTGGTTTTCAGCTTGCCGTTTTccacc                                                                                                                                                                      >  1:142970/1‑251 (MQ=255)
                                                                            tcgattcCCCCTCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTACCATCCCTTTCATCGGGTTGGTTTTCAGCTTGCCGTTTTCCACCAGCt                                                                                                                                                                  <  2:89455/250‑1 (MQ=255)
                                                                              attcccccTCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATGCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTACCATCCCTTTCATCGGGTTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTg                                                                                                                                                              <  1:151864/247‑1 (MQ=255)
                                                                                             ccAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTACCATCCCTTTCATCGGGTTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCtgcc                                                                                                                                                 <  1:192165/251‑4 (MQ=255)
                                                                                              cagcgCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTACCATCCCTTTCATCGGGTTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCa                                                                                                                                                <  1:436470/247‑1 (MQ=255)
                                                                                                         aggtGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTACCATCCCTTTCATCGGGTTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACGc                                                                                                                                       >  1:242193/4‑249 (MQ=255)
                                                                                                              ggttgCTCTTGAAGGGTGAACGGCCTGCCTTCTGCACTCCCCTGCCCGTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTACCATCCCTTTCATCGGGTTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACGCAGGCAcc                                                                                                                                <  2:75426/247‑1 (MQ=255)
                                                                                                              aGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTCCCATCCCTTTCATCGGGTTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACGCAGGCAcc                                                                                                                                <  1:439689/251‑1 (MQ=255)
                                                                                                                     ctTCGTGGGTGCACGGCTCGCCTTCTGCCGCCCCCTGCACTTCAATGATGCGCCGGGCTTCGGTCACCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTAGTCCAGATCGCAAATGGCTTCACCGTTTACAATACCGACAGAAACTGCGTCTACCATCCCTTTCATCGGGTTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACGCAGGCACCCGTaa                                                                                                                           <  2:403232/249‑1 (MQ=255)
                                                                                                                      ttCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTACCATCCCTTTCATCGGGTTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACGCAGGCACCCGTAATCg                                                                                                                        >  1:457362/1‑251 (MQ=255)
                                                                                                                         aTGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTACCATCCCTTTCATCGGGTTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACGCAGGCACCCGTAATCGAcg                                                                                                                     >  1:144360/1‑251 (MQ=255)
                                                                                                                           gggTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTACCATCCCTTTCATCGGGTTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACGCAGGCACCCGTAATCGAcgcg                                                                                                                   >  1:371356/1‑251 (MQ=255)
                                                                                                                               gAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTACCATCCCTTTCATCGGGTTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACGCAGGCACCCGTAATCGACGCGGTAc                                                                                                               >  1:75060/1‑251 (MQ=255)
                                                                                                                                aaCGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACGACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTACCATCCCTTTCATCGGGTTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACGCAGGCACCCGTAATCGACGCGGTAcg                                                                                                              <  1:197792/251‑1 (MQ=255)
                                                                                                                                   ggCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTACTTCACCGTTTCCAATACCGACAGAAACTGCGGCTACCATCCCTTTCATCGGGTTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCCCGCAGGCCCCCGTAATCGACGCGGTACTCGt                                                                                                           >  2:136320/1‑251 (MQ=255)
                                                                                                                                    gCTCGCCTTCTGCCGACCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGCGTCTTCAACGTATTCCAGATCGCAAATTTCTTCACCGTTTACAATACCGACAGAAACTGCGGCTACCATCACTTTCATCGGGTTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGAGCATCTGCCAGCGCCACGCAGGCACCCGTAATCGCCGCGGAACGCGTg                                                                                                          >  2:280443/1‑251 (MQ=255)
                                                                                                                                       cGCCTTCTGCCGTCCCCTGCACTTCAATGANGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTACCATCCCTTTCATCGGGTTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACGCAGGCACCCGTAATCGACGCGGTACGCGTGcc                                                                                                        <  2:209214/250‑1 (MQ=255)
                                                                                                                                             cTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTACCATCCCTTTCATCGGGTTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACGCAGGCACCCGTAATCGACGCGGTACGCGTGCCACCATc                                                                                                  >  1:213509/1‑250 (MQ=255)
                                                                                                                                                 cGTCCCCTGCACTTCAATGATGCGCCGGTCTACGGTCATCACTACGTTGTGGTATGCCCCTGCGGCTCAGTCTTCAACGTATTCCAGATCGCAAATTGCTACACCGTTTACAATTCAGACAGAAACTGCGGCTCCCATCCCTTTCATCGGGTTGGTTTTCAGCTTGCCGTTTTCCACCAGCTGCTGTAGCGCATCTGCCAGCGCCACGCAGGCACCCGTAATCGACGCGGTACGCGTGCCACCATCAGCCt                                                                                             <  2:253718/251‑1 (MQ=255)
                                                                                                                                                   tCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTACCATCCCTTTCATCGGGTTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACGCAGGCACCCGTAATCGACGCGGTACGCGTGCCACCATCAGCCTg                                                                                            >  2:337636/1‑250 (MQ=255)
                                                                                                                                                    ccccTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATGCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGGTTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTg                                                                                                                                                              <  1:302447/183‑1 (MQ=255)
                                                                                                                                                    ccccTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTACCATCCCTTTCATCGGGTTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTg                                                                                                                                                              >  2:302447/1‑183 (MQ=255)
                                                                                                                                                    ccccTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTACCATCCCTTTCATCGGGTTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACGCAGGCACCCGTAATCGACGCGGTACGCGTGCCACCATCAGCCTGa                                                                                           >  1:414840/1‑250 (MQ=255)
                                                                                                                                                         gCACTTCAATGATGCGCCCTTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTACCAGCCCTTTCATCGGGTTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCTCATCTGCCAGCGCCACGCAGGCACCCGTAATCGACGCGGTACGCGTGCCACCATCAGCCTGAAGCa                                                                                       <  2:405135/249‑1 (MQ=255)
                                                                                                                                                          cccTTCAATGATGCGCCCGTCTTTGGTCATCACTACGTTCATGTCTGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTACACCGTTTACAATTCCGACAGAAACTGCGTCTACCATCCCTTTCATCGGGTTGGTTTTCAGCTTGCCGTTTTCCACCAGCTGCTGTAGCGCATCTGCCAGCGCCACGCATGCACCCGTAATCGACGCGGGACGCTTGCCACCATCAGCCTGAAGCACGt                                                                                    <  2:31129/249‑1 (MQ=255)
                                                                                                                                                               cAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTACCATCCCTTTCATCGGGTTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACGCAGGCACCCGTAATCGACGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAg                                                                               <  2:42203/251‑1 (MQ=255)
                                                                                                                                                                           cGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTACCATCCCTTTCATCGGGTTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACGCAGGCACCCGTAATCGACGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATg                                                                   <  1:392940/251‑1 (MQ=255)
                                                                                                                                                                            gTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCCCCGTTTACAATACCGACAGAAACTGCGGCTACCATCCCTTTCATCGGGTTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACGCAGGCACCCGTAATCGACGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATgg                                                                  <  1:102178/251‑1 (MQ=255)
                                                                                                                                                                            gTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTACCATCCCTTTCATCGGGTTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACGCAGGCACCCGTAATCGACGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATgg                                                                  >  1:166733/1‑251 (MQ=255)
                                                                                                                                                                            gTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTACCATCCCTTTCATCGGGTTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACGCAGGCACCCGTAATCGACGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATgg                                                                  <  1:122259/251‑1 (MQ=255)
                                                                                                                                                                               ttCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAg                                                                                                                                                                                                                              <  1:385859/92‑1 (MQ=255)
                                                                                                                                                                               ttCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTACCATCCCTTTCATCGGGTTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACGCAGGCACCCGTAATCGACGCGGTACGCGTGCCACCATCAGCNTGAAGCACGTCGCAGTCCAGCGTAATGGTg                                                                >  2:231849/1‑250 (MQ=255)
                                                                                                                                                                                 cGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTACCATCCCTTTCATCGGGTTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACGCAGGCACCCGTAATCGACGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAAc                                                             >  2:177701/1‑251 (MQ=255)
                                                                                                                                                                                  ggTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTACCATCCCTTTCATCGGGTTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACGCAGGCACCCGTAATCGACGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACt                                                            >  1:137371/1‑251 (MQ=255)
                                                                                                                                                                                     catcaCTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTACCATCCCTTTCATCGGGTTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACGCAGGCACCCGTAATCGACGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCa                                                          <  1:324415/250‑1 (MQ=255)
                                                                                                                                                                                        caCTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTACCATCCCTTTCATCGGGATGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACGCAGGCACCCGTAATCGACGCGGTACGCGTGCCAACATCAGCCTGAAGCACGACGCAGTCCAGCGTAATGGTGAACTCcccc                                                       >  1:333762/1‑250 (MQ=255)
                                                                                                                                                                                                 cATGTCGGTATCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCCCCGTTTACAATACCGACAGAAACTGCGGCTACCATCCCTTTCATCGGGTTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACGCAGGCACCCGTAATCGACGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTCATGGTGAACTCACCCAGCGCTTTcc                                             >  1:72493/1‑250 (MQ=255)
                                                                                                                                                                                                  aTGTCGGTCTCTGCGGCTGAGTCTTCAACGTCTTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGCCAGAAACTGCGGCTACCTTCCCTTTCCTCTGGTTGGTTTTCAGCTTGCCGTTTTCCCCCAGCTTCTGTAGCGCATCTGCCAGCGCCCCGCGGGCACCCGTAATCGACGCGGAACGCGTGCCACCATCAGCCTGACGCACGTCGCAGTCCAGCGTAATGGTGAACTCACCCAGCTCttcca                                             >  2:270489/1‑247 (MQ=255)
                                                                                                                                                                                                   tGTCGGTCTCTGCGGCTGAGACTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTACCAGCCCTTTCATCGGGTTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAACGCCACGCAGGCACCCGTAATCGACGCGGTACGCGTGCCACCATCCGCCTGAAGAACGTCGCAGTCCCGCGTAATGGTGAACTCACCCAGCGCTGTCaac                                           >  2:37734/1‑250 (MQ=255)
                                                                                                                                                                                                        gTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTACCATCCCTTTCATCGGGTTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACGCAGGCACCCGTAATCGACGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTAc                                      >  2:191466/1‑251 (MQ=255)
                                                                                                                                                                                                            ctGCGGCTGAGTCTTCACCGGATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACATAAACTGCGGCTACCATCCCTTTCATCGGGTTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACGCAGGCACCCGTAATCGACGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTAcc                                     <  2:101675/248‑2 (MQ=255)
                                                                                                                                                                                                              gCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTACCATCCCTTTCATCGGGTTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACGCAGGCACCCGTAATCGACGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCcgc                                >  1:13992/1‑251 (MQ=255)
                                                                                                                                                                                                                       tCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTACCATCCCTTTCATCGGGTTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACGCAGGCACCCGTAATCGACGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCCGCGCGAAGAGc                       >  2:34213/1‑251 (MQ=255)
                                                                                                                                                                                                                          tCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTACCATCCCTTTCATCGGGTTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACGCAGGCACCCGTAATCGACGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCCGCGCGAAGAGCa                      <  1:128501/249‑1 (MQ=255)
                                                                                                                                                                                                                            aaCGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTACCATCCCTTTCATCGGGTTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACGCAGGCACCCGTAATCGACGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTAc                                      <  2:121711/231‑1 (MQ=255)
                                                                                                                                                                                                                            aaCGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTACCATCCCTTTCATCGGGTTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACGCAGGCACCCGTAATCGACGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTAc                                      >  1:121711/1‑231 (MQ=255)
                                                                                                                                                                                                                                   tCCATATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTACCATCCCTTTCATCGGGTTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACGCAGGCACCCGTAATCGACGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCCGCGCGAAGAGCACGGGCGATCa            <  1:314612/250‑1 (MQ=255)
                                                                                                                                                                                                                                   tCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTACCATCCCTTTCATCGGGTTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACGCAGGCACCCGTAATCGACGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCCGCGCGAAGAGCACGGGCGATCa            <  1:5418/250‑1 (MQ=255)
                                                                                                                                                                                                                                     cAGATCTCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTCCCATCCCTTTCATCGGGTTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACGCAGGCACCCGTAATCGACGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCTTGCAAATCTACTGCCGCGCGAAGAGCACGGGCGATCAGAc         <  2:142970/251‑1 (MQ=255)
                                                                                                                                                                                                                                        aTCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTACCATCCCTTTCATCGGGTTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACGCAGGCACCCGTAATCGACGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTAACCCAGCGCTTTCAAATCTACTGCCGCGCGAAGAGCACGGGCGATCAGACGCt      >  2:392711/1‑251 (MQ=255)
                                                                                                                                                                                                                                        aTCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGAGGCTACCATCCCTTTCATCGGGTTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACGCAGGCACCCGTAATCGACGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCACATCTACTGCCGCGCGAAGAGCACGGGCGATCAGACGCt      >  2:270009/1‑251 (MQ=255)
                                                                                                                                                                                                                                            cAAACTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTACCATCCCTTGCATCGGGTTGGTTTTCATCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACGCAGGCACCCGTAATCGACGCGGTACGCGTGCCACCATCAGCCTGAAGCACTTCGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAAGCTACTGCCGCGCGAAGAGCACGGGCGATCAGACGCTGGAt  <  2:48915/251‑1 (MQ=255)
                                                                                                                                                                                                                                                ||                                                                                                                                                                                                                                                     
GACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTACCAGCGCCACGCAGGCACCCGTAATCGACGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCCGCGCGAAGAGCACGGGCGATCAGACGCTGGAT  >  NC_000913/3815798‑3816282

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