Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,093,736 T→G 21.0% V90G (GTG→GGG)  hisB → histidinol‑phosphatase and imidazoleglycerol‑phosphate dehydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,093,7360TG21.0% 134.3 / 4.6 75V90G (GTG→GGG) hisBhistidinol‑phosphatase and imidazoleglycerol‑phosphate dehydratase
Reads supporting (aligned to +/- strand):  ref base T (39/16);  new base G (10/7);  total (52/23)
Fisher's exact test for biased strand distribution p-value = 3.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCCTGATTAGCGAACCGCCGAGTGATTTTCAGGTGGACCGTTTTGATAAACTCGCCTTTGAACCGGGCGTGATCCCGGAACTGCTGAAGCTGCAAAAAGCGGGCTACAAGCTGGTGATGATCACTAATCAGGATGGTCTTGGAACACAAAGTTTCCCACAGGCGGATTTCGATGGCCCGCACAACCTGATGATGCAGATCTTCACCTCGCAAGGCGTACAGTTTGATGAAGTGCTGATTTGTCCGCACCTGCCCGCCGATGAGTGCGACTGCCGTAAGCCGAAAGTAAAACTGGTGGAACGTTATCTGGCTGAGCAAGCGATGGATCGCGCTAACAGTTATGTGATTGGCGATCGCGCGACCGACATTCAACTGGCGGAAAACATGGGCATTACTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCT  >  NC_000913/2093505‑2093974
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                              
cccTGATTAGCGACCGGGCGAGTGATGTGCAGGTGGACCGTTTTGATAAACGCGCCTTTGAACCGGGCGTGATCCCGGAACTGCTGAAGCTGCAAAAAGCGGGCTACAAGCTGGTGATGATCACTAATCAGGATGGTCTTGGAACACAAAGTTTCCCACAGGCGGATTTCGATGGCCCGCACAACCTGATGATGCAGATCTTCACCTCGCAAGGCGTACAGTTTGATGAAGTGCTGAt                                                                                                                                                                                                                                          <  1:386632/238‑1 (MQ=255)
     aTTAGCGAACCGCCGAGTGATTTTCAGGTGGACCGTTTTGATAAACTCGCCTTTGAACCGGGCGTGATCCCGGAACTGCTGAAGCTGCAAAAAGCGGGCTACAAGCTGGTGATGATCACTAATCAGGATGGTCTTGGAACACAAAGTTTCCCACAGGCGGATTTCGATGGCCCGCCCAACCTGATGATGCAGATCTTCACCTCGCAAGGCGTACCGTTTGATTAAGGGCTGATTTGTc                                                                                                                                                                                                                                     >  1:899924/1‑238 (MQ=255)
     aTTAGCGAACCGCCGAGTGATTTTCAGGTGGACCGTTTTGATAAACTCGCCTTTGAACCGGGCGTGATCCCGGAACTGCTGAAGCTGCAAAAAGCGGGCTACAAGCTGGTGATGATCACTAATCAGGATGGTCTTGGAACACAAAGTTTCCCACAGGCGGATTTCGATGGCCCGCCCAACCTGATGATGCAGATCTTCACCTCCCAAGGCGTACAGTTTGATTAAGGGCTGATTTTTc                                                                                                                                                                                                                                     >  1:890641/1‑238 (MQ=255)
     aTTAGCGAACCGCCGAGTGATTTTCAGGTGGACCGTTTTGATAAACTCGCCTTTGAACCGGGCGTGATCCCGGAACTGCTGAAGCTGCAAAAAGCGGGCTACAAGCTGGCGATGATCACTAATCAGGATGGTCTTGGAACACAAAGTTTCCCACAGGCGGATTTCGATGGCCCGCACAACCTGATGATGCAGACCTTCACCTCGCAAGGCGTACAGTTTGCAGAACTGCTGATTTGTcc                                                                                                                                                                                                                                    >  1:704564/1‑239 (MQ=255)
      ttgccgaccccgcgggtgtTTTTCAGGTGGACCGTTTGGATAAACTCGCCTTTGAACCGGGCGTGATCCCGGAACTGCTGAAGCTGCAAAAAGCGGGCTACAAGCTGGTGATGATCACTAATCAGGATGGTCTTGGAACACAAAGTTTCCCACAGGCGGATTTCGATGGCCCGCACAACCTGATGATGCAGATCTTCACCTCGCAAGGCGTACAGTTTGATGAAGTGCTGATTTGTCCg                                                                                                                                                                                                                                   <  2:385404/224‑1 (MQ=255)
          cGAACCGCCGAGTGATTTTCAGGTGGACCGTTTTGATAAACTCGCCTTTGAACCGGGCGTGATCCCGGAACTGCTGAAGCTGCAAAAAGCGGGCTACAAGCTGGTGATGATCACTAATCAGGATGGTCTTGGAACACAAAGTTTCCCACAGGCGGATTTCGATGGCCCGCCCAACCTGATGATGCAGATCTTCACCTCGCAAGGCGTACAGTTTGATGAAGGGCTGATTTTGCCGccc                                                                                                                                                                                                                                >  2:258114/1‑236 (MQ=255)
          cGAACCGCCGAGTGATTTTCAGGTGGACCGTTTTGATAAACTCGCCTTTGAACCGGGCGTGATCCCGGAACTGCTGAAGCTGCAAAAAGCGGGCTACAAGCTGGTGATGATCACTAATCAGGATGGTCTTGGAACACAAAGTTTCCCACAGGCGGATTTCGATGGCCCGCACAACCTGATGCTTCAGATCTTCACCTCGCAAGGCGTACAGTTTGATGAAGGGCTGATTTGTCCGCAc                                                                                                                                                                                                                                >  1:1038548/1‑238 (MQ=255)
            aaCCGCCGAGTGATTTTCAGGTGGACCGTTTTGATAAACTCGCCTTTGAACCGGGCGTGATCCCGGAACTGCTGAAGCTGCAAAAAGCGGGCTACAAGCTGGTGATGATCACTAATCAGGATGGACTTGGAACACAAAGTTTCCCACAGGCGGATTTCGATGGCCCGCACAACCTGATGATGCAGATCTTCACCTCGCAAGGCGTACAGTTTGATGAAGGGCGGATTTGTCCGCGccc                                                                                                                                                                                                                              >  1:126556/1‑237 (MQ=255)
               cgccgAGTGATTTTCAGGTGGACCGTTTTGATAAACTCGCCTTTGAACCGGGCGTGATCCCGGAACTGCTGAAGCTGCAAAAAGCGGGCTACAAGCTGGTGATGATCACTAATCAGGATGGTCTTGGAACACAAAGTTTCCCCCAGGCGGATTTCGATGGCCCGCCCACCCTGATGTTGCAGATCTTCACCTCGCAAGGCTTACAGTTTGATGAAGGGCTGGTTTGTCCGCACCTGcc                                                                                                                                                                                                                           >  1:151797/1‑238 (MQ=255)
               cgccgAGTGATTTTCAGGTGGACCGTTTTGATAAACTCGCCTTTGAACCGGGCGTGATCCCGGAACTGCTGAAGCTGCAAAAAGCGGGCTACAAGCTGGTGATGATCACTAATCAGGATGGTCTTGGAACACAAAGTTTCCCACAGGCGGATTTCGATGGCCCGCACAACCTGATGATGCAGATCTTCACCTCGCAAGGCGTACAGTTTGATGAAGTGATGATTTTTCCGCACCTGcc                                                                                                                                                                                                                           >  1:927401/1‑238 (MQ=255)
               cgccgAGTGATTTTCAGGTGGACCGTTTTGATAAACTCGCCTTTGAACCGGGCGTGATCCCGGAACTGCTGAAGCTGCAAAAAGCGGGCTACAAGCTGGTGATGATCACTAATCAGGATGGTCTTGGAACACAAAGTTTCCCACAGGCGGATTTCGATGGCCCCCACAACCTGATGATGCAGATCTTCACCTCGCAAGGCGTACAGTTTGATGAAGAGCTGATTTGTCCGCACCTGcc                                                                                                                                                                                                                           >  1:754464/1‑238 (MQ=255)
                         ttttCAGGTGGACCGTTTTGATAAACTCGCCTTTGAACCGGGCGTGATCCCGGAACTGCTGAAGCTGCAAAAAGCGGGCTACAAGCTGGTGATGATCACTAATCAGGATGGTCTTGGAACACAAAGTTTCCCACAGGCGGATTTCGATGGCCCGCACAACCTGATGATGCAGATCTTCACCTCGCAAGGCGTACAGTTTGATGAAGTGCTGATTTGTCCGCACCTGCCCGCCGATGAGt                                                                                                                                                                                                                >  1:796380/1‑239 (MQ=255)
                             ggggtggccttttttGAAAACCCCGCTTTGAACCCGGGGGTGACCCCGGAATGGCAGAAGCGACAAAAAGGGGGCTAAAAGCTGGGGATGATCACTAATGAGGTGGGTCTGGGAAAAAAAAGTTTCCCACAGGCGGATTTCGAGGGCCCGCACAACCTGATGATGAAGATCTCCACCCCGCAAGGCGTACAGTTTGATGAAGTGCTGTTTTGTCCGCACCGGCCCGCCGATGAGTGCGa                                                                                                                                                                                                            <  2:880110/228‑1 (MQ=255)
                                   gACCGTTTTGATAAACTCGCCTTTGAACCGGGCGTGATCCCGGAACTGCTGAAGCTGCAAAAAGCGGGCTACAAGCTGGTGATGATCACTAATCAGGATGGTCTTGGAACACAAAGTTTCCCACAGGCGGATTTCGATGGCCCGCACAACCTGATGATGCAGATCTTCACCTCGCAAGGCGTACAGTTTGATGAAGGGCTGATTTGTCCGCACCTGCCCCCCGATGAGGGCGACTGCCg                                                                                                                                                                                                      >  2:810774/1‑239 (MQ=255)
                                        ttttGATAAACTCGCCTTTGAACCGGGCGTGATCCCGGAACTGCTGAAGCTGCAAAAAGCGGGCTACAAGCTGGTGATGATCACTAATCAGGATGGTCTTGGAACACAAAGTTTCCCACAGGCGGATTTCGATGGCCCGCACAACCTGATGATGCAGATCTTCACCTCGCAAGGCGTACAGTTTGATGAAGTGCTGATTTGTCCGCACCTGCCCGCCGATGAGTGCGACTGCCGTAAg                                                                                                                                                                                                  >  1:328238/1‑238 (MQ=255)
                                                aaCTCGCCTTTGAACCGGGCGTGATCCCGGAACTGCTGAAGCTGCAAAAAGCGGGCTACAAGCTGGTGATGATCACTAATCAGGATGGTCTTGGAACACAAAGTTTCCCACAGGCGGATTTCGATGGCCCGCCCAACCTGATGATGCAGATCTTCACCTCGCAAGGCGTACCGTTTGCTGAAGGGCGGATTTGTCCCCCCCCGCCCCCCGAGTAGGGGGACTGCCGAAACCCGaaaaaa                                                                                                                                                                                         >  2:950137/1‑236 (MQ=255)
                                                  cTCGCCTTTGAACCGGGCGTGATCCCGGAACTGCTGAAGCTGCAAAAAGCGGGCTACAAGCTGGTGATGATCACTAATCAGGATGGTCTTGGAACACAAAGTTTCCCACAGGCGGATTTCGATGGCCCGCACAACCTGATGATGCAGATCTTCACCTCCCAAGGCGTAACGGTTGATGAAGTGCTGATTTGTCCGCCCCTGCCCGCCGCTGAGGGCGACTGCAGGAAACCGaaaagaaa                                                                                                                                                                                       >  2:602954/1‑234 (MQ=255)
                                                       cTTTGAACCGGGCGTGATCCCGGAACTGCTGAAGCTGCAAAAAGCGGGCTACAAGCTGGTGATGATCACTAATCAGGATGGTCTTGGAACACAAAGTTTCCCACAGGCGGATTTCGATGGCCCGCACAACCTGATGATGCAGATCTTCACCTCGCAAGGCGTACAGTTTGATGAAGTGCTGATTTGTCCGCACCTGCCCGCCGATGAGTGCGACTGCCGTAAGCCGAAAGTAAAACtgg                                                                                                                                                                                  >  1:603309/1‑239 (MQ=255)
                                                              ccccggggttgttccGGAAATGGTGAAGGCTGAAAAAGGCGGTCTACAAGCGGGGGATGACACCTAAGCAGTATGGTCTGGAAACACAAGTTTCCCCACGGGCGGATTTCGATGGCCCGCACACCCTGAGGAGGCAGCTCTTCCCCTCGAAAGGGGTACATTTTGATGAAGTGCTGTTTTCCCCGCCCCTGCCCGCCGATGAGTGCGACTGCCGTAAGCCGAAAGTAAAACTGGTGGaa                                                                                                                                                                             <  2:1038548/226‑1 (MQ=255)
                                                              cccGGGCGTAACCCCGGAACTGCGAAAGCTGAAAAAGGCGGGTAAAAAGCTGGTAATGATCATTAATCAGGATGGTCTTGGAACACAAAGTTTCCCACGGGCGGATTCCGATGGCCCGCACACCCTGATGATGCAGATCTTCACCTCGCAAGGCGTACAGTTTGATGAAGTGCTGATTTGTCCGCACCTGCCCGCCGATGAGTGCGACTGCCGTAAGCCGAAAGTAAAACTGGTGGaa                                                                                                                                                                             <  1:357124/237‑1 (MQ=255)
                                                              acgggggGTGACCCGGGACGTGCTGAAGCGGAAAAAGGGGGCAGAAGAGCGGGTGATGACCAATAAGCGGGATGGTCTTGAAACAAAAATTTCCCCAAAGGCGAATTCGGATGGCCGGCACACCGTGATGATGCAGATCTTCCCCGGCAAGGGGGTACAGTTTGATGAAGGGCTGTTTTGTCCGCCCCTGCCCGCCGATGAGTGCGACTGCCGTAAGCCGAAAGTAAAACTGGTGGaa                                                                                                                                                                             <  2:867741/235‑1 (MQ=255)
                                                                  ggCGTGATCCCGGAACTGCTGAAGCTGCAAAAAGCGGGCTACAAGCTGGTGATGATCACTAATCAGGATGGTCTTGGAACACAAAGTTTCCCACAGGCGGATTTCGATGGCCCGCACAACCTGATGATGCAGATCTTCACCTCGCAAGGCGTACAGTTTGATGAAGTGCTGCTTTGTCCGCACCTGCCCGCCGATGAGTGCGACTGCCGTAAGCCGAAAGTAAAACTGGTGGAACGTTa                                                                                                                                                                        >  1:377543/1‑239 (MQ=255)
                                                                  ggCGTGATCCCGGAACTGCTGAAGCTGCAAAAAGCGGGCTACAAGCTGGTGATGATCACTAATCAGGATGGTCTTGGAACACAAAGTTTCCCACAGGCGGATTTCGATGGCCCGCACAACCTGATGATGCAGATCTTCACCTCGCAAGGCGTACAGTTTGATGAAGTGCTGATTTGTCCGCACCTGCCCGCCGATGAGTGCGACTGCCGTAAGCCGAAAGTAAAACTGGTGGAACGTTa                                                                                                                                                                        >  1:927538/1‑239 (MQ=255)
                                                                  ggCGTGATCCCGGAACTGCTGAAGCTGCAAAAAGCGGGCTACAAGCTGGTGATGATCACTAATCAGGATGGTCTTGGAACACAAAGTTTCCCACAGGCGGATTTCGATGGCCCGCACAACCTGATGATGCAGATCTTCACCTCGCAAGGCGTACAGTTTGATGAAGTGCTGATTTGTCCGCACCTGCCCGCCCCTGAGTGCGACTGCCGTAAGCCCAAAGAaaaaacggggggaagtta                                                                                                                                                                        >  2:299228/1‑225 (MQ=255)
                                                                  ggCGTGATCCCGGAACTGCTGAAGCTGCAAAAAGCGGGCTACAAGCTGGTGATGATCACTAATCAGGATGGACTTGGAACACAAAGTTTCCCACAGGCGGATTTCGATGGCCCGCACAACCTGATGATGCAGATCTTCACCTCGCAAGGCGTACAGTTTGCTGAAGGGCTGATTTGTCCGCACCTGACCGCCGCTGAGTGCGACGGCCGTAAGCCGAAAGGAAAACTGGTGGAAcagta                                                                                                                                                                        >  1:519446/1‑235 (MQ=255)
                                                                       gATCCCGGAACTGCTGAAGCTGCAAAAAGCGGGCTACAAGCTGGTGATGATCACTAATCAGGATGGTCTTGGAACACAAAGTTTCCCACAGGCGGATTTCGATGGCCCGCACAACCTGATGATGCAGATCTTCACCTCGCAAGGCGTACAGTTTGATGAAGTGCTGATTTGTCCGCACCTGCCCGCCGATGAGTGCGACTGCCGTAAGCCGAAAGTAAAACTGGTGGAACGTTATCTg                                                                                                                                                                    >  1:292547/1‑238 (MQ=255)
                                                                              gaaccggctgacggcgaaaaaaGGGGGCTAAAAGTTGGGTATGCCCACTAAGAGGGAGGGCCTGGGAACAAAAAGTTCCCCACGGGCGGATTCGGATGCCCCGCACAACTGGAGGATGCAGATCTTCCCCCCGAAGGGCGTACATTTGAAGGAAGTGTGGTTTTGCCCCCCCCGCCCCGCCGATGAGTGCGACTGCCGTAAGCCGAAAGTAAAACTGGTGGAACGTTATCTGGCTGAGc                                                                                                                                                             <  2:899924/222‑1 (MQ=255)
                                                                                    ctgAAGCTGCAAAAAGCGGGCTACAAGCTGGTGATGATCACTAATCAGGATGGTCTTGGAACACAAAGTTTCCCACAGGCGGATTTCGATGGCCCGCACAACCTGCTGATGCAGATCTTCCCCTCGCAAGGCGTACAGTTTGATGAAGTGCTGATTTGTCCGCACCTGCCCGCCGCTCAGTGCGACTGCCGTAAGCCGAAAGCAAAACTGGGGGAACGTTATCTGGCGGAGCAAACgca                                                                                                                                                      >  2:382430/1‑237 (MQ=255)
                                                                                          gCTGCAAAAAGCGGGTTAAAAGCGGGGTATGATCACTAATCAGGATGGTCTGGGAACACAAAGTTTCCCACAGGCGGATTTCGATGGCCCGCACAACCTGATGATGCAGATCTTCACCGCGCAAGGCGTACAGTTTGATGAAGTGCTGATTTGTCCGCCCCTGCCCGCCGATGAGTGCGACTGCCGTAAGCCGAAAGTAAAACTGGTGGAACGTTATCTGGCTGAGCAAGCGATGGAt                                                                                                                                                  <  1:950137/238‑1 (MQ=255)
                                                                                              cAAAAAGCGGGCTACAAGCTGGTGATGATCACTAATCAGGATGGTCTTGGAACACAAAGTTTCCCACAGGCGGATTTCGATGGCCCGCACAACCTGATGATGCAGATCTTCACCTCGCAAGGCGTACAGTTTGATGAAGTGCTGATTTGTCCGCACCTGCCCGCCGATGAGTGCGACTGCCGTAAGCCGAAAGTAAAACTGGCGGAACGTTATCTGGCCGAGCAAGCGATGGGTcgc                                                                                                                                               >  1:583886/1‑237 (MQ=255)
                                                                                                 aaaGCGGGCTAAAAGCGGGTGGTGAACAATAACCAGGATGGTCTTGAAACAAAAAGTTCCCCAAGGGCGGATTTCGATGGCCCGCACAACCTGATGATGCAGATCTTCACCTCGCAAGGCGTACAGTTTGATGAAGTGCTGATTTGTCCGCACCTGCCCGCCGATGAGTGCGACTGCCGTAAGCCGAAAGTAAAACTGGTGGAACGTTATCTGGCTGAGCAAGCGATGGATCGCGCTaa                                                                                                                                          <  1:382430/239‑1 (MQ=255)
                                                                                                  taaaaagacaagaaGGTGGTGATGATCAATAAGCAGGATGGTCTTGGAACAAAAAGTTCCCCACAGGCGGATTGCGATGGCCCGCACAACCTGATGATGCAGATCTTCACCTCGCAAGGCGTACAGTTTGATGAAGTGCTGATTTGTCCGCACCTGCCCGCCGATGAGTGCGACTGCCGTAAGCCGAAAGTAAAACTGGTGGAACGTTATCTGGCTGAGCAAGCGATGGATCGCGCTaa                                                                                                                                          <  1:740663/227‑1 (MQ=255)
                                                                                                    agacagttaaaaGTTGTTGATGATCACTAACCGGGGTGTCTTGGAAAAAAAAGTTTTCCACAGGGCGGATTGCGTGGGCCCGCAAACCGGGATGTGGAAGTTCTTCCCCCCGAAGGGCGTACGTTTGGATGAAGTGCTGATTTGCCCGCCCCTGCCCGCCGATGAGTGCGACTGCCGTAAGCCGAAAGTAAAACTGGTGGAACGTTATCTGGCTGAGCAAGCGATGGATCGCGCTaa                                                                                                                                          <  2:196591/227‑1 (MQ=255)
                                                                                                                   ggTGATGATCACTAATCAGGATGGTCTTGGAACACAAAGTTTCCCACAGGCGGATTTCGATGGCCCGCACAACCTGATGATGCAGATCTTCACCTCGCAAGGCGTACAGTTTGATGAAGTGCTGATTTGTCCGCACCTGCCCGCCGATGAGTGCGACTGCCGTAAGCCGAAAGGAAAACTGGTGGAACGTTATCTGGCTGAGCCAGGCAGGGAACCGGCCAAAAGTTAtggtgatgggg                                                                                                                         >  2:522222/1‑230 (MQ=255)
                                                                                                                    gTGATGATCACTAATCAGGATGGTCTTGGAACACAAAGTTTCCCACAGGCGGATTTCGATGGCCCGCACAACCTGATGATGCAGATCTTCACCTCGCAAGGCGTACAGTTTGATGAAGGGCGGATTTGTCCGCACCTGCCCGCCGATGAGTGCGACTGCCGTAAGCCGAAAGGAAAACTGTGGGAACCGTTCCTGGCTGGGAAGGGGGGGGAAGGGCCTAACACTTTAGTGATTGGCg                                                                                                                         >  2:914027/1‑238 (MQ=255)
                                                                                                                     tgatgaTCACTAATCAGGATGGTCTTGGAACACAAAGTTTCCCACAGGCGGATTTCGATGGCCCGCACAACCTGATGATGCAGATCTTCACCTCGCAAGGCGTACAGTTTGATGAAGTGCTGATTTGTCCGCACCTGCCCGCCGATGAGTGCGACTGCCGTAAGCCGAAAGTAAAACTGGTGGAACGTTATCTGGCTGAGCAAGCGATGGATCGCGCTAACAGTTAtgtg                                                                                                                                >  1:196591/1‑230 (MQ=255)
                                                                                                                         gatCACTAATCAGGATGGTCTTGGAACACAAAGTTTCCCACAGGCGGATTTCGATGGCCCGCACAACCTGATGATGCAGATCTTCACCTCGCAAGGCGTACAGTTTGATGAAGTGCTGATTTGTCCGCACCTGCCCGCCGATGAGTGCGACTGCCGTAAGCCGAAAGTAAAACTGGGGGAACGTTATCTGGCTGAGCAAGCGATGGAACCGGCTAACAGTTAttgg                                                                                                                                >  2:740663/1‑223 (MQ=255)
                                                                                                                          atCACTAATCAGGATGGTCTTGGAACACAAAGTTTCCCACAGGCGGATTTCGATGGCCCGCACAACCTGATGATGCAGATCTTCACCTCGCAAGGCGTACAGTTTGATGAAGTGCTGATTTGTCCGCACCTGCCCGCCGATGAGTGCGACTGCCGTAAGCCGAAAGTAAAACTGGTGGAACGTTATCTGGCTGAGCAAGCGATGGATCGCGCTAACAGTTATGTGATTGGCGATcgcgc                                                                                                                  >  1:17217/1‑239 (MQ=255)
                                                                                                                              cTAATCAGGATGGTCTTGGCACACAACGTTTCACACAGGCGGATTTAGATGGCCCGCGCCACCTGATGATGCAGATCTTCACCTCGCAAGGCGTACAGTTTGATGAAGTGCTGATTTGTCCGCACCTGCCCGCCGATGAGTGCGACTGCCGGCAGCCGAAAGAAAAACTGGCGGGCGGTTATCTGCCTGCGCAAACGATGGCGCGCGCTCACGATTATGAGATTGGGGATcgcgcaccc                                                                                                              >  2:120718/1‑235 (MQ=255)
                                                                                                                                 aTCAGGATGGTCTTGGAACACAAAGTTTCCCACAGGCGGATTTCGATGGCCCGCACAACCTGATGATGCAGATCTTCACCTCGCAAGGCGTACAGTTTGATGAAGTGCTGATTTGTCCGCACCTGCCCGCCGATGAGTGCGACTGCCGTAAGCCGAAAGTAAAACTGGTGGAACGTTATCTGGCTGAGCAAGCGATGGATCGCGCTAACAGTTAGTGGATTGGCGATGCGGcgaccgac                                                                                                           >  1:1033833/1‑239 (MQ=255)
                                                                                                                                         gtttttgaaaacaaaaTTTTCCCCCGGGCGTTTTCCGATGCCCCACCAACCCGTATGTTGCATGTTTCCCCCCCGCAAGGCTTACAGTTTGATGAAGTGCGAATTTGCCCGCCCCTGCCCGCCGATGAGTGCGACTGCCGTAAGCCGAAAGTAAAACTGGTGGAACGTTATCTGGCTGAGCAAGCGATGGATCGCGCTAACAGTTATGTGATTGGCGATCGCGCGACCGACATTCAAc                                                                                                    <  2:867300/225‑1 (MQ=255)
                                                                                                                                           tCTTGGAACACAAAGTTTCCCACAGGCGGATTTCGATGGCCCGCACAACCTGATGATGCAGATCTTCACCTCGCAAGGCGTACAGTTTGATGAAGTGCTGATTTGTCCGCACCTGCCCGCCGATGAGTGCGACTGCCGTAAGCCGAAAGGAAAACTGGGGGAACGTTATCTGGCTGAGCAAGGGAAGGATCGCGCTACAAGTATTGTGTTTGGGGATCGCGCGACCGAAATTCAACTgg                                                                                                 >  2:792285/1‑239 (MQ=255)
                                                                                                                                            cTTGGAACACAAAGTTTCCCACAGGCGGATTTCGATGGCCCGCACAACCTGGTGATGCAGATCTTCACCTCGCAAGGCGTACAGTTTGATGAAGTGCTGATTTGTCCGCCCCTGCCCGCCGATGAGTGCGACTGCCGTAAGCCGAAAGTAAAACTGGGGGAACGTCTTATGGCTGAGAACGCGGTGGACCGCGCAAACCGCTATTGGATTGGGGTACCGGcgaccgaacatcacctgcg                                                                                                >  2:638537/1‑227 (MQ=255)
                                                                                                                                             ttGGAACACAAAGTTTCCCACAGGCGGATTTCGATGGCCCGCACAACCTGATGATGCAGATCTTCACCTCGCAAGGCGTACAGTTTGATGAAGTGCTGATTTGTCCGCACCTGCCCGCCGATGAGTGCGACTGCCGTAAGCCGAAAGTAAAACTGGTGGAACGTTATCTGGCTGAGCAAGCGATGGATCGCGCTAACAGTTATGGGATTGGCGATCGCGCGACCGACATTCAACTGGCg                                                                                               >  1:518261/1‑239 (MQ=255)
                                                                                                                                             ttGGAACACAAAGTTTCCCACAGGCGGATTTCGATGGCCCGCACAACCTGATGATGCAGATCTTCACCTCGCAAGGCGTACAGTTTGATGAAGTGCTGATTTGTCCGCACCTGCCCGCCGATGAGTGCGACTGCCGTAAGCCGAAAGTAAAACTGGCGGAACGTTATCTGGCTGGGCAAGCGATGGATCGCGCTAACAGTTATGTGATTGGCGATCGCGCGACCGCCATTCAACTGGCg                                                                                               >  1:608361/1‑239 (MQ=255)
                                                                                                                                             ttGGAACACAAAGTTTCCCACAGGCGGATTTCGATGGCCCGCACAACCTGATGATGCAGATCTTCACCTCGCAAGGCGAACAGTTTGATGAAGTGCTGATTTGTCCGCACCCGCCCGCCGCTGAGTGCGACTGCCGTAAGCCGAAAGTAAAACTGGGGGAACGTTATCTGGCTGAGCAAGCGATGGAACGCGCTAAACGTTATGGGATTGGAGATCCCGCGACCGACATTAAACAgggg                                                                                               >  2:182927/1‑237 (MQ=255)
                                                                                                                                                gAACACAAAGTTTCCCACAGGCGGATTTCGATGGCCCGCACAACCTGATGATGCAGATCTTCACCTCGCAAGGCGTACAGTTTGATGAAGTg                                                                                                                                                                                                                                               >  1:61212/1‑92 (MQ=255)
                                                                                                                                                gAACACAAAGTTTCCCACAGGCGGATTTCGATGGCCCGCACAACCTGATGATGCAGATCTTCACCTCGCAAGGCGTACAGTTTGATGAAGTg                                                                                                                                                                                                                                               <  2:61212/92‑1 (MQ=255)
                                                                                                                                                    acaAAGTTTCCCACAGGCGGATTTCGATGGCCCGCACAACCTGATGATGCAGATCTTCACCTCGCAAGGCGTACAGTTTGATGAAGTGCTGATTTGTCCGCACCTGCCCGCCGATGAGTGCGACTGCCGTAAGCCGAAAGTAAAACTGGTGGAACGTTATCTGGCTGAGCAAGCGATGGATCGCGCTAACAGTTATGTGATTGGCGATCGCGCGACCGACATTCAACTGGCGGaaaa                                                                                          >  1:867300/1‑237 (MQ=255)
                                                                                                                                                     cagaGTTTCCCACAGGCGGATTTCGATGGCCCGCACAACCTGATGATGCAGATCTTCACCTCGCAAGGCGTACAGTTTGATGAAGTGCTGATTTGTCCGCACCTGCCCGCCGATGAGTGCGACTGCCGTAAGCCGAAAGTAAAACTGGTGGAACGTTATCTGGCTGAGCAAGCGATGGATCGCGCTAACAGTTATGGGGTTGGGGGACGCGCCACCGCAATTAAACGGGGGGaaaaaag                                                                                       >  2:913915/4‑236 (MQ=255)
                                                                                                                                                     caAAGTTTCCCACAGGCGGATTTCGATGGCCCGCACAACCTGATGATGCAGATCTTCACCTCGCACGGCGTACAGTTTGATGAAAAGCTGATTTGTCCGCCCCTGCCCGCCGATTAGTGCGCATGCCGGAAGACGAAAGTAAAACTGGTGGAACGTTAGCTGGCGGAGCAAGCGATGGAACGCGCTAACGGTTAAGTGCTTAGCGATCGCGCGGCCCACAACCAAATCGCGGGAAACAt                                                                                       >  2:130162/1‑239 (MQ=255)
                                                                                                                                                     caAAGTTTCCCACAGGCGGATTTCGATGGCCCGCACAACCTGATGATGCAGATCTTCACCTCGCAAGGCGTACAGTTTGATGAAGTGCTGATTTGTCCGCACCTGCCCGCCGATGAGTGCGACTGCCGTAAGCCGAAAGTAAAACTGGTGGAACGTTATCTGGCTGAGCAAGCGATGGATCGCGCTAACa                                                                                                                                        >  1:1038866/1‑190 (MQ=255)
                                                                                                                                                     caAAGTTTCCCACAGGCGGATTTCGATGGCCCGCACAACCTGATGATGCAGATCTTCACCTCGCAAGGCGTACAGTTTGATGAAGTGCTGATTTGTCCGCACCTGCCCGCCGATGAGTGCGACTGCCGTAAGCCGAAAGTAAAACTGGGGGAACGTTATCTGGCTGAGCACGCGATGGAACCGCCAAACAATTATGGGATTGGGGTTGGCGGGGCCGCAATTCAAATGGCGGGAAAAAt                                                                                       >  2:740481/1‑239 (MQ=255)
                                                                                                                                                             cccAGAGGCGGTTTCGGATGGCCCCAAAAACCTGTTGATGAAAATCTTCACCCCGAAAGGCGTACATTTGGATGAAGTGCTGATTTGCCCGCCCCTGCCCGCCGATGAGTGCGACTGCCGTAAGCCGAAAGTAAAACTGGTGGAACGTTATCTGGCTGAGCAAGCGATGGATCGCGCTAACAGTTATGTGATTGGCGATCGCGCGACCGACATTCAACTGGCGGAAAACATGGGCATTa                                                                               <  1:258114/239‑1 (MQ=255)
                                                                                                                                                             cccACAGGCGGATTTCGATGGCCCGCACAACCTGATGATGCAGATCTTCACCTCGCAAGGCGTACAGTTTGATGAAGTGCTGATTTGTCCGCACCTGCCCGCCGATGAGTGCGACTGCCGTAAGCCGAAAGTAAAACTGGTGGAACGTTATCTGGCTGAGCAAGCGATGGATCGCGCTAACAGTTATGTGATTGGCGATCGCGCGACCGACATTcttgctcttatacacatctccgag                                                                                >  1:717773/1‑215 (MQ=255)
                                                                                                                                                             cccACAGGCGGATTTCGATGGCCCGCACAACCTGATGATGCAGATCTTCACCTCGCAAGGCGTACAGTTTGATGAAGTGCTGATTTGTCCGCACCTGCCCGCCGATGAGTGCGACTGCCGTAAGCCGAAAGTAAAACTGGTGGAACGTTATCTGGCTGAGCAAGCGAAGGATCGCGCTAACAGTTATGTGATTGGCGATCGCGCGACCGACATTCAACTGGCGCAAAACATGGGCAtt                                                                                >  1:864551/1‑238 (MQ=255)
                                                                                                                                                                 caGGCGGTTTTCGATGGCCCGCACAGCCTGATTGTGCAGATCCTCCCCTCGCAAGGCGTACAGTTTGATGAATTTCTGATTTTCCCGCCCCTGCCCGCCGATGAGTGCGACTGCCGTAAGCCGAAAGTAAAACTGGTGGAACGTTATCTGGCTGAGCAAGCGATGGATCGCGCTAACAGTTATGTGATTGGCGATCGCGCGACCGACATTCAACTGGCGGAAAACATGGGCATTACTgg                                                                           <  1:330176/239‑1 (MQ=255)
                                                                                                                                                                    gggggttttGGAGGGCCCGCCAAACCTGATGAGGCAGCTCTTCACCTCGAAAGGCGTACATTTTGATGAAGGGCTGTTTTGTCCGCACCGGCCCGCCGATGAGTGCGACTGCCGTAAGCCGAAAGTAAAACTGGTGGAACGTTATCTGGCTGAGCAAGCGATGGATCGCGCTAACAGTTATGTGATTGGCGATCGCGCGACCGACATTCAACTGGCGGAAAACATGGGCATTACTGGtt                                                                         <  1:336096/233‑1 (MQ=255)
                                                                                                                                                                     gCGGATTTCGATGGCACCCCCAACCTGATGATGCAGATCTTCACCTCGCAAGGCGTAGAGTTTGGTGAAGCGCTTGATTGTCCTCCACTGCCCGCCGCTGAGTGCGACTGCCGTAGGCCGAAAGTCAAACTGGGGGGAAGTTATCTGGTCGAGCAAGCGAAGGATCGCGTTACATGTTATGCAGTGGGCGACCGCGGGCCCGCCGTACAACAGGCGGGACACATgggactaacggataa                                                                        >  2:267951/1‑227 (MQ=255)
                                                                                                                                                                                  gcccgacaacccTGGTGATGAAGTCTTTCCCCGCGCAGGGCGTCCGGTTTGATGAAGTGTTGTTTTGTCCGCACTGGCCCGCCGATGAGTGCGACTGCCGTAAGCCGAAAGTAAAACTGGTGGAACGTTATCTGGCTGAGCAAGCGATGGATCGCGCt                                                                                                                                            <  2:61078/148‑1 (MQ=255)
                                                                                                                                                                                   cccGCACAACCTGATGATGCAGATCTTCACCTCGCAAGGCGTACAGTTTGATGAAGTGCTGATTTGTCCGCACCTGCCCGCCGATGAGTGCGACTGCCGTAAGCCGAAAGTAAAACTGGTGGAACGTTATCTGGCTGAGCAAGCGATGGATCGCGCTAACAGTTATGTGATTGGCGATCGGGCGACCGACATTCAACCTGCGGAAAACAAGGGCATTCCGGGTTTACGccagaccccca                                                          >  2:400433/1‑229 (MQ=255)
                                                                                                                                                                                    ctccacaCCCTGATGATGCAGATCTTCCCCTGGCAGGGGGTCGTGTTGGATGAAGGTCTGTTTTGTCCGCCCCTGCCCGCCGATGAGTGCGACTGCCGTAAGCCGAAAGTAAAACTGGTGGAACGTTATCTGGCTGAGCAAGCGATGGATCGCGCTAACAGTTATGTGATTGGCGATCGCGCGACCGACATTCAACTGGCGGAAAACATGGGCATTACTGGTTTACGCTACGACcgcg                                                          <  1:400433/235‑1 (MQ=255)
                                                                                                                                                                                    cggaacaACCGGAGGAGGCAGATCTTCACCCGGAAAGGCGTACAGTTTGATGAAGGGCTAATTGGCCCGCACCTGACCGCCGATGAGTGCGACTGCCGTAAGCCGAAAGTAAAACTGGTGGAACGTTATCTGGCTGAGCAAGCGATGGATCGCGCTAACAGTTATGTGATTGGCGATCGCGCGACCGACATTCAACTGGCGGAAAACATGGGCATTACTGGTTTACGCTACGACcgcg                                                          <  1:947396/234‑1 (MQ=255)
                                                                                                                                                                                     ccgccaaacttgatgtgcaaGATTTTCACCCCGCAGGGCTTCAAGTTTGATGAATGGCTGATTTGCCCGCACCGGCCCGCCGATGAGTGCGACTGCCGTAAGCCGAAAGTAAAACTGGTGGAACGTTATCTGGCTGAGCAAGCGATGGATCGCGCTAACAGTTATGTGATTGGCGATCGCGCGACCGACATTCAACTGGCGGAAAACATGGGCATTACTGGTTTACGCTACGACcgcg                                                          <  1:702412/219‑1 (MQ=255)
                                                                                                                                                                                            ccTGATGATGCAGATCTTCACCTCGCAAGGCGTACAGTTTGATGAAGTGCTGATTTGTCCGCACCTGCCCGCCGATGAGTGCGACTGCCGTAAGCCGAAAGTAAAACTGGTGGAACGTTATCTGGCTGAGCAAGCGATGGATCGCGCTAACAGTTAtg                                                                                                                                  >  1:61078/1‑158 (MQ=255)
                                                                                                                                                                                             gtgatgaGGAATATCTTCCACTCGCAGGGGGTTAAGTTGGAGTAGGGGGTGATTTGCCCGAACCTGCCCGCCGATGAGTGCGACTGCCGTAAGCCGAAAGTAAAACTGGTGGAACGTTATCTGGCTGAGCAAGCGATGGATCGCGCTAACAGTTATGTGATTGGCGATCGCGCGACCGACATTCAACTGGCGGAAAACATGGGCATTACTGGTTTACGCTACGACCGCGAAACCCTGaa                                                <  1:522222/238‑1 (MQ=255)
                                                                                                                                                                                                       cAGATCTTCACCTCGCAAGGCGTACAGTTTGATGAAGTGCTGATTTGTCCGCACCTGCCCGCCGATGAGTGCGACTGCCGTAAGCCGAAAGTAAAACTGGTGGAACGTTATCTGGCTGAGCAAGCGATGGATCGCGCTAACAGTTATGTGATTGGCGATCGCGCGACCGACATTCAACTGGCGGAAAACATGGGCATTACTGGTTTACGCTACGACCGCGAAACCCTGACCTGGCCAAt                                       >  1:109041/1‑239 (MQ=255)
                                                                                                                                                                                                       cAGATCTTCACCTCGCAAGGCGTACAGTTTGATGAAGTGCTGATTTGTCCGCACCTGCCCGCCGATGAGTGCGACTGCCGTAAGCCGAAAGTAAAACTGGTGGAACGTTATCTGGCTGAGCAAGCGATGGATCGCGCTAACAGTTATGTGATTGGCGATCGCGCGACCGACATTCAACTGGCGGAAAACATGGGCATTACTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAAt                                       >  1:1008830/1‑239 (MQ=255)
                                                                                                                                                                                                            cTTCACCTCGCAAGGCGTACAGTTTGATGAAGTGCTGATTTGTCCGCACCTGCCCGCCGATGAGTGCGACTGCCGTAAGCCGAAAGTAAAACTGGTGGAACGTTATCTGGCTGAGCAAGCGATGGATCGCGCTAACAGTTATGTGATTGGCGATCGCGCGACCGACATTCAACTGGCGGAAAACATGGGCATTACGGGTTTACGCCACGACCCCGAAACCCTGAACGGGCCAATAATTg                                  >  2:854947/1‑239 (MQ=255)
                                                                                                                                                                                                                    cGCAAGGCGTACAGTTTGATGAAGTGCTGATTTGTCCGCACCTGCCCGCCGATGAGTGCGACTGCCGTAAGCCGAAAGTAAAACTGGTGGAACGTTATCTGGCTGAGCAAGCGATGGATCGCGCTAACAGTTATGTGATTGGCGATCGCGCGACCGACATTCAACTGGCGGAAAACATGGGCATTACTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGACCa                           >  1:531058/1‑238 (MQ=255)
                                                                                                                                                                                                                    cGCAAGGCGTACAGTTTGATGAAGTGCTGATTTGTCCGCACCTGCCCGCCGATGAGTGCGACTGCCGTAAGCCGAAAGTAAAACTGGTGGAACGTTATCTGGCTGAGCAAGCGATGGATCGCGCTAACAGTTATGTGATTGCCGATCGGGCGACGGACATTCAACTGGCGGAAAACATGGGCCTTACTGGTTTACGCTACGACCGCGAAACCCTGACCTGGCCAATGATggggcgtaac                          >  2:621668/1‑232 (MQ=255)
                                                                                                                                                                                                                     cgcggggtctACTGTTTGATGGTGGGTTGTTTGGCCCGCCCCTGCCCGCCGATGAGGGCGACTGCCGTAACCCGAAAGTAAAACTGGTGGAACGTTATCTGGCTGAGCAAGCGATGGATCGCGCTAACAGTTATGTGATTGGCGATCGCGCGACCGACATTCAACTGGCGGAAAACATGGGCATTACTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCaa                          <  2:183308/230‑1 (MQ=255)
                                                                                                                                                                                                                     ccaagggCGTACGTTTGGATGAAGTTCTGATTTGCCCGCCCCTGCCCGCCGATGAGTGCGACTGCCGTAAGCCGAAAGTAAAACTGGTGGAACGTTATCTGGCTGAGCAAGCGATGGATCGCGCTAACAGTTATGTGATTGGCGATCGCGCGACCGACATTCAACTGGCGGAAAACATGGGCATTACTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCaa                          <  1:257019/234‑1 (MQ=255)
                                                                                                                                                                                                                      cAAGGCGTACAGTTTGATGAAGTGCTGATTTGTCCGCACCTGCCCGCCGATGAGTGCGACTGCCGTAAGCCGAAAGTAAAACTGGTGGAACGTTATCTGGCTGAGCAAGCGATGGATCGCGCTAACAGTTATGTGATTGGCGATCGCGCGACCGACATTCAACTGGCGGAAAACATGGGCATTACTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATAATTGGCGGACAAcc                        >  2:257019/1‑238 (MQ=255)
                                                                                                                                                                                                                                     tgatgaAGTGCTGATTTGTCCGCACCTGCCCGCCGATGAGTGCGACTGCCGTAAGCCGAAAGTAAAACTGGTGGAACGTTATCTGGCTGAGCAAGCGATGGATCGCGCTAACAGTTATGTGATTGGCGATCGCGCGACCGACATTCAACTGGCGGAAAACATGGGCATTACTGGTTTACGCTACGACCCCGAAACCCTGAACTGCCCAATGATTGGGGAGAAAACCACCAGACGTGccc          >  2:799919/1‑239 (MQ=255)
                                                                                                                                                                                                                                             tGCTGATTTGTCCGCACCTGCCCGCCGATGAGTGCGACTGCCGTAAGCCGAAAGTAAAACTGGTGGAACGTTATCTGGCTGAGCAAGCGATGGATCGCGCTAACAGTTATGTGATTGGCGATCGCGCGACCGACATTCAACTGGCGGAAAACATGGGCATTACTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGGGGGCAACTCACCAGACGGGACCGTTaaccc  >  2:865198/1‑235 (MQ=255)
                                                                                                                                                                                                                                             tGCTGATTTGTCCGCACCTGCCCGCCGATGAGTGCGACTGCCGTAAGCCGAAAGTAAAACTGGTGGAACGTTATCTGGCTGAGCAAGCGATGGATCGCGCTAACAGTTATGTGATTGGCGATCGCGCGACCGACATTCAACTGGCGGAAAACATGGGCATTACTGGTTTACGCCACGACCGCGAAACCCTGAACTGGCCAATGATTGGGGCGCAACTACCCAGACGGAACCGTTACGc   >  2:737814/1‑238 (MQ=255)
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                              
CCCTGATTAGCGAACCGCCGAGTGATTTTCAGGTGGACCGTTTTGATAAACTCGCCTTTGAACCGGGCGTGATCCCGGAACTGCTGAAGCTGCAAAAAGCGGGCTACAAGCTGGTGATGATCACTAATCAGGATGGTCTTGGAACACAAAGTTTCCCACAGGCGGATTTCGATGGCCCGCACAACCTGATGATGCAGATCTTCACCTCGCAAGGCGTACAGTTTGATGAAGTGCTGATTTGTCCGCACCTGCCCGCCGATGAGTGCGACTGCCGTAAGCCGAAAGTAAAACTGGTGGAACGTTATCTGGCTGAGCAAGCGATGGATCGCGCTAACAGTTATGTGATTGGCGATCGCGCGACCGACATTCAACTGGCGGAAAACATGGGCATTACTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCT  >  NC_000913/2093505‑2093974

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