Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 81,102 G→A S286L (TCG→TTG)  leuB ← 3‑isopropylmalate dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091381,1020GA97.8% 138.0 / NA 46S286L (TCG→TTG) leuB3‑isopropylmalate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base G (1/0);  major base A (26/18);  minor base C (1/0);  total (28/18)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.84e-01

ATTACACCCCTTCTGCTACATAGCGGGCAATGATATCGCCCATTTCATCGGTACTAACGGCGGCAGCGCCACGGGCTAAATCCCCGGTGCGAATGCCTTCTTCTAATGCGCGGTTAATGGCGCGTTCAATGGCGCAAGCCGCATCATCGGCATCCAGGCTGTAACGCAGCAGCAGTGCCAGCGAAAGGATTTGTGCAATCGGGTTGGCGATGTTTTTGCCTGCGATATCTGGTGCCGAGCCGCCCGCCGGTTCATACAGTCCAAAACCTTGCTCGTTCAGGCTGGCGGAAGGCAACATCCCCATCGAGCCAGTGATCATTGCGCACTCGTCAGACAGAATGTCGCCAAACAGGTTGGAGCACAGCAGAACGTCAAACTGTGATGGATCTTTAATCAGCTGCATGGTGGCGTTGTCGATGTACATATGCGCCAGTTCGACATCCGGGTATTCCGTGGCGATCTCGTT  >  NC_000913/80866‑81331
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                     
aTTACACCCCTTCTGCTACATAGCGGGCAATGATATCGCCCATTTCATCGGTACTAACGGCGGCAGCGCCCCGGGCTAAATCCCCGGTGCGAATGCCTTCTTCTAATGCGCGGTTAATGGCGCGTTCAATGGCGCAAGCCGCATCATCGGCATCCAGGCTGTAACGCAGCAGCAGTGCCAGCGAAAGGATTTGGGCAATCGGGTTGGCGATGTTTTTGCCTGCGGTATCTGGTGCCAAGCCGCCCGCCggg                                                                                                                                                                                                                         >  2:329947/1‑250 (MQ=255)
   acacCCCTTCTGCTACATAGCGGGCAATGATATCGCCCATTTCATCGGTACTAACGGCGGCAGCGCCACGGGCTAAATCCCCGGTGCGAATGCCTTCTTCTAATGCGCGGTTAATGGCGCGTTCAATGGCGCAAGCCGCATCATCGGCATCCAGGCTGTAACGCAGCAGCAGTGCCAGCGAAAGGATTTGTGCAATCGGGTTGGCGATGTTTTTGCCTGCGATATCTGGTGCCAAGCCGCCCGCCGGTTc                                                                                                                                                                                                                       >  1:220110/1‑250 (MQ=255)
           tCTGCTACATCGCGGGCAATGATATCGACCATTTCATCGGTACTAACGGCGGCAGCGCCACGGGCTAACTACCCGGTGCGAATGCCTTCTTCTAATGCGCGGTTAATGGCGCGTTCAATGGCGCAAGCCGCCTCATCGGCATCCAGGCTGTAACGCAGCAGCCGTGCCAGCGAAAGGATTTGGGCAATCGGGTTGGCGGTGTTTTTGCCTGCGATATCTGGTGCCAAGCCGCCCGCAGGGTCCTACAGTcc                                                                                                                                                                                                              >  2:279681/1‑251 (MQ=255)
           tCTGCTACATAGCGGGCAATGATATCGCCCATTTCATCGGTACTAACGTCGGCAGCGCCACGGGCTAAATCCCCGGTGCGAATGCCTTCTTCTAATGCGCGGTTAATGGCGCGTTCAATGGCGCAAGCCGCATCATCGGCATCCAGGCTGTAACGCAGCAGCAGTGCCAGCGAAAGGATTTGTGCAATCGGGTTGGCGATGTTTTTGCCTGCGATATCTGGTGCCAAGCCGCCCGCCGGTTCATACAGTcc                                                                                                                                                                                                              >  1:410067/1‑251 (MQ=255)
                          gCAATGATATCGCCCATTTCATCGGTACTAACGGCGGCAGCGCCACGGGCTAAATCCCCGGTGCGAATGCCTTCTTCTAATGCGCGGTTAATGGCGCGTTCAATGGCGCAAGCCGCATCATCGGCATCCAGGCTGTAACGCAGCAGCAGTGCCAGCGAAAGGATTTGTGCAATCGGGTTGGCGATGTTTTTGCCTGCGATATCTGGTGCCAAGCCGCCCGCCGGTTCATACAGTCCAAAACCTTGCTCGtt                                                                                                                                                                                               >  2:260425/1‑251 (MQ=255)
                                    gccccATTCCATCGGTACTACCGGCGGCAGCGCCCCGGGCTAAATCCCCGGTGCGAATGCCTTCTTCTAATGCGCGGTTAAGGGCGCGTTCAAGGGCGCAAGCCGCATCATCGGCATCCGGGGTGTAACGCAGCAGCAGTGCCAGCGAAAGGATTTGTGCAATCGGGTTGGCGATTTTTTTGCCTGCGATATCTGGTGCCAAGCCGCCCGCCGGTTCATACAGTCCAAAACCTTGCTCGTTCAGGCTGGc                                                                                                                                                                                      <  2:175228/248‑1 (MQ=255)
                                       ccATTTCATCGGTACTAACGGCGGCAGCGCCACGGGCTAAATCCCCGGTGCGAATGCCTTCTTCTAATGCGCGGTTAATGGCGCGTTCAATGGCGCAAGCCGCATCATCGGCATCCAGGCTGTAACGCCGCAGCCGTGCCAGCGAAAGGATTTGTGCCATCGGGGTGGCGATGGTTTTGCCTGCGATCTCTGGGGCCAAGCCGCCCCCCGGGTTATACAGTCCAAAACCTTGCTTGGTCAGGGTGGCGGa                                                                                                                                                                                   >  2:320542/1‑250 (MQ=255)
                                           ttCATCGGTACTAACGGCGGCAGCGCCACGGGCTAAATCCCCGGTGCGAATGCCTTCTTCTAATGCGCGGTTAATGGCGCGTTCAATGGCGCAAGCCGCATCATCGGCATCCAGGCTGTAACGCAGCAGCAGTGCCAGCGAAAGGATTTGTGCAATCGGGTTGGCGATGTTTTTGCCTGCGATATCTGGTGCCAAGCCGCCCGCCGGTTCATACAGTCCAAAACCTTGCTCGTTCAGGCTGGCGGAAGGc                                                                                                                                                                               >  1:96339/1‑250 (MQ=255)
                                             cATCGGTACTAACGGCGGCAGCGCCACGGGCTAAATCCCCGGTGCGAATGCCTTCTTCTAATGCGCGGTTAATGGCGCGTTCAATGGCGCAAGCCGCATCATCGGCATCCAGGCTGTAACGCAGCAGCAGTGCCAGCGAAAGGATTTGTGCAATCGGGTTGGCGATGTTTTTGCCTGCGATATCTGGTGCCCAGCCGCCCGCCGGTTCATACAGTCCAAAAACTTGCTCGTTCAGGCTGGCGGAAGGCAAc                                                                                                                                                                            >  2:10905/1‑251 (MQ=255)
                                                  gTACTAACGGCGGCAGCGCCACGGGCTAAATCCCCGGGGCGAATGCCTTCTTCTAATGCGCGGTTTAGGGCGCGTTCGATGGCGCAAGCCGCATCATCGGCATCCAGGCTGTAACGCAGCAGCCGTGCCAGCGAAATGCTTTGGGCAATCGGGCTGGCGATGTTTTTGCCTACGCTATCTGGGGCCACGCCTCCCGCCGGTTCATACCTTCCCAAACCGTGCTCTTTCCGGCGGGCCGACGGCAACCTccc                                                                                                                                                                       >  2:151564/1‑251 (MQ=255)
                                                  gTACTAACGCCGCCAGCGCCCCGGCCTAAACCCCGGGTGCGAAGGCCTTCTTCAAATGCGGGGTTAATGGCGCGTTAATTGGCGAAGGCCGCATCATCGGCCGCCAGGCTGTAACGCAGCAGCATGGCCAGCGAAAGGATTTGTGCAATCGGGTTGGCGTTGTTTTTGCCGGCGATATCTGGTGCCAAGCCGCCCGCCGGTTCATACAGTCCAAAACCTTGCTCGTTCAGGCTGGCGGAAGGCAACATccc                                                                                                                                                                       <  2:252958/251‑1 (MQ=255)
                                                                 gcgcCTCGGGTTAAACCCCCGGTGCGAATGTCTCTTTCTAATGCGCGTTTAATGGCGCGTCAAAGGGCGCAACCCGCATCATCGGCATACAGGCTGTAACGCAGCAGCAGTGCCAGCAAAAGGATTTGTGCAATCGGGTTGGCGATGTTTTTTCCTGCGTTATCTGGTGCCAAGCCGCCCGCCGGTTCATACAGTCCAAAACCTTGCTCGTTCAGGCTGGCGGAAGGCAACATCCCCATCGAGCCAGTGAt                                                                                                                                                        <  1:11918/251‑1 (MQ=255)
                                                                  cgcCACGGGCTAAATCCCCGGTGCGAATGCCGTTTTCTAATGCGCGTTTAAGGGCGCTTTCAATGGCGCAAGCCGCATCATCGGCATCCAGGCTGTACCGCAGCAGCAGTGCCGGCGAAAGGATTTGTGCAATCGGGTTGGCGATGTTTTTGCCTGCGATATCGGGTGCCAAGCCGCCCGCCGGTTCATACAGTCCAAAACCTTGCTCGTTCAGGCTGGCGGAAGGCAACATCCCCATCGAGCCAGTGATc                                                                                                                                                       <  1:428402/251‑1 (MQ=255)
                                                                   gcCACGGACTAAATCCCCGGTGCGAATGCCTTCTTCTAATGCGCGGTTAATAGCGCGTTCAATGGCGCAAGCCGCATCATCGGCATCCAGGCTGTAACGCAGCAGCAGTGCCAGCGAAAGGATTTGTGCAATCGGGTTGGCGATGTTTTTGCCTGCGATATCTGGTGCCAAGCCGCCCGCCGGTTCATACAGTCCAAAACCTTGCTCGTTCAGGCTGGCGGAAGGCAACATCACCATCGAGCCAGTGATc                                                                                                                                                       >  1:175960/1‑250 (MQ=255)
                                                                     cACGGGCTAAATCCCCGGTGCGAATGCCTTCTTCTAATGCGCGGTTAATGGCGCGTTCAATGGCGCAAGCCGCATCATCGGCATCCAGGCTGTAACGCAGCAGCAGTGCCAGCGAAAGGATTTGTGCAATCGGGTTGGCGATGTTTTTGCCTGCGATATCTGGTGCCAAGCCGCCCGCCGGTTCATACAGTCCAAAACCTTGCTCGTTCAGGCTGGCGGAAGGCAACATCCCCATCGAGCCAGTGATCAtt                                                                                                                                                    >  2:28128/1‑251 (MQ=255)
                                                                                cccccGGTGGGAATGCCTTCTTCTAATGCGCGGTTAAGGGCGCGTTCAAGGGCGAAAGCCGCATCATCGGCATCCAGGCTGTAACGCAGCAGCAGTGCCAGCGAAAGGATTTGTGCAATCGGGTTGGCGATTTTTTTGCCTGCGATATCTGGTGCCAAGCCGCCCGCCGGTTCATACAGTCCAAAACCTTGCTCGTTCAGGCTGGCGGAAGGCAACATCCCCATCGAGCCAGTGATCATTGCGCACTCGTc                                                                                                                                         <  1:159829/250‑1 (MQ=255)
                                                                                     ggTGCGAATGCCTTCTTCTAATGCGCGGTTAATGGCGAGTTCAATGGCGCAAGCCGCATCATCGGCATCCAGGCTGTAACGCAGCAGCAGTGCCAGCGAAAGGATTTGTGCAATCGGGTTGGCGATGTTTTTGCCTGCGATATCTGGTGCCAAGCCGCCCGACGGTTCATACAGTCCAAAACCTTGCTCGTTCAGGCTGGCGGAAGGCAACATCCCCATCGAGCCAGTGATCCTTGCGCACTCGTcagac                                                                                                                                     >  1:384959/1‑250 (MQ=255)
                                                                                      gTGCGAATGCCTTCTTCTAATGCGCGGTTAATGGCGCGTTCAATGGCGCAAGCCGCATCATCGGCATCCAGGCTGTAACGCAGCAGCAGTGCCAGCGAAAGGATTTGTGCAATCGGGTTGGCGATGTTTTTGCCTGCGATATCTGGTGCCAAGCCGCCCGCCGGTTCATACAGTCCAAAACCTTGCTCGTTCAGGCTGGCGGAAGGCAACATCCCCATCGAGCCAGTGATCATTGCGCACTCGTcagac                                                                                                                                     >  2:340955/1‑249 (MQ=255)
                                                                                                    ttttaaggcgcGTTTAAGGGCGCGTTGACTGCGCAAAGCCGCATCATGGGCTCCCAGGGTGTAACGCAGCAGTGTGGCGAGCAAAAGTATTGGTGCAACCGGTGTGGCTTTGTTTTTGCCTGCGATTTCGGGTGCCACCCCGCCCGCCGGTTCATACAGTCCAAAACCTTGCTCGTTCAGGCTGGCGGAAGGCAACATCCCCATCGAGCCAGTGATCATTGCGCACTCGTCAGACAGAATGTCGCCAAAc                                                                                                                      <  2:4367/243‑1 (MQ=255)
                                                                                                     tctAATGCGCGGTTAATGGCGCGTTCAATGGCGCAAGCCGCATCATCGGCATCCAGGCTGTAACGCAGCAGCAGTGCCAGCGAAAGGATTTGTGCAATCGGGTTGGCGATGTTTTTGCCTGCGATATCTGGTGCCAAGCCGCCCGCCGGTTCATACAGTCCAAAACCTTGCTCGTTCAGGCTGGCGGAAGGCAACATCCCCATCGAGCCAGTGATCATTGCGCACTCGTCAGACAGAATGTCGCCAAAc                                                                                                                      >  1:72068/1‑249 (MQ=255)
                                                                                                         aggcgcGGTTAAGGGCGCGTTAAAGGGCGCAACCCGCATCATCGGCATCCAGGCTGTAACGCAGCAGCAGTGCCAGCGAAAGGTTTTGTGCAATCGGGTTGGCGATGTTTTTGCCTGCGATATCTGGTGCCAAGCCGCCCGCCGGTTCATACAGTCCAAAACCTTGCTCGTTCAGGCTGGCGGAAGGCAACATCCCCATCGAGCCAGTGATCATTGCGCACTCGTCAGACAGAATGTCGCCAAACAGGTTg                                                                                                                <  2:47769/249‑1 (MQ=255)
                                                                                                                gTTAATGGCGCGTTCAATGGCGCAAGCCGCATCATCGGCATCCAGGCTGTAACGCAGCAGCAGTGCCAGCGAAAGGATTTGTGCAATCGGGATGGCGATGTTTTTGCCTGCGATATCTGGTGCCAAGCCGCCCGCCGGTTCATACAGTCCAAAACCTTGCTCGTTCAGGCTGGCGGAAGGCAACATCCCCATCGAGCCAGTGATCATTGCGCACTCGTCAGACAGAATGTCGCCAAACAGGTTGGAGcac                                                                                                          >  1:38925/1‑250 (MQ=255)
                                                                                                                 ttAATGGCGCGTTCAATGGCGCAAGCCGCATCATCGGCCTACAGGCTGTAACGCAGAAGCAGTGCCAGCGAAAGGATTTGTGCACTCGGGTTGGCGATGTTTTTGCCTGCGATATCTGGTGCCAAGCCGCCCGCCGGTTCATACAGTCCAAAACCTTGCTCGTTCAGGCTGGCGGAAGGCAACATCCCCATCGAGCCAGTGATCATTGCGCACTCGTACGACAGAAAGTCGCCAAACAGGGTggcgaac                                                                                                          >  2:436824/1‑244 (MQ=255)
                                                                                                                    aagggcgcgtcaatgggcgaaggCCGCATCAGCGTCTTCCAGGCGGTAACGCAGCGGCAGGGCCGGCGAAGGGATTTGTGCAACCGGGTTGGCGATGTTTTTGCCTGCGATATCTGGTGCCAAGCCGCCCGCCGGTTCATACAGGCCAAAACCTTGCTCGTTCAGGCTGGCGGAAGGCAACATCCCCATCGAGCCAGTGATCATTGCGCACTCGTCAGACAGAATGTCGCCAAACAGGTTGGAGCAcagca                                                                                                      <  2:410067/229‑1 (MQ=255)
                                                                                                                    aTGGCGCGTTCAATGGCGCAAGCCGCATCATCGGCATCCAGGCTGTAACGCAGCAGCAGTGCCAGCGAAAGGATTTGTGCAATCGGGTTGGCGATGTTTTTGCCTGCGATATCTGGTGCCAAGCCGCCCGCCGGTTCATACAGTCCAAAACCTTGCTCGTTCAGGCTGGCGGAAGGCAACATCCCCATCGAGCCAGTGATCATTGCGCACTCGTCAGACAGAATGTCGCCAAACAGGTTGGAGCGcagcag                                                                                                     <  1:126857/251‑1 (MQ=255)
                                                                                                                               aaTGGCGCAAGCCGCATCATCGGCATCCAGGCTGTAACGCAGCAGCAGTGCCAGCGCAAGGATTTGTGCAATCGGGTTGGCGATGTTTTTGCCTGCGATATCTGGTGCCAAGCCGCCCGCCGGTTCATACAGTCCAAAACCTTGCTCGTTCAGGCTGGCGGAAGGCAACATCCCCATCGAGCCAGTGATCATTGCGCACTCGTCAGACAGAATGTCGCCAAACAGGTTGGAGCACAGCAGAACGTCAAACt                                                                                          >  1:385993/1‑251 (MQ=255)
                                                                                                                                  ggCGCAAGCCGCATCATCGGCATCCAGGCTGTAACGCAGCAGCAGTGCCAGCGAAAGGATTTGTGCAATCGGGTTGGCGATGTTTTTGCCTGCGATATCTGGTGCCAAGCCGCCCGCCGGTTCATACAGTCCAAAACCTTGCTCGTTCAGGCTGGCGGAAGGCAACATCCCCATCGAGCCAGTGATCATTGCGCACTCGTCAGACAGAATGTCGCCAAACAGGTTGGAGCACAGCAGAACGTCAAACtgtg                                                                                       >  1:94924/1‑251 (MQ=255)
                                                                                                                                       aaGCCGCATCATCGGCATCCAGGCTGTAACGAAGCAGCAGTACAAGCGACAGGCTTTGTGCAATCGGGTTGGCGATGTTTTTGCCTGCGATATCTGGTGCCGAGCAGCACGCCGGCTAATAAAGTCCAAAACCTTGATCGTTCAGGCTGGCGGAAGGAAACATCCCCATCGCGCTAGTGATCACTGCGCACTCGTTAGACAGAATGGCGCCAAACAGGGTGGAGCACAGCAGAAc                                                                                                  >  1:449433/1‑235 (MQ=255)
                                                                                                                                                    cgccATCCAGGCTGTACCGATGCAGCAGGGCCAGCGAAAGGTTTTGTGCAATGGGGTTGGCGATGTTTTTGCCTGCGTTTTCTGTTGCCAAGCGGCCCGCCGGTTCATACAGTCCAAAACCTTGCTCGTTCAGGCTGGCGGAAGGCAACATCCCCATCGAGCCAGTGATCATTGCGCACTCGTCAGACAGAATGTCGCCAAACAGGTTGGAGCACAGCAGAACGTCAAACTGTGATGGATCTTTAATCAGc                                                                      <  1:28128/248‑1 (MQ=255)
                                                                                                                                                                 gTAACGCAGCAGCAGTGCCAGCGAAAGGATTTGTGCAATCGGGTTGGCGATGTTTTTGCCTGCGATATCTGGTGCCAAGCCGCCCGCCGGTTCATACAGTCCAAAACCTTGCTCGTTCAGGCTGGCGGAAGGCAACATCCCCATCGAGCCAGTGATCATTGCGCACTCGTCAGACAGAATGTCGCCAAACCGGTTGAAGCACAGCAGAACGTCCAACTGTGATGGATCTTTAATCAGCTGCATGGTGGCGt                                                         >  2:50946/1‑251 (MQ=255)
                                                                                                                                                                   aaCGCAGCAGCAGTGCCAGCGAAAGGCTTTGTGCAATCGGGTTGGCGATGTTTTTGCCTGCGATATCTGGTGCCAAGCCGCCCGCCGGTTCATACAGTCCAAAACCTTGCTCGTTCAGGCTGGCGGAAGGCAACATCCCCATCGAGCCAGTGATCATTGCGCACTCGTCAGACAGAATGTCGCCAAACAGGTTGGAGCACAGCAGAACGTCAAACTGTGATGGATCTTTAATCAGCTGCATGGTGGCGTTg                                                       >  2:430223/1‑251 (MQ=255)
                                                                                                                                                                     ccgctgcgtggtgcccAGGGAAAGGTTTTGTGCAAGCGGGTGGGCGATTTTTGTGCCTGCGATAGCTGGTGCCAAGCCGCCCGCCGGTTCATTCAGTCCAAACCCTTGCTCGTTCAGGCTGGCGGAAGGCAACATCCCCATCGAGCCAGTGATCATTGCGCACTCGTCAGACAGAATGTCGCCAAACAGGTTGGAGCACAGCAGAACGTCAAACTGTGATGGAACTTTAATCAGCTGCATGGTGGCGTTGt                                                      <  2:29875/237‑1 (MQ=255)
                                                                                                                                                                     cGCAGCAGCAGTGCCAGCGAAAGGATTTGTGCAATCGGGTTGGCGATGTTTTTGCCTGCGATATCTGGTGCCAAGCCGCCCGCCGGTTCATACAGTCCAAAACCTTGCTCGTTCAGGCTGGCGGAAGGCAACATCCCCATCGAGCCAGTGATCATTGCGCACTCGTCAGACAGAATGTCGCCAAACAGGTTGGAGCACAGCAGAACGTCAAACTGTGATGGATCTTTAATCAGCTGCATGGTGGCGTTGTc                                                     >  1:185978/1‑251 (MQ=255)
                                                                                                                                                                     cGCAGCAGCAGTGCCAGCGAAAGGATTTGTGAAAGGGGGTTGGCGATTGTTTGGCTCGCGATATCGGGTGCCAAGCCGCCCGCCGGTTCATACAGTCCAAAACCTTGCTCGTTCAGGCTGGCGGAAGGCAACATCCCCATCGAGCCAGTGATCATTGCGCACTCGTCAGACAGAATGTCGCCAAACAGGTTGGAGCACAGCAGAACGTCAAACTGTGATGGATCTTTAATCAGCTGCATGGTGGCGTTGTc                                                     <  1:332921/251‑1 (MQ=255)
                                                                                                                                                                     cGCAGCAGCAGGGCCAGCGAAAGGATTTGTGCAATCGGGTTGGCGATGTTTTTGCCTGCGATATCTGGTGCCAAGCCGCCCGCCGGTTCATACAGTCCAAAACCTTGCTCGTTCAGGCTGGCGGAAGGCAACATCCCCATCGAGCCAGTGATCATTGCGCACTCGTCAGACAGAATGTCGCCAAACAGGTTGGAGCACAGCAGAACGTCAAACTGTGATGGATCTTTAATCAGCTGCATGGTGGCGTTGTc                                                     <  1:279681/251‑1 (MQ=255)
                                                                                                                                                                         gcagcagTGCCAGCGAAAGGATTTGTGCAATCGGGTTGGCGATGTTTTTGCCTGCGATATCTGTTGCCAAGCCGCCCGCCGGTTCATACAGTCCAAAACCTTGCTCGTTCAGGCTGGCGGAAGGCAACATCCCCATCGAGCCAGTGATCATTGCGCACTCGTCAGACAGAATGTCGCCAAACAGGTTGGAGCACAGCAGAACGTCAAACTGTGATGGATCTTTAATCAGCTGCATGGTGGCGTTGTCGATg                                                 <  1:66839/251‑1 (MQ=255)
                                                                                                                                                                          ccgctattgtcAGCGAAAGGGTTTGTGCAATCGGGTTTGCTATGTTTTTGCCTGCAATATGTGGTGCCACGCCTCCCGCCGTTTCATACAGTCCAAAACCTTGCTCGTTCAGGCTGGCGGAAGGCAACATCCCCATCGAGCCAGTGATCATTGCGCACTCGTCAGACAGAATGTCGCCAAACAGGTTGGAGCACAGCTGAACGTCAAACTGTGATGGATCTTTAATCAGCTGCATGGTTGCGTTGTCGAtt                                                 <  2:67793/241‑2 (MQ=255)
                                                                                                                                                                           agcagTGCCAGCGAAAGGATTTGTGCAATCGGGTTGGCGATGTTTTTGCCTGCGATATCTGGTGCCAAGCCGCCCGCCGGTTCATACAGTCCAAAACCTTGCTCGTTCAGGCTGGCGGAAGGCAACATCCCCATCGAGCCAGTGATCATTGCGCACTCGTCAGACAGAATGTCGCCAAACAGGTTGGAGCACAGCAGAACGTCAAACTGTGATGGATCTTTAATCAGCTGCATGGTGGCGTTGTCGATGt                                                >  2:76987/1‑250 (MQ=255)
                                                                                                                                                                                tGCCAGCGAAAGGATTTGTGCAATCGGGTTGGCGATGTTTTTGCCTGCGATATCTGGTGCCAAGCCGCCCGCCGGTTCATACAGTCCAAAACCTTGCTCGTTCAGGCTGGCGGAAGGCAACATCCCCATCGAGCCAGTGATCATTGCGCACTCGTCAGACAGAATGTCGCCAAACAGGTTGGAGCACAGCAGAACGTCAAACTGTGATGGATCTTTAATCAGCTGCATGGTGGCGTTGTCGATGTACatat                                          >  1:292326/1‑251 (MQ=255)
                                                                                                                                                                                      gCGAAAGGATTTGTGCAATCGGGTTGGCGATGTTTTTGCCTGCGATATCTGGTGCCAAGCCGCCCGCCGGTTCATACAGTCCAAAACCTTGCTCGTTCAGGCTGGCGGAAGGCAACATCCCCATCGAGCCAGTGTTCATTGCGCACTCGTCAGACAGAATGTCGCCAAACAGGTTGGAGCACCGCAGAACGTCAACCTGGGATGGGTCTTTAATCAGCTGCATGGTGGCGTTGTCGATGGACATATGCGcc                                     >  2:112624/1‑251 (MQ=255)
                                                                                                                                                                                      gCGAAAGGATTTGTGCAATCGGGTTGGCGAGGTTTTGGCCTGCGATTTCTGGTGCCAAGCCGCCCGCCGGTTCATACAGTCCAAAACCTTGCTCGTTCAGGCTGGCGGAAGGCAACATCCCCATCGAGCCAGTGATCATTGCGCACTCGTCAGACAGAATGTCGCCAAACAGGTTGGAGCACAGCAGAACGTCAAACTGTGATGGATCTTTAATCAGCTGCATGGTGGCGTTGTCGATGTACATATGCGcc                                     <  1:430223/251‑1 (MQ=255)
                                                                                                                                                                                          aaGGATTTGTGCAATCGGGTTGGCGATGTTTTTGCCTGCGATATCTGGTGCCAAGCCGCCCGCCGGTTCATACAGTCCAAAACCTTGCTCGTTCAGGCTGGCGGAAGGCAACATCCCCATCGAGCCAGTGATCATTGCGCACTCGTCAGACAGAATGTCGCCAAACAGGTTGGAGCACAGCAGAACGTCAAACTGTGATGGATCTTTAATCAGCTGCATGGTGGCGTTGTCGATGTACATATGCGCCAGtt                                 <  1:260425/251‑1 (MQ=255)
                                                                                                                                                                                             gATTTGTGCAATCGGGTTGGCGATGTTTTTGCCTGCGATATCTGGTGCCAAGCCGCCCGCCGGTTCATACAGTCCAAAACCTTGCTCGTTCAGGCTGGCGGAAGGCAACATCCCCATCGAGCCAGTGATCATTGCGCACTCGTCAGACAGAATGTCGCCAAACAGGTTGGAGCACAGCAGAACGTCAAACTGTGATGGATCTTTAATCAGCTGCATGGTGGCGTTGTCGATGTACATATGCGCCAGTTCg                               >  1:306700/1‑250 (MQ=255)
                                                                                                                                                                                                               ggCGATGTTTTTGCCTGCGATATCTGGTGCCAAGCCGCCCGCCGGTTCATACAGTACAAAACCTTGCTCGTTCAGGCTGGCGGAAGGCAACATCCCCATCGAGCCCGTGATCATTGCGCACTCGTCAGACAGAATGTCGCCAAACAGGTTGGAGCACCGCAGAACGTCAAACTGTGATGGATCTTTAATCAGCTGCATGGTGGCGTTGTCGATGTACATATGCGCCAGTTCGACATCCGGGTATTCCGTgg            >  1:181107/1‑251 (MQ=255)
                                                                                                                                                                                                                       ttttGCTTGCGATATCTGGTGCCAAGCCGCCCGCCGGTTCATACAGTCCAAAACCTTGCTCGTTCAGGCTGGCGGAAGGCAACATCCCCATCGAGCCAGTGATCATTGCGCACTCGTCAGACAGAATGTCGCCAAACAGGTTGGAGCACAGCAGAACGTCAAACTGTGATGGATCTTTAATCAGCTGCATGGTGGCGTTGTCGATGTACATATGCGCCAGTTCGACATCCGGGTATTCCGTGGCGAtctc     <  1:383651/250‑1 (MQ=255)
                                                                                                                                                                                                                          tGCCTGCGATATCTGGTGCCAAGCCGCCCGCCGGTTCATACAGTCCAAAACCTTGCTCGTTCAGGCTGGCGGAAGGCAACATCCCCATCGAGCCAGTGATCATTGCGCACTCGTCAGACAGAATGTCGCCAAACAGGTTGGAGCACAGCAGAACGTCAAACTGTGATGGATCTTTAATCAGCTGCATGGTGGCGTTGTCGATGTACATATGCGCCAGTTCGACATCCGGGTATTCCGTGGCGATCTCGtt  >  1:395987/1‑250 (MQ=255)
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                     
ATTACACCCCTTCTGCTACATAGCGGGCAATGATATCGCCCATTTCATCGGTACTAACGGCGGCAGCGCCACGGGCTAAATCCCCGGTGCGAATGCCTTCTTCTAATGCGCGGTTAATGGCGCGTTCAATGGCGCAAGCCGCATCATCGGCATCCAGGCTGTAACGCAGCAGCAGTGCCAGCGAAAGGATTTGTGCAATCGGGTTGGCGATGTTTTTGCCTGCGATATCTGGTGCCGAGCCGCCCGCCGGTTCATACAGTCCAAAACCTTGCTCGTTCAGGCTGGCGGAAGGCAACATCCCCATCGAGCCAGTGATCATTGCGCACTCGTCAGACAGAATGTCGCCAAACAGGTTGGAGCACAGCAGAACGTCAAACTGTGATGGATCTTTAATCAGCTGCATGGTGGCGTTGTCGATGTACATATGCGCCAGTTCGACATCCGGGTATTCCGTGGCGATCTCGTT  >  NC_000913/80866‑81331

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