Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 82,957 C→T V191V (GTG→GTA)  leuA ← 2‑isopropylmalate synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091382,9570CT95.2% 166.1 / NA 66V191V (GTG→GTA) leuA2‑isopropylmalate synthase
Reads supporting (aligned to +/- strand):  ref base C (0/1);  major base T (34/26);  minor base G (3/0);  total (39/27)
Fisher's exact test for biased strand distribution p-value = 2.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGCTGGTGCGCCATATCTCCTGGTGATTAATGGCGGTGTGGACGTTGAGAATATCCTTACGAACTTTG‑ATCGCCATGATGACTTCTTCCAGGGAACAGTTTCCGGCACGCTCGCCGATCCCGTTCATT‑GCGCCTTCCACCTGGCGTGCACCGGCATGTACCGCCGCCAGTGAGTTTCCGACCGCCAGGCCCAAATCGTCGTGGGTATGTACGGAGATAATGGCTTTGTCGATGTTAGGCACGCGTTCATACAGGCCGCTGATGATTCCGGCGAACTCAAACGGCATGGTGTAGCCCACGG‑TGTCCGGAATGTTGATGGTGGTGGCACCGGCATTAATCGCCGCTTCGACCACTCGCGCCAGATCGGCAATGGGTGTACGCCCGGCATCTTCGCAAGAAAATTCAACATCATCGGTGTAATTACGGGCGCGTTTCACCATATAGATAG‑CGCGTTCGATCACCTCGTCCAGCGTGCTGCGCAGCTTGGTGG  >  NC_000913/82719‑83207
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                            
ggctggTGCGCCATATCTCCTGGTGATTAATGGCGGTGTGGACGTTGAGAATATCCTTACGAACTTTG‑ATCGCCATGATGACTTCTTCCAGGGAACAGTTTCCGGCACGCTCGCCGATCCCGTTCATT‑GCGCCTTCCACCTGGCGTGCACCGGCATGTACCGCCGCCAGTGAGTTGCCGACCGCCAGGCCTAAATCGTCGTGGGTATGTACGGAAATAATGGCTTTGTCGATGTTAGGTACGCGTTCAt                                                                                                                                                                                                                                                    <  1:264602/249‑1 (MQ=255)
ggcgtgTGCGCCATACCTCCTGGTTATTAATGGGGGGGTGGGCGTTGAGAATATCCTTACGATCTTTG‑ACCGCCAAGATGACTTCTTCTAGGGACCAGTTTCCGGCACGCTCGCCGATCCCGTTCATG‑GCGCCTTCCCCCTGGCGTGCACCGGCATGTACCGCCGCCAGTGAGTTGCCGACCGCCAGGCCTAAATCGTCGTGGGTATGTACGGAAATAATGGCTTTGTCTATGTTAGGTACGCGTTCAt                                                                                                                                                                                                                                                    <  2:149375/244‑1 (MQ=255)
  ctggTGCGCCATATCTCCTGGTGATTAATGGCGGTGTGGAAGTTGAGAATATCCTTACGAACTTTT‑ATCGCCATGATGCCTTCTTCCAGGGAACAGTTTCCGGCACGCTCGCCGATCCCGGTCATT‑GCGCCTTCCACCTGGCGTGCACCGGCATGTACCGCCGCCACTGAGGTGCCGACCGCCAGGCCTAAATCCTCGTGGGTCTGTTCGGAAAAAATGGCTTTTTCGATGTTGGGGACGCGTTCATACa                                                                                                                                                                                                                                                 >  1:327290/1‑250 (MQ=255)
         gcCATATCTCCTGGTGATTAATGGCGGTGGGGACGTTGAGAATATCCTTACGAACTTTG‑ATCGCCATGATGACTTCTTCCAGGGAACAGTTTCCGGCACGCTCGCCGATACCGTTCATT‑GCGCCTTCCACCTGGCGTGCCCCGGCATGTAACGCCGCCAGTGAGTTGCCGACCGCCCGGCCTAAATCGGCGTGGGTATGTCCGGAAAAAATGGCTTTGTCGATGGTAGGGACGCGTTCATACAGGCCGCtg                                                                                                                                                                                                                                         >  1:592660/1‑251 (MQ=255)
           cATATCTCCTGGTGATTAATGGCGGTGTGGACGTTGAGAATATCCTTACGAACTTTG‑ATCGCCATGATGACTTCTTCCAGGGAACAGTTTCCGGCACGCTCGCCGATCCCGTTCATT‑GCGCCTTCCACCTGGCGTGCACCGGCATGTACCGCCGCCAGTGAGTTGCCGACCGCCAGGCCTAAATCGTCGTGGGTATGTACGGAAATAATGGCTTTGTCGATGTTAGGTACGCGTTCATACAGGCCGCtgat                                                                                                                                                                                                                                       <  1:484385/251‑1 (MQ=255)
            atatCTCCTGGTGATTAATGGCGGTGTGGACGTTGAGAATATCCTTACGAACTTTG‑ATCGCCATGATGACTTCTTCCAGGGAACAGTTTCCGGCACGCTCGCCGATCCCGTTCATT‑GCGCCTTCCACCTGGCGTGCACCGGCATGTACCGCCGCCAGTGAGTTGCCGACCGCCAGGCCTAAATCGTCGTGGGTATGTACGGAAATAATGGCTTTGTCGATGTTAGGTACGCGTTCATACAGGCCGCtgatg                                                                                                                                                                                                                                      <  1:72922/251‑1 (MQ=255)
            atatCTCCTGGTGATTAATGGCGGTGTGGAAGTTGAGAATATCCTTACGGACTTTG‑ATCGCCATGATGACTTCTTCCAGGGAACAGTTTCCGGCACGCTCGCCGATCCCGTTCATT‑GAGCCTTCCACCTGGGGTGCACCGGCATGTAACGCCGCCCGGGAGGTGCCGACCGCCCGGCCTAAATTGTCGTGGGGATGTAAGGGAATAATGGGTTTGTTCATGTTAGGTACGCGGTCATACAgggcggtgggg                                                                                                                                                                                                                                      >  1:647212/1‑242 (MQ=255)
               tctcCTGGTGATTAATGGCGGTGTGGACGTTGAGAAAATCCTTACGCACTTTG‑ATCGCCATGATGACTTTTTCCAGGGAACCGGTTCCGGCACGCTCGCCGATCCCGTTCATT‑GCGCCTTCCACCTGGCGTGCACCGGCATGTACCGCCGCCCGTGAGGTGGCGACCGCCAGGCCCAAATCGGCGGGGGTATGTACGGGAATAATGGCTTTTTCGGTGGTAGGTACGCGGTCATAAAGGCCGCtgatgttt                                                                                                                                                                                                                                   >  1:168196/1‑248 (MQ=255)
                          taaaTGGCGGGGTGGACGTTGAAAATTTCCTTCCGACCTTTG‑ATCCCCATGATGTCTCCTTCCAGGGAACAGTTTCCGGCACGCGCGCCGATCCCGTTCATT‑GCGCCTTCCACCTGGCGTGCACCGGCATGTACCGCCGCCAGTGAGTTTCCGACCGCCAGGCCCAAATCGTCGTGGGTATGTACGGAGATAATGGCTTTGTCGATGTTAGGCACGCGTTCATACAGGCCGCTGATGATTCCGGCGAACt                                                                                                                                                                                                                         <  2:506947/248‑1 (MQ=255)
                                gtgttgtgtaagttaagaATATCCATACGAA‑TGTGTATCGCCATGATGACTTCTTCGAGGTAACTGTTTCCGGCACGCCCGCCGATCCCGTTCACT‑GCGCCTTCCACCTGGCGTGCACCGTCATGTACCGCCGCCAGTGAGTTGCCGACCGCCAGGCCTAAATCGTCGTGGGTATGTACGGAAATAATGGCTTTGTCGATGTTAGGTACGCGTTCATACAGGCCGCTGATGATTCCGGCGAACTCAAACgg                                                                                                                                                                                                                  <  2:164090/236‑1 (MQ=255)
                                  ggTGTGGACGTTGAGAATATCCTTACGAACTTTG‑ATCGCCATGATGACTTCTTCCAGGGAACAGTTTCCGGCACGCTCGCCGATCCCGTTCATT‑GCGCCTTCCACCTGGCGTGCACCGGCATGTACCGCCGCCAATGAGTTGCCGACCGCCAGGCCTAAATCGTCGTGGGTATGTACGGAAATAATGGCTTTGTCGATGTTAGGTACGCGTTCATACAGGCCGCTGATGATTCCGGCGAACTCAAACGGc                                                                                                                                                                                                                 >  1:680482/1‑250 (MQ=255)
                                        gACGTTGAGAATATCCTTACGAACTTTG‑ACCGCCATGATGACTCCTTCCAGGGAACAGTTTCCGGCACGCCCGCCGATCCCGTTCTTT‑GCGCCTTCCACCTGGCGTGCACCGGCATGTACCGCCGCCAGTGAGTTGCCGACCGCCAGGCCTAAATCGTCGTGGGTATGTACGGAAATAATGGCTTTGTCGATGTTAGGTACGCGTTCATACAGGCCGCTGATGATTCCGGCGAACTCAAACGGCATGgtgt                                                                                                                                                                                                          <  1:572654/251‑1 (MQ=255)
                                          cGTTGAGAATATCCTTACGAACTTTG‑ATCGCCATGACGACTTCTTCCAGGGAACAGTTTCCGGCACGCTCGCCGATCCCGTTCATT‑GCGCCTTCCACCTGGCGTGCACCGGCATGTACCGCCGCCAGTGATTTGCCGACCGCCAGGCCTAAATCGTCGTGGGTATGTACGGAAATAATGGCTTTGTCGATGTGAGGTACGCGTTCATACAGGCCGCTGATGATTCCGGCGAACTCAAACGGCATGGTGTAg                                                                                                                                                                                                        <  1:283031/251‑1 (MQ=255)
                                                gaATATCCTTACGAACTTTG‑ATCGCCATGATGACTTCTTCCAGGGAACAGTTTCCGGCACGCTCGCCGATCCCGTTCATT‑GCGCCTTCCACCTGGCGTGCACCGGCATGTACCGCCGCCAGTGAGTTGCCGACCGCCAGGCCTAAATCGTCGTGGGTATGTACGGAAATAATGGCTTTGTCGATGTTAGGTACGCGTTCATACAGGCCGCTGATGATTCCGGCGAACTCAAACGGCATGGTGTAGCCCACg                                                                                                                                                                                                  <  1:147707/251‑1 (MQ=255)
                                                    atCCTTACGAACTTTG‑ATCGCCATGATGACTTCTTCCAGGGACCAGTTTCCGGCACACTCGCCGATCCCGTTCATT‑GCGCCTTCCACCTGGCGTGCCCCGGCATGTACCGCCGCCAGTGAGTTGCCGACCGCCAGGCCTAAATCGTCGTGGGTATGTACGGAAATAATGGCTTTGTCGGTGTTAGGTACGCGTTCATACAGGCCGCTGATGATTCCGGCGAAATCAAACGGCATGGTGTAGCCCACggagt                                                                                                                                                                                              >  2:424831/1‑248 (MQ=255)
                                                      ccTTACGAACTTTG‑ATCGCCATGATGACTTCTTCCAGGGAACAGTTTCCGGCACGCTCGCCGATCCCGTTCATT‑GCGCCTTCCACCTGGCGTGCACCGGCATGTACCGCCGCCAGTGAGTTGCCGACCGCCAGGCCTAAATCGTCGTGGGTATGTACGGAAATAATGGCTTTGTCGATGTTAGGTACGCGTTCATACAGGCCGCTGATGATTCCGGCGAACTCAAACGGCATGGTGTAGCCCACGG‑TGTcc                                                                                                                                                                                           <  1:233896/251‑1 (MQ=255)
                                                                tttG‑ATCGCCATGATGACTTCTTCCAGGGAACAGTTTCCGGCACGCTCGCCGATCCCGTTCATT‑GCGCCTTCCACCTGGCGTGCACCGGCATGTACCGCCGCCAGTGAGTTGCCGACCGCCAGGCCTAAATCGTCGTGGGTCTGTACGGAAATAATGGCTTTGTCGATGTTAGGTACGCGGTCATACAGGCCGCTGATGATTCCGGCGAACTCAAACGGCATCGTGTAGCCCACGC‑GGTCCGGCATGTTg                                                                                                                                                                                  >  2:327453/1‑250 (MQ=255)
                                                                     aTCGCCATGATGACTTCTTCCAGGGAACAGTTTCCGGCACGCTCGCCGATCCCGTTCATT‑GCGCCTTCCACCTGGCGTGCACCGGCATGTACCGCCGCCAGTGAGTTGCCGACCGCCAGGCCTAAATCGTCGTGGGTATGTACGGAAATAATGGCTTTGTCGATGTTAGGTACGCGTTCATACAGGCCGCTGATGATTCCGGCGAACTCAAACGGCATGGTGTAGCCCACGG‑TGTCCGGAATGTTGAtggg                                                                                                                                                                             >  1:662928/1‑250 (MQ=255)
                                                                         ccATGATGACTTCTTCCAGGGAACAGTTTCCGGCACGCTCGCCGATCCCGTTCATT‑GCGCCTTCCACCTGGCGTGCACCGGCATGTACCGCCGCCAGTGAGTTGCCGACCGCCAGGCCTAAATCGTCGTGGGTATGTACGGAAATAATGGCTTTGTCGATGTTAGGTACGCGTTCATACAGGCCGCTGATGATTCCGGCGAACTCAAACGGCATGGTGTAGCCCACGG‑TGTCCGGAATGTTGAtggtggtg                                                                                                                                                                         >  1:41399/1‑251 (MQ=255)
                                                                          cATGATGACTTCTTCCAGGGAACAGTTTCCGGCACGCTCGCCGATCCCGTTCATT‑GCGACTTCCACCTGGCGTGCACCGGCATGTACCGCCGCCAGTGAGTTGCCGACCGCCAGGCCTAAATCGTCGTGGGTATGTACGGAAATAATGGCTTTGTCGATGTTAGGTACGCGTTCATACAGGCCGCTGATGATTCCGGCGAACTCAAACGGCATGGTGTAGCCCACGG‑TGTCCGGAATGTTGAtggtggtgg                                                                                                                                                                        >  1:298919/1‑251 (MQ=255)
                                                                                   ttcttcCAGGGAACAGTTTCCGGCACGCTCGCCGATCCCGTTCATT‑GCGCCTTCCACCTGGCGTGCACCGGCATGTACCGCCGCCAGTGAGTTGCCGACCGCCAGGCCTAAATCGTCGTGGGTATGTACGGAAATAATGGCTTTGTCGATGTTAGGTACGCGTTCATACAGGCCGCTGATGATTCCGGCGAACTCAAACGGCATGGTGTAGCCCACGG‑TGTCCGGAATGTTGATGGTGGTGGCACCGGCAt                                                                                                                                                               >  1:37263/1‑251 (MQ=255)
                                                                                     cttcCAGGGAACAGTTTCCGGCACGCTCGCCGATCCCGTTCATT‑GCGCCTTCCACCTGGCGTGCACCGGCATGTACCGCCGCCCGTGAGTTGCCGACCGCCAGGCCTAAATCGTCGTGGGTATGTACGGAAATAATGGCTTTGTCGATGTTAGGTACGCGTTCATACAGGCCGCTGATGATTACGGCGAACTCAAACGGCATGGTGTAGCCCACGG‑TGTCCGGAATGTTGATGGTGGTGGCACCGGCAtt                                                                                                                                                              >  1:441902/1‑250 (MQ=255)
                                                                                            ggAACAGTTTCCGGCACGCTCGCCGATCCCGTTCATT‑GCGCCTTCCACCTGGCGTGCACCGGCATGTACCGCCGCCAGTGAGTTGCCGACCGCCAGGCCTACATCGTCGTGGGTATGTACGGAAATAATGGCTTTGTCGCTGTTAGGTACGCGTTCATACAGGCCGCTGATGATTCCGGCGAACTCAAACGGCATTGTGTAGCCCACGG‑TGTCCGGAATGTTGCTGGTGGTGGCACCGGCATTAAccgccg                                                                                                                                                      >  2:349086/1‑251 (MQ=255)
                                                                                             gAACAGTTTCCGGCACGCTCGCCGATCCCGTTCATT‑GCGCCTTCCACCTGGCGTGCACCGGCATGTACCGCCGCCCGTGAGTTGCCAACCGCCAGGCCTAAATCGTCGTGGGTATGTACGGAACTAATAGCTTTGTCGCTGTTAGGTACGCGTTCATACAGGCAGCTGATGATTCCGTCGAACTCAAACGGCATGGTGTAGCCCACGG‑TGTCCGGAATGTTGTTGGTGGTGGCACCGGCATTAATcgccgc                                                                                                                                                     >  1:420152/1‑251 (MQ=255)
                                                                                                 ggttcccggccggccctccGATCCCGT‑CATTTGCGCCTTCCCCGGGGCGGGCCCCGGGATGTACCGCCGCCAGTGAGTTTCCGACCGCCATGCCTAAATCGTCGTGTGTATGTACGGAAATAATGGCTTTGTCGATGTTAGGTACGCGTTCATACAGGCCGCTGATGATTCCGGCGAACTCAAACGGCATGGTGTAGCCCACGG‑TGTCCGGAATGTTGATGGTGGTGGCACCGGCATTAATCGCCGCTTCg                                                                                                                                                 <  2:669302/234‑1 (MQ=255)
                                                                                                    ttCCGGCACGCTCGCCGATCCCGTTCATT‑GCGCCTTCCACCTGGCGTGCACCGGCATGTACCGCCGCCAGTGAGTTGCCGACCGCCAGGCCTAAATCGTCGTGGGCATGTACGGAAATAATGGCTTTGTCGATGTTAGGGACGCGTTCATACAGGCCGCTGCTGATTCCGGCGAACTCAAACGGCATGGTGTAGCCCACGG‑TGTCCGGAATGTTGATGGTGGTGGCCCCGGCATTAATCGCCGCTTCGAcc                                                                                                                                              >  2:572364/1‑251 (MQ=255)
                                                                                                        ggCACGCTCGCCGACCCCGTTCATT‑GCGCCTTCCACCTGGCGTGCACCGGCATGTACCGCCGCCAGTGAGTTGCCGACCGCCAGGCCTAAATCGTCGTGGGTATGTACGGAAATAATGGCTTTGTCGATGTTAGGTACGCGTTCATACAGGCCGCTGATGATTCCGGCGAACTCAAACGGCATGGTGTAGCCCACGG‑TGTCCGGAATGTTGATGGTGGTGGCACCGGCATTAATCGCCGCTTCGACCACTc                                                                                                                                          <  1:675963/251‑1 (MQ=255)
                                                                                                          cACGCTCGCCGATCCCGTTCATT‑GCGCCTTCCACCTGGCGTGCACCGGCCTGTACCGCCGCCAGTGAGTTGCCGACCGCCAGGCCTAAATCGTCGTAGGTATGTACGGAAATAATGGCTTTGTCGATGTTAGGTACGCGTTCATACAGGCCGCTGATGATTCCGGCGAACTCAAACGGAATGGTGTAGCCCACGG‑TGTCCGGCATGTTGATGGCGGTGGCACCGGCATTAATCGACGCTTCGGCCACTcgc                                                                                                                                        >  2:604572/1‑251 (MQ=255)
                                                                                                             gCTCGCCGATCCCGTTCATT‑GCGCCTTCCACCTGGCGTGCACCGACATGTACCGCCGCCAGTGAGTTGCCGGCCGCCAGGACTAAATCGTCGTGGGTATGGACGGACATAATGGCCTTGTCGATGTTAGGTACGCGTTCCTACAGGCCACTGATGATTCCGGCGATCACAAACGGCTTGGAGTAGCCCACGGCGGTCCGGG‑TGGTGATGGTGGTGGCGCCGGCATTaaccaccgactcgaacaatcacccc                                                                                                                                     >  2:688845/1‑228 (MQ=255)
                                                                                                                  ccGATCCCGTTCATT‑GCGCCTTCCACCTGGCGTGCACCGGCATGTACCGCCGCCAGTGAGTTGCCGACCGCCAGGCCTAAATCGTCGTGGGTATGTACGGAAATAATGGCTTTGTCGATGTTAGGTACGCGTTCATACAGGCCGCTGATGATTCCGGCGAACTCAAACGGCATGGTGTAGCCCACGG‑TGTCCGGAATGTTGATGGTGGTGGCACCGGCATTAATCGCCGCTTCGACCACTCGCGCCAGATc                                                                                                                                <  1:173255/251‑1 (MQ=255)
                                                                                                                  ccGATCCCGTTCATT‑GCGCCTCCCTCCTGGCGTGCACCGGCATGTACCGCGGCCAGTGAGTTGCCGACCGCCAGGCCTAAATCGTCGTGGGTATGTACGGAAATAATGGCTTTGTCGATGTTAGGTACGCGTTCATACAGGCCGCTGATGATTCCGGCGAACGCAAACGGCATGGTGTAGCCCACGG‑TGCCCGTAATGTTGATGGTGGTGGCACCGGCATTAATCGCCGCTTCGACCACTCGCGCCAGATc                                                                                                                                <  1:710583/251‑1 (MQ=255)
                                                                                                                   cGATCCCGTTCATT‑GCGCCTTCCACCTGGCGTGCACCGGCATGTACCGCCGCCAGTGAGTTGCCGACCGCCAGGCCTAAATCGTCGTGGGTATGTACGGAAATAATGGCTTTGTCGATGTTAGGTACGCGTTCATACAGGCCGCTGATGATTCCGGCGAACTCAAACGGCATGGTGTAGCCCCCGG‑TGTCCGGAATGTTGATGGTGGTGGCACCGGCATTCATCGCCGCTTCGACCACTCGCGCCAGGTCg                                                                                                                               >  2:21121/1‑251 (MQ=255)
                                                                                                                    cgttcccttgcgtggcgcCTTCCCCCTGGCGTGCCCCGGCATGTACCGCCGCCAGTGAGTTGCCGACCGCCAGGCCTAAATCGTCGTGGGTATGTACGGAAATAATGGCTTTGTCGATGTTAGGTACGCGTTCATACAGGCCGCTGATGATTCCGGCGAACTCAAACGGCATGGTGTAGCCCACGG‑TGTCCGGAATTTTGATGGTGGTGGCACCGGCATTAATCGCCGCTTCGACCACTCGCGCCAGATCg                                                                                                                               <  1:651896/237‑1 (MQ=255)
                                                                                                                          gTTCATT‑GCGCCTTCCACCTGGCGTGCACCGGCATGTACCGCCGCCAGTGAGTTGCCGACCGCCAGGCCTAAATCGTCGTGGGTATGTACGGAAATAATGGCTTTGTCGATGTTAGGTACGCGTTCATACAGGCCGCTGATGATTCCGGCGAACTCAAACGGCATGGTGTAGCCCACGG‑TGTCCGGAATGTTGATGGTGGTGGCACCGGCATTAATCGCCGCTTCGACCACTCGCGCCAGATCGGCAATgg                                                                                                                        >  1:39286/1‑251 (MQ=255)
                                                                                                                          gTTCATT‑GCGCCTTCCACCTGGCGTGCACCGGCATGTACCGCCGCCAGTGAGTTGCCGACCGCCAGGCCTAAATCGTCGTGGGTATGTACGGAAATAATGGCTTTGTCGATGTTAGGTACGCGTTCATACAGGCCGCTGATGATTCCGGCGAACTCAAACGGCATGGTGTAGCCCACGG‑TGTCCGGAATGTTGATGGTGGTGGCACCGGCATTAATCGCCGCTTCGACCACTCGCGCCAGATCGGCAATgg                                                                                                                        >  1:239218/1‑251 (MQ=255)
                                                                                                                             cATT‑GCGCCTTCCACCGGGCGTGCACCGGCATGTACCGCCGCCAGTGAGTTGCCGACCGCCAGGCCTAAATCGTCGTGGGTATGTACGGAAATAATGGCTTTGTCGATGTTAGGTACGCGTTCATACAGGCCGCTGATGATTCCGGCGAACTCAAACGGCATGGTGTAGCCCACGG‑TGTCCGGAATGTTGATGGTGGTGGCACCGGCATTAATCGCCGCTTCGACCACTCGCGCCAGATCGGCAATGGgt                                                                                                                      <  1:278022/250‑1 (MQ=255)
                                                                                                                                                 cGTGCACCGGCATGTACCGCCGCCAGTGAGTTGCCGACCGCCAGGCCTAAATCGTCGTGGGTATGTACGGAAATAATGGCTTTGTCGATGTTAGGTACGCGTTCATACAGGCCGCTGATGATTCCGGCGAACTCAAACGGCATGGTGTAGCCCACCG‑TGTCCGGTATGTTGATGGTGGTGGCACCGGCATTAATCGCCGCTTCGACCACTCGCGCCAGGTCGGCAATGGGTGTACGCCCGGCTACTTCGcc                                                                                                  >  2:174577/1‑250 (MQ=255)
                                                                                                                                                  gTGCACCGGCATGTACCGCCGCCAGTGAGTTGCCGACCGCCAGGCCTAAATCGTCGTGGGTATGTACGGAAATAATGGCTTTGTCGATGTTAGGTACGCGTTCATACAGGCCGCTGATGATTCCGGCGAACTCAAACGGCATGGTGTAGCCCACGG‑TGTCCGGAATGTTGATGGTGGTGGCACCGGCATTAATCGCCGCTTCGACCACTCGCGCCa                                                                                                                                    <  1:328279/216‑1 (MQ=255)
                                                                                                                                                  gTGCACCGGCATGTACCGCCGCCAGTGAGTTGCCGACCGCCAGGCCTAAATCGTCGTGGGTATGTACGGAAATAATGGCTTTGTCGATGTTAGGTACGCGTTCATACAGGCCGCTGATGATTCCGGCGAACTCAAACGGCATGGGGTAGCCCACGG‑TGTCCGGAATGTTGATGGTGGTGGCACCGGCATTAATCGCCGCTTCGACCAACCGCGCCa                                                                                                                                    >  2:328279/1‑216 (MQ=255)
                                                                                                                                                     cacccgcaggcaCCGCCGCCAGTGAGTTGCCGACCGCCAGGCCTAAATCGTCGTTGGTATGTACGGAAAGAATTGTTTTGTCGATGTTAGGTACGCGTTCATACAGGCCGCTGATGATTCCGGCGAACTCAAACGGCATGGTGTAGCCCACGG‑TGTCCGGAATGTTGATGGTGGTGGCACCGGCATTAATCGCCGCTTCGATCACTCGCGCCAGATCGGCAATGGGTGTACGCCCGGCATCTTCGCAAGaa                                                                                              <  2:359843/240‑1 (MQ=255)
                                                                                                                                                        ggccttgtcccgccgacgtggAGTTGCCGACCGCCAGGCCTAAATCGTCGTGGGTATGTACGGAAATAATGGCTTTGTCGATGTTAGGTACTCGTTCATACAGGCCGCTGATGATTCCGGCGAACGCAAACGGCATGGGGTAGCCCACGG‑TGTCCGGAATGGTGATGGTGGTGGCACCGGCATTAATCGCCGCTACGACCACACGCTACAGAGCGGCAATGGGCGTACGCCCGGCATCTTTGCAAGAAAAt                                                                                           <  2:73885/231‑1 (MQ=255)
                                                                                                                                                         ggcctGTACCGCCGCCAGTGAGTTGCCGACCGCCAGGCCTAAATCGTCGTGGGTATGTACGGAAATAATGGCTTTGTCGATGTTAGGAACGCGTTCATACAGGCCGCTGATGATTCCGGCGAACTCAAACGGCATGGTGTAGCCCACGG‑TGTCCGGAATGTTGATGGTGGTGGCACCGGCATTAATCGCCGCTTCGACCACTCGCGCCAGATCGGCAATGGGTGTACGCCCGGCATCTTCGCACGAAAAtt                                                                                          >  1:283726/5‑251 (MQ=255)
                                                                                                                                                           caggTACCGCCGCCAGTGAGGTGCCGACCGCCAGGCCTAAATCGTCGTGGGTTAGTACGGAAATAATGGCTTTGTCGATGGTAGGTACGCGTTCATACAGGCCGCTGATGATTCCGGAGAACACAAACGGCATGGTGTAGCCCACGG‑TGGCAGGTATGTTGATGGTGGTGGCACCGGCATTAATCGCCGCGTCGACAACTCGCGCCAGATCGGAACTTGGTGAAAGCCCAGCATATTCGCAAGACAAtacc                                                                                        >  2:259906/4‑248 (MQ=255)
                                                                                                                                                                  cgccgctattgagtggccGACCGTCAGGGCTAAATCTTCGGGGGTATGTACGGAAATAATGTCTTTGTCGATGTTAGTTTCGCGTTCATACAGGCCGCTGATGATTCCGGCGAACTCAAACGGCATGGTGGAGCCCACGG‑TGTCCGGAATGTTGATGGTGGTGGCACCGGCATTAATCGCCGCTTCGACCACTCGCGCCAGATCGGCAATGGGTGTACGCCCGTCATCTTCGCAAGAAAATTCAAcatcat                                                                                 <  2:664650/235‑1 (MQ=255)
                                                                                                                                                                              gttgccGACCGCCAGGCCTAAATCGTCGTGGGTATGTACGGAAATAATGGCTTTGTCGATGTTAGGTACGCGTTCATACAGGCCGCTGATGATTCCGGCGAACTCAAACGGCATGGTGTAGCCCACGG‑TGTCCGGAATGTTGATGGTGGTGGCACCGGCATTAATCGCCGCTTCGACCACTCGCGCCAGATCGGCAATGGGTGTACGCCCGGCATCTTCGCAAGAAAATTCAACATCATCGGTGTAATTAc                                                                     <  1:683944/247‑1 (MQ=255)
                                                                                                                                                                              gttgccGACCGCCAGGCCTAAATCGTCGTGGGTATGTACGGAAATAATGGCTTTGTCGATGTTAGGTACGCGTTCATACAGGCCGCTGATGATTCCGGCGAACTCAAACGGCATGGTGTAGCCCACGG‑TGTCCGGAATGTTGATGGTGGTGGCACCGGCATTAATCGCCGCTTCGACCACTCGCGCCAGATCGGCAATGGGTGTACGCCCGGCATCTTCGCAAGAAAATTCAACATCATCGGTGTAATTAc                                                                     <  1:21121/247‑1 (MQ=255)
                                                                                                                                                                                 tCCGACCGCCAGGCCCAAATCGTCGTGGGTATGTACGGAGATAATGGCTTTGTCGATGTTAGGAACGCGTTCATACAGGCCGCTGATGATTCCGGCGAACTCAAACGGCATGGCGAAGCCCACGG‑TGTCCGGAATGTTGATGGTGGTGGCACCGGCATTAATCGCCGCTTCGACCACTCGCGCCAGATCGGCACTGGGTGTACGCCCGGCATCTTCGCTAGAAAATTCacaataatcggtcaaattaaggg                                                                  >  2:361362/1‑229 (MQ=255)
                                                                                                                                                                                          cAGGCCTAAATCGTCGTGGGTATGTACGGAAATAATGGCTTTGTCGATGTTAGGTACGCGTTCATACAGGCCGCTGATGATTCCGGCGAACTCAAACGGCATGGTGTAGCCCACGG‑TGTCCGGAATGTTGATGGTGGTGGCACCGGCATTAATCGCCGCTTCGACCACTCGCGCCAGATCGGCAATGGGTGTACGCCCGGCATCTTCGCAAGAAAATTCAACATCATCGGTGTAATTACGGGCGCGTTTc                                                          >  1:96507/1‑250 (MQ=255)
                                                                                                                                                                                          cAGGCCTAAATCGTCGTGGGTATGTACGGAAATAATGGCTTTGTCGATGTTAGGTACGCGTTCATACAGGCCGCTGATGATTCCGGCGAACACAAACGGCATGGTGTAGCCCACGG‑TGTCCGGCATGTTGATGGTGGTGGCACCGGCATTAATCGCCGCTTCGACCACTCGCGCCAGATCGGCAATGGGTGTACGCCCGGCATCTTCGCAAGAAAATTCAACATCATCGGTGTAATTACGGGCGCGTTTc                                                          >  2:274896/1‑250 (MQ=255)
                                                                                                                                                                                           aGGCCTAAATCGTCGTGGGTATGTACGGAAATAATGGCTTTGTCGATGTTAGGTACGCGTTCATACAGGCCGCTGATGATTCCGGCGAACTCAAACGGCATGGTGGAGCCCACGG‑TGTCCGGAATGTTGATGGTGGTGGCACCGGCATTAATCGCCGCTTCGACCACTCGCGCCAGATCGGCAATGGGTGTACGCCCGGCATCTTCGCAAGAAAATTCAACATCATCGGTGTAATTACGGGCGCGTTTCAc                                                        <  1:42867/251‑1 (MQ=255)
                                                                                                                                                                                                            ggTATGTACGGAAATAATGGCTTTGTCGATGTTAGGTACGCGTTCATACAGGCCGCTGATGATTCCGGCGAACTCAAACGGCATGGTGTAGCCCACGG‑TGTCCGGAATGTTGATGGTGGTGGCACCGGCATTAATCGCCGCTTCGACCACTCGCGCCAGATCGGCAATGGGTGTACGCCCGGCCTCTTCGCAAGAAAATTCAACATCATCGGTGTAATTACGGGCGCGTTTCACCATATAGATAG‑CGCGtt                                      >  2:358840/1‑251 (MQ=255)
                                                                                                                                                                                                            ggTATGTACGGAAATAATGGCTTTGTCGATGTTAGGTACGCGTTCATACAGGCCGCTGATGATTCCGGCGAACTCAAACGGCATGGTGTAGCCCACGG‑TGTCCGGAATGGTGATGGTGGTGGCACCGGCATTAATCGCCGCTTCGACCACTCGCGCCAGATCGGCAATGGGTGTACGCCCGGCATCTTCGCAAGAAAATTCAACATCATCGGTGTAATTACGGGCGCGTTTCACCATATAGATAG‑CGCGtt                                      <  1:603805/251‑1 (MQ=255)
                                                                                                                                                                                                             gTATGTACGGAAATAATGGCTTTGTCGATGTTAGGTACGCGTTCATACAGGCCGCTGATGATTCCGGCGAACTCAAACGGCATGGTGTAGCCCACGG‑TGTCCGGAATGTTGATGGTGGTGGCACCGGCATTaa                                                                                                                                                            <  2:236824/133‑1 (MQ=255)
                                                                                                                                                                                                             gTATGTACGGAAATAATGGCTTTGTCGATGTTAGGTACGCGTTCATACAGGCCGCTGATGATTCCGGCGAACTCAAACGGCATGGTGTAGCCCACGG‑TGTCCGGAATGTTGATGGTGGTGGCACCGGCATTaa                                                                                                                                                            >  1:236824/1‑133 (MQ=255)
                                                                                                                                                                                                               aTGTACGGAAATAATGGCTTTGTCGATGTTAGGTACGCGTTCATACAGGCCGCTGATGATTCCGGCGAACTCAAACGGCATGGTGTAGCCCACGG‑TGTCCGGAATGTTGATGGTGGTGGCACCGGAATTAATCGCCGCTTCGACCACTCGCGCCAGATCGGCAATGGGTGTACGCCCGGCATCTTCGCAAGAAAATTCAACATCATCGGTGTAATTACGGGCGCGTTTCACCATATATACAG‑CGCGTTCg                                    >  1:702109/1‑250 (MQ=255)
                                                                                                                                                                                                                tGTACGGAAATAATGGCTTTGTCGATGTTAGGTACGCGTTCATACAGGCCGCTGATGATTCCGGCGAACTCAAACGGCATGGTGTAGCCCACGG‑TGTCCGGAATGTTGATGGTGGTGGCACCGGCATTAATCGCCGCTTCGACCACTCGCGCCAGATCGGCAATGGGTGTACGCCCGGCATCTTCGCAAGAAAATTCACCATCATCGGTGTAATTACGGGCGCGTTTCACCATATAGATAG‑CGCGTTCGAt                                  >  2:81305/1‑251 (MQ=255)
                                                                                                                                                                                                                   aCGGAAATAATGGCTTTGTCGATGTTAGGTACGCGTTCATACAGGCCGCTGATGATTCCGGCGAACTCAAACGGCATGGTGTAGCCCACGG‑TGTCCGGAATGTTGATGGTGGTGGCACCGGCATTAATCGCCGCTTCGACCACTCGCGCCAGATCGGCAATGGGTGTACGCCCGGCATCTTCGCAAGAAAATTCAACATCATCGGTGTAATTACGGGCGCGTTTCACCATATAGATAG‑CGCGTTCGATCAc                               >  1:289356/1‑251 (MQ=255)
                                                                                                                                                                                                                    cGGAAATAATGGCTTTGTCGATGGTAGGTACGCGTTCATACAGGCCGCTGATGATTCCGGCGAACTCAAACGGCATGGTGTAGCCCACGG‑TGTCCGGGATGTTGATGGTGGTGGCACCGGCATTAATCGCAGCTTCGACCACTCGCGCCAGATCGGCAATGGGTGTACGCCCGGCATCTTCGCCAGAAAATTCAACATCCTCGGTGTAATTACGGGCGCGTTTCACCATATAG‑TAGGCGCGTTctctcaat                              >  2:452708/1‑244 (MQ=255)
                                                                                                                                                                                                                     ggcaaTAATGGCTTTGTCGATGTTAGGTCCGAGTTCATACAGGCAGCTGATGATTCCGGCGAACTCAAACGGCATGGTGTAGCCCACGG‑AGTCCGGAATGTTGATGGTGGTGGCACCGGCATTAATCGCCGCTTCGACCACTCGCGCCAGATCGGCAATGGGTGTACGCCCGGCATCTTCGCAAGAAAATTCAACATCATCGGGGTACTTACGGGCGCGTTTCACAATATAGATCG‑CGGGTTCGATCActt                             >  2:689066/5‑249 (MQ=255)
                                                                                                                                                                                                                     ggaaaTAATGGCTTTGTCGATGTTAGGTACGCGTTCATACAGGCCGCTGATGATTCCGGCGAACTCAAACGGCATGGTGTAGCCCACGG‑TGTCCGGAATGTTGATGGTGGTGGCACCGGCATTAATCGCCGCTTCGACCACTCGCGCCAGATCGGCAATGGGTGTACGCCCGGCATCTTCGCAAGAAAATTCAACATCATCGGTGTAATTACGGGCGCGTTTCACCATATAGATAG‑CGCGTTCGATCaatt                             >  2:18316/5‑248 (MQ=255)
                                                                                                                                                                                                                      gaaaTAATGGCTTTGTCGATGTTAGGTACGCGTTCATACAGGCCGCTGATGATTCCGCCGAACTCAAACGGCATGGTGTAGCCCACGG‑TGTCCGGAATGTTGATGGTGGTGGCACCGGCATTAATCGCCGCTTCGACCACTCGCGCCAGATCGGCAATGGGTGTACGCCCGGCATCTTCGCAAGAAAATTCAACATCATCGGTGTAATTACGGGCGCGTTTCACCATATAGATAG‑CTCGTTCGCTCActtc                            >  1:408502/4‑248 (MQ=255)
                                                                                                                                                                                                                           aaTGGCTTTGTCGATGTTAGGTACGCGTTCATACAGGCCGCTGATGATTCCGGCGAACTCAAACGGCATGGTGTAGCCCACGG‑TGTCCGGAATGTTGATGGTGGTGGCACCGGCATTAATCGCCGCTTCGACCACTCGCGCCAGATCGGCAATGGGTGTACGCCCGGCATCTTCGCAAGAAAATTCAACATCATCGGTGTAATTACGGGCGCGTTTCACCATATAGATAG‑CGCGTTCGATCActtcatcc                        <  1:666727/250‑8 (MQ=255)
                                                                                                                                                                                                                               gCTTTGTCGATGTTAGGTCCGCGTTCATACAGGCCGCTGATGATTCCGGCGCACTCAAACGGCATGGTGTAGCCCACGG‑TGTCCGGAATGTTGATGGAGGTGGCACCGGCATTAATCGCCGCTTCGACCACTCGCGCCCGATCGGCAATGGGTGTACGCCCGGCATCTTCGCACGAAACTTCAACATCATCGGTGTCATTACGGGCGCGTTTCACCATCTAGATAG‑CGCGTTCGATCACTTCATCCAGCGt                   >  2:69574/1‑251 (MQ=255)
                                                                                                                                                                                                                                cTTTGTCGATGTTAGGTCCGCGTCCATACAGGCCGCTGATGATTCCGGCCAACTCAACCGGAATGCGGTAGCCCACGG‑CGTCCGGAATGTTGATGGTGGTGGCACCGGCATTAATCGCCGCATCGCCCACTCGCGCCAGATCGGCAATGGGTACACGCCCGGCTTCTTCGCCAGAAAATTCAACATCATCGGTGTAGTCACGGGCGAGTTTCACCATATCGCAAG‑CACGTTCGATCACGTACTCCAGCGtg                  >  2:601103/1‑251 (MQ=255)
                                                                                                                                                                                                                                cTTTGTCGATGTTAGGTACGCGTTCATACAGGCCGCTGATGATTCCGGCGAACTCAAACGGCATGGTGTAGCCCACGG‑TGTCCGGAATGTTGATGGTGGTGGCACCGGCATTAATCGCCGCTTCGACCACTCGCGCCAGATCGGCAATGGGTGTACGCCCGGCATCTTCGCAAGAAAATTCAACATCATCGGTGTAATTACGGGCGCGTTTCACCATATAGATAG‑CGCGTTCGATCAATTCATCCAGCGtg                  >  2:463404/1‑251 (MQ=255)
                                                                                                                                                                                                                                    gTCGATGTTAGGTACGCGTTCATACAGGCCGCTGATGATTCCGGCGAACTCAAACGGCATGGTGTAGCCCACGG‑TGTCCGGAATGTTGATGGTGGTGGCACCGGCATTAATCGCCGCTTCGACCACTCGCGCCAGATCGGCAATGGGTGTACGCCCGGCATCTTCGCAAGAAAATTCAACATCATCGGTGTAATTACGGGCGCGTTTCACCATATAGATAG‑CGCGTTCGATCACTTCATCCAGCGtgctgc              <  1:297145/251‑1 (MQ=255)
                                                                                                                                                                                                                                                cACGCGTTCATACAGGCCGCTGATGATTCCGGCGAACTCAAACGGCATGGTGTAGCCCACGG‑TGTCCGGAATGTTGATGGTGGTGGCACCGGCATTAATCGCCGCTTCGACCACTCGCGCCAGATCGGCAATGGGTGTACGCCCGGCATCTTCGCAAGAAAATTCAACATCATCGGTGTCATTACGGGCGCGTTTCACCATATAGATAG‑CGCGTTCGATCACCTCGTCCAGCGTGCTGCGCAGCTtggtgg  >  2:125975/1‑251 (MQ=255)
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                            
GGCTGGTGCGCCATATCTCCTGGTGATTAATGGCGGTGTGGACGTTGAGAATATCCTTACGAACTTTG‑ATCGCCATGATGACTTCTTCCAGGGAACAGTTTCCGGCACGCTCGCCGATCCCGTTCATT‑GCGCCTTCCACCTGGCGTGCACCGGCATGTACCGCCGCCAGTGAGTTTCCGACCGCCAGGCCCAAATCGTCGTGGGTATGTACGGAGATAATGGCTTTGTCGATGTTAGGCACGCGTTCATACAGGCCGCTGATGATTCCGGCGAACTCAAACGGCATGGTGTAGCCCACGG‑TGTCCGGAATGTTGATGGTGGTGGCACCGGCATTAATCGCCGCTTCGACCACTCGCGCCAGATCGGCAATGGGTGTACGCCCGGCATCTTCGCAAGAAAATTCAACATCATCGGTGTAATTACGGGCGCGTTTCACCATATAGATAG‑CGCGTTCGATCACCTCGTCCAGCGTGCTGCGCAGCTTGGTGG  >  NC_000913/82719‑83207

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