Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 83,371 G→A G53G (GGC→GGT leuA ← 2‑isopropylmalate synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091383,3710GA94.4% 208.5 / 0.8 75G53G (GGC→GGTleuA2‑isopropylmalate synthase
Reads supporting (aligned to +/- strand):  ref base G (1/3);  new base A (29/39);  total (30/45)
Fisher's exact test for biased strand distribution p-value = 6.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.92e-01

TACGGGCGCGTTTCACCATATAGATAGCGCGTTCGATCACCTCGTCCAGCGTGCTGCGCAGCTTGGTGGCGATGTGCATTGGCGAAGTGGCAATAAAGGTATGAATACGGAAGGCTTCGGCGACTTTCAGGGATTCGGCCGCCACGTCGATATCTTTTTCCACGCAGCGAGCTAACGCACATACGCGGCTGTTTTTAACCTGGCGGGCGATGGTTTGCACCGATTCAAAATCGCCCGGCGAAGAGACGGGGAAACCGACTTCCATCACGTCAACACCCATACGCTCAAGGGCCAGCGCAATTTGCAGTTTTTCTTTCACACTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTG  >  NC_000913/83139‑83608
                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                             
tACGGGCGCGTTTCACCATATAGATAGCGCGTTCGATCACTTCATCCAGCGTGCTGCGCAGCTTGGTGGCGATGTGCATTGGCGAAGTGGCAATAAAGGTATGAATACGGAAGGCTTCGGCGACTTTCAGGGATTCGGCCGCCACGTCGATATCTTTCTCCACGCAGCGAGCTAACGCGCATACGCGGCTGTTTTTAACCTGGCGAGCGATGGTTTGCACTGATTCAAAATCACCCGGCGAAGAAACggg                                                                                                                                                                                                                              <  1:349086/250‑1 (MQ=255)
 aCGGGCGCTTTTCACCATATAGATAGCGCGTTCGATCACGTCACACAGCGTGCTGCGCAGCTTGGGGGCGATGTGCATTGGCGAAGTGGCAATAAAGGTATGATTACGGAAGGCTTCGGCGACTTTCAGGGATTCGGCCGCCACGTCGATATCTTTCTCCACGCAGCGAGCTAACGCGCATACGCGGCTGTTTTTAACCTGGCGAGCGATGGTTTGCACTGATTCAAAATCACCCGGCGAAGAAACgggg                                                                                                                                                                                                                             <  2:77596/250‑1 (MQ=255)
 aCGGGCGCGTTTCACCATATAGATAGCGCGTTCGATCACTTCATCCAGCGTGCTGCGCAGCTTGGTGGCGATGTGCATTGGCGAAGTGGCAATAAAGGTATGAATACGGAAGGCTTCGGCGACTTTCAGGGATTCGGCCGCCACGTCGATATCTTTCTCCACGCAGCGAGCTAACGCGCATACGCGGCTGTTTTTAACCTGGCGAGCGATGGTTTGCACTGATTCAAAATCACCCGGCGAAGAAACgggg                                                                                                                                                                                                                             <  1:81305/250‑1 (MQ=255)
 aCGGCCGCGTTTCACCATATAGATAGCGCTTTCGATAACTTCATCCAGCGGGCGGCGCAGCTTGGTGGCGATGTGCATTGGCGAAGTGGCAATAAAGGTATGAATCCGGAAGGCTTCGGCGACTTTCAGGGATTCGGCCGCCACGTCGATATCTTTCTCCACGCAGCTAGCTAACGCGCATACTCGGCTGTTTTTAACCTGGCGATCGATGGTTTGCACTGATTCAAAATCACCCGTCGAAGAAACGGGGa                                                                                                                                                                                                                            <  2:647212/251‑1 (MQ=255)
          tttCACCATATAGATAGCGCGTTCGATCACCTCGTCCAGCGTGCTGCGCAGCTTGGGGGCGATGTGCATTGGCGAAGTGGCAATAATGGTATGAATACGGAAGGCTTCGGCGACTTTCAGGGATTCGGCCGCCACGTCGATATCTTTTTCCACGCAGCGAGCTAACGCACATACGCGGCTGTTTTTAACCTGGCGGGCGATGGTTTGCACCGGTTCAAAATCGCCCGGCGAAGAGACGGGGAAACCGACtt                                                                                                                                                                                                                   >  1:55223/1‑251 (MQ=255)
             cACCATATAGATAGCGCGTTCGATCACTTCATCCAGCGTGCTGCGCAGCTTGGTGGCGATGTGCATTGGCGAAGTGGCAATAAAGGTATGAATACGGAAGGCTTCGGCGACTTTCAGGGATTCGGCCGCCACGTCGATATCTTTCTCCACGCAGCGAGCTAACGCGCATACGCGGCTGTTTTTAACCTGGCGAGCGATGGTTTGCACTGATTCAAAATCACCCGGCGAAGAAACGGGGAAACCGACTTcc                                                                                                                                                                                                                 <  1:185919/250‑1 (MQ=255)
                     agatagCGCGTTCGATCACTTCATCCAGCGTGCTGCGCAGCTTGGTGGCGATGTGCATTGGCGAAGTGGCAATAAAGGTATGAATACGGAAGGCTTCGGCGACTTTCAGGGATTCGGCCGCCACGTCGATATCTTTCTCCACGCAGCGAGCTAACGCGCATACGCGGCTGTTTTTAACCTGGCGAGCGATGGTTTGCACTGATTCAAAATCACCCGGCGAAGAAACGGGGAAACCGACTTCCATCACGTc                                                                                                                                                                                                         >  1:447966/1‑250 (MQ=255)
                      gatagCGCGTTCGATCACTTCATCCAGCGTGCTGCGCAGCTTGGTGGCGATGTGCATTGGCGAAGTGGCAATAAAGGTATGAATACGGAAGGCTTCGGCGACTTTCAGGGATTCGGCCGCCACGTCGATATCTTTCTCCACGCAGCGAGCTAACGCGCATACGCGGCTGTTTTTAACCTGGCGAGCGATGGTTTGCACTGATTCCAAATCACCCGGGGAAGCAACGGGGAAACCGCCTTCCATCACGTCa                                                                                                                                                                                                        >  2:21517/1‑250 (MQ=255)
                       atagCGCGTTCGATCACTTCATCCAGCGTGCTGCGCAGCTTGGTGGCGATGTGCATTGGCGAAGTGGCAATAAAGGTATGAATACGGAAGGCTTCGGCGACTTTCAGGGATTCGGCCGCCACGTCGATATCTTTCTCCACGCAGCGAGCTAACGCGCATACGCGGCTGTTTTTAACCTGGCGAGCGATGGTTTGCACTGATTCAAAATCACCCGGCGAAGAAACGGGGAAACCGACTTCCATCACGTCAac                                                                                                                                                                                                      >  1:241202/1‑251 (MQ=255)
                           cgcgTTCGATCACTTCATCCAGCGTGCTGCGCAGCTTGGTGGCGATGTGCATTGGCGAAGTGGCAATAAAGGTATGAATACGGAAGGCTTCGGCGACTTTCAGGGATTCGGCCGCCACGTCGATATCTTTCTCCACGCAGCGAGCTAACGCGCATACGCGGCTGTTTTTAACCTGGCGAGCGATGGTTTGCACTGATTCAAAATCACCCGGCGAAGAAACGGGGAAACCGACTTCCATCACGTCAACAccc                                                                                                                                                                                                  >  1:185003/1‑251 (MQ=255)
                           cgcgTTCGATCACTTCATCCAGCGTGCTGCGCAGCTTGGAGGCGTTGTGAATTGGCGAAGTGGCAATAAAGGTATGAATACGGAAGGCTTCGGCGACTTTCAGGGATTCGGCCGCCACGTCGATATCTTTCTCCACGCAGCGAGCTACCGCGCATACGCGGCTGTTTTTAACCTGGCGAGCGATGGTTTGCACTGATTCAAAATCACCCGGCGAAGAAACGGGGAAACCGACTTCCATCACGTCAACAccc                                                                                                                                                                                                  >  1:557466/1‑251 (MQ=255)
                             cgTTCGATCACTTCATCCAGCGTGCTGCGCAGCTTGGTGGCGATGTGCATTGGCGAAGTGGCAATAAAGGTATGAATACGGAAGGCTTCGGCGACTTTCAGGGATTCGGCCGCCACGTCGATATCTTTCTCCACGCAGCGAGCTAACGCGCATACGCGGCTGTTTTTAACCTGGCGAGCGATGGTTTGCACTGATTCAAAATCACCCGGCGAAGAAACGGGGAAACCGACTTCCATCACGTCAACACCCAt                                                                                                                                                                                                >  1:337034/1‑251 (MQ=255)
                              gTTCGATCACTTCATCCAGCGTGCTGCGCAGCTTGGTGGCGATGTGCATTGGCGAAGTGGCAATAAAGGTATGAATACGGAAGGCTTCGGCGACTTTCAGGGATTCGGCCGCCACGTCGATATCTTTCTCCACGCAGCGAGCTAACGCGCATACGCGGCTGTTTTTAACCTGGCGAGCGATGGTTTGCACTGATTCAAAATCACCCGGCGAAGAAACGGGGAAACCGACTTCCATCACGTCAACACCCATa                                                                                                                                                                                               <  1:604572/251‑1 (MQ=255)
                                tCGATCACTTCATCCAGCGTGCTGCGCAGCTTGGTGGCGATGTGCAGTGGCGAAGTGGCAATAAAGGTATGAATACGGAAGGCTTCGGCGACTTTCAGGGATTCGGCCGCCACGTCGATATCTTTCTCCACGCAGCGAGCTAACGCGCATACGCGGCTGTTTTTAACCTGGCGAGCGATGGTTTGCACTGATTCAAAATCACCCGGCGAAGAAACGGGGAAACCGACTTCCATCACGTCAACACCCATAc                                                                                                                                                                                              <  1:452708/250‑1 (MQ=255)
                                   aTCACTTCATCCAGCGTGCTGCGCAGCTTGGTGGCGATGTGCATTGGCGAAGTGGCAATAAAGGTATGAATACGGAAGGCTTCGGCGACTTTCAGGGATTCGGCCGCCACGTCGATATCTTTCTCCACGCAGCGAGCTAACGCGCATACGCGGCTGTTTTTAACCTGGCGAGCGATGGTTTGCACTGATTCAAAATCACCCGGCGAAGAAACGGGGAAACCGACTTCCATCACGTCAACACCCATACGCTc                                                                                                                                                                                          >  1:366018/1‑251 (MQ=255)
                                          catCCAGCGTGCTGCGCAGCTTGGTGGCGATGTGCATTGGCGAAGTGGCAATAAAGGTATGAATACGGAAGGCTTCGGCGACTTTCAGGGATTCGGCCGCCACGTCGATATCTTTCTCCACGCAGCGAGCTAACGCGCATACGCGGCTGTTTTTAACCTGGCGAGCGATGGTTTGCACTGATTCAAAATCACCCGGCGAAGAAACGGGGAAACCGACTTCCATCACGTCAACACCCATACGCTCAAGGGcc                                                                                                                                                                                   <  1:114382/249‑1 (MQ=255)
                                                 cGTGCTGCGCAGCTTGGTGGCGATGTGCATTGGCGAAGTGGCAATAAAGGTATGAATACGGAAGGCTTCGGCGACTTTCAGGGATTCGGCCGCCACGTCGATATCTTTCTCCACGCAGCGAGCTAACGCGCATACGCGGCTGTTTTTAACCTGGCGAGCGATGGTTTGCACTGATTCAAAATCACCCGGCGAAGAAACGGGGAAACCGACTTCCATCACGTCAACACCCATACGCTCAAGGGCCAGCGCaa                                                                                                                                                                            <  1:21517/251‑1 (MQ=255)
                                                 cGTGCTGCGCAGCTTGGTGGCGATGTGCATTGGCGAAGTGGCAATAAAGGTATGAATACGGAAGGCTTCGGCGACTTTCAGGGATTCGGCCGCCACGTCGATATCTTTCTCCACGCAGCGAGCTAACGCGCATACGCGGCTGTTTTTAACCTGGCGAGCGATGGTTTGCACTGATTCAAAATCACCCGGCGAAGAAACGGGGAAACCGACTTCCATCACGTCAACACCCATACGCTCAAGGGCCAGCGCaa                                                                                                                                                                            <  1:424831/251‑1 (MQ=255)
                                                     ctgcGCAGCTTGGTGGCGCTGTGCATTGGCGAAGTGGCAATAAAGGTATGAATACGGAAGGCTTCGGCGACTTTCAGGGATTCGGCCGCCACGTCGATATCTTTCTCCACGCAGCGAGCTAACGCGCATACGCGGCTGTTTTTAACCTGGCGAGCGATGGTTTGCACTGATTCAAAATCACCCGGCGAAGAAACGGGGAAACCGACTTCCATCACGTCAACACCCATCCGCTCAAGGGCCAGCGCAATTTg                                                                                                                                                                        >  2:110206/1‑251 (MQ=255)
                                                                ggtggCGATGGGCATGGGCGAAGTGGCAATAAAGGTATGAATACTGAAGGCTGCGGCGACTTTCAGGGATTCGGCCGCCACGTCGATATCTTTCTCCACGCAGCGAGCTAACGCGCATACGCGGCTGTTTTTAACCTGGCGAGCGATGGTTTGCACTGATTCAAAATCACCCGGCGAAGAAACGGGGAAACCGACTTCCATCACGTCAACACCCATACGCTCAAGGGCCAGCGCAATTTGCAGTTtttctt                                                                                                                                                             <  1:104456/251‑1 (MQ=255)
                                                                 gtggCGATGTGCATTGGCGAAGTGGCAATAAAGGTATGAATACGGAAGGCTTCGGCGACTTTCAGGGATTCGGCCGCCACGTCGATATCTTTCTCCACGCAGCGAGCTAACGCGCATACGCGGCTGTTTTTAACCTGGCGAGCGATGGTTTGCACTGATTCAAAATCACCCGGCGAAGAAACGGGGAAACCGACTTCCATCACGTCAACACCCATACGCTCAAGGGCCAGCGCAATTTGCAGTTtttcttt                                                                                                                                                            <  1:358840/251‑1 (MQ=255)
                                                                   ggCGATGTGCATTGGCGAAGTGGCAATAAAGGTATGAATACGGAAGGCTTCGGCGACTTTCAGGGATTCGGCCGCCACGTCGATATCTTTCTCCACGCAGCGAGCTAACGCGCATACGCGGCTGTTTTTAACCTGGCGAGCGATGGTTTGCACTGATTCAAAATCACCCGGCGAAGAAACGGGGAAACCGACTTCCATCACGTCAACACCCATACGCTCAAGGGCCAGCGCAATTTGCAGTTtttctttc                                                                                                                                                           <  1:202776/250‑1 (MQ=255)
                                                                               tGGCGAAGTGGCAATAAAGGTATGAATACGGAAGGCTTCGGCGACTTTCAGGGATTCGGCCGCCACGTCGATATCTTTCTCCACGCAGCGAGCTAACGCGCATACGCGGCTGTTTTTAACCTGGCGAGCGATGGTTTGCACTGATTCAAAATCACCCGGCGAAGAAACGGGGCAACCGACTTCCATCACGTCAACACCCATACGCTCAAGGGCCAGCGCAATTTGCAGTTTTTCTTTCACACTCACGCTTg                                                                                                                                              >  2:509201/1‑251 (MQ=255)
                                                                                        ggCAATAAAGGTATGAATACGGAAGGCTTCGGCGACTTTCAGGGATTCGGCCGCCACGTCGATATCTTTCTCCACGCAGCGAGCTAACGCGCATACGCGGCTGTTTTTAACCTGGCGAGCGATGGTTTGCACTGATTCAAAATCACCCGGCGAAGAAACGGGGAAACCGACTTCCATCACGTCAACACCCATACGCTCAAGGGCCAGCGCAATTTGCAGTTTTTCTTTCACACTCAAGCTTGCCTGTAACg                                                                                                                                     >  1:154546/1‑251 (MQ=255)
                                                                                              aaagggaggAATACGGAAGGCTTCGGCGACTTTCAGGGATTCGGCCGCCACGTTGATATCTTTCTCCACGCAGCGAGCTAACGCGCATACGCGGCTGTTTTTAACCTGGCGAGCGATGGTGTGCACTGATTCAAAAGCACCCGGCGAAGAAACGGGGAAACCGACTTTCATCACGTCAACACCCATACGCTCAAGGGCCAGCGCAATTTGCAGTTTTTCTTTCCCACTCAAGCTTGCCTGGAACGCCTGtt                                                                                                                               <  2:298919/243‑1 (MQ=255)
                                                                                               aaGGTATGAATACGGAAGGCTTCGGCGACTTTCAGGGATTCGGCCGCCACGTCGATATCTTTCTCCACGCAGCGAGCTAACGCGCATACGCGGCTGTTTTTAACCTGGCGAGCGATGGTTTGCACTGATTCAAAATCACCCGGCGAAGAAACGGGGAAACCGACTTCCATCACGTCAACACCCATACGCTCAAGGGCCAGCGCAATTTGCAGTTTTTCTTTCACACTCAAGCTTGCCTGTAACGCCTGTTc                                                                                                                              <  1:110206/251‑1 (MQ=255)
                                                                                                      gAATACGGAAGGCTTCGGCGACTTTCAGGGATTCGGCCGCCACGTCGATATCTTTCTCCACGCAGCGAGCTAACGCGCATACGCGGCTGTTTTTAACCTGGCGAGCGATGGTTTGCACTGATTCAAAATCACCCGGCGAAGAAACGGGGAAACCGACTTCCATCACGTCAACACCCATACGCTCAAGGGCCAGCGCAATTTGCAGTTTTTCTTTCACACTCAAGCTTGCCTGTAACGCCTGTTCACCGTc                                                                                                                        <  1:145222/250‑1 (MQ=255)
                                                                                                       aaTACGGAAGGCTTCGGCGACTTTCAGGGATTCGGCCGCCACGTCGATATCTTTCTCCACGCAGCGAGCTAACGCGCATACGCGGCTGTTTTTAACCTGGCGAGCGATGGTTTGCACTGATTCAAAATCACCCGGCGAAGAAACGGGGAAACCGACTTCCATCACGTCAACACCCATACGCTCAAGGGCCAGCGCAATTTGCAGTTTTTCTTTCACACTCAAGCTTGCCTGTAACGCCTGTTCACCGTcgc                                                                                                                      >  1:204901/1‑251 (MQ=255)
                                                                                                                ggCTTCGGCGCCTTTCAGGGATTCGGCCGCCACGTCGATATCTTTCTCCACGCAGCGAGCTAACGCGCATACGCGGCTGTTTTTAACCTGGCGAGCGATGGTTTGCACTGATTCACAATCACCCGGCGAAGAAACGGGGAAACCGACTTCCATCACGTCAACACCCATACGCTCAAGGGCCAGCGCACTTTGCAGTTTTTCTTTCACACTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTgg                                                                                                             >  2:697299/1‑251 (MQ=255)
                                                                                                                ggCTTCGGCGACTTTCAGGGATTCGGCCGCCACGTCGATATCTTTCTCCACGCAGCGAGCTAACGCGCATACGCGGCTGTTTTTAACCTGGCGAGCGATGGTTTGCACTGATTCAAAATCACCCGGCGAAGAAACGGGGAAACCGACTTCCATCACGTCAACACCCATACGCTCAAGGGCCAGCGCAATTTGCAGTTTTTCTTTCACACTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTgg                                                                                                             >  1:289220/1‑251 (MQ=255)
                                                                                                                  cTTCGGCGACTTTCAGGGATTCGGCCGCCACGTCGATATCTTTCTCCACGCAGCGAGCTAACGCGCATACGCGGCTGTTTTTAACCTGGCGAGCGATGGTTTGCACTGATTCAAAATCACCCGGCGAAGAAACGGGGAAACCGACTTCCATCACGTCAACACCCATACGCTCAAGGGCCAGCGCAATTTGCAGTTTTTCTTTCACACTCAAGCTTGCCTGTAACGCCTGTTCACCGTcgcg                                                                                                                     >  1:24779/1‑241 (MQ=255)
                                                                                                                  cTTCGGCGACTTTCAGGGATTCGGCCGCCACGGCGATATCTTTCTCCACGCAGCGAGCTAACGCGCATACGCGGCTGTTTGTAACCTGGCGAGCGATGGTTTGCACTGATTCAAAATCACCCGGCGAAGAAACGGGGAAACCGACTTCCATCACGTCAACACCCATACGCTCAAGGGCCAGCGCAATTTGCAGTTTTTCTTTCACACTCAAGCTTGCCTGTAACGCCTGTTCACCGTcgcg                                                                                                                     <  2:24779/241‑1 (MQ=255)
                                                                                                                     gctcgtctgtCAGGTATGCGGCCGCCACGGCGATATCTTACTCCACGCAGCGAGCTAACGCGCATACGCGTCTGTTTTTAACCTGGCGAGCGATTGTTTGCACTGATTCAAAATCACCCGGCGAAGAAACGGGGACACCGACTTCCATCACGTCAACCCCCATACGCTCAAGGGCCAGCGCAATTTGCAGTTTTTCTTTCACACTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCg                                                                                                        <  2:574614/242‑1 (MQ=255)
                                                                                                                     cgtcgccTTTCAGGAGTTCGGCCGCCACGTCGATAGCGTTCTCCACGCAGCGAGCTAACGCGCATACGCGGCTGTTTTTAACCTGGCGAGCGATGGTTTGCACTGATGCAAAATCACCCGGCGAAGAAACGGGGAAACCGACTTCCATCACGGCAACACCCATACGCTCAAGGGCCAGCGCAATTTGCAGTTTTTCTTTAACTCTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCg                                                                                                        <  2:283726/245‑1 (MQ=255)
                                                                                                                     cGGCGACTTTCAGGGATTCGGCCGCCACGTCGATATCTTTCTCCACGCAGCGAGCTAACGCGCATACGCGGCTGTTTTTAACCTGGCGAGCGATGGTTTGCACTGATTCAAAATCACCCGGCGAAGAAACGGGGAAACCGACTTCCATCACGTCCACACCCATACGCTCAAGGGCCAGCGCAATTTGCAGTTTTTCTTTCACACTCAAGCTTGCCTGTAACGCATGTTCACCGTCGCGCacagtgggatcg                                                                                                        >  2:477107/1‑240 (MQ=255)
                                                                                                                         gACTTTCAGTGATTCGGCCGCCACGTCGATATCTGTCTCCACGCAGCGAGCTAACGCGCATACGCGGCTGGTTTTAACCTGGCGAGCGATGGTTTGCACTGATTCAAAATCACCCGGCGAAGAAACGGGGAAACCGACTTCCATCACGTCAACACCCATACGCTCAAGGGCCAGCGCAATTTGCAGTTTTTCTTTCACACTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTTGTATCGaaaa                                                                                                    <  1:132866/251‑1 (MQ=255)
                                                                                                                         gACTTGCAGGGACTCGGCAGCCACGTCGATATCTTTCCCCACGCAGCGAGCTAACGCGAATACGCGGCTGTTTTTAACCTGGCGAGCGATGGTTTGCACTGATTCAAAATCTCCCGGCGAAGAAACGGGGAAACCGACTTCCATCACGTCAACACCCATACGCTCAAGGGCCAGCGCAAGTTGCAGTTTTTCTTTCACACTCAAGCTTGCCTGTAACGCCTGTTCACCGGCGCGCAATGTGGTATCGaaaa                                                                                                    <  1:11198/251‑1 (MQ=255)
                                                                                                                             ttCAGGGATTCGGCCGCCACGTCGATATCTTTCTCCACGCAGCGAGCTAACGCGCATACGCGGCTGTTTTTAACCTGGCGAGCGATGGTTTGCACTGATTCAAAATCACCCGGCGAAGAAACGGGGAAACCGACTTCCAt                                                                                                                                                                                                               >  2:482442/1‑140 (MQ=255)
                                                                                                                             ttCAGGGATTCGGCCGCCACGTCGATATCTTTCTCCACGCAGCGAGCTAACGCGCATACGCGGCTGTTTTTAACCTGGCGAGCGATGGTTTGCACTGATTCAAAATCACCCGGCGAAGAAACGGGGAAACCGACTTCCAt                                                                                                                                                                                                               <  1:482442/140‑1 (MQ=255)
                                                                                                                                aGGGATTCGGCCGCCACGTCGATATCTTTCTCCACGCAGCGAGCTAACGCGCATACGCGGCTGTTTTTAACCTGGCGAGCGATGGTTTGCACTGATTCAAAATCACCCGGCGAAGAAACGGGGAAACCGACTTCCATCACGTCAACACCCATACGCTCAAGGGCCAGCGCAATTTGCAGTTTTTCTTTCACACTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGAc                                                                                             >  1:355734/1‑251 (MQ=255)
                                                                                                                                  ggATTCGGCCGCCACGTCGATATCTTTCTCCACGCAGCGAGCTAACGCGCATACGCGGCTGTTTTTAACCTGGCGAGCGATGGTTTGCACTGATTCAAAATCACCCGGCGAAGAAACGGGGAAACCGACTTCCATCACGTCAACACCCATACGCTCAAGGGCCAGCGCAATTTGCAGTTTTTCTTTCACACTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACtt                                                                                           >  1:686313/1‑251 (MQ=255)
                                                                                                                                    aTTCGGCCGCCACGTCGATATCTTTCTCCACGCAGCGAGCTAACGCGCATACGCGGCTGTTTTTAACCTGGCGAGCGATGGTTTGCACTGATTCAAAATCACCCGGCGAAGAAACGGGGAAACCGACTTCCATCACGTCAACACCCATACGCTCAAGGGCCAGCGCAAttt                                                                                                                                                                         >  2:616978/1‑171 (MQ=255)
                                                                                                                                    aTTCGGCCGCCACGTCGATATCTTTCTCCACGCAGCGAGCTAACGCGCATACGCGGCTGTTTTTAACCTGGCGAGCGATGGTTTGCACTGATTCAAAATCACCCGGCGAAGAAACGGGGAAACCGACTTCCATCACGTCAACACCCATACGCTCAAGGGCCAGCGCAAttt                                                                                                                                                                         <  1:616978/171‑1 (MQ=255)
                                                                                                                                        ggCCGCCACGTCGATATCTTTCTCCACGCAGCGAGCTAACGCGCATACGCGGCTGTTTTTAACCTGGCGAGCGATGGTTTGCACTGATTCAAAATCACCCGTCGAAGAAACGGGGAAACCGACTTCCATCACGTCAACACCCATACGCTCAAGGGCCAGCGCAATTTGCAGTTTTTCTTTCACACTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGATAATAATGACTTGCTGGc                                                                                     >  1:540324/1‑251 (MQ=255)
                                                                                                                                         gccgccACGTCGATATCTTTCTCCACGCAGCGAGCTAACGCGCATACGCGGCTGTTTTTAACCTGGCGAGCGATGGTTTGCACTGATTCAAAATCACCCGGCGAAGAAACGGGGAAACCGACTTCCATCACGTCAACACCCATACGCTCAAGGGCCAGCGCAATTTGCAGTTTTTCTTTCACACTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCt                                                                                    >  1:506307/1‑251 (MQ=255)
                                                                                                                                         gccgccACGTCGATATCTTTCTCCACGCAGCGAGCTAACGCGCATACGCGGCTGTTTTTAACCTGGCGAGCGATGGTTTGCACTGATTCAAAATCACCCGGCGAAGAAACGGGGAAACCGACTTCCATCACGTCAACACCCATACGCTCAAGGGCCAGCGCAATTTGCAGTTTTTCTTTCACACTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGGATCGAAAATAATGACTTGCTGGCt                                                                                    <  1:514274/251‑1 (MQ=255)
                                                                                                                                              cACGTCGATATCTTTCTCCACGCAGCGAGCTAACGCGCATACGCGGCTGTTTTTAACCTGGCGAGCGATGGTTTGCACTTATTCAAAATCACCCGGCGAAGAAACGGGGAAACCGACTTCCATCACGTCAACACCCATACGCTCAAGGGCCAGCGCAATTTGCAGTTTTTCTTTCACACGCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATgg                                                                               <  2:355734/251‑1 (MQ=255)
                                                                                                                                               aCGTCGATATCTTTCTCCACGCAGCGAGCTAACGCGCATACGCGGCTGTTTTTAACCTGGCGAGCGATGGTTTGCACTGATTCAAAATCACCCGGCGAAGAAACGGGGAAACCGACTTCCATCACGTCAACACCCATCCGCTCAAGGGCCCGCGCAATTTGCAGTTTTTCTTTCACACTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATACCTTGCTGGCCCATGGt                                                                              >  2:429486/1‑251 (MQ=255)
                                                                                                                                               aCGTCGATATCTTTCTCCACGCAGCGAGCTAACGCGCATACGCGGCTGTTTTTAACCTGGCGAGCGATGGTTTGCACTGATTCAAAATCACCCGGCGAAGAAACGGGGAAACCGACTTCCATCACGTCAACACCCATACGCTCAAGGGCCAGCGCAATTTGCAGTTTTTCTTTCACACTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGt                                                                              >  1:386515/1‑251 (MQ=255)
                                                                                                                                                cGTCGATATCTTTCTCCACGCAGCGAGCTAACGCGCATACGCGGCTGTTTTTAACCTGGCGAGCGATGGTTTGCACTGATTCAAAATCACCCGGCGAAGAAACGGGGAAACCGACTTCCATCACGTCAACACCCATACGCTCAAGGGCCAGCGCAATTTGCAGTTTTTCTTTCACACTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGtt                                                                             <  1:343315/251‑1 (MQ=255)
                                                                                                                                                   cGATATCTTTTTCCACGCAGCGAGCTAACGCACATACGCGGCTGTTTTTAACCTGGCGGGCGATGGTTTGCACCGATTCAAAATCGCCCGGCGAAGAGACGGGGAAACCGACTTCCATCACGTCAACACCCATACGCTCAAGGGCCAGCGCAATTTGCAGTTTTTCTTTCACACTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTgg                                                                          <  1:361362/251‑1 (MQ=255)
                                                                                                                                                   cGATAGCTTTTTCCACGCAGCGAGCTAACGCACATACGCGGCTGTTTTTAACCTGGCGGGCGATGGTTTGCACCGATTCAAAATCGCCCGGCGAAGAGACGGGGAAACCGACTTCCATCACGTCAACACCCATACGCTCAAGGGCCAGCGCAAGTTGCAGTTTTTCTTTCACACTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTg                                                                           <  2:149244/250‑1 (MQ=255)
                                                                                                                                                         ttttctCCACGCCGCTAGCTAACGCGCATACGCGGCTGTTTTTAACCTGGCGAGCGATGGTTTGCACTGATTCAAAATCACCCGGCTAAGAAACGGGGAAACCGACTTCCATCACGTCAACACCCATACGCTCAAGGGCCAGCGCAATTTGCAGTTTTTCTTTCACACTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCTCAATGTGGTATCGAAAATAATGACTTGCTGGCGCATGGTTTGCGTCCt                                                                     <  2:150905/245‑1 (MQ=255)
                                                                                                                                                               cctcGCAGCGAGCTAACGCGCATACGCGGCTGCTGTGAACCTGGCGAGCGATGGTTTGCACTGATTCAAAATCACCCGGCGAAGAAACGGGGAAAACGACTTCCATCACGTCAACACCCATACGCTCAAGGGCCAGAGCAATTTGCAGTTTTTCTTTCACACTCAATCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGtc                                                                 <  2:334020/245‑1 (MQ=255)
                                                                                                                                                                  cGCAGCGAGCTAACGCGCATACGCGGCTGTTTTTAACCTGGCGAGCGATGGTTTGCACTGATTCAAAATCACCCGGCGAAGAAACGGGGAAACCGACTTCCATCACGTCAACACCCATACGCTCAAGGGCCAGCGCAATTTGCAGTTTTTCTTTCACACTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCtttt                                                           <  1:689066/251‑1 (MQ=255)
                                                                                                                                                                   gCAGCGAGCTAACGCGCATACGCGGCTGTTTTTAACCTGGCGAGCGATGGTTTGCACTGATTCAAAATCACCCGGCGAAGAAACGGGGAAACCGACTTCCATCACGTCAACACCCATACGCTCAAGGGCCAGCGCAATTTGCAGTTTTTCTTTCACACTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTa                                                          <  1:301160/251‑1 (MQ=255)
                                                                                                                                                                          gatAACGCGCATACGCGGCTGTTTTTAACCTGGCGAGCGATGGTTTGCACTGATTCAAAATCTCCCGGCGAAGAAACGGGGAAACCGACTTCCATCACGTCAACACCCATACGCTCAAGGGCCAGCGCAAGTTGCAGTTTTTCTTGCACACTCAAGCTGGTCTGTAACGCCGGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAgcgtc                                                   <  2:225277/249‑3 (MQ=255)
                                                                                                                                                                          gCTAACGCGCATACGCGGCTGTTTTTAACCTGGCGAGCGATGGTTTGCACTGATTCAAAATCACCCGGCGAAGAAACGGGGAAACCGACTTCCATCACGTCAACACCCATACGCTCAAGGGCCAGCGCAATTTGCAGTTTTTCTTTCACACTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGcc                                                   <  1:40821/251‑1 (MQ=255)
                                                                                                                                                                            tAACGCGCATACGCGGCTGTTTTTAACCTGGCGAGCGATGGTTTGCACTGATTCAAAATCACCCGGCGAAGAAACGGGGAAACCGACTTCCATCACGTCAACACCCATACGCTCAAGGGCCAGCGCAATTTGCAGTTTTTCTTTCACACTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTc                                                 >  1:92054/1‑251 (MQ=255)
                                                                                                                                                                            tAACGCGCAAACGCGGTTGTTTTTGACCTGACGAACGATGGTTTGCACTGATTCAAAATCACCCGGCGACGAAACGGGGAAACCGACTTACATCACGTAAACACCCATACGCTCAAGGGCCAGCGCAATTTGCAGTTTTTCTTTCACCCTCACGCTTGCCTGTACCGCCTGTTCACCGTCTCGCTATGTGGTTTTGACAATAATTCCTTGCTGGCTAATGGTTTGGGTCATTGTCTCTTTTattgagcctc                                                 >  2:611936/1‑242 (MQ=255)
                                                                                                                                                                            caaCGCGCAGACGCGGCTGGTTTTAACCTGGCGAGCGCTGGTGTGCACTTAATCAAAATCCCCCGGCGAGGAAACGGGGAAACCGACTTCCATCACGTCAACCCCCATACGCCCAAGGGCCAGCGCAAGTTGCAGTTGTTCTTTCACACTCAAGCTTGCCTGTAACGCCTGTTCTCCGGCGCGCAATGTGGTATCGAAAATAATGCCTTGCTGGCTCATGGTTCGGGTTCTTGTATTTTTGAGAGCGcccc                                                 <  2:408502/250‑3 (MQ=255)
                                                                                                                                                                             aaCGCGCATACGCGGCTGTTTTTAACCTGGCGAGCGATGGTTTGCACTGATTCAAAATCACCCGGCGAAGAAACGGGGAAACCGACTTCCATCACGTCAACACCCATACGCTCAAGGGCCAGCGCAATTTGCAGTTTTTCTTTCACACTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTc                                                 <  2:289220/250‑1 (MQ=255)
                                                                                                                                                                               cgcgaaTACTCGGCTGTTTTTAACCTGGCGAGCGATGGTTTGCACTGATTCAAAATCACCCGGCGAAGAAACGGGGAAACCGACTTCCATCACGTCAACACCCATACGCTCAAGGGCCAGCGCAATTTGCAGTTTTTCTTTCACACTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCt                                              <  2:507375/246‑1 (MQ=255)
                                                                                                                                                                                gcgcaTACGCGGATGTTTTTAACCTGGCGAGCGATGGTTTGCACTGATTCAAAATCACCCGGCGAAGAAACGGGGAAACCGACTTCAATCACGTCAACACCCATACGCTCAAGGGCCAGCGCAATTTGCAGTTTTTCTTTCCAACTCAAGCTTGCCTGTAACGCCTGTTCAACATCGCGCAATGTGTAATCGAAAATAAAGACTTGCTGTCTACTGTTTTTGGTCTTTGTCTCTTTTAGAGCGCCTCTCtt                                             >  2:245715/4‑251 (MQ=255)
                                                                                                                                                                                 cgcctACGCGGGTGTTTGTAACCTGGCGAGCTATGGTTTGCACTGATTCAAAATCACCCGGCGAAGAAACGGGGAAACCGACTTCCATCACGTCAACACCCATACGCTCAAGGGCCAGCGCAATTTGCAGTTTTTCTTTCACACTCAAGCTTGCCTGTAACGCCTGTTCACCGGCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTc                                            <  2:460180/247‑1 (MQ=255)
                                                                                                                                                                                   caTACGCGGCTGTTTTTAACCTGGCGAGCGATTGTTTGCACTGATTCAAAATCACCCGGCGAAGAACCGGGGAAACCGACTTCCATCACGTCAACACCCATACGCTCAAGGGCCAGCGCAATTTGCAGTTTTTCTTTCACACTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCgg                                          <  2:383829/251‑1 (MQ=255)
                                                                                                                                                                                   caTACGCGGCTGTTTTTAACCTGGCGAGCGATGGTTTGCACTGATTCAAAATCACCCGGCGAAGAAACGGGGAAACCGACTTCCATCACGTCAACACCCATACGCTCAAGGGCCAGCGCAATTTGCAGTTTTTCTTTCACACTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCgg                                          <  1:548360/251‑1 (MQ=255)
                                                                                                                                                                                         cgGCTGTTTTTAACCTGGCGAGCGATGGTTTGCACTGATTCAAAATCACCCGGCGAAGAAACGGGGAAACCGACTTCCATCACGTCAACACCCATACGCTCAAGGGCCAGCGCAATTTGCAGTTTTTCTTTCACACTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAAGAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATaa                                    <  2:366018/251‑1 (MQ=255)
                                                                                                                                                                                           gCTGTTTTTAACCTGGCGAGCGATGGTTTGCACTGATTCAAAATCACCCGGCGAAGAAACGGGGAAACCGACTTCCATCACGTCAACACCCATACGCTCAAGGGCCAGCGCAATTTGCAGTTTTTCTTTCACACTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATaaa                                   >  1:59948/1‑250 (MQ=255)
                                                                                                                                                                                                     cccTGGCGAGCGATGGTTTGCACTGATTCAAAATCACCCGGCGAAGAAACGGGGAAACCGACTTCCATCACGTCAACACCCATACGCTCAAGGGCCAGCGCAATTTGCAGTTTTTCTTTCACACTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAACTAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCTTACAAAAACCCgcgc                        >  2:141605/2‑251 (MQ=255)
                                                                                                                                                                                                       cTGGCGAGCGATGGTTTGCACTGATGCAAAATCACCCGGCGAAGAAACGGGGAAACCGCCTTCCATCCCGTCAACACCCATACGATCAAGGGCCAGCGCACTTTGCAGTTTTTCTTTCACACTCAAGCTTGCCTGTAACGCCTGTTAACCGTCGCGAAATGTGTTATCGAAAATAATGACTTGCTGGCTCAGGGTTTGGGACCTCGTCTCTTTTAGAGCGCCTGGCATCAGACATAACAAAACCCGCGCa                       >  2:474821/1‑250 (MQ=255)
                                                                                                                                                                                                           cgtgcgttgggttGCACTGATTCAAAATCACCCGGCGAAGAAACTGGGAAACCGACTGCCATCACGCCAACACCCATATGCTCAAGGGCCAGCGCAATTGGCAGTTTTCCTTGCACACCCAAGCTTGCCTGTAACGCCTGTTCACCGACGCGCAATTTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTGTAGAGCGCCTCTCTTCGGGCATAAAAAACCCCGCTCAATGgc                  <  2:93486/240‑1 (MQ=255)
                                                                                                                                                                                                            gggcgATGGTTTGCACCGAGTCAAAATCGCCCGGGGAAGAGACGGGGAAACCGACTTCCATCTCGTCAACACCCATACGCTCAAGGGCCAGCGCAATTTGCAGTTTTTCTTTCTCAATCAATCTTGCCTGGAACGTCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTGGGGACCTTGTCTCTTTCAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGgcg                 <  2:55223/251‑1 (MQ=255)
                                                                                                                                                                                                                      ttGCACTGATTCAAAATCACCCGGCGAAGAAACGGGGAAACCGACTTCCATCACGTCAACACCCATACGCTCAAGGGCCAGCGCAATTTGCAGTTTTTCTTTCACACTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGttttt        <  2:289356/250‑1 (MQ=255)
                                                                                                                                                                                                                             gaatCAAAATCACCCGGCGAAGAAACGGGGAAACCGACTTCCATCACGTCAACACCCATACGCTCAAGGGCCAGCGCAATTTGCAGTTTTTCTTTCACACTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCTCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTtttgtttg  <  2:154546/246‑1 (MQ=255)
                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                             
TACGGGCGCGTTTCACCATATAGATAGCGCGTTCGATCACCTCGTCCAGCGTGCTGCGCAGCTTGGTGGCGATGTGCATTGGCGAAGTGGCAATAAAGGTATGAATACGGAAGGCTTCGGCGACTTTCAGGGATTCGGCCGCCACGTCGATATCTTTTTCCACGCAGCGAGCTAACGCACATACGCGGCTGTTTTTAACCTGGCGGGCGATGGTTTGCACCGATTCAAAATCGCCCGGCGAAGAGACGGGGAAACCGACTTCCATCACGTCAACACCCATACGCTCAAGGGCCAGCGCAATTTGCAGTTTTTCTTTCACACTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTG  >  NC_000913/83139‑83608

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