Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,688,337 A→G 5.1% L14P (CTC→CCC)  fumA ← fumarase A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,688,3370AG 5.1% 278.4 / 3.4 84L14P (CTC→CCC) fumAfumarase A
Reads supporting (aligned to +/- strand):  ref base A (42/37);  new base G (2/2);  total (45/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.01e-01

TTGCAGCGCCACATATTTATCATTTTCGCTGGCCTCCGGGTCACGCAGAATGTCGGCCACCTGTTGTTGGTGCGCCGGACGCAGCATGAACGACGCATCATGAAACGCCTGGCGCGCCAACAGAGTTAACGCTTCGGGTGCGACTTTCAAAATCTCCTGCCCTTCAAATTCAGATACGCTAACGTGTTCGCTGGTTAGCAGGTAATACTCAGTATCATCTTTTTTGAGTGGAAAAGGAGCCTGATAATGAAAGGGTTTGTTTGACATTGTTCTCTCACTTACTGCCTGGTTTGGTTATGCTCTGGGCGGGTGTTCCGTTGCCCTGTTAAAAGCGAGTAACAATATCCTACACACTTTTTTAACAAAAACTGAGACTAGTACGAC‑TTTTTGCGGCTCCAGGTTACTTCCCGTAGGATTCTTGCTTTAATAGTGGGATT‑AATTTCCACATTAAAACAGGGATTGATCATGCAAAAACTCATTAACTC  >  NC_000913/1688111‑1688595
                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                                    
tggcacggcaacaTTTTTTTCATTTTGGCTGGCCTCCGGGTCACGCAGAATGTCGGCCACCTGTTGTTGGTGCGCCGGATGCAGCATGAACGACGCATCATGAAACGCCTGGCGCGCCAACAGAGTTAACGCTTCGGGTGCGACTTTCAAAATCTCCTGCCCTTCAAATTCAGATACGCTAACGTGTTCGCTGGTTAGCAGGTAATACTCAGTATCATCTTTTTTGAGTGGAAAAGGAGCCTGATAATGAAAGGGTTTGTTTGACATTGTtc                                                                                                                                                                                                                         <  2:227796/262‑1 (MQ=255)
      cgcCACATATTTATCATTTTCGCTGGCCTCCGGGTCACGCAGAATGTCGGCCACCTGTTGTTGGTGCGCCGGACGCAGCATGAACGACGCATCATGAAACGCCTGGCGCGCCAACAGAGTTAACGCTTCGGGTGCGACTTTCAAAATCTCCTGCCCTTCAAATTCAGATACGCTAACGTGTTCGCTGGTTAGCAGGTAATACTCAGTATCATCTTTTTTGAGTGGAAAAGGAGCCTGATAATGAAAGGGTTTGTTTGACATTGTTTTTTCAc                                                                                                                                                                                                                   >  1:1109041/1‑272 (MQ=255)
          acaTATTTATCATTTTCGCTGGCCTCCGGGTCACGCAGAATGTCGGCCACCTGTTGTTGGTGCGCCGGACGCAGCATGAACGACGCATCATGAAACGCCTGGCGCGCCAACAGAGTTAACGCTTCGGGGGCGACTTTCAAAATCTCCTGCCCTTCAAATTCAGATACGCTAACGTGTTCGCTGGTTAGCAGGTAATACTCAGTATCATCTTTTTTGAGTGGAAAAGGAGCCTGATAATGAAAGGGTTTGTTTGAAATTTTTCTCCCACTTAc                                                                                                                                                                                                               >  1:1114360/1‑272 (MQ=255)
                     aTTTTCGCTGGCCTCCGGGTCACGCAGAATGTCGGCCACCTGTTGTTGGTGCGCCGGACGCAGCATGAACGACGCATCATGAAACGCCTGGCGCGCCAACAGAGTTAACGCTTCGGGTGCGACTTTCAAAATCTCCTGCCCTTCAAATTCAGATACGCTAACGTGTTCGCTGGTTAGCAGGTAATACTCAGTATCATCTTTTTTTAGTGGAAAAGGAGCCTGATAATGAAAGGGTTTGTTTGACATTGTTCTCTCACTTACTGccctgtct                                                                                                                                                                                                     >  1:250676/1‑265 (MQ=255)
                           gCTGGCCTCCGGGTCACGCAGAATGTCGGCCACCTGTTGTTGGTGCGCCGGACGCAGCATGAACGACGCATCATGAAACGCCTGGCGCGCCAACAGAGTTAACGCTTCGGGTGCGACTTTCAAAATCTCCTGCCCTTCAAATTCAGATACGCTAACGTGTTCGCTGGTTAGCAGGTAATACTCAGTATCATCTTTTTTGAGTGGAAAAGGAGCCTGATAATGAAAGGGTTTGTTTGACATTGTTCTCTCACTTACTGCCTGGTTTGGTTa                                                                                                                                                                                                >  2:189605/1‑270 (MQ=255)
                            cTGGCCTCCGGGTCACGCAGAATGTCGGCCACCTGTTGTTGGTGCGCCGGACGCAGCATGAACGACGCATCATGAAACGCCTGGCGCGCCAACAGAGTTAACGCTTCGGGTGCGACTTTCAAAATCTCCTGCCCTTCAAATTCAGATACGCTAACGTGTTCGCTGGTTAGCAGGTAATACTCAGTATCATCTTTTTTGAGTGGAAAAGGAGCCTGATAATGAAAGGGTTTGTTTGACATTGTTCTCTCACTTACTGCCTGGTTTGGTTAt                                                                                                                                                                                               <  1:189605/270‑1 (MQ=255)
                                             cAGAATGTCGGCCACCTGTTGTTGGTGCGCCGGACGCAGCATGAACGACGCATCATGAAACGCCTGGCGCGCCAACAGAGTTAACGCTTCGGGTGCGACTTTCAAAATCTCCTGCCCTTCAAATTCAGATACGCTAACGTGTTCGCTGGTTAGCAGGTAATACTCAGTATCATCTTTTTTGa                                                                                                                                                                                                                                                                      <  1:128460/182‑1 (MQ=255)
                                                             tgttgtTGGTGCGCCGGACGCAGCATGAACGACGCATCATGAAACGCCTGGCGCGCCAACAGAGTTAACGCTTCGGGTGCGACTTTCAAAATCTCCTGCCCTTCAAATTCAGATACGCTAACGTGTTCGCTGGTTAGCAGGTAATACTCAGTATCATCTTTTTTGAGTGGAAAAGGAGCCTGATAATGAAAGGGTTTGTTTGACATTGTTCTCTCACTTACTGCCt                                                                                                                                                                                                          >  2:248485/1‑226 (MQ=255)
                                                              gttgttGGTGCGCCGGACGCAGCATGAACGACGCATCATGAAACGCCTGGCGCGCCAACAGAGTTAACGCTTCGGGTGCGACTTTCAAAATCTCCTGCCCTTCAAATTCAGATACGCTAACGTGTTCGCTGGTTAGCAGGTAATACTCAGTATCATCTTTTTTGAGTGGAAAAGGAGCCTGATAATGAAAGGGTTTGTTTGACATTGTTCTCTCACTTACTGCCtg                                                                                                                                                                                                         <  1:248485/226‑1 (MQ=255)
                                                                     gTGCGCCGGACGCAGCATGAACGACGCATCATGAAACGCCTGGCGCGCCAACAGAGTTAACGCTTCGGGTGCGACTTTCAAAATCTCCTGCCCTTCAAATTCAGATACGCTAACGTGTTCGCTGGTTAGCAGGTAATACTCAGTATCATCTTTTTTGAGTg                                                                                                                                                                                                                                                                   >  2:651354/1‑161 (MQ=255)
                                                                     gTGCGCCGGACGCAGCATGAACGACGCATCATGAAACGCCTGGCGCGCCAACAGAGTTAACGCTTCGGGTGCGACTTTCAAAATCTCCTGCCCTTCAAATTCAGATACGCTAACGTGTTCGCTGGTTAGCAGGTAATACTCAGTATCATCTTTTTTGAGTGGAAAAGGAGCCTGATAATGAAAGGGTTTGTTTGACATTGTTCTCTCACTTACTGCCTGGTTTGGTTATGctct                                                                                                                                                                                          >  1:834179/1‑234 (MQ=255)
                                                                      tGCGCCGGACGCAGCATGAACGACGCATCATGAAACGCCTGGCGCGCCAACAGAGTTAACGCTTCGGGTGCGACTTTCAAAATCTCCTGCCCTTCAAATTCAGATACGCTAACGTGTTCGCTGGTTAGCAGGTAATACTCAGTATCATCTTTTTTGAGTgg                                                                                                                                                                                                                                                                  <  1:651355/161‑1 (MQ=255)
                                                                        cgcCGGACGCAGCATGAACGACGCATCATGAAACGCCTGGCGCGCCAACAGAGTTAACGCTTCGGGTGCGACTTTCAAAATCTCCTGCCCTTCAAATTCAGATACGCTAACGTGTTCGCTGGTTAGCAGGTAATACTCAGTATCATCTTTTTTGAGTGGAAAAGGAGCCTGATAATGAAAGGGtttgtttg                                                                                                                                                                                                                                  >  2:1092336/1‑191 (MQ=255)
                                                                         gcCGGACGCAGCATGAACGACGCATCATGAAACGCCTGGCGCGCCAACAGAGTTAACGCTTCGGGTGCGACTTTCAAAATCTCCTGCCCTTCAAATTCAGATACGCTAACGTGTTCGCTGGTTAGCAGGTAATACTCAGTATCATCTTTTTTGAGTGGAAAAGGAGCCTGATAATGAAAGGGTTTGTTTGACATTGTTCTCTCACTTACTGCCtggtttg                                                                                                                                                                                                    >  2:753309/1‑220 (MQ=255)
                                                                          ccGGACGCAGCATGAACGACGCATCATGAAACGCCTGGCGCGCCAACAGAGTTAACGCTTCGGGTGCGACTTTCAAAATCTCCTGCCCTTCAAATTCAGATACGCTAACGTGTTCGCTGGTTAGCAGGTAATACTCAGTATCATCTTTTTTGAGTGGAAAAGGAGCCTGATAATGAAAGGGTTTGTTTGACATTGTTCTCTCACTTACTGCCtggtttgg                                                                                                                                                                                                   <  1:753310/220‑1 (MQ=255)
                                                                            ggACGCAGCATGAACGACGCATCATGAAACGCCTGGCGCGCCAACAGAGTTAACGCTTCGGGTGCGACTTTCAAAATCTCCTGCCCTTCAAATTCAGATACGCTAACGTGTTCGCTGGTTAGCAGGTAATACTCAGTATCATCTTTTTTGAGTGGAAAAGGAGCCTGATAATGAAAGGgtt                                                                                                                                                                                                                                        >  1:449477/1‑181 (MQ=255)
                                                                             gACGCAGCATGAACGACGCATCATGAAACGCCTGGCGCGCCAACAGAGTTAACGCTTCGGGTGCGACTTTCAAAATCTCCTGCCCTTCAAATTCAGATACGCTAACGTGTTCGCTGGTTAGCAGGTAATACTCAGTATCATCTTTTTTGAGTGGAAAAGGAGCCTGATAATGAAAGGgttt                                                                                                                                                                                                                                       <  2:449476/181‑1 (MQ=255)
                                                                             gACGCAGCATGAACGACGCATCATGAAACGCCTGGCGCGCCAACAGAGTTAACGCTTCGGGTGCGACTTTCAAAATCTCCTGCCCTTCAAATTCAGATACGCTAACGTGTTCGCTGGTTAGCAGGTAATACTCAGTATCATC‑TTTTTGAGTGGAAAAGGAGCCTGATAATGAAAGGGTTTGTTTGACATTGTTCTCTCACTTACTGCCtg                                                                                                                                                                                                         <  1:594881/210‑1 (MQ=255)
                                                                                               caccgataaaCGCCGGCCGCCCCAACAGATTTAAGCCTTCGGGGGCGACTTTCAAAATCTCCTGCCCTTCAAATTCAGATACGCTAACGTGTTCGCTGGTTAGCAGGTAATACTCAGTATCATCTTTTTTGAGTGGAAAAGGAGCCTGATAATGAAAGGGTTTGTTTGACATTGTTCTCTCACTTACTGCCTGGTTTGGTTATGCTCTGGGCGGGTGTTCCGTTGCCCTGTTAAAAGCGAGTAACAATATCCTACACACTTTTTTAACaaaa                                                                                                                          <  2:258521/265‑1 (MQ=255)
                                                                                                 tcatGAAACGCCTGGCGCGCCAACAGAGTTAACGCTTCGGGTGCGACTTTCAAAATCTCCTGCCCTTCAAATTCAGATACGCTAACGTGTTCGCTGGTTAGCAGGTAATACTCAGTATCATCTTTTTTGAGTGGAAAAGGAGCCTGATAATGAAAGGGTTTGTTTGACATTg                                                                                                                                                                                                                            >  2:581792/1‑172 (MQ=255)
                                                                                                  catGAAACGCCTGGCGCGCCAACAGAGTTAACGCTTCGGGTGCGACTTTCAAAATCTCCTGCCCTTCAAATTCAGATACGCTAACGTGTTCGCTGGTTAGCAGGTAATACTCAGTATCATCTTTTTTGAGTGGAAAAGGAGCCTGATAATGAAAGGGTTTGTTTGACATTGt                                                                                                                                                                                                                           <  1:581793/172‑1 (MQ=255)
                                                                                                     gAAACGCCTGGCGCGCCAACAGAGTTAACGCTTCGGGTGCGACTTTCAAAATCTCCTGCCCTTCAAATTCAGATACGCTAACGTGTTCGCTGGTTAGCAGGTAATACTCAGTATCATC‑TTTTTGAGTGGAAAAGGAGCCTGATAATGAAAGGGTTTGTTTGACATTGTTCTCTCACTTACTGCCtggtttg                                                                                                                                                                                                    >  2:640371/1‑191 (MQ=255)
                                                                                                      aaaCGCCTGGCGCGCCAACATAGTTAACGCTTCGGGTGCGACTTTCAAAATCTCCTGCCCTTCAAATTCAGATACGCTAACGTGTTCGCTGGTTAGCAGGTAATACTCAGTATCATC‑TTTTTGAGTGGAAAAGGAGCCTGATAATGAAAGGGTTTGTTTGACATTGTTCTCTCACTTACTGCCtggtttgg                                                                                                                                                                                                   <  1:640372/191‑1 (MQ=255)
                                                                                                          gCCTGGCGCGCCAACAGAGTTAACGCTTCGGGTGCGACTTTCAAAATCTCCTGCCCTTCAAATTCAGATACGCTAACGTGTTCGCTGGTTAGCAGGTAATACTCAGTATCATCTTTTTTGAGTGGAAAAGGAGCCTGATAATGAAAGGGTTTGTTTGACATTGTTCTCTCACTTACTGCCTGGTTTGGTTATGCTCTGGGCGGGTGTTCCGTTGCCCTGTTAAAAGCGAGTAACAATATCCTACACACTTTTTTAACaaa                                                                                                                           >  1:258521/1‑260 (MQ=255)
                                                                                                                 gcgcCAACAGAGTTAACGCTTCGGGTGCGACTTTCAAAATCTCCTGCCCTTCAAATTCAGATACGCTAACGTGTTCGCTGGTTAGCAGGTAATACTCAGTATCATCTTTTTTGAGTGGAAAAGGAGCCTGATAATGAAAGGGTTTGTTTGACATTGTTCTCTCACTTACTGCCTGGTTTGGTTa                                                                                                                                                                                                >  1:524679/1‑184 (MQ=255)
                                                                                                                 gcgcCAACAGAGTTAACGCTTCGGGTGCGACTTTCAAAATCTCCTGCCCTTCAAATTCAGATACGCTAACGTGTTCGCTGGTTAGCAGGTAATACTCAGTATCATCTTTTTTGAGTGGAAAAGGAGCCTGATAATGAAAGGGTTTGTTTGACATTGTTCTCTCACTTACTGCCTGGTTTGGTTATGCTCTGGGCGGGTGTTCCGTTGCCCTGTTAAAAGCGAGTAACAATATCCTAACCACTTTTTTAAAAAAAACTGAGACTAGTAcaaac                                                                                                        >  2:950468/1‑268 (MQ=255)
                                                                                                                  cgcCAACAGAGTTAACGCTTCGGGTGCGACTTTCAAAATCTCCTGCCCTTCAAATTCAGATACGCTAACGTGTTCGCTGGTTAGCAGGTAATACTCAGTATCATCTTTTTTGAGTGGAAAAGGAGCCTGATAATGAAAGGGTTTGTTTGACATTGTTCTCTCACTTACTGCCTGGTTTGGTTAt                                                                                                                                                                                               <  2:524678/184‑1 (MQ=255)
                                                                                                                   aatggacaggTTAACGCTTCGGGTGCGACTTTCAAAATCTCCTGCCCTTCAAATTCAGATACGCTAACGTGTTCGCTGGTTAGCAGGTAATACTCAGTATCATCTTTTTTGAGTGGAAAAGGAGCCTGATAATGAAAGGGTTTGTTTGACATTGTTCTCTCACTTACTGCCTGGTTTGGTCATGCTCTGTGCGTGTGTTCCGTTGCCCTGTTAAAAGCGAGTAACAATATCCTACACACTTTTTTAACAAAAACTGAGACTAGTACGAC‑ttt                                                                                                     <  1:1051106/263‑1 (MQ=255)
                                                                                                                                aCGCTTCGGGTGCGACTTTCAAAATCTCCTGCCCTTCAAATTCAGATACGCTAACGTGTTCGCTGGTTAGCAGGTAATACTCAGTATCATCTTTTTTGAGTGGAAAAGGAGCCTGATAATGAAAGGGTTTGTTTGACATTGTTCTCTCACTTACTGCCt                                                                                                                                                                                                          >  2:1024172/1‑159 (MQ=255)
                                                                                                                                aCGCTTCGGGTGCGACTTTCAAAATCTCCTGCCCTTCAAATTCAGATACGCTAACGTGTTCGCTGGTTAGCAGGTAATACTCAGTATCATCTTTTTTGAGTGGAAAAGGAGCCTGATAATGAAAGGGTTTGTTTGACATTGTTCTCTCACTTACTGCCTGGTTTGGTTATGCTCTGGGCg                                                                                                                                                                                     >  1:57221/1‑180 (MQ=255)
                                                                                                                                 cGCTTCGGGTGCGACTTTCAAAATCTCCTGCCCTTCAAATTCAGATACGCTAACGTGTTCGCTGGTTAGCAGGTAATACTCAGTATCATCTTTTTTGAGTGGAAAAGGAGCCTGATAATGAAAGGGTTTGTTTGACATTGTTCTCTCACTTACTGCCtg                                                                                                                                                                                                         <  1:1024173/159‑1 (MQ=255)
                                                                                                                                 cGCTTCGGGTGCGACTTTCAAAATCTCCTGCCCTTCAAATTCAGATACGCTAACGTGTTCGCTGGTTAGCAGGTAATACTCAGTATCATCTTTTTTGAGTGGAAAAGGAGCCTGATAATGAAAGGGTTTGTTTGACATTGTTCTCTCACTTACTGCCTGGTTTGGTTATGCTCTGGGCgg                                                                                                                                                                                    <  2:57221/180‑1 (MQ=255)
                                                                                                                                        ggTGCGACTTTCAAAATCTCCTGCCCTTCAAATTCAGATACGCTAACGTGTTCGCTGGTTAGCAGGTAATACTCAGTATCATCTTTTTTGAGTGGAAAAGGAGCCTGATAATGAAAGGGTTTGTTTGACATTGTTCTCTCACTTACTGCCTGGTTTGGTTATGc                                                                                                                                                                                             >  2:780137/1‑164 (MQ=255)
                                                                                                                                         gTGCGACTTTCAAAATCTCCTGCCCTTCAAATTCAGATACGCTAACGTGTTCGCTGGTTAGCAGGTAATACTCAGTATCATCTTTTTTGAGTGGAAAAGGAGCCTGATAATGAAAGGGTTTGTTTGACATTGTTCTCTCACTTACTGCCTGGTTTGGTTATGct                                                                                                                                                                                            <  1:780138/164‑1 (MQ=255)
                                                                                                                                          tGCGACTTTCAAAATCTCCTGCCCTTCAAATTCAGATACGCTAACGTGTTCGCTGGTTAGCAGGTAATACTCAGTATCATCTTTTTTGAGTGGAAAAGGAGCCTGATAATGAAAGGGTTTGTTTGACATTGTTCTCTCACTTACTGCCt                                                                                                                                                                                                          >  2:1118764/1‑149 (MQ=255)
                                                                                                                                          tGCGACTTTCAAAATCTCCTGCCCTTCAAATTCAGATACGCTAACGTGTTCGCTGGTTAGCAGGTAATACTCAGTATCATCTTTTTTGAGTGGAAAAGGAGCCTGATAATGAAAGGGTTTGTTTGACATTGTTCTCTCACTTACTGCCTGGTTTGGTTATGCTCTGGGCGGGTGTTCCGTTGCCCTGTTAAAAGCGAGTAACAATATCCTACACACTTTTTTAACAAAAACTGAGACTAATACGCCATTTTTTGGGCCCCAgggtggttttc                                                                               >  2:1051105/1‑263 (MQ=255)
                                                                                                                                           gCGACTTTCAAAATCTCCTGCCCTTCAAATTCAGATACGCTAACGTGTTCGCTGGTTAGCAGGTAATACTCAGTATCATCTTTTTTGAGTGGAAAAGGAGCCTGATAATGAAAGGGTTTGTTTGACATTGTTCTCTCACTTACTGCCtg                                                                                                                                                                                                         <  1:1118765/149‑1 (MQ=255)
                                                                                                                                                tttCAAAATCTCCTGCCCTTCAAATTCAGATACGCTAACGTGTTCGCTGGGTAGCAGGTAATACTCAGTATCATCTTTTTTGAGTGGAAAAGGAGCCTGATAATGAAAGGGTTTGTTTGACATTGTTCTCTCACTTACTGCCTGGTTTGGTTATGCTCTGGGCGGGTGTTCCGTTGCCCTGTTAAAAGCGAGTAACAATATCCTACACACTTTTTTAACAAAAACTGAGACTAGTACGAC‑TTTTTGCGGCTCCAGGTTACTTCCCGTAGGAt                                                                        <  1:950469/272‑1 (MQ=255)
                                                                                                                                                            cTGCCCTTCAAATTCAGATACGCTAACGTGTTCGCTGGTTAGCAGGTAAAACTCAGTCTCATCTTTTTTGAGTCGAAAAGTAGCCTGATAATTAAAGGGTTTGTTTGACATTGTTCTCTCACTTACTTCCTGGTTTTTTTATTCTCTGTGCGGGTTGTCCGTTTCCCTGTTaac                                                                                                                                                               >  1:974144/1‑173 (MQ=255)
                                                                                                                                                             ggCCCTTCAAATGCAGATACGCTAACGTGTGCGCTGGTTAGCAGGTAATACTCAGTATCATCTTTTGTGAGTGGAAAAGGAGCCTGATAATGAAAGGGTTTGTTTGACATTGTTCTCTCACTTACTGCCTGGTTTGGTTATGCTCTGGGCGGGTGTTCCGTTGCCCTGTTaaaa                                                                                                                                                              <  2:974143/173‑1 (MQ=255)
                                                                                                                                                               cccGTCAAATTCAGATACGCTAACGTGTTCGCTGGTTAGCAGGTAATACTCAGTATCATCTTTTTTGAGTGGAAAAGGAGCCTGATAATGAAAGGGTTTGTTTGACATTGTTCTCTCACTTACTGCCTGGTTTGGTTATGCTCTGGGCg                                                                                                                                                                                     >  2:426144/1‑149 (MQ=255)
                                                                                                                                                                ccgtCAAATTCAGATACGCTAACGTGTTCGCTGGTTAGCAGGTAATACTCAGTATCATCTTTTTTGAGTGGAAAAGGAGCCTGATAATGAAAGGGTTTGTTTGACATTGTTCTCTCACTTACTGCCTGGTTTGGTTATGCTCTGGGCgg                                                                                                                                                                                    <  1:426145/146‑1 (MQ=255)
                                                                                                                                                                       aTTCAGATACGCTAACGTGTTCGCTGGTTAGCAGGTAATACTCAGTATCATCTTTTTTGTGTGGAAAAGGAGCCTGATAATGaaa                                                                                                                                                                                                                                             >  1:376672/1‑85 (MQ=255)
                                                                                                                                                                       aTTCAGATACGCTAACGTGTTCGCTGGTTAGCAGGTAATACTCAGTATCATCTTTTTTGAGTGGAAAAGGAGCCTGATAATGAAAGGGTTTGTTTGACATTGTTCTCTCACTTACTGCCTGGTTTGGTTATGCTCTGGGCGGGTGTTCCGTTGCCCTGTTAAAAGCGAGTAACAATATCCTACACACTTTTTTAACAAAAACTGAGACTAGTACGAC‑TTTTTGCg                                                                                                <  2:347091/225‑1 (MQ=255)
                                                                                                                                                                       aTTCAGATACGCTAACGTGTTCGCTGGTTAGCAGGTAATACTCAGTATCATCTTTTTTGAGTGGAAAAGGAGCCTGATAATGAAAGGGTTTGTTTGACATTGTTCTCTCACTTACTGCCTGGTTTGGTTATGCTCTGGGCGGGTGTTCCGTTGCCCTGTTAAAAGCGAGTAACAATATCCTACACACTTTTTTAACAAAAACTGAGACTAGTACGAC‑TTTTTGCg                                                                                                >  2:915497/1‑225 (MQ=255)
                                                                                                                                                                       aTTCAGATACGCTAACGTGTTCGCTGGTTAGCAGGTAATACTCAGTATCATCTTTTTTGAGTGGAAAAGGAGCCTGATAATGAAAGGGTTTGTTTGACATTGTTCTCTCACTTACTGCCTGGTTTGGTTATGCTCTGGGCGGGTGTTCCGTTGCCCTGTTAAAAGCGAGTAACAATATCCTACACACTTTTTTAACAAAAACTGAGACTAGTACGAC‑TTTTTGCg                                                                                                >  1:347092/1‑225 (MQ=255)
                                                                                                                                                                        ttCAGATACGCTAACGTGTTCGCTGGTTAGCAGGTAATACTCAGTATCATCTTTTTTGAGTGGAAAAGGAGCCTGATAATGAAAg                                                                                                                                                                                                                                            <  2:376671/85‑1 (MQ=255)
                                                                                                                                                                        ttCAGATACGCTAACGTGTTCGCTGGTTAGCAGGTAATACTCAGTATCATCTTTTTTGAGTGGAAAAGGAGCCTGATAATGAAAGGGTTTGTTTGACATTGTTCTCTCACTTACTGCCTGGTTTGGTTATGCTCTGGGCGGGTGTTCCGTTGCCCTGTTAAAAGCGAGTAACAATATCCTACACACTTTTTTAACAAAAACTGAGACTAGTACGAC‑TTTTTGCgg                                                                                               <  1:915498/225‑1 (MQ=255)
                                                                                                                                                                                             gCTGGTTAGCAGGTAATACTCAGTATCATCTTTTTTGAGTGGAAAAGGAGCCTGATAATGAAAGGGTTTGTTTGACATTGTTCTCTCACTTACTGCCTGGTTTGGTTATGCTCTGGGCGGGTGTTCCGTTGCCCTGTTAAAAGCGAGTAACAATATCCTACACACTTTTTTAACAAAAACTGAGACTAGTACGAC‑TTTTTGCGGCTCCAGGTTACTTCCCGTAGGATTCTTGCTTTAATAGTGGGATT‑AATTTACACATTTAAAAAGGGAtt                          >  2:773104/1‑272 (MQ=255)
                                                                                                                                                                                              cTGGTTAGCAGGTAATACTCAGTATCATCTTTTTTGAGTGGAAAAGGAGCCTGATAATGAAAGGGTTTGTTTGACATTGTTCTCTCACTTACTGCCTGGTTTGGTTATGCTCTGGGCGGGTGTTCCGTTGCCCTGTTAAAAg                                                                                                                                                             >  2:200251/1‑142 (MQ=255)
                                                                                                                                                                                              cTGGTTAGCAGGTAATACTCAGTATCATCTTTTTTGAGTGGAAAAGGAGCCTGATAATGAAAGGGTTTGTTTGACATTGTTCTCTCACTTACTGCCTGGTTTGGTTATGCTCTGGGCGGGTGTTCCGTTGCCCTGTTAAAAGCGAGt                                                                                                                                                        >  2:761614/1‑147 (MQ=255)
                                                                                                                                                                                               tGGTTAGCAGGTAATACTCAGTATCATCTTTTTTGAGTGGAAAAGGAGCCTGATAATGAAAGGGTTTGTTTGACATTGTTCTCTCACTTACTGCCTGGTTTGGTTATGCTCTGGGCGGGTGTTCCGTTGCCCTGTTAAAAGc                                                                                                                                                            <  1:200251/142‑1 (MQ=255)
                                                                                                                                                                                               tGGTTAGCAGGTAATACTCAGTATCATCTTTTTTGAGTGGAAAAGGAGCCTGATAATGAAAGGGTTTGTTTGACATTGTTCTCTCACTTACTGCCTGGTTTGGTTATGCTCTGGGCGGGTGTTCCGTTGCCCTGTTAAAAGCGAGTa                                                                                                                                                       <  1:761615/147‑1 (MQ=255)
                                                                                                                                                                                               tGGTTAGCAGGTAATACTCAGTATCATCTTTTTTGAGTGGAAAAGGAGCCTGATAATGAAAGGGTTTGTTTGACATTGTTCTCTCACTTACTGCCTGGTTTGGTTATGCTCTGGGCGGGTGTTCCGTTGCCCTGTTAAAAGCGAGTAACAATATCCTACACACTTTTTTAACAAAAACTGAGACTAGTACGAC‑TTTTTGCGGCTCCAGGTTACTTCCCGTAGGATTCCTGCTTTAATTGTTGGATT‑AATTTCCAAATTaaaaacgggttata                        >  2:306072/1‑262 (MQ=255)
                                                                                                                                                                                                  ttAGCAGGTAATACTCAGTATCATCTTTTTTGAGTGGAAAAGGAGCCTGATAATGAAAGGgtttgttt                                                                                                                                                                                                                                   >  2:1098680/1‑68 (MQ=255)
                                                                                                                                                                                                   tAGCAGGTAATACTCAGTATCATCTTTTTTGAGTGGAAAAGGAGCCTGATAATGAAAGGGtttgtttg                                                                                                                                                                                                                                  <  1:1098681/68‑1 (MQ=255)
                                                                                                                                                                                                    aGCAGGTAATACTCAGTATCATCTTTTTTGAGTGGAAAAGGAGCCTGATAATGAAAGGGTTTGTTTGACATTGTTCTCTCACTTACTGCCTGGTTTGGTTa                                                                                                                                                                                                >  2:382116/1‑101 (MQ=255)
                                                                                                                                                                                                    aGCAGGTAATACTCAGTATCATCTTTTTTGAGTGGAAAAGGAGCCTGATAATGAAAGGGTTTGTTTGACATTGTTCTCTCACTTACTGCCTGGTTTGGTTATGCTCTGGGCGGGTGTTCCGTTGCCCTGTTAAAAGCGAGTAACAATATCCTACACACTTTTTTAACAAAAACTGAGCCTAGTACGAC‑TTTTTGCGGCCCCAGGTTACTTCCAGTAGGATACTTGATTTAATAGTGGGATTAAATATCCACATAAAAACaagattggatc                      >  2:1011924/1‑260 (MQ=255)
                                                                                                                                                                                                     gCAGGTAATACTCAGTATCATCTTTTTTGAGTGGAAAAGGAGCCTGATAATGAAAGGGTTTGTTTGACATTGTTCTCTCACTTACTGCCTGGTTTGGTTAt                                                                                                                                                                                               <  1:382117/101‑1 (MQ=255)
                                                                                                                                                                                                       aGGTATTATTTATTATTATCTTTTTTGAGTGGAAAAGGAGCCTGATAATGAAAGGTTTTGTTTGACTTTTTTCTTCCACTTACTGCCTGGTTTGGTTATGCTCTGGGCGGGTTTTCCGTTGCCCTGTTAAAAGCGAGTAAAAATATCCTACACACTTTTTTAACAAAAACTGAGACTAGTACGAC‑TTTTTGCGGCTCCAGGTTACTTCCCGTAGGATTCTTGCTTTAATAGTGGGATT‑AATTTCCACATTAAAACAGGGATTGATCATGCaa                <  1:773105/272‑1 (MQ=255)
                                                                                                                                                                                                              aCTCAGTATCATCTTTTTTGAGTGGAAAAGGAGCCTGATAATGAAAGGgtttgttt                                                                                                                                                                                                                                   >  1:904278/1‑56 (MQ=255)
                                                                                                                                                                                                              aCTCAGTATCATCTTTTTTGAGTGGAAAAGGAGCCTGATAATGAAAGGGTTTGTTTGACATTGTTCTCTCACTTACTGCCTGGTTTGGTTATGCTCTGGGCGGGTGTTCCGTTGCCCTGTTAAAAGCGAGTAACAATATCCTACACACTTTTTTAACAAAAACTGAGACTAGTACGAC‑TTTTTGCGGCTCCAGGTTACTTCCCGTAGGATTCTTGCTTTAATAGTGGGATT‑AATTTCCACAtt                                      >  2:281870/1‑243 (MQ=255)
                                                                                                                                                                                                               cTCAGTATCATCTTTTTTGAGTGGAAAAGGAGCCTGATAATGAAAGGGtttgtttg                                                                                                                                                                                                                                  <  2:904277/56‑1 (MQ=255)
                                                                                                                                                                                                               cTCAGTATCATCTTTTTTGAGTGGAAAAGGAGCCTGATAATGAAAGGGTTTGTTTGACATTGTTCTCTCACTTACTGCCTGGTTTGGTTATGCTCTGGGCGGGTGTTCCGTTGCCCTGTTAAAAGCGAGTAACAATATCCTACACACTTTTTTAACAAAAACTGAGACTAGTACGAC‑TTTTTGCGGCTCCAGGTTACTTCCCGTAGGATTCTTGCTTTAATAGTGGGATT‑AATTTCCACATTa                                     <  1:281871/243‑1 (MQ=255)
                                                                                                                                                                                                                     atcatcTTTTTTGAGTGGAAAAGGAGCCTGATAATGAAAGGGTTTGTTTGACATTg                                                                                                                                                                                                                            >  2:868243/1‑56 (MQ=255)
                                                                                                                                                                                                                     atcatcTTTTTTGAGTGGAAAAGGAGCCTGATAATGAAAGGGTTTGTTTGACATTGTTCTCTCACTTACTGCCTGGTTTGGTTATGCTCTGGGCGGGTGTTCCGTTGCCCTGTTAAAAGCGAGTAACAATATCCTACACACTTTTTTAACAAAAACTGAGACTAGTACGAC‑tt                                                                                                      >  2:37255/1‑173 (MQ=255)
                                                                                                                                                                                                                      tcatcTTTTTTGGGTGGAAAAGGAGCCTGATAATGAAAGGGTTTGTTTGACATTGTTCTCTCACTTACTGCCTGGTTTGGTTATGCTCTGGGCGGGTGTTCCGTTGCCCTGTTAAAAGCGAGTAACAATATCCTACACACTTTTTTAACAAAAACTGAGACTAGTACGAC‑tt                                                                                                      >  1:799008/1‑172 (MQ=255)
                                                                                                                                                                                                                      tcatcTTTTTTGAGTGGAAAAGGAGCCTGATAATGAAAGGGTTTGTTTTACATTGTTCTCTCACTTACTGCCTGGTTTGGTTATGCTCTGGGCGGGTGTTCCGTTGCCCTGTTAAAAGCGAGTAACAATATCCTACACACTTTTTTAACAAAAACTGAGACTAGTACGAC‑TTTTTGCGTCTCCAGGTTACTTCCCGTAGGATTCTTGCTTTAATAGTGTGATT‑AATTTCCACATTAAAACAGGGATTGATCATGCAAAAACTCATTAACTc  >  1:542827/1‑271 (MQ=255)
                                                                                                                                                                                                                      tcatcTTTTTTGAGTGGAAAAGGAGCCTGATAATGAAAGGGTTTGTTTGACATTGt                                                                                                                                                                                                                           <  1:868244/56‑1 (MQ=255)
                                                                                                                                                                                                                      tcatcTTTTTTGAGTGGAAAAGGAGCCTGATAATGAAAGGGTTTGTTTGACATTGTTCTCTCACTTACTGCCTGGTTTGGTTATGctct                                                                                                                                                                                          >  2:241319/1‑89 (MQ=255)
                                                                                                                                                                                                                      tcatcTTTTTTGAGTGGAAAAGGAGCCTGATAATGAAAGGGTTTGTTTGACATTGTTCTCTCACTTACTGCCTGGTTTGGTTATGCTCTGGGCGGGTGTTCCGTTGCCCTGTTAAAAGCGAGTAACAATATCCTACACACTTTTTTAACAAAAACTGAGACTAGTACGAC‑ttt                                                                                                     <  1:37255/173‑1 (MQ=255)
                                                                                                                                                                                                                       catcTTTTTTGGGTGGAAAAGGAGCCTGATAATGAAAGGGTTTGTTTGACATTGTTCTCTCACTTACTGCCTGGTTTGGTTATGCTCTGGGCGGGTGTTCCGTTGCCCTGTTAAAAGCGAGTAACAATATCCTACACACTTTTTTAACAAAAACTGAGACTAGTACGAC‑ttt                                                                                                     <  2:799007/172‑1 (MQ=255)
                                                                                                                                                                                                                       catcTTTTTTGGGTGGAAAAGGAGCCTGATAATGAAAGGGTTTGTTTGACATTGTTCTCTCACTTACTGCCTGGTTTGGTTATGCTCTGGGCGGGTGTTCCGTTGCCCTGTTAAAAGCGAGTAACAATATCCTACACACTTTTTTAACAAAAACTGAGACTAGTACGAC‑tt                                                                                                      >  1:798980/1‑171 (MQ=255)
                                                                                                                                                                                                                       catcTTTTTTGAGTGGCAAAGGAGCCTGCTAATGAAAGGGTTTGTTTGACATTGTTCTCTCACTTACTGCCTGGTTTGGTTATGctct                                                                                                                                                                                          >  1:771384/1‑88 (MQ=255)
                                                                                                                                                                                                                       catcTTTTTTGAGTGGAAAAGGAGCCTGATAATGAAAGGGTTTGTTTGACATTGTTCTCTCACTTACTGCCTGGTTTGGTTATGCTCTg                                                                                                                                                                                         <  1:241319/89‑1 (MQ=255)
                                                                                                                                                                                                                       catcTTTTTTGAGTGGAAAAGGAGCCTGATAATGAAAGGGTTTGTTTGACATTGTTCTCTCACTTACTGCCTGGTTTGGTTATGCTCTGGGCg                                                                                                                                                                                     >  2:1009590/1‑93 (MQ=255)
                                                                                                                                                                                                                       catcTTTTTTGAGTGGAAAAGGAGCCTGATAATGAAAGGGTTTGTTTGACATTGTTCTCTCACTTACTGCCTGGTTTGGTTATGCTCTGGGCGGGTGTTCCGTTGCCCTGTTAAAAg                                                                                                                                                             >  2:127332/1‑117 (MQ=255)
                                                                                                                                                                                                                        atcTTTTTTGGGTGGAAAAGGAGCCTGATAATGAAAGGGTTTGTTTGACATTGTTCTCTCACTTACTGCCTGGTTTGGTTATGCTCTGGGCGGGTGTTCCGTTGCCCTGTTAAAAGCGAGTAACAATATCCTACACACTTTTTTAACAAAAACTGAGACTAGTACGAC‑ttt                                                                                                     <  2:798979/171‑1 (MQ=255)
                                                                                                                                                                                                                        atcTTTTTTGAGTGGAAAAGGAGCCTGATAATGAAAGGGTTTGTTTGACATTGTTCTCTCACTTACTGCCTGGTTTGGTTATGCTCTg                                                                                                                                                                                         <  2:771383/88‑1 (MQ=255)
                                                                                                                                                                                                                        atcTTTTTTGAGTGGAAAAGGAGCCTGATAATGAAAGGGTTTGTTTGACATTGTTCTCTCACTTACTGCCTGGTTTGGTTATGCTCTGGGCgg                                                                                                                                                                                    <  1:1009591/93‑1 (MQ=255)
                                                                                                                                                                                                                        atcTTTTTTGAGTGGAAAAGGAGCCTGATAATGAAAGGGTTTGTTTGACATTGTTCTCTCACTTACTGCCTGGTTTGGTTATGCTCTGGGCGGGTGTTCCGTTGCCCTGTTAAAAGc                                                                                                                                                            <  1:127332/117‑1 (MQ=255)
                                                                                                                                                                                                                               ttGAGTGGAAAAGGAGCCTGATAATGAAAGGGTTTGTTTGACATTGTTCTCTCACTTACTGCCTGGTTTGGTTATGCTCTGGGCGGGTGTTCCGt                                                                                                                                                                           >  1:779726/1‑95 (MQ=255)
                                                                                                                                                                                                                                tGAGTGGAAAAGGAGCCTGATAATGAAAGGGTTTGTTTGACATTGTTCTCTCACTTACTGCCTGGTTTGGTTATGCTCTGGGCGGGTGTTCCGtt                                                                                                                                                                          <  2:779725/95‑1 (MQ=255)
                                                                                                                                                                                                                                 gAGTGGAAAAGGAGCCTGATAATGAAAGGGTTTGTTTGACATTGTTCTCTCACTTACTGCCTGGTTTGGTTATGCTCTGGGCGGGTGTTCCGTTGCCCTGTTAAAAGCGAGTAACAATATCCTACACACTTTTTTAACAAAAACTGAGACTAGTACGAC‑TTTTTGCGGCTCCAGGTTACTTCCCGTAgg                                                                          >  2:1058587/1‑189 (MQ=255)
                                                                                                                                                                                                                                 gAGTGGAAAAGGAGCCTGATAATGAAAGGGTTTGTTTGACATTGTTCTCTCACTTACTGCCTGGTTTGGTTATGCTCTGGGCGGGTGTTCCGTTGCCCTGTTAAAAGCGAGTAACAATATCCTACACACTTTTTTAACAAAAACTGAGACTAGTACGAC‑TTTTTGCGGCTCCAGGTTACTTCCCGTAgg                                                                          <  1:1058588/189‑1 (MQ=255)
                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                                    
TTGCAGCGCCACATATTTATCATTTTCGCTGGCCTCCGGGTCACGCAGAATGTCGGCCACCTGTTGTTGGTGCGCCGGACGCAGCATGAACGACGCATCATGAAACGCCTGGCGCGCCAACAGAGTTAACGCTTCGGGTGCGACTTTCAAAATCTCCTGCCCTTCAAATTCAGATACGCTAACGTGTTCGCTGGTTAGCAGGTAATACTCAGTATCATCTTTTTTGAGTGGAAAAGGAGCCTGATAATGAAAGGGTTTGTTTGACATTGTTCTCTCACTTACTGCCTGGTTTGGTTATGCTCTGGGCGGGTGTTCCGTTGCCCTGTTAAAAGCGAGTAACAATATCCTACACACTTTTTTAACAAAAACTGAGACTAGTACGAC‑TTTTTGCGGCTCCAGGTTACTTCCCGTAGGATTCTTGCTTTAATAGTGGGATT‑AATTTCCACATTAAAACAGGGATTGATCATGCAAAAACTCATTAACTC  >  NC_000913/1688111‑1688595

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