Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 82,894 T→G G212G (GGA→GGC leuA ← 2‑isopropylmalate synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091382,8940TG92.5% 101.8 / ‑4.5 40G212G (GGA→GGCleuA2‑isopropylmalate synthase
Reads supporting (aligned to +/- strand):  ref base T (3/0);  new base G (16/21);  total (19/21)
Fisher's exact test for biased strand distribution p-value = 9.81e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAATGCGCCGCTGCCAACAATGGCTTTGTTTGCCGGGATCGGCATATTACAAATCTGGCTAACTAACTGGCTGGTGCGCCATATCTCCTGGTGATTAATGGCGGTGTGGACGTTGAGAATATCCTTACGAACTTTGATCGCCATGATGACTTCTTCCAGGGAACAGTTTCCGGCACGCTCGCCGATCCCGTTCATTGCGCCTTCCACCTGGCGTGCACCGGCATGTACCGCCGCCAGTGAGTTTCCGACCGCCAGGCCCAAATCGTCGTGGGTATGTACGGAGATAATGGCTTTGTCGATGTTAGGCACGCGTTCATACAGGCCGCTGATGATTCCGGCGAACTCAAACGGCATGGTGTAGCCCACGGTGTCCGGAATGTTGATGGTGGTGGCACCGGCATTAATCGCCGCTTCGACCACTCGCGCCAGATCGGCAATGGGTGTACGCCCGGCATCTTCGCAAGAAAATTCAACATCATCGGTGTAAT  >  NC_000913/82651‑83138
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                       
gAATGCGCCGCTGCCAACAATGGCTTTGTTTGCCGGGATCGGCATATTACAAATCTGGCTAACTAACTGGCTGGTGCGCCATATCTCCTGGTGATTAATGGCGGTGTGGACGTTGAGCATATCCTTACGAACTTTGATCGCCCTGATGACTTCTTCCAGGGGAAAGTTTCCGGCACGCTCGCCCCTCCCGTTCATTTCGCCCTCCACCTGGCGTGCACCGGCATGTACCCCCGCCCGTGAGGTTcccagcc                                                                                                                                                                                                                                                  >  2:344818/1‑246 (MQ=255)
            gCCAACAATGGCTTTGTTTGCCGGGATCGGCATATTACAAATCTGGCTAACTAACTGGCTGGTGCGCCATATCTCCTGGTGATTAATGGCGGTGTGGACGTTGAGAATATCCTTACGAACTTTGATCGCCATGATGACTTCTTCCAGGGAACAGTTTCCGGCACGCTCGCCGATCCCGTTCATTGCGCCTTCCACCTGGCGTGCACCGGCATGTACCGCCGCCAGTGAGTTGCCGACCGCCAGGcctaa                                                                                                                                                                                                                                       <  1:366690/249‑4 (MQ=255)
                   aTGGCTTTGTTTGCCGGGATCGGCATATTACAAATCTGGCTAACTAACTGGCTGGTGCGCCATATCTCCTGGTGATTAATGGCGGTGTGGACGTTGAGAATATCCTTACGAACTTTGATCGCCATGATGACTTCTTCCAGGGAACAGTTTCCGGCACGCTCGCCGATCCCGTTCATTGCGCCTTCCACCTGGCGTGCACCGGCATGTACCGCCGCCAGTGAGTTGCCGACCGCCAGGCCTAAATCGTCGTg                                                                                                                                                                                                                             <  1:49814/251‑1 (MQ=255)
                     ggCTTTGTTTGCCGGGATCGGCATATTACAAATCTGGCTAACTAACTGGCTGGTGCGCCATATCTCCTGGTGATTAATGGCGGTGTGGACGTTGAGAATATCCTTACGACCTTTGATCGCCATGATGACTTCTTCCAGGGAACAGTTTCCGGCACGCTCGCCGATCCCGTTCATTGCGCCTTCCACCTGGCGTGCACCGGCATGTACCGCCGCCAGTGAGTTGCCGACCGCCAGGCCTAAATCGTCGTggg                                                                                                                                                                                                                           <  1:81736/251‑1 (MQ=255)
                     ggCTTTGTTTGCCGGGATCGGCATATTACAAATCTGGCTAACTAACTGGCTGGTGCGCCATATCTCCTGGTGATTAATGGCGGTGTGGACGTTGAGAATATCCTTACGAACTTTGATCGCCATGATGACTTCTTCCAGGGAACAGTTTCCGGCACGCTCGCCGATCCCGTTCATTGCGCCTTCCACCTGGCGTGCACCGGCATGTACCGCCGCCAGTGAGTTGCCGACCGCCAGGCCTAAATCGTCGTggg                                                                                                                                                                                                                           <  1:370843/251‑1 (MQ=255)
                                         gCATATTACAAATCTGGCTAACTAACTGGCTGGTGCGCCATATCTCCGGGTGATTAATGGCGGTGTGGACGTTGAGAATATCCTTACGAACTTTGATCGCCATGATGACTTCTTCCAGGGAACAGTTTCCGGCACGCTCGCCGATCCCGTTCATTGCGCCTTCCACCTGGCGTGCACCGGCATGTACCGCCGCCAGTGAGTTGCCGACCGCCAGGCCTAAATCGTCGTGGGTATGTACGGAAATAATGGCt                                                                                                                                                                                                      <  2:165703/251‑1 (MQ=255)
                                            tatTACAAATCTGGCTAACTAACTGGCTGGTGCGCCATATCTCCTGGTGATTAATGGCGGCGTGGACGTTGAGAATATCCTTACGAACTTTGATCGCCATGATGGCTTCTTCCAGGGAACAGTTTCCGGCACTCTCGCCGATCCCGTTCATTGCGCCTTCCACCTGGCGTGCACCGGCCTGTACCCCCGCCAGTGAGTTGCCGACCGCCCGGGCTAAATCGACGTGGGTATTTACGGAAATAAAGGCtttt                                                                                                                                                                                                   >  2:199346/1‑250 (MQ=255)
                                                  aaaTCTGGCTAACTAACTGGCTGGTGCGCCATATCTCCTGGTGATTAATGGCGGTGTGGACGTTGAGAATATCCTTACGAAATTTGATCGCCATGATGACTTCTTCCAGGGAACAGTTTCCGGCCCGCTCGCCGATCCCGTTCATTGCGCCTTCCACCTGGCGGGCACCGGCATGTACCCCCCCCAGTGAGTTGCCGACCGCCCGGCCTCAAACGTCGTGGGTATGTACGGAAATAATGGCTTTGTCgctg                                                                                                                                                                                             >  2:83007/1‑248 (MQ=255)
                                                   aaTCTGGCTAACTAACTGGCTGGTGCGCCATATCTCCTGGTGATTAATGGCGGTGTGGACGTTGAAAAAATCCTTACGAACTTTGATCGCCATGATGACTTCTTCCCGGGAACAGTTTCCGGCACGCTCGCCGATCCCGTTCATTGCGCCTTCCACCTGGCGGGCACCGGCCTGTACCGGCGCCAGTGGGTTTCCGACCGCCAGGGCTAAATCGTCGCGGGTATGTACGgaaaaaaaggtcttgttgatgt                                                                                                                                                                                             >  1:190223/1‑231 (MQ=255)
                                                                       tggTGCGCCATATCTCCTGGTGATTAATGGCGGTGTGGACGTTGAGAATATCCTTACGAACTTTGATCGCCATGATGACTTCTTTCCGGGAACAGTTTCCGGCACGCTCGCCGATCCCGTTCATTTCGCCTTCCCCCTGGCGTGCACCGGCATGTACCGCCGCCAGTGAGGTGCCGACCGCCAGGCCTAAATCGTCGGGGGTATGTACGGAAATAATGGCTTTGTCGGTGTTTGGTACGCGTTCATACAgg                                                                                                                                                                        >  1:164126/1‑251 (MQ=255)
                                                                        ggTGCGCCATATCTCCTGGTGATTAATGGCGGTGTGGACGTTGAGAATATCCTTACGAACTTTGATCGCCATGATGACTTCTTCCAGGGAACAGTTTCCGGCACGCTCGCCGATCCCGTTCATTGCGCCTTCCACCTGGCGTGCACCGGCATGTACCGCCGCCAGTGAGTTGCCGACCGCCAGGCCTAAATCGTCGTGGGTATGTACGGAAATAATGGCTTTGTCGATGTTAGGTACGCGTTCATACAg                                                                                                                                                                         <  1:36958/249‑1 (MQ=255)
                                                                            cgcCATATCTCCTGGTGATTAATGGCGGTGTGGTCGTTGGGAATATCCCTACGCACTTTGATCGCCATGATGACTTCTTTCCGGGGAACGGTTCCCGCACGCCCGCCGATCCCGTTCATTGCGCCTTCCACCTGGCGTGCCCCGGCATGTACCGCCGCCCGTGAGTTGCCCACCGCCCGGCCTAAATCGCCGTGGGGCTGTACGGAAATCATGGCCTTGTCGGTGTTAAGGACGCGGTCTTAAAGGCCGcc                                                                                                                                                                   >  2:25845/1‑250 (MQ=255)
                                                                             gcCATATCTCCTGGTGATTAATGGCGGTGTGGACGTTGAGAATATCCTTACGAACTTTGATCGCCATGATGACTTCTTCCAGGGAACAGTTTCCGGCACGCTCGCCGATCCCGTTCATTGCGCCTTCCACCTGGCGTGCACCGGCATGTACCGCCGCCAGTGAGTTGCCGACCGCCCGGCCTAAATCGTCGTGGGGATGTACGGAAATAATGGCTTTTTCGATGTTAGGTACGCGTTCATACAGGCCGCtt                                                                                                                                                                  >  1:373588/1‑250 (MQ=255)
                                                                              ccATATCTCCTGGTGATTAATGGCGGTGTGGACGTTGAGAATATCCTTACGCACTTTGATCGCCATGATGACTTCTTCCAGGGAACAGTTTCCGGCACGCTCGCCGATCCCGTTCATTGCGCCTTCCACCTGGCGTGCACCGGCATGTACCGCCGCCCGTGAGTTGCCGACCGCCCGGCCTAAATCGTCGTGGGTATGTACGGAAATAATGGCTTTGTCGATGTTAGGTACGCGTTCATACAGGCCGCtgg                                                                                                                                                                 >  2:70970/1‑250 (MQ=255)
                                                                                atatCTCCTGGTGATTAATGGCGGTGTGGACGTTGAGAATATCCTTACGAACTTTGATCGCCATGATGACTTCTTCCAGGGGACAGGTTCCGGCACGCTCGCCGATCCCGTTCATTGCGCCTTCCACCTGGCGTGCACCGGACTGTACCGCCGCCAGTGAGGTGGCGACCGCCCGGCCTAAATCGTCGTGGGCATGTACGGGAAACATTGCTTTGTCGATGTTAGGGACGCGTTCATACAGGCCGCtgatg                                                                                                                                                               >  1:370728/1‑251 (MQ=255)
                                                                                  atCTCCTGGTGATTAATGGCGGTGTGGACGTTGAGAATATCCTTACGAACTTTGATCGCCATGATGACTTCTTCCAGGGAACAGTTTCCGGCACGCTCGCCGATCCCGTTCATTGCGCCTTCCACCTGGCGTGCACCGGCATGTACCGCCGCCAGTGAGTTGCCGACCGCCAGGCCTAAATCGTCGTGGGTATGTACGGAAATAATGGCTTTGTCGATGTTAGGTACGCGTTCATACAGGCCGCtgatga                                                                                                                                                              <  1:70970/250‑1 (MQ=255)
                                                                                      ccTGGTGATTAATGGCGGTGGGGACGTTGAGAATATCCTTACGCACTTTGATCGCCATGATGACTTCTTCCAGGGAACAGTTTCCGGCCCGCTCGCCGATCCCGTTCCTTGCGCCCTCCCCCTGGCGTGCACCCGCATGTACCCCCGCCAGGGAGGTTCCGACCCCCACGGCTCAATCCCCGGGGGTTTTTACGGGAATAAAGGCTTTGTCGGTGTTTGGTACGCGTTCAtcccggcccctcgtggtttcc                                                                                                                                                         >  2:309971/1‑231 (MQ=255)
                                                                                         ggTGATTAATGGCGGTGTGGACGTTGAGAATATCCTTACGAACTTTGATCGCCATGATGACTTCTTCCAGGGAACAGTTTCCGGCACGCTCGCCGATCCCGTTCATTGCGCCTTCCACCTGGCGTGCACCGGCATGTACCGCCGCCAGTGAGTTGCCGACCGCCAGGCCTAAATCGTCGTGGGGATGTACGGAAATAATGGCTTTGTCGATGTTAGGTACGCGTTCATAAAGGCCGCTGATGATTCCGGCg                                                                                                                                                      >  1:105114/1‑251 (MQ=255)
                                                                                                     cGGTGTGGACGTTGAGAATATCCTTACGAACTTTGATCGCCATGATGACTTCTTCCAGGGAACAGTTTCCGGCACGCTCGCCGATCCCGTTCATTGCGCCTTCCACCTGGCGTGCACCGGCATGTACCGCCGCCAGTGAGTTGCCGACCGCCAGGCCTAAATCGTCGTGGGTATGTACGGAAATAATGGCTTTGTCGATGTTAGGTACGCGTTCATACAGGCCGCTGATGATTCCGGCGAACTCAAACGGc                                                                                                                                          >  1:366625/1‑251 (MQ=255)
                                                                                                         gtgGACGTTGAGAATATCCTTACGAACTTTGATCGCCATGATGACTTCTTCCAGGGAACAGTTTCCGGCACGCTCGCCGATCCCGTTCATTGCGCCTTCCACCTGGCGTGCACCGGCATGTACCGCCGCCAGTGAGTTGCCGACCGCCAGGCCTAAATCGTCGTGGGTATGTACGGAAATAATGGCTTTGTCGATGTTAGGTACGCGTTCATACAGGCCGCTGATGATTCCGGCGAACTCAAACGGCATgg                                                                                                                                      <  1:344818/251‑1 (MQ=255)
                                                                                                                      atatCCTTACGAACTTTGATCGCCATGATGACTGCTTCCAGGGAACAGTTTCCGGCACGCTCGCCGATCCCGTTCATTGCGCCGTCCACCTGGCGGGCACCGGCATGTACCGCCGCCAGTGAGTTGCCGACCGCCAGGCCTAAATCGTCGTGGGTATGTACGGAAATAATGGCTTTGTCGATGTTAGGTACGCGTTCATACAGGCCGCTGATGATTCCGGCGAACTCAAACGGCATGGTGTAGCCCACGgt                                                                                                                         <  1:343265/251‑1 (MQ=255)
                                                                                                                                        aTCGCCATGATGACTTCTTCCAGGGAACAGTTTCCGGCACGCTCGCCGATCCCGTTCATTGCGCCTTCCACCTGGCGTGCACCGGCATGTACCGCCGCCAGTGAGTGGCCGACCGCCAGGCCTAAATCGTCGTGGGTATGTACGGAAATAATGGCTTTGTCGATGTTAGGTACGCGTTCATACAGGCCGCTGATGATTCCGGCGAACTCAAACGGCATGGTGTAGCCCACGGTGTCCGGAATGTTGAtggt                                                                                                       >  2:45009/1‑251 (MQ=255)
                                                                                                                                               tgatgaCTTCTTCCAGGGAACAGTTTCCGGCACGCTCGCCGATCCCGTTCATTGCGCCTTCCACCTGGCGTGCACCGGCATGTACCGCCGCCAGTGAGTTGCCGACCGCCAGGCCTAAATCGTCGTGGGTATGTACGGAAATAATGGCTTTGTCGATGTTAGGTACGCGTTCATACAGGCCGCTGATGATTCCGGCGAACTCAAACGGCATGGTGTAGCCCACGGTGTCCGGAATGTTGATGGTGGTGGc                                                                                                 >  1:280987/1‑250 (MQ=255)
                                                                                                                                                   gaCTTCTTCCAGGGAACAGTTTCCGGCACGCTCGCCGATCCCGTTCATTGCGCCTTCCACCTGGCGTGCACCGGCATGTACCGCCGCCAGTGAGTTGCCGACCGCCAGGCCTAAATCGTCGTGGGTATGTACGGAAATAATGGCTTTGTCGATGTTAGGTACGCGTTCATACAGGCCGCTGATGATTCCGGCGAACTCAAACGGCATGGTGTAGCCCACGGTGTCCGGAATGTTGATGGTGGTGGCACCg                                                                                             <  1:136203/250‑1 (MQ=255)
                                                                                                                                                            ccAGGGAACAGTTTCCGGCACGCACGCCGATCCCGTTCATGGCGCCTTCCGCCTGGCGTGCACCGGCATGTACCGCCGCCAGTGAGTTGCCGACCGCCAGGCCTAAATCGTCGTGGGTATGTACGGAAATAATGGCTTTGTCGATGTTAGGTACGCGTTCATACAGGCCGCTGATGATTCCGGCGAACTCAAACGGCATGGTGTAGCCCACGGTGTCCGGAATGTTGATGGTGGTGGCACCGGCATTAATc                                                                                    <  2:282119/251‑1 (MQ=255)
                                                                                                                                                               gggAACAGTTTCCGGAACGCTCGCCGATCCCGTTCATTGCGCCTTCCCCCTGGAGTGCACCGGCATGTACCGCCGCCAGTGAGTTGCCGACCGCCAGGCCTAAATCGTCGTGGGTATGTACGGAAATAATGGCTTTGTCGATGTTAGGTACGCGTTCATACAGGCCGCTGATGATTCCGGCGAACTCAAACGGAATGGTGTAGCCCACGGTGTCCGGAATGTTTATGGTGGTGGCACCGGCATTAATcgcc                                                                                 >  2:182242/1‑251 (MQ=255)
                                                                                                                                                                                 gCTCGCCGATCCCGTTCATTGCGCCTTCCACCTGGCGTGCACCGGCATGTACCGCCGCCAGTGAGTTGCCGACCGCCAGGCCTAAATCGTCGTGGGTATGTACGGAAATAATGGCTTTGTCGATGTTAGGTACGCGTTCATACAGGCCGCTGATGATTCCGGCGAACTCAAACGGCATGGTGTAGCCCACGGTGTCCGGAATGTTGATGGTGGTGGCACCGGCATTAATCGCCGCTTCGACCACTCGCGcc                                                               >  2:310796/1‑251 (MQ=255)
                                                                                                                                                                                     cccgctCCCGTTCATTGCGCCTTCCACCTGGCGTGCACCGGCAGGTCCCGCCGCCAGTGAGTTGCCGACCGCCAGGCCTAAAGCGTCGTGGGTATGTACGGAAATAATGGCTTTGTCGATGTTAGGTACGCGTTCATACAGGCCGCTGATGATTCCGGCGAACTCAAACGGCATGGTGTAGCCCACGGTGTCCGGAATGTTGATGGTGGTGGCACCGGCATTAATCGCCGCTTCGACCACTCGCGCCAGAt                                                           <  2:320886/246‑1 (MQ=255)
                                                                                                                                                                                        gATCCCGTTCATTGCGCCTTCCACCTGGCGTGCACCGGCATGTACCGCCGCCAGTGAGTTGCCGACCGCCAGGCCTCAATCGTCGTGGGTATGTACGGAAATAATGGCTTTGTCGATGTTAGGTACGCGTTCATACAGGCCGCTGATGATTCCGGCGAACTCAAACGGCATGGTGTAGCCCACGGTGTCCGGAATGTTGATGGTGGTGGCACCGGCATTAATCGCCGCTTCGACCACTCGCGcc                                                               >  2:184638/1‑244 (MQ=255)
                                                                                                                                                                                                 cATTGCGCCTTCCACCGGGCGTGCACCGGCATGTACCGCCGCCAGTGAGTTGCCGACCGCCAGGCCTAAATCGTCGTGGGTATGTACGGAAATAATGGCTTTGTCGATGTTAGGTACGCGTTCATACAGGCCGCTGATGATTCCGGCGAACTCAAACGGCATGGTGTAGCCCACGGTGTCCGGAATGTTGATGGTGGTGGCACCGGCATTAATCGCCGCTTCGACCACTCGCGCCAGATCGGCAATGGgt                                                <  2:107245/250‑1 (MQ=255)
                                                                                                                                                                                                         cTTCCATCTGGCGTGCACCGGCATGTACCGCCGCCAGTGAGTTGCAGACCGCAAGGCCTAAATCGTCGTGGGTATGTACGGAAATAATGGCTTTGTCGATGTTAGGTACGCGTTCATACAGGCCGCTGATGATTCCGGCGAACTCAAACGGCATGGTGTAGCCCACGGTGTCCGGAATGTTGATGGTGGTGGCACCGGCATTAATCGCCGCTTCGACCACTcgcg                                                                 >  1:306076/1‑225 (MQ=255)
                                                                                                                                                                                                         cTTCCACCTGGCGTGCACCGGCATGGACCGCCGCCAGTGAGTTGCCGACCGCCAGGCCTAAATCGTCGTGGGTATGTACGGAAATAATGGCTTTGTCGATGTTAGGTACGCGTTCATACAGGCCGCTGATGATTCCGGCGAACTCAAACGGCATGGTGTAGCCCACGGTGTCCGGAATGTTGATGGTGGTGGCACCGGCATTAATCGCCGCTTCGACCACTcgcg                                                                 <  2:306076/225‑1 (MQ=255)
                                                                                                                                                                                                           tCCACCTGGCGTGCACCGGCATGTACCGCCGCCAGTGAGTTGCCGACCGCCAGGCCTAAATCGTCGTGGGTATGTACGGAAATAATGGCTTTGTCGATGTTAGGTACGCGTTCATACAGGCCGCTGATGATTCCGGCGAACTCAAACGGCATGGTGTAGCCCACGGTGTCCGGAATGTTGATGGTGGTGGCACCGGCATTAATCGCCGCTTCGACCACTCGCGCCAGATCGGCAATGGGTGTACGCCCgg                                      <  2:373588/250‑1 (MQ=255)
                                                                                                                                                                                                                       gCACCGGCATGTACCGCCGCCAGTGAGTTGCCGACCGCCAGGCCTAAATCGTCGTGGGTATGTACGGAAATAATGGCTTTGTCGATGTTAGGTACGCGTTCATACAGGCCGCTGATGATTCCGGCGAACTCAAACGGCATGGTGTAGCCCACGGTGTCCGGAATGTTGATGGTGGTGGCACCGGCATTAATCGCCGCTTCGACCACTCGCGCCAGATCGGCAATGGGTGTACGCCCGGCATCTTCGCAAGa                         <  1:294818/251‑1 (MQ=255)
                                                                                                                                                                                                                          ccGGCATGTACCGCCGCCAGTGAGTTGCCGACCGCCAGGCCTAAATCGTCGTGGGTATGTACGGAAATAATGGCTTTGTCGATGTTAGGTACGCGTTCATACAGGCCGCTGATGATTCCGGCGAACTCAAACGGCATGGTGTAGCCCACGGTGTCCGGAATGTTGATGGTGGTGGCACCGGCATTAATCGCCGCTTCGACCACTCGCGCCAGATCGGCAATGGGTGTACGCCCGGCATCTTCGCAAGaaaa                      <  1:309971/251‑1 (MQ=255)
                                                                                                                                                                                                                               aTGTACCGCCGCCAGTGAGTTGCCGACCGCCAGGCCTAAATCGTCGTGGGTATGTACGGAAATAATGGCTTTGTCGATGTTAGGTACGCGTTCATACAGGCCGCTGATGATTCCGGCGAACTCAAACGGCATGGTGTAGCCCACGGTGTCCGGAATGTTGATGGTGGTGGCACCGGCATTAATCGCCGCTTCGACCACTCGCGCCAGATCGGCAATGGGTGTACGCCCGGCATCTTCGCAAGAAAATTCaa                 <  1:91341/251‑1 (MQ=255)
                                                                                                                                                                                                                                        cgccAGTGAGTTGCCGACCGCCAGGCCTAAATCGTCGTGGGTATGTACGGAAATAATGGCGTTGTCGATGTTAGGTACGCGTTCATACAGGCCGCTGATGATTCCGGCGAACTCAAACGGCATGGTGTAGCCCACGGTGTCCGGAATGTTGATGGTGGTGGCACCGGCATTAATCGCCGCTTCGACCACTCGCGCCAGATCGGCAATGGGTGTACGCCCGGCATCTTCGCAAGAAAATTCAACATCATCgg        <  2:72167/251‑1 (MQ=255)
                                                                                                                                                                                                                                         gccAGTGAGTTGCCGACCGCCAGGCCTAAATCGTCGTGGGTATGTACGGAAATAATGGCTTTGTCGATGTTAGGTACGCGTTCATACAGGCCGCTGATGATTCCGGCGAACTCAAACGGCATGGTGTAGCCCACGGTGTCCGGAATGTTGATGGTGGTGGCACCGGCATTAATCGCCGCTTCGACCACTCGCGCCAGATCGGCAATGGGTGTACGCCCGGCATCTTCGCAAGAAAATTCAACATCATCGgt       <  1:56540/251‑1 (MQ=255)
                                                                                                                                                                                                                                            aGTGAGTTGCCGACCGCCAGGCCTAAATCGTCGTGGGTATGTACGGAAATAATGGCTTTGTCGATGTTAGGTACGCGTTCATACAGGCCGCTGATGATTCCGGCGAACTCAAACGGCATGGTGTAGCCCACGGTGTCCGGAATGTTGATGGTGGTGGCACCGGCATTAATCGCCGCTTCGACCACTCGCGCCAGATCGGCAATGGGTGTACGCCCGGCATCTTCGCAAGAAAATTCAACATCATCGGTGTa    <  1:184638/251‑1 (MQ=255)
                                                                                                                                                                                                                                              tGAGTTGCCGACCGCCAGGCCTAAATCGTCGTGGGTATGTACGGAAATAATGGCTTTGTCGATGTTAGGTACGCGTTCATACAGGCCGCTGATGATTCCGGCGAACTCAAACGGCATGGTGTAGCCCACGGTGTCCGGAATGTTGATGGTGGTGGCACCGGCATTAATCGCCGCTTCGACCACTCGCGCCAGATCGGCAATGGGTGTACGCCCGGCATCTTCGCAAGAAAATTCAACATCATCGGTGTAAt  >  1:113182/1‑251 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                       
GAATGCGCCGCTGCCAACAATGGCTTTGTTTGCCGGGATCGGCATATTACAAATCTGGCTAACTAACTGGCTGGTGCGCCATATCTCCTGGTGATTAATGGCGGTGTGGACGTTGAGAATATCCTTACGAACTTTGATCGCCATGATGACTTCTTCCAGGGAACAGTTTCCGGCACGCTCGCCGATCCCGTTCATTGCGCCTTCCACCTGGCGTGCACCGGCATGTACCGCCGCCAGTGAGTTTCCGACCGCCAGGCCCAAATCGTCGTGGGTATGTACGGAGATAATGGCTTTGTCGATGTTAGGCACGCGTTCATACAGGCCGCTGATGATTCCGGCGAACTCAAACGGCATGGTGTAGCCCACGGTGTCCGGAATGTTGATGGTGGTGGCACCGGCATTAATCGCCGCTTCGACCACTCGCGCCAGATCGGCAATGGGTGTACGCCCGGCATCTTCGCAAGAAAATTCAACATCATCGGTGTAAT  >  NC_000913/82651‑83138

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