Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 82,933 G→A I199I (ATC→ATT leuA ← 2‑isopropylmalate synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091382,9330GA89.8% 124.7 / 2.7 49I199I (ATC→ATTleuA2‑isopropylmalate synthase
Reads supporting (aligned to +/- strand):  ref base G (4/1);  new base A (24/20);  total (28/21)
Fisher's exact test for biased strand distribution p-value = 3.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.10e-01

GCATATTACAAATCTGGCTAACTAACTGGCTGGTGCGCCATATCTCCTGGTGATTAATGGCGGTGTGGACGTTGAGAATATCCTTACGAACTTTGATCGCCATGATGACTTCTTCCAGGGAACAGTTTCCGGCACGCTCGCCGATCCCGTTCATTGCGCCTTCCACCTGGCGTGCACCGGCATGTACCGCCGCCAGTGAGTTTCCGACCGCCAGGCCCAAATCGTCGTGGGTATGTACGGAGATAATGGCTTTGTCGATGTTAGGCACGCGTTCATACAGGCCGCTGATGATTCCGGCGAACTCAAACGGCATGGTGTAGCCCACGGTGTCCGGAATGTTGATGGTGGTGGCACCGGCATTAATCGCCGCTTCGACCACTCGCGCCAGATCGGCAATGGGTGTACGCCCGGCATCTTCGCAAGAAAATTCAACATCATCGGTGTAATTACGGGCGCGTTTCACCATATAGATAGCGCGTTCGA  >  NC_000913/82692‑83174
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                 
gCATATTACAAATCTGGCTAACTAACTGGCTGGTGCGCCATATCTCCTGGTGATTAATGGCGGTGTGGGCGCTGAGAATATCCTTACGAACTTTGATCGCCATGATGACTTCTTCCAGGGAACAGGTTCCGGCACGCTCGCCGATCCCGTTCATTGCGCCTTCCCCCCTGCGTGCACCGGCATGTACCGCCGCCCGTGAGTTTCCGACCGCCACGGCTAAATCGTCGAGGGTATTGACGGAAATAATGGCt                                                                                                                                                                                                                                          >  1:37020/1‑251 (MQ=255)
   tatTACAAATCTGGCTAACTAACTGGCTGGTGCGCCATATCTCCTGGTGATTAATGGCGGTGTGGACGGTGGGAAAATCCTTACGAACTTTGATCGCCATGATGACTTCTTCCAGGGAACAGTTTCCGGCACGCTCGCCGATCCCGTTCATTGCGCCTTCCACCTGGCGTGCACCGGCATGGACCGCCGCCAGTGAGTTGCCGACCGCCAGGCCTAAATCGGCGGGGGGATGTACGGGAATAATgggcttg                                                                                                                                                                                                                                       >  1:132201/1‑246 (MQ=255)
        cAAATCTGGCTAACTAACTGGCTGGTGCGCCATATCTCCTGGTGATTAATGGCGGTGTGGGCGTTGAGAATATCCTTACGAACTTTGATCGCCATGATGACTTCTTCCAGGGAACAATTTCCGGCACGCTCGCCGGTCCCGTTCATTGCGCCTTCCACCTGGCGGGCACCGGCCTGTACCGCCGCCAGTGAGTTGCCGCACGCCAGGCCTAAATCGACGGGGGTATGTACGGAAATAATGGCTTTGTCggt                                                                                                                                                                                                                                  >  1:139856/1‑249 (MQ=255)
               ggCTAACTAACGGGCTGGGCCGCCTTATCCCCTGGTGATTAAGGGCGGTTGGGACGTTGAGAATATCCTTACGAACTTTGATCGCCATGATGACTTCTTCCAGGGAACATTTTCCGGCCCGCTCGCCGATCCCGTTCATTGCGCCTTCCCCCTGGCGTGCACCGGCATGTACCGCCGCCAGTGAGTTGCCGACCGCCAGGCCTAAATCGTCGTGGGTATGTACGCAAATAATGGCTTTGTCGATGtt                                                                                                                                                                                                                               <  2:278203/247‑1 (MQ=255)
                gCTAACTAGCTGGCTGGTGGGCCATATCTCCTGGTGAGTAATGGCTGTGTGGACGTTGAGAATATCCTTACGAACTTTGATCGCCATGATGACTTCTTCCGGGGAACAGTTTCCGGCTCTCTCGCCGATCCCGTTCATTGCGCTTTCCACCTGGCGTGCTCCGGCATGTACCGCCGCCAGTGAGTTGCCTACCGCCAGGCCTAAATCGTCGTGGGTATGTACGGAAATAATGGCTTTGTCGATGTTAggta                                                                                                                                                                                                                          <  1:205993/251‑3 (MQ=255)
                                              cTGGTGATTAATGGCGGTGTGGACGTTGAGAATATCCTTACGAACTTTTATCGCCATGATGACTTCTTCCAGGGAACAGTTTCCCGCACGCTCGCCGCTCCCGTTCCTTGCGCCTTCCACCTGGCGTGCACCGGCCTGTACCCCCGCCAGGGAGTTTCCCACCCCCCGGCCCCAAACCTCGGGGGGTTGTACGGAGATAAAGGCGTTGTCGATGTTTGGCACGCGTTCATAccgagcgctgttgaatgccg                                                                                                                                                                                            >  2:392107/1‑232 (MQ=255)
                                                ggTGATTAATGGCGGGGTGGGCGGTGGGAATATCCTTACGGACTTTTATTGCCATGGTGGCTTCTTCCCGGGGACAGGTTCCGGCACGCTCGCCGATCCCGTTTATTGCGCCTTCCACCTGGCGTGCACCGGCATTGACCGCCCCCCGTGAGTTGCCCGCCGCCCGGGCCAAAACGTCCGGGGGATGTACGGAAATAACGGCTTTGTGGATGTTAGGGTCGCGGTCATACAGGGCGCtggtgtattcggcg                                                                                                                                                                                          >  2:242247/1‑239 (MQ=255)
                                                 gTGATTAATGGCGGTGTGGACTTTGAGAATATCCTTACGAACTTTGATCGCCATGATGACTTCTTCCAGGGAACAGTTTCCGGCACGCTCGCCGATCCCGTTCATTGCGCCTTCCACCTGGCGTGCACCGGCATGTACCGCCGCCAGTGAGTTTCCGACCGCCAGGCCCAAATCGTCGTGGGTATGTACGGAGATAATGGCTTTGTCGATGTTAGGCACGCGTTCATACAGGCCGCTGATGATTCCGGCGa                                                                                                                                                                                         <  1:392107/251‑1 (MQ=255)
                                                                    aCGTTGAGAATATCCTTACGAACTTTGATCGCCATGATGACTTCTTCCAGGGAACAGTTTCCGGCACGCTCGCCGATCCCGTTCATTGCGCCTTCCACCTGGCGTGCACCGGCATGTACCGCCGCCAGTGAGTTGCCGACCGCCAGGCCTAAATCGTCGTGGGTATGTACGGAAATAATGGCTTTGTCGATGTTAGGTACGCGTTCATACAGGCCGCTGATGATTCCGGCGAACTCAAACGGCATGGTGTa                                                                                                                                                                      <  1:391869/251‑1 (MQ=255)
                                                                    aCGTTGAGAATATCCTTACGAACTTTGATCGCCATGATGACTTCTTCAAGGGACCAGTTTCCGGCACGCTCTCCGATCCCGTTCATTGCGCCTTCCACCTGGCGTGCACCGGCATGTACCGCCGCCAGTGAGTTGCCTACCGCCATGCCTAAATCGACGTGGGTCTGTACGGAAATAATGGCTTTGTAGATGTTAGGTACGCGTTCATACAGGCCGCAGCTGCTTCTGTCGACCTCAAACGGCATGgtgt                                                                                                                                                                       >  1:387873/1‑250 (MQ=255)
                                                                        tGAGAATATCCTTACGAACTTTGATCGCCATGATGACTTCTTCCAGGGAACAGTTTCCGGCACGCTCGCCGATCCCGTTCATTGCGCCTTCCACCTGGCGTGCACCGGCATGTACCGCCGCCAGTGAGTTGCCGACCGCCAGGCCTAAATCGTCGTGGGTATGTACGGAAATAATGGCTTTGTCGATGTTAGGTACGCGTTCATACAGGCCGCTGATGATTCCGGCGAACTCCACCGGCATGGTGTAGccc                                                                                                                                                                  >  2:179173/1‑251 (MQ=255)
                                                                               atCCTTACGAACTTTGATCGCCATGATGACTTCTTCCAGGGAACAGTTTCCGGCACGCTCGCCGATCCCGTTCATTGCGCCTTCCACCTGGCGTGCACCGGCATGTACCGCCGCCAGTGAGTTGCCGACCGCCAGGCCTAAATCGTCGTGGGTATGTACGGAAATCATGGCTTTGTCGATGTTAGGTACGCGTTCATACAGGCCGCTGATGATTCCGGCGAACTCAAACGGCATGGTGTAGCCCACGgtgt                                                                                                                                                           >  1:310927/1‑251 (MQ=255)
                                                                                     cccccctcccgctcccttgttgtcttcttCCGGGGAACCGTTCCCGGCCCGCTCGCCGTTCCCGTTCATTGCGCCTTCCCCCTGGCGTGCACCGGCATGTCCCGCCGCCAGTGAGTTGCCGACCGCCAGGCCTAAATCGTCGTGGGTATGTACGGAAATAATGGCTTTGTCGATGTTAGGTACGCGTTCATACAGGCCGCTGATGATTCCGGCGAACTCAAACGGCATGGTGTAGCCCACGGTGTCCGGaa                                                                                                                                                     <  2:309978/228‑1 (MQ=255)
                                                                                                 cGCCATGATGACTTCTTCCAGGGAACAGTTTCCGGCACGCTCGCCGATCCCGTTCATTGAGCCTTCCACCTGGAGTGCACCGGCATGTACCGCCGCCAGTGAGTTGCCGACCGCCAGGCCTAAATCGTCGTGGGTATGTACGGAAATAATGACTTTATCGATGTTAGGTACGAGTTCGTAAAGGCCGATGATGATTCAGGCGAACTCAAACGGCAAGGGGTAGccccaggtggccagcatgtttgagggtg                                                                                                                                         >  2:445738/1‑226 (MQ=255)
                                                                                                  gCCATGCTGACATCTTCACGGAAACAGTTTCCGGCACGCTCGCCGATCCCGTTCATTGCGCCTTCCACCTGGCGTGCACCGGCATGTACCGCCGCAAGTGAGTTGCCGACCGCCAGGCCTAAATCGTCGTGGGTATGGACGGAAATAATGGCTTTGTCGATGTTAGGTCCGCGTTCATACAGGCCGCTGATGATTCCGGCGAAATCAAACGGCATGGGGTAGCCCACGGTGTCCGGGATGTtggtggtggg                                                                                                                                        >  2:429399/1‑250 (MQ=255)
                                                                                                  cccATGATGACTTCTTCCAGGGAACAGTTTCCGGCACGCTCGCCGATCCCGTTCATTGCGCCTTCCACCTGGCGTGCACCGGCATGTACCGCCGCCAGTGAGTTGCCGACCGCCAGGCCTAAATCGTCGTGGGTATGTACGGAAATAATGGCTTTGTCGATGTTAGGTACGCGTTCATACAGGCCGCTGATGATTCCGGCGAACTCAAACGGCATGGTGTAGCCCACGGTGTCCGGAATGTTGAtggtggt                                                                                                                                        <  1:242247/250‑1 (MQ=255)
                                                                                                  cccATGATGACTTCTCCCAGGTAACATTGCCCGGCACGCTCGCCGATCCCTTTCATTGCGCCTCCCACCTGGCGTGCCCCGGCATGTACCGCCGCCAGTGAGTTGCCGACCGCCAGGCCTAAATCGTCGTGGGTATGTACGGAAATAATGGCTTTGTCGATGTTAGGTACGCGTTCATACAGGCCGCTGATGATTCCGGCGAACTCAAACGGCATGGTGTAGCCCACGGTGTCCGGAATGTTGAtggtggt                                                                                                                                        <  2:132201/250‑1 (MQ=255)
                                                                                                    cATGATGACTTCTTCCAGGGAACAGTTTCCGGCACGCTCGCCGATCCCGTTCATTGCGCCTTCCACCTGGCGTGCACCGGCATGTACCGCCGCCAGTGAGTTGCCGACCGCCAGGCCTAAATCGTCGTGGGTATGTACGGAAATAATGGCTTTGTCGATGTTAGGTACGCGTTCATACAGGCCGCTGATGATTCCGGCGAACTCAAACGGCATGGTGTAGCCCACGGTGTCCGGaa                                                                                                                                                     >  1:309978/1‑236 (MQ=255)
                                                                                                        atgaCTTCTTCCAGGGAACCGTTTCCGGCACGCTCGCCGATCCCGTTCATTGCGCCTTCCACCTGGCGTGCACCGGCATGTACCGCCGCCAGTGAGTTGCCGACCGCCAGGCCTAAATCGTCGTGGGTATGTACGGAAATACTGGCTTTGTCGATGTTAGGTACGCGtt                                                                                                                                                                                                                    >  2:270866/1‑169 (MQ=255)
                                                                                                        atgaCTTCTTCCAGGGAACAGTTTCCGGCACGCTCGCCGATCCCGTTCATTGCGCCTTCCACCTGGCGTGCACCGGCATGTACCGCCGCCAGTGAGTTGCCGACCGCCAGGCCTAAATCGTCGTGGGTATGTACGGAAATAATGGCTTTGTCGATGTTAGGTACGCGtt                                                                                                                                                                                                                    <  1:270866/169‑1 (MQ=255)
                                                                                                         tgaCTTCTTCCAGGGAACAGTTTCCGGCACGCTCGCCGATCCCGTTCATTGCGCCTTCCACCTGGCGTGCACCGGCATGTACCGCCGCCAGTGAGTTGCCGACCGCCAGGCCTAAATCGTCGTGGGTATGTACGGAAATAATGGCTTTGTCGATGTTAGGTACGCGTTCATACAGGCCGCTGATGATTCCGGCGAACTCAAACGGCATGGTGTAGCCCACGGGGTCCGGAATGTTGATGGGGGTGGCACCg                                                                                                                                 >  1:149662/1‑251 (MQ=255)
                                                                                                                 tcCAGGGAACAGTTTCCGGCACGCTCGCCGATCCCGTTCATTGCGCCTTCCCCCTGGCGTGCCCCGGCATGTACCGCCGCCAGTGAGTTGCCGACCGCCAGGCCTAAATCGTCGTGGGTATGTACGGAAATAATGGCTTTGTCGATGTTAGGTACGCGTTCATACAGGCCGCTGATGATTCCGGCGAACTCAAACGGCATGGTGTAGCCCACGGTGTCCGGAATGTTGATGGTGGTGGCACCGGCATTa                                                                                                                           <  1:17695/249‑1 (MQ=255)
                                                                                                                    aGGGAACAGTTTCCGGCACGCTCGCCGATCCCGTTCATTGCGCCTTCCACCTGGCGTGCACCGGCATGTACCGCCGCCAGTGAGTTGCCGACCGCCAGGCCTAAATCGTCGTGGGTATGTACGGAAATAATGGCTTTGTCGATGTTAGGTACGCGTTCATACAGGCCGCTGATGATTCCGGCGAACTCAAACGGCATGGTGTAGCCCACGGTGTCCGGAATGTTGATGGTGGTGGCACCGGCATTAATcgc                                                                                                                      >  2:255779/1‑251 (MQ=255)
                                                                                                                     gggAACAGTTTCCGGCACGCTCGCCGATCCCGTTCATTGCGCCTTCCACCTGGCGTGCACCGGCATGTACCGCCGCCAGAGAGTTTCCAACCCCCAGGCCCAAATCGTCGTGGGTATGTCCGGAGATAATGGCTTTGTCGATGTTAGGCAAGCGTTCATAAAGGCCGCTGATGATTCCGGGGAACTCACACGGCCACGTGTAGCCCACCGGGTCCGGAATGTTGATGGTGGTGGCACCGGCATCAAtaccc                                                                                                                     >  2:38693/1‑247 (MQ=255)
                                                                                                                          cagttccgggcacccgcGCCGTTCCCGTGCATTGTGCCTCCCCCCGGGCGGGCACCGGCATGTACCGCCGCCAGTGAGTTGCCGACCGCCAGGCCTAAATCGTCGTGGGTATGTACGGAAATAATGGCTTTGTCGATGTTAGGTACGCGTTCATACAGGCCGCTGATGATTCCGGCGAACTCAAACGGCATGGTGTAGCCCACGGTGTCCGGAATGTTGATGGTGGTGGCACCGGCATTAATCGCCGCTTc                                                                                                                <  2:121077/235‑1 (MQ=255)
                                                                                                                                 cGGCCCGCCCGCCGATCCCGTTCATTGCGCCTGCCACCTGGCTTGCACCGGCATGTACCTCCGCCAGTGAGTTGCCGACCGCCAGGCCTAAATCGTCGTGGGTATGTACGGAAATAATGGCTTTGTCGATGTTAGGTACGCGTTCATACAGGCCGCTGATGATTCCGGCGAACTCAAACGGCATGGTGTAGCCCACGGTGTCCGGAATGTTGATGGTGGTGGCACCGGCATTAATCGCCGCTTCGACCACt                                                                                                         <  1:336200/251‑1 (MQ=255)
                                                                                                                                       acTCGCCGAACCCGTTCCTTCCGCCTTACACCTGGCGTGCACCGGCATGTACCGCCGCCAGTGAGTTGCCGGCCGCCAGGCCTCAATCGTCGTGGGTATGTACGGAAATAATGGCTTTATCGATTTTAGGTACGCGTTCATCCAGGCCGCTGATGGTTCCGGCGAACTCAACCCTCACTGTGTAGCCCACGGCTTCCGGCATGTTATTGGTGGAGGCCCCGGCATTacttgccccttctccccaccgcgcc                                                                                                   >  2:234619/2‑227 (MQ=255)
                                                                                                                                          cGCCGATCCCGTTCATTGCGCCTTCCACCTGGCGTGCACCGGCATGTACCGCCGCCAGTGAGTTGCCGACCGCCAGGCCTAAATCGTCGTGGGGATGTACGGAAATAATGGCTTTGTCGATGTTAGGTACGCGTTCATACAGGCCGCTGATGATTCCGGCGAACTCAAACGGCATGGTGTAGCCCACGGTGTCCGGAATGTTGATGGTGGTGGCACCGGCATTAATCGCCGCTTCGACCACTCGCGCCAgc                                                                                                >  2:60044/1‑250 (MQ=255)
                                                                                                                                            ccGATCCCGTTCATTGCGCCTTCCACCTGGCGTGCACCGGCATGTACCGCCGCCAGTGAGTTGCCGACCGCCAGGCCTAAATCGTCGTGGGTATGTACGGAAATAATGGCTTTGTCGATGTTAGGTACGCGTTCATACAGGCCGCTGATGATTCCGGCGAACTCAAACGGCATGGTGTAGcc                                                                                                                                                                   >  1:400171/1‑182 (MQ=255)
                                                                                                                                                        cgctccagtatccccgggCGTGCCCCGGAATGTACCGCCGCCAGTGAGTTGCCGCCCGCCAGGCCTAAATCGTCGTGGGTATGTACGGAAATAATGGCTTTGTCGATGTTAGGTACGCGTTCATACAGGCCGCTGATGATTCCGGCGAACTCAAACGGCATGGTGTAGCCCACGGTGTCAGGAATGTTGATGGTGGTGGCACCGGCATTAATCGCCGCTTCGACCACTCGCGCCAGATCGGCAATGGgt                                                                                    <  2:261220/233‑1 (MQ=255)
                                                                                                                                                            cgcCTTCCCCCTGGCGTGCACCGGCATGTACCGCCGACAGTGAGTTGCCGACCGCCAGGCCTAAATCGTCGTGGGTATGTACGGAAATAATGGCTTTGTCGATGTTAGGTACGCGTTCATACAGGCCGCTGATGATTCCGGCGAACTCAAACGGCATGGTGTAGCCCACGGTGTCCGGAATGTTGATGGTGGTGGCACCGGCATTAATCGCCGCTTCGACCACTCGCGCCAGATCGGCAATGGGTGTACGc                                                                              <  1:120663/251‑1 (MQ=255)
                                                                                                                                                             cccTTCCACCTGGCGTGCACCGGCATGTACCGCCGCCAGTGAGTTGCCGACCGCCAGGCCTAAATCGTCGTGGGTATGTACGGAAATAATGGCTTTGTCGATGTTAGGTACGCGTTCATACAGGCCGCTGATGATTCCGGCGAACTCAAACGGCATGGTGTAGCCCACGGTGTCCGGAATGTTGATGGTGGTGGCACCGGCATTAATCGCCGCTTCGACCACTCGCGCCAGATCGGCAATGGGTGTACGcc                                                                             <  1:419732/250‑1 (MQ=255)
                                                                                                                                                                 tCCACCTGGCGTGCACCGGCATGTACCGCCGCCAGTGAGTTGCCGACCGCCAGGCCTAAATCGTCGTGGGTATGTACGGAAATAATGGCTTTGTCGATGTTAGGTACGCGTTCATACAGGCCGCTGATGATTCCGGCGAACTCAAACGGCATGGTGTAGCCCACGGTGTCCGGAATGTTGATGGTGGTGGCACCGGCATTAATCGCCGCTTCGACCACTCGCGCCAGATCGGCAATGGGTGTACGCCCgg                                                                          <  1:109261/250‑1 (MQ=255)
                                                                                                                                                                    aCCTGGCGTGCACCGGCATGTACCGCCGCCAGTGAGTTGCCGACCGCCAGGCCTAAATCGTCGTGGGTATGTACGGAAATAATGGCTTTGTCGATGTTAGGTACGCGTTCATACAGGCCGCTGATGATTCCGGCGAACTCAAACGGCATGGTGTAGCCCACGGTGTCCGGAATGTTGATGGTGGTGGCACCGGCATTAATCGCCGCTTCGACCACTCGCGCCAGATCGGCAATGGgt                                                                                    >  1:261220/1‑237 (MQ=255)
                                                                                                                                                                         gCGTGCACCGGCATGTACCGCCGCCAGTGAGTTGCCGACCGCCAGGCCTAAATCGTCGTGGGTATGTACGGAAATAATGGCTTTGTCGATGTTAGGTACGCGTTCATACAGGCCGCTGATGATTCCGGCGAACTCAAACGGCATGGTGTAGCCCACGGTGTCCGGAATGTTGATGGTGGTGGCACCGGCATTAATCGCCGCTTCGACCACTCGCGCCAGATCGGCAATGGGTGTACGCCCGGCATCTTCGc                                                                 <  1:129832/251‑1 (MQ=255)
                                                                                                                                                                                 cGGCATGTACCGCCGCCAGTGAGTTGCCGACCGCCAGGCCTAAATCGTCGTGGGTATGTACGGAAATAATGGCTTTGTCGATGTTAGGTACGCGTTCATACAGGCCGCTGATGATTCCGGCGAACTCAAACGGCATGGTGTAGCCCACGGTGTCCGGAATGTTGATGGTGGTGGCACCGGCATTAATCGCCGCTTCGACCACTCGCGCCAGATCGGCAATGGTTGTTCGCCCGGCATCTTCGCAAGAAAAt                                                         <  1:181450/251‑1 (MQ=255)
                                                                                                                                                                                    cATGTACCGCCGCCAGTGAGTTGCCGACCGCCAGGCCTAAATCGTCGTGGGTATGTACGGAAATAATGGCTTTGTCGATGTTAGGTACGCGTTCATACAGGCCGCTGCTGATTCCGGCGAACTCAAACGGCATGGTGTAGCCCACGGTGTCCGGAATGTTGATGGTGGTGGCACCGGCATTAATCGCCGCTTCGACCACTCGCGCCAGATCGGCAATGGGTGTACGCCCGGCATCTTCGCAAGAAAATTc                                                       >  1:270910/1‑250 (MQ=255)
                                                                                                                                                                                    cATGTACCGCCGCCAGTGAGTTGCCGACCGCCAGGCCTAAATCGTCGTGGGTATGTACGGAAATAATGGCTTTGTCGATGTTAGGTACGCGTTCATACAGGCCGCTGATGATTCCGGCGAACTCAAACGGCATGGTGTAGCCCACGGTGTCCGGAATGTTGATGGTGGTGGCACCGGCATTAATCGCCGCTTCGACCACTCGCGCCAGATCGGCAATGGGTGTACGCCCGGCATCTTCGCAAGAAAATTc                                                       >  1:130258/1‑250 (MQ=255)
                                                                                                                                                                                    aaTGTACCGCCGACAGTGAGTTGCCGACCGCCAGGTCTAACTCGTAGTGGGTATGTACGGAGATACTGGCTTTGTCGCGGTTAGGTAAGCGTTCATACAGGCCGCTGATGATTCCGGCGAACTCAAAAGGCATGGTGTAGCACACGGACTCCGGAATGTTGCTGGAGGAGGAACCGGCATTAATCGACGCTTTGACCACTCGCGCCAGGTAGGCAATCGGTGTCCGCCCCGCATCTTCGCAAggaatttc                                                       >  2:231796/2‑243 (MQ=255)
                                                                                                                                                                                           cgccgcCAGTGAGTTGCCGACCGCCAGGCCTAAATCGTCTTGGGTATGTACGGAAATAATGGCTTTGTCGATGTTAGGTACGCGTTCATACAGGCCGCTGATGATTCCGGCGAACTCAAACGGCATGGTGTAGCCCACGGTGTCCGGAATGTTGATGGTGGTGGCACCGGCATTAATCGCCGCTTCGACCACTCGCGCCAGATCGGCAATGGGTGTACGCCCGGCATCTTCGCAAGAAAATTCAAcatcat                                               <  1:260268/251‑1 (MQ=255)
                                                                                                                                                                                                           ccGACCGCCAGGCCTAAATCGTCGTGGGTATGTACGGAAATAATGGCTTTGTCGATGTTAGGTACGCGTTCATACAGGCCGCTGATGATTCCGGCGAACTCAAACGGCATGGTGTAGCCCACGGTGTCCGGAATGTTGATGGTGGTGGCACCGGCATTAATCGCCGCTTCGACCACTCGCGCCAGATCGGCAATGGGTGTACGCCCGGCATCTTCGCAAGAAAATTCAACATCATCGGTGTAATTACGGgc                               >  1:275203/1‑251 (MQ=255)
                                                                                                                                                                                                           ccGACCGCCAGGACTAAATCGTCGAGGGTATGTACGGAAATAATGGCTTTGTCGATGTTAGGTACGCGTTCATACAGGCCGCTGATGATTCCGGCGAACTCAAACGGCATGGTGTAGCCCACGGTGTCAGGAATGTTGATGGTGGTGGCCCCGGCATTAATAGCCGCTTCGACCACTCGCGCCAGATCGGCAATGGGCGTACACCCGGCATCTTCGCAAGAAAATTCAACATCCACGGTGTAAATCCGGgc                               >  2:224246/1‑251 (MQ=255)
                                                                                                                                                                                                               ccaccAGGCCTAAATCGTCGTGGGTATGTACGGAAATAATGGCTTTGTCGATGTTAGGTACGCGTTCATACAGGCCGCTGATGATTCCGGCGAACTCAAACGGCATGGTGTAGCCCACGGTGTCCGGAATGTTGATGGTGGTGGCACCGGCATTAATCGCCGCTTCGACCACTCGCGCCAGATCGGCAATGGGTGTACGCCCGGCATCTTCGCAAGAAAATTCAACATCATCGGTGTAATTCCGGGCGCGt                           >  2:10241/4‑251 (MQ=255)
                                                                                                                                                                                                                          aaaTCGTCGTGGGTATGTACGGAAATAATGGCTTTTTCGTTGTTAGGTACGCGTTCATACAGGCCGCTGATGATTCCGGCGAACTCAAACGGCATGGTGTAGCCCACGGTGTCCGTAATGTTGATGGTGGTGGCACCGGCATTAATCGCCGCTTCGACCACTCGCGCCAGATCGGCAATGGGTGTACGCCCGGCATCTGCGCAAGAAAATTCAACATCATCGGTGTAATTACGGGCGCGTTTCACCatat                 <  1:454140/250‑1 (MQ=255)
                                                                                                                                                                                                                              cgtcgtGGGTATGTACGGAAATAATGGCTTTGTCGATGTTAGGTACGCGTTCATACAGGCCGCTGATGATTCCGGCGAACTCAAACGGCATGGTGTAGCCCACGGTGTCCGGAATGTTGATGGTGGTGGCACCGGCATTAATCGCCGCTTCGACCACTCGCGCCAGATCGGCAATGGGTGTACGCCCGGCATCTTCGCAAGAAAATTCAACATCATCGGTGTAATTACGGGCGCGTTTCACCATatagat             <  1:348370/250‑1 (MQ=255)
                                                                                                                                                                                                                                    gggTATGTACGGAAATAATGGCTTTGTCGATGTTAGGTACGCGTTCATACAGGCCGCTGATGATTCCGGCGAACTCAAACGGCATGGTGTAGCCCACGGTGTCCGGAATGTTGATGGTGGTGGCACCGGCATTAATCGCCGCTTCGACCACTCGCGCCAGATCGGCAATGGGTGTACGCCCGGCATCTTCGCAAGAAAATTCAACATCATCGGTGTAATTACGGGCGCGTTTCACCATATAGATAGCGCGt      <  1:179173/251‑1 (MQ=255)
                                                                                                                                                                                                                                     ggTATGTACGGAGATAATGGCTTTGTCGATGTTAGGCACGCGTTCATACAGGCCGCTGATGATTCCGGCGAACTCAAACGGCATGGTGTAGCCCACGGTGTCCGGAATGTTGATGGTGGTGGCACCGGCATTAATCGCCGCTTCGACCACTCGCGCCAGATCGGCAATGGGTGTACGCCCGGCATCTTCGCAAGAAAATTCAACATCATCGGTGTAATTACGGGCGCGTTTCACCATATAGATAGCGCGtt     >  1:229067/1‑251 (MQ=255)
                                                                                                                                                                                                                                        aTGTACGGAAATAATGGCTTTGTCGATGTTAGGTACGCGTTCATACAGGCCGCTGATGATTCCGGCGAACTCAAACGGCATGGTGTAGCCCACGGTGTCCGGAATGTTGATGGTGGTGGCACCGGCATTAATCGCCGCTTCGACCACTCGCGCCAGATCGGCAATGGGTGTACGCCCGGCATCTTCGCAAGAAAATTCAACATCATAGGTGTAATTACGGGCGCGTTTCACCATATAGATAGCGCGTTCg   >  1:88233/1‑250 (MQ=255)
                                                                                                                                                                                                                                         tGTACGGAAATAATGGCTTTGTCGATGTTAGGTACGCGTTCATACCAGCCGCTGATGAGTCCGGCGAAATCAAACGGCATGATGTAGCCCACGGTGGCCGGAGTGTTGATGGGGGGGGCAGCGGCATTAATCGCCGCGTCGACCACTCGCGCCAGATCGGCAATGGGTGAACGCCCGGCATCTTCGCAAGAAAATTCAACATCATCGGTGTAATTACGGGCGCGTTTCCCCATATAGATAGCGCGGTCgc  >  1:382/1‑249 (MQ=255)
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                 
GCATATTACAAATCTGGCTAACTAACTGGCTGGTGCGCCATATCTCCTGGTGATTAATGGCGGTGTGGACGTTGAGAATATCCTTACGAACTTTGATCGCCATGATGACTTCTTCCAGGGAACAGTTTCCGGCACGCTCGCCGATCCCGTTCATTGCGCCTTCCACCTGGCGTGCACCGGCATGTACCGCCGCCAGTGAGTTTCCGACCGCCAGGCCCAAATCGTCGTGGGTATGTACGGAGATAATGGCTTTGTCGATGTTAGGCACGCGTTCATACAGGCCGCTGATGATTCCGGCGAACTCAAACGGCATGGTGTAGCCCACGGTGTCCGGAATGTTGATGGTGGTGGCACCGGCATTAATCGCCGCTTCGACCACTCGCGCCAGATCGGCAATGGGTGTACGCCCGGCATCTTCGCAAGAAAATTCAACATCATCGGTGTAATTACGGGCGCGTTTCACCATATAGATAGCGCGTTCGA  >  NC_000913/82692‑83174

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