Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,817,090 G→A D139N (GAC→AAC)  yicC → putative RNase adaptor protein YicC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,817,0900GA98.4% 214.1 / ‑4.7 65D139N (GAC→AAC) yicCputative RNase adaptor protein YicC
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (29/34);  total (30/35)
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.61e-01

GCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGCCAGATGTTAGCGCGCAAGGTGAGCTGATCCTCAACGAAAAACTGGCTAAACAGCTGGTAACTGCCGCGAACTGGGTAAAAATGCAGAGTGACGAAGGGGAAATCAACCCGGTTGATATTCTACGCTGGCCGGGCGTGATGGCAGCCCAGGAGCAGGATCTTGACGCCATTGCCGCTGAAATTCTCGCGGCGCTGGATGGTACGCTGGACGACTTTATTGTCGCGCGCGAAACCGAAGGTCAGGCACTGAAAGCATTGATCGAGCAGCGTCTGGAAGGCGTCACCGCCGAAGTGGTCAAAGTCCGCTCCCATATGCCGGAAATCCTGCAATGGCAGCGTGAGCGTCTGGTCGCGAAGCTGGAAGATGCTCAGGTGCAACTGGAAAACAACCGTCTGGAGCAGGAACTGGTTCTGCTGGCACAACGAATTGACGTTGCCGAAGAACTGGATCGCC  >  NC_000913/3816849‑3817336
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                       
gcgcgGTAAAGTGGAATGTACCCTGCGCTATGAGCCAGATGTTAGCGCGCAAGGTGAGCTGATCCTCAACGAAAAACTGGCTAAACAGCTGGTAACTGCCGCGAACTGGGTAAAAATGCAGAGTGACGAAGGGGAAATCAACCCGGTTGATATTCTACGCTGGCCGGGCGTGATGGCAGCCCAGGAGCAGGATCTTGGCGCCATTGCCGCTGAAATTCTCGCGGCGCTGGATGGTACGCTGAACGACTTTa                                                                                                                                                                                                                                               <  2:429259/251‑1 (MQ=255)
                 gTACCCTGCGCTATGAGCCAGATGTTAGCGCGCAAGGTGAGCTGATCCTCAACGAAAAACTGGCTAAACAGCTGGAAACTGCCGCGAACTGGGTAAAAATGCAGAGTGACGAAGGGGAAATCAACCCGGTTGATATTCTACGCTGGCCGGGCGTGATGGCAGCCCAGGAGCAGGATCTTGACGCCATTGCCGCTGAAATTCTCGCGGCGCTGGATGGTACGCTGAACGACTTTATTGTCGCGCGCGAAAcc                                                                                                                                                                                                                              <  2:375188/251‑1 (MQ=255)
                    cccTGCGCTATGAGCCAGATGTTAGCGCGCAAGGTGAGCTGATCCTCAACGAAAAACTGGCTAAACAGCTGGTAACTGCCGCGAACTGGGTAAAAAGGCAGAGTGACGAAGGGGAAATCAACCCGGTTGATATTCTACGCTGGCCGGGCGTGATGGCAGCCCAGGAGCAGGATCTTGACGTCATTGCCGCTGAAATTCTCTCGGCGCTGGATGGTACGCTGAAATTCTTTATTGTCGCGTGCGAAACCGaa                                                                                                                                                                                                                           <  2:57293/251‑1 (MQ=255)
                        gcgcTATGAGCCAGATGTTAGCGCGCAAGGTGAGCTGATCCTCAACGAAAAACTGGCTAAACAGCTGGTAACTGCCGCGAACTGGGTAAAAATGCAGAGTGACGAAGGGGAAATCAACCCGGTTGATATTCTACGCTGGCCGGGCGTGATGGCAGCCCAGGAGCAGGATCTTGACGCCATTGCCGCTGAAATTCTCGCGGCGCTGGATGGTACGCTGAACGACTTTATTGTCGCGCGCGAAACCGAAGGt                                                                                                                                                                                                                        <  1:4846/250‑1 (MQ=255)
                                      aTGTTAGCGCGCAAGGTGAGCTGATCCTCAACGAAAAACTGGCTAAACAGCTGGTAACTGCCGCGAACTGGGTAAAAATGCAGAGTGACGAAGGGGAAATCAACCCGGTTGATATTCTACGCTGGCCGGGCGTGATGGCAGCCCAGGAGCAGGATCTTGACGCCATTGCCGCTGAAATTCTCGCGGCGCTGGATGGTACGCTGAACGACTTTATTGTCGCGCGCGAAACCGAAGGTCAGGCACTGAAAGc                                                                                                                                                                                                          >  1:108216/1‑250 (MQ=255)
                                       tGTTAGCGCGCAAGGTGAGCTGATCCTCAACGAAAAACTGGCTAAACAGCTGGTAACTGCCGCGAACTGGGTAAAAATGCAGAGTGACGAAGGGGAAATCAACCCGGTTGATATTCTACGCTGGCCGGGCGTGATGGCAGCCCAGGAGCAGGATCTTGACGCCATTGCCGCTGAAATTCTCGCGGCGCTGGATGGTACGCTGAACGACTTTATTGTCGCGCGCGAAACCGAAGGTCAGGCACTGAAAGCAt                                                                                                                                                                                                        >  1:11617/1‑251 (MQ=255)
                                             caagcacggTGAGCTGATCCTCAACGAAAAACTGGCTAAACAGCTGGTAACTGCCGCGAACTGGGTAAAAATGAAGAGAGACGAAGGGGAAATCAACCCGGTTGATATTCTACGCTGGCCGGGCGTGATGGCAGCCCAGGAGCAGGATCTTGACGCCATTGCCGCTGAAATTCTCGCGGCGCTGGATGGTACGCTGAACGACTTTATTGTCGCGCGCGAAACCGACGGTCAGGCACTGACAGCATTGATCg                                                                                                                                                                                                  >  2:56867/8‑251 (MQ=255)
                                                       gAGCTGATCCTCAACGAAAAACTGGCTAAACAGCTGGTAACTGCCGCGAACTGGGTAAAAATGCAGAGTGACGAAGGGGAAATCAACCCGGTTGATATTCTACGCTGGCCGGGCGTGATGGCAGCCCAGGAGCAGGATCTTGACGCCATTGCCGCTGAAATTCTCGCGGCGCTGGATGGTACGCTGAACGACTTTATTGTCGCGCGCGAAACCGAAGGTCAGGCACTGAAAGCATTGATCGAGCAGCGTCt                                                                                                                                                                                        <  1:239521/251‑1 (MQ=255)
                                                            gATCCGCAACGAAAAACCGGCTAAACAGCTGGTAACTGCCGCGAACTGTGTAAAAATGCAGAGTGACGAAGGGGAAATCAACCCGGTTGATATTCTACGCTGGCCGGGCGTGATGGCAGCCCAGGAGCAGGATCTTGACGCCATTGCCGCTGAAATTCTCGCGGCGCTGGATGGTACGCTGAACGACTTTATTGTCGCGCGCGAAACCGAAGGTCAGGCACTGAAAGCATTGATCGAGCAGCGTCTGGaa                                                                                                                                                                                    <  2:129834/250‑1 (MQ=255)
                                                                    aCGAAAAACTGGCTAAACAGCTGGTAACTGCCGCGAACTGGGTAAAAATGCAGAGTGACGAAGGGGAAATCAACCCGGTTGATATTCTACGCTGGCCGGGCGTGATGGCAGCCCAGGAGCAGGATCTTGACGCCATTGCCGCTGAAATTCTCGCGGCGCTGGATGGTACGCTGAACGACTTTATTGTCGCGCGCGAAACCGAAGGTCAGGCACTGAAAGCATTGATCGAGCAGCGTCTTGAAGGCGTCAcc                                                                                                                                                                           <  1:88211/251‑1 (MQ=255)
                                                                         aaaCTGGCTAAACAGCTGGTAACTGCCGCGAACTGGGTAAAAATGCAGAGTGACGAAGGGGAAATCAACCCGGTTGATATTCTACGCTGGCCGGGCGTGATGGCAGCCCAGGAGCAGGATCTTGACGCCATTGCCGCTGAAATTCTCGCGGCGCTGGATGGTACGCTGAACGACTTTATTGTCGCGCGCGAAACCGAAGGTCAGGCACTGAAAGCATTGATCGAGCAGCGTCTGGAAGGCGTCAccgccg                                                                                                                                                                       >  1:330193/1‑250 (MQ=255)
                                                                            cTGGCTAAACAGCTGGTAACTGCCGCGAACTGGGTAAACATGCAGAGTGACGAAGGGGAAATCAACCCGGTTGATATTCTACGCTGGCCGGGCGTGATGGCAGCCCAGGAGCAGGATCTTGACGCCATTGCCGCTGAAATTCTCGCGGCGCTGGATGGTACGCTGAACGCCTTTATTGTCGCGCGCGAAACCGAAGGTCAGGCACTGAAAGCATTGCTCGAGCAGCGTCTGGCAGGCGTCACCGCCGAAgg                                                                                                                                                                   >  2:119670/1‑250 (MQ=255)
                                                                                  aaaCAGCTGGTACCTGCCGCGAACTGGGTAAAAATGCAGAGTGGCGAAGGGGAAATCAACCCGGTTGATATTCTACGCTGGCCGGGCGTGATGGCAGCCCAGGAGCAGGATCTGGACGCCATTGCCGCTGAAATTCTCGCGGCGCTGGAGGGTACGCTGAACGACTTTATTGTCGCGCGCGAAACCGAAGGTCAGGCACTGAAAGCATTGATCGAGCAGCGTCTGGAAGGCGTCACCGCCGAAGTGGTCa                                                                                                                                                              <  2:11617/250‑1 (MQ=255)
                                                                                    aCAGCTGGTAACTGCCGCGAACTGGGTAAAAATGCAGAGTGACGAAGGGGAAATCAACCCGGTTGATATTCTACGCTGGCCGGGCGTGATGGCAGCCCAGGAGCAGGATCTTGACGCCATTGCCGCTGAAATTCTCGCGGCGCTGGATGGTACGCTGAACGACTTTATTGTCGCGCGCGAAACCGAAGGTCAGGCACTGAAAGCATTGATCGAGCAGCGTCTGGAAGGCGTCACCGCCGAAGTGGTCAAAg                                                                                                                                                           <  1:413145/251‑1 (MQ=255)
                                                                                     cAGCTGGTAACTGCCGCGAACTGGGTAAAAATGCAGAGTGACGAAGGGGAAATCAACCCGGTTGATATTCTACGCTGGCCGGGCGTGATGGCAGCCCAGGAGCAGGATCTTGACGCCATTGCCGCTGAAATTCTCGCGGCGCTGGATGGTACGCTGAACGACTTTATTGTCGCGCGCGAAACCGAAGGTCAGGCACTGAAAGCATTGATCGAGCAGCg                                                                                                                                                                                           >  1:407509/1‑218 (MQ=255)
                                                                                     cAGCTGGTAACTGCCGCGAACTGGGTAAAAATGCAGAGTGACGAAGGGGAAATCAACCCGGTTGATATTCTACGCTGGCCGGGCGTGATGGCAGCCCAGGAGCAGGATCTTGACGCCATTGCCGCTGAAATTCTCGCGGCGCTGGATGGTACGCTGAACGACTTTATTGTCGCGCGCGAAACCGAAGGTCAGGCACTGAAAGCATTGATCGAGCAGCg                                                                                                                                                                                           <  2:407509/218‑1 (MQ=255)
                                                                                           gTAACTGCCGCGAACTGGGTAAAAATGCAGAGTGACGAAGGGGAAATCAACCCGGTTGATATTCTACGCTGGCCGGGCGTGATGGCAGCCCAGGAGCAGGATCTTGACGCCATTGCCGCTGAAATTCTCGCGGCGCTGGATGGTACGCTGAACGACTTTATTGTCGCGCGCGAAACCGAAGGTCAGGCACTGAAAGCATTGATCGAGCAGCGTCTGGAAGGCGTCACCGCCGAAGTGGTCAAAGTCCGCTc                                                                                                                                                    >  2:314783/1‑251 (MQ=255)
                                                                                              aCTGCCGCGAACTGGGTAAAAATGCAGAGTGACGAAGGGGAACTCAACCCGGTTGATATTCTACGCTGGCCGGGCGTGATGGCAGCCCAGGAGCAGGATCTTGACGCCATTGCCGCTGAAATTCTCGCGGCGCTGGATGGTACGCTGAACGACTTTATTGTCGCGCGCGAAACCGAAGGTCAGGCACTGAAAGCATTGATCGAGCAGCGTCTGGAAGGCGTCACCGCCGCAGTGGTCAAAGGCCGCTccc                                                                                                                                                  >  1:355443/1‑250 (MQ=255)
                                                                                                    gcgAACTGGGTAAAAATGCAGAGTGACGACGGGGAAATCAACACGGTTGATATTCTACGCTGGCCGGGCGCGATGGCCGCCAAGGAGCAGGATCTTGACGCCATCGCCGCTGAAATTATCGAGGCGCTGGATGGTACGCTGAACGACTTTATTGTCGCGCGCGAAACCGACGGTCAGGCACTGAAAGCATTGATCGAGCAACGTCTGGAAGGCGTCACCGCCGAGGTGGTCAAAGTCCGCTCCCATATGcc                                                                                                                                           >  1:206677/1‑251 (MQ=255)
                                                                                                         cTGGGTAAAAATGCAGAGTGACGAAAGGGAAATCAACCCGGTTGATATTCTACGCTGGCAGGGCGTGATGGCCGCCCAGGCGCAGGATCTTGACGCCATTGCCGCTGAAATTCTCGCGGCGCTGGATGGTACGCTGAACGACTTTATTGTCGCGCGCGAAACCGACGGTCAGGCACTGAAAGCATTGATCGAGCAGCGTCTGGAAGGCGTCACCGCCGAAGTGGTCAAAGTCCGCTCCCATCTGCCGGaa                                                                                                                                       >  2:368457/1‑250 (MQ=255)
                                                                                                             gTAAAAATGCAGAGTGACGAAGGGGAAATCAACCCGGTTGATATTCTACGCTGGCCGGGCGTGATGGCAGCCCAGGAGCAGGATCTTGACGCCATTGCCGCTGAAATTCTCGCGGCGCTGGATGGTACGCTGAACGACTTTATTGTCGCGCGCGAAACCGAAGGTCAGGCACTGAAAGCATTGATCGAGCAGCGTCTGGAAGGCGTCACCGCCGAAGTGGTCAAAGTCCGCTCCCATATGCCGGAAATCCt                                                                                                                                  <  2:330828/251‑1 (MQ=255)
                                                                                                                aaaaTGCAGAGTGACGAAGGGGAAATCAACCCGGTTGATATTCTACGCTGGCCGGGCGTGATGGCAGCCCAGGAGCAGGATCTTGACGCCATTGCCGCTGAAATTCTCGCGGCGCTGGATGGTACGCTGAACGACTTTATTGTCGCGCGCGAAACCGAAGGTCAGGCACTGAAAGCATTGATCGAGCAGCGTCTGGAAGGCGTCACCGCCGAAGTGGTCAAAGTCCGCTCCCATATGCCGGAAATCCTGCa                                                                                                                               <  1:185230/251‑1 (MQ=255)
                                                                                                                 aaaTGCAGAGGGACGACGGGGAAATCAACCCGGTTGATATTCTACGCTGGCCGGGCGTGATGGCAGCCCAGGAGCAGGATCTTGACGCCATTGCCGCTGAAATTCTCGCGGCGCTGGATGGTACGCTGAACGACTTTATTGTCGCGCGCGAAACCGAAGGTCAGGCACTGAAAGCATTGATCGAGCAGCGTCTGGAAGGCGTCACCGCCGAAGTGGTCAAAGTCCGCTCCCATATGCCGGAAATCCTGCa                                                                                                                               >  2:418963/1‑250 (MQ=255)
                                                                                                                       agagTGACGAAGGGGAAATCAACCCGGTTGATATTCTACGCTGGCCGGGCGTGATGGCAGCCCAGGAGCAGGATCTTGACGCCATTGCCGCTGAAATTCTCGCGGCGCTGGATGGTACGCTGAACGACTTTATTGTCGCGCGCGAAACCGAAGGTCAGGCACTGAAAGCATTGATCGAGCAGCGTCTGGAAGGCGTCACCGCCGAAGTGGTCAAAGTCCGCTCCCATATGCCGGAAATCCTGCACTGGc                                                                                                                          >  1:173208/1‑249 (MQ=255)
                                                                                                                         agTGACGAAGGGGAAATCAACCCGGTTGATATTCTACGCTGGCCGGGCGTGATGGCAGCCCAGGAGCAGGATCTTGACGCCATTGCCGCTGAAATTCTCGCGGCGCTGGATGGTACGCTGAACGACTTTATTGTCGCGCGCGAAACCGAAGGTCAGGCACTGAAAGCATTGATCGAGCAGCGTCTGGAAGGCGTCACCGCCGAAGTGGTCAAAGTCCGCTCCCATATGCCGGAAATCCTGCAATGGCAGCg                                                                                                                      >  1:228324/1‑251 (MQ=255)
                                                                                                                            gACGAAGGGGAAACCAACCCGGTTGATATTCTACGCGGGCCGGGCGTGATGGCAGCCCAGGAGCCGGATCTTGACGCCATTGCCGCTGAAAGTCTCGCGGCGCTGGATGGTACGCTGAACGACTTTATTGTCGCGCGCGAAACCGAAGGTCAGGCACTGAAAGCATTGATCGAGCAGCGTCTGGAAGGCGTCACCGCCGAAGTGGTCAAAGTCCGCTCCCATATGCCGGAAATCCTGCAATGGCAGCGTGa                                                                                                                   <  2:335984/251‑1 (MQ=255)
                                                                                                                              cGAAGGGGAAATCAACCCGGTTGATATTCTACGCTGGCCGGGCGTGATGGCAGCCCAGGAGCAGGATCTTGACGCCATTGCCGCTGAAATTCTCGCGGCGCTGGATGGTACGCTGAACGACTTTATTGTCGCGCGCGAAACCGAAGGTCAGGCACTGAAAGCATTGATCGAGCAGCGTCTGGAAGGCGTCACCGCCGAAGTGGTCAAAGTCCGCTCCCATATGCCGGAAATCCTGCAATGGCAGCGTGaac                                                                                                                 <  2:330193/251‑3 (MQ=255)
                                                                                                                              cGAAGGGGAAATCAACCCGGTTGATATTCTACGCTGGCCGGGCGTGATGGCAGCCCAGGAGCAGGATCTTGACGCCATTGCCGCTGAAATTCTCGCGGCGCTGGATGGTACGCTGAACGACTTTATTGTCGCGCGCGAAACCGAAGGTCAGGCACTGAAAGCATTGATCGAGCAGCGTCTGGAAGGCGTCACCGCCGAAGTGGTCAAAGTCCGCTCCCATATGCCGGAAATCCTGCAATGGCAGCGTGaac                                                                                                                 >  1:358158/1‑249 (MQ=255)
                                                                                                                                 aGGGGAAATCAACCCGGTTGATATTCTACGCTGGCCGGGCGTTATGGCAGCCCAGGAGCAGGATCTTGACGCCATTGCCGCTGAAATTCTCGCGGCGCTGGATGGTACGCTGAACGACTTTATTGTCGCGCGCGAAACCGAAGGTCAGGCACTGAAAGCATTGATCGAGCAGCGGCTGGAAGGCGTCACCGCCGAAGTGGTCAAAGTCCGCTCCCATATGCCGGAAATCCTGCAATGGCAGCGTGaacgct                                                                                                              <  2:228324/251‑6 (MQ=255)
                                                                                                                                     gAAATCAACCCGGTTGATATTCTACGCTGGCCGGGCGTGATGGCAGCCCAGGAGCAGGATCTTGACGCCATTGCCGCTGAAATTCTCGCGGCGCTGGATGGTACGCTGAACGACTTTATTGTCGCGCGCGAAACCGAAGGTCAGGCACTGAAAGCATTGATCGAGCAGCGTCTGGAAGGCGTCACCGCCGAAGTGGTCAAAGTCCGCTCCCATATGCCGGAAATCCTGCAATGGCAGCGTGaacgcttggt                                                                                                          >  1:147233/1‑242 (MQ=255)
                                                                                                                                           agcccGGTTGATATTCTACGCGGGCCGGGCGTGATGGCAGCCCAGGAGCAGGATCTTGACGCCATTGCCGCTGAAATTCTCGCGGCGCTGGATGGTACGCTGAACGACTTTATTGTCGCGCGCGAAACCGAAGGTCAGGCACTGAAAGCATTGATCGAGCAGCGCCTGGAAGGCGTCACCGCCGAAGTGGTCAAAGTCCGCTCCCATATGCCGGAAATCCTGCAATGGCAGCGTGAACGCTTGGTCGCGaa                                                                                                    <  2:256546/249‑1 (MQ=255)
                                                                                                                                              ccGGTTGATATTCTACGCTGGCCGGGCGTGATGGCAGCCCAGGAGCAGGATCTTGACGCCATTGCCGCTGAAATTCTCGCGGCGCTGGATGGTACGCTGAACGACTTTATTGTCGCGCGCGAAACCGAAGGTCAGGCACTGAAAGCATTGATCGAGCAGCGTCTGGAAGGCGTCACCGCCGAAGTGGTCAAAGTCCGCTCCCATATGCCGGAAATCCTGCAATGGCAGCGTGAACGCTTGGTCGCGAAGCt                                                                                                 >  1:292365/1‑251 (MQ=255)
                                                                                                                                               cGGTTGATATTCTACGCTGGCCGGGCGTGATGGCAGCCCAGGAGCAGGATCTTGACGCCCTTGCCGCTGAAATTCTCGAGGCGCTGGATGGTCCGCTGCACGACTTTATTGTCGCGCGCGAAACCGAAGGTCAGGCACTGAAAGCATTGATCGAGCAGCGTCTGGAAGGCGTCACCGCCGAAGTGGTCAAAGTCCGCTCCCATATGCCGGAAATCCTGCAATGGCAGCgggaacgcttggacgcgtagctg                                                                                                >  2:159501/1‑229 (MQ=255)
                                                                                                                                                   agATATTCTACGCTGGCCGGGCGTGCTGGCAGCCCAGGAGCAGGATCTTGAAGCCATTGCCGCTGAAATTCTCGCGGCGCTGGATGGTACGCTGAACGACTTTATTGTCGCGCGCGAAACCGAAGGTCAGGCACTGAAAGCATTGATCGAGCAGCGTCTGGAAGGCGTCACCGCCGAAGTGGTCAAAGTCCGCTCCCATATGCCGGAAATACTGCAATGGCAGCGTGAACGCTTGGTAGCGAAGCTGGacg                                                                                            >  2:58230/2‑249 (MQ=255)
                                                                                                                                                    gATATTCTACGCTGGCCGGGCGTGATGGCAGCCCAGGAGCAGGATCTTGACGCCATTGCCGCTGAAATTCTCGCGGCGCTGGATGGTACGCTGAACGACTTTATTGTCGCGCGCGAAACCGAAGGTCAGGCACTGAAAGCATTGATCGAGCAGCGTCTGGAAGGCGTCACCGCCGAAGTGGTCAAAGTCCGCTCCCATATGCCGGAAATCCTGCAATGGCAGCGTGAACGCTTGGTCGCGAAGCTGGAAGa                                                                                           <  2:83501/251‑1 (MQ=255)
                                                                                                                                                    gATATTCTACGCGGGCCGGGCGGGATGGCAGCCCAGGAGCAGGATCTGGACGCCAGTGCCGCTGAAATTCTCGCGGCGCTGGATGGTACGCTGAACGACTTTAGTGGCGCGCGCGAAACCGAAGGTCAGGCACTGAAAGCATTGATCGAGCAGCGTCTGGAAGGCGTCACCGCCGAAGTGGTCAAAGTCCGCTCCCATATGCCGGAAATCCTGCAATGGCAGCGTGAACGCTTGGTAGCGAAGCTGGAAg                                                                                            <  2:173208/250‑1 (MQ=255)
                                                                                                                                                      tatTCTACGCTGGCCGGGCGTGATGGCAGCCCAGGAGCAGGATCTTGACGCCATTGCCGCTGAAATTCTCGCGGCGCTGGATGGTACGCTGAACGACTTTATTGTCGCGCGCGAAACCGAAGGTCAGGCACTGAAAGCATTGATCGAGCAGCGTCTGGAAGGCGTCACCGCCGAAGTGGTCAAAGTCCGCTCCCATATGCCGGAAATCCTGCAATGGCAGCGTGAACGCTTGGTCGCGAAGCTGGAAGATg                                                                                         >  2:68698/1‑251 (MQ=255)
                                                                                                                                                      tatTCTACGCTGGCCGGGCGTGATGGCAGCCCAGGAGCAGGATCTTGACGCCATTGCCGCTGAAATTCTCGCGGCGCTGGATGGTACGCTGAACGACTTTATTGTCGCGCGCGAAACCGAAGGTCAGGCACTGAAAGCATTGATCGAGCAGCGTCTGGAAGGCGTCACCGCCGAAGTGGTAAAAGTCCGCTCCCATATGCCGGAAATCCTGCAATGGCAGCGTGAACGCTTGGTCGCGAAGCTGGACGATg                                                                                         >  2:107278/1‑251 (MQ=255)
                                                                                                                                                            aCGCTGGCCGGGCGTGATGGCAGCCCAGGAGCAGGATCTTGACGCCATTGCCGCTGAAATTCTCGCGGCGCTGGATGGTACGCTGAACGACTTTATTGTCGCGCGCGAAACCGAAGGTCAGGCACTGAAAGCATTGATCGAGCAGCGTCTGGAAGGCGTCACCGCCGAAGTGGTCAAAGTCCGCTCCCATATGCCGGAAATCCTGCAATGGCAGCGTGACCGCTTGGTCGCGAAGCTGGAAGATGCTCa                                                                                     <  2:126746/249‑1 (MQ=255)
                                                                                                                                                            aCGCTGGCCGGGCGTGATGGCAGCCCAGGAGCAGGATCTTGACGCCATTGCCGCTGAAATGCTCGCGGCGCTGGATGGTACGCTGAACGACTTTATTGTCGCGCGCGAAACCGAAGGTCAGGCACTGAAAGCATTGATCGAGCAGCGTCTGGAAGGCGTCACCGCCGAAGTGGTCAAAGTCCGCTCCCATATGCCGGAAATCCTGCAATGGCAGCGTGAACGCTTGGTCGCGAAGCTGGAAGATGCTCa                                                                                     <  2:376341/249‑1 (MQ=255)
                                                                                                                                                                 ggCCGGGCGTGATGGCAGCCCAGGAGCAGGATCTTGACGCCATTGCCGCTGGAATTCTCGCGGCGCTGGATGGTACGCTGAACGACTTTATTGTCGCGCGCGAAACCGAAGGTCAGGCACTGAAAGCATTGATCGAGCAGCGTCTGGAAGGCGTCACCGCCGAAGTGGTCAAAGTCCGCTCCCATATGCCGGAAATCCTGCAATGGCAGCGTGAACGCTTGGTCGCGAAGCTGGAAGATGCTCAGGttcag                                                                              <  1:147541/251‑5 (MQ=255)
                                                                                                                                                                   ccGGGCGTGATGGCAGCCCAGGAGCAGGATCTTGACGCCATTGCCGCTGAAATTCTCGCGGCGCTGGATGGTACGCTGAACGTCTTTATTGTCGCGCGCGAAACCGAAGGTCAGGCACTGAAAGCATTGATCGAGCAGCGTCTGGAAGGCTTCACCGCCGAAGTGGTCAAAGTCCGCTCCCATATGCCGGAAATCCTGCAATGGCAGCGTGAACGCTTGGTCGCGAAGCTGGAAGATGCTCAGGttcagtt                                                                            <  1:29951/251‑7 (MQ=255)
                                                                                                                                                                     gggCGTGATGGCAGCCCAGGAGCAGGATCTTGACGCCATTGCCGCTGAAATTCTCGCGGCGCTGGATGGTACGCTGAACGACTTTATTGTCGCGCGCGAAACCGAAGGTCAGGCACTGAAAGCATTGATCGAGCAGCGTCTGGAAGGCGTCACCGCCGAAGTGGTCAAAGTCCGCTCCCATATGCCGGAAATCCTGCAATGGCAGCGTGAACGCTTGGTCGCGAAGCTGGAAGATGCTCAgggtcagttgg                                                                          >  2:101866/1‑242 (MQ=255)
                                                                                                                                                                      ggCGTGATGGCAGCCCAGGAGCAGGATCTTGACGCCATTGCCGCTGAAATTCTCGCGGCGCTGGATGGTACGCTGAACGACTTTATTGTCGCGCGCGAAACCGAAGGTCAGGCACTGAAAGCATTGATCGAGCAGCGTCTGGAAGGCGTCACCGCCGAAGTGGTCAAAGTCCGCTCCCATATGCCGGAAATCCTGCAATGGCAGCGTGAACGCTTGGTCGCGAAGCTGGAAGATGCTCAGGttcagttgga                                                                         <  2:147233/251‑10 (MQ=255)
                                                                                                                                                                      ggCGTGATGGCAGCCCAGGAGCAGGATCTTGACGCCATTGCCGCTGAAATTCTCGCGGCGCTGGATGGTACGCTGAACGACTTTATTGTCGCGCGCGAAACCGAAGGTCAGGCACTGAAAGCATTGATCGAGCAGCGTCTGGAAGGCGTCACCGCCGAAGTGGTCAAAGTCCGCTCCCATATGCCGGAAATCCTGCAATGGCAGCGTGAACGCTTGGTCGCGAAGCTGGAAGATGCTCAGGttcagttgg                                                                          >  2:21702/1‑242 (MQ=255)
                                                                                                                                                                        cgagATGGCCGCCCAGGAGCAGGATCTTGACGCCATTGCCGCTGAAATTCTCGCGGCGCTGGATGGTGCGCTGAACGACTTTATTGTCGCGCGCGAAACCGAAGGTCAGGCACTGCAAGCATTGATCGAGCAGCGTCTGGAAGGCGTCACCGCCGAAGTGGTCAAAGTCCGCTCCCATATGCCGGAAATCCTGCAATGGCAGCGTGAACGCTTGGTCGCGGAGCTGGACGATGCTCAGGttcagttggaa                                                                        >  2:54922/4‑240 (MQ=255)
                                                                                                                                                                                     cAGGAGCAGGATCTTGACGCCATTGCCGCTGAAATTCTCGCGGCGCTGGATGGTACGCTGAACGACTTTATTGTCGCGCGCGAAACCGAAGGTCAGGCACTGAAAGCATTGATCGAGCAGCGTCTGGAAGGCGTCACCGCCGAAGTGGTCAAAGTCCGCTCCCATATGCCGGAAATCCTGCAATGGCAGCGTGAACGCTTGGTCGCGAAGCTGGAAGATGCTCAGGTTCAGTTGGAAAACAACCGTCTgg                                                           >  1:309861/1‑250 (MQ=255)
                                                                                                                                                                                      aGGAGCAGGATCTTGACGCCATTGCCGCTGAAATTCTCGCGGCGCTGGATGGTACGCTGAACGACTTTATTGTCGCGCGCGAAACCGAAGGTCAGGCACTGAAAGCATTGATCGAGCAGCGTCTGGAAGGCGTCACCGCCGAAGTGGTCAAAGTCCGCTCCCATATGCCGGAAATCCTGCAATGGCAGCGTGAACGCTTGGTCGCGAAGCTGGAAGATGCTCAGGTTCAGTTGGAAAACAACCGTCTgg                                                           >  1:52493/1‑249 (MQ=255)
                                                                                                                                                                                       ggAGCAGGATCTTGACGCCATTGCCGCTGAAATTCTCGCGGCGCTGGATGGTACGCTGAATGACTTTATTGTCGCGCGCGAAACCGAAGGTCAGGCACTGAAAGCATTGATCGAGCAGCGTCTGGAAGGCGTCACCGCCGAAGTGGTCAAAGTCCGCTCCCATATGCCGGAAATCCTGCAATGGCAGCGTGAACGCTTGGTCGCGAAGCTGGAAGATGCTCAGGTTCAGTTGGAAAACAACCGTCTGGAGc                                                        <  1:38142/251‑1 (MQ=255)
                                                                                                                                                                                          gCAGGATCTTGACGCCATTGCCGCTGAAATTCTCGCGGCGCTGGATGGTACGCTGAACGACTTTATTGTCGCGCGCGAAACCGAAGGTCAGGCACTGAAAGCATTGATCGAGCAGCGTCTGGAAGGCGTCACCGCCGAAGTGGTCAAAGTCCGCTCCCATATGCCGGAAATCCTGCAATGGCAGCGTGAACGCTTGGTCGCGAAGCTg                                                                                                >  1:99614/1‑208 (MQ=255)
                                                                                                                                                                                          gCAGGATCTTGACGCCATGGCCGCTGAAATTCTCGCGGCGCTGGATGGTACGCTGAACGACTTTATTGTCGCGCGCGAAACCGAAGGTCAGGCACTGAAAGCATTGATCGAGCAGCGTCTGGAAGGCGTCACCGCCGAAGTGGTCAAAGTCCGCTCCCATATGCCGGAAATCCTGCAATGGCAGCGTGAACGCTTGGTCGCGAAGCTg                                                                                                <  2:99614/208‑1 (MQ=255)
                                                                                                                                                                                             ggATCTTGACGCCATTGCCGCTGAAATTCTCGCGGCGCTGGATGGTACGCTGAACGACTTTATTGTCGCGCGCGAAACCGAAGGTCAGGCACTGAAAGCATTGATCGAGCAGCGTCGGGAAGGCGTCACCGCCGAAGTGGTCAAAGTCCGCTCCCATATGCCGGAAATCCTGCAATGGCAGCGTGAACGCTTGGTCGCGAAGCTGGAAGATGCTCAGGTTCAGTTGGAAAACAACCGTCTGGAGCAGGagt                                                  <  2:108216/251‑3 (MQ=255)
                                                                                                                                                                                                   ggACGCCAGTGCCGCTGAAATGCTCGCGGCGCTGGATGGTACGCTGAACGACTTTATTGTCGCGCGCGAAACCGAAGGTCAGGCACTGAAAGCATTGATCGAGCAGCGTCTGGAAGGCGTCACCGCCGAAGTGGTCAAAGTCCGCTCCCATATGCCGGAAATCCTGCAAGGGCAGCGTGAACGCTTGATCGCGAAGCTGGAAGATGCTCAGGTTCAGTTGGAAAACAACCGTCTGGAGCAGGagttggtt                                             <  2:125693/249‑8 (MQ=255)
                                                                                                                                                                                                             gCCGCTGAAATGCTCGCGGCGCTGGATGGTACGCGGAGCTCCTTTATTGTCGCGCGCGAAACCGAAGGTCAGGCACTGAAAGCAGTGTTCGCGCAGCGTCTGGAAGGCGTCACCGCCGAAGTGGTCAAAGTCCGCGCCCATATGCCGGAAATCCTGCAATGGCAGCGTGAACGCTTGGTCGCGAAGCTGGAAGATGCTCAGGTTCAGTTGGAAAACAACCGTCTGGAGCAGGAGTTGGTTCTGCTGGcaca                                  <  2:99004/251‑1 (MQ=255)
                                                                                                                                                                                                              ccGCTGAAATTCTCGCGGCGCTGGATGGTACGCTGAACGACTTTATTGTCGCGCGCGAAACCGAAGGTCAGGCACTGAAAGCATTGATCGAGCAGCGTCTGGAAGGCGTCACCGCCGAAGTGGTCAAAGTCCGCTCCCATATGCCGGAAATCCTGCAATGGCAGCGTGAACGCTTGGTCGCGAAGCTGGAAGATGCTCAGGTTCAGTTGGAAAACAACCGTCTGGAGCAGGagtt                                                 <  2:368359/235‑4 (MQ=255)
                                                                                                                                                                                                              ccGCTGAAATTCTCGCGGCGCTGGATGGTACGCTGAACGACTTTATTGTCGCGCGCGAAACCGAAGGTCAGGCACTGAAAGCATTGATCGAGCAGCGTCTGGAAGGCGTCACCGCCGAAGTGGTCAAAGTCCGCTCCCATATGCCGGAAATCCTGCAATGGCAGCGTGAACGCTTGGTCGCGAAGCTGGAAGATGCTCAGGTTCAGTTGGAAAACAACCGTCTGGAGCAGGagtt                                                 >  1:368359/1‑232 (MQ=255)
                                                                                                                                                                                                                   gAAATTCTCGCGGCGCTGGATGGTACGCTGAACGACTTTATTGTCGCGCGCGAAACCGAAGGTCAGGCACTGAAAGCATTGATCGAGCCGCGTCTGGAAGGCGTAACCGCCGAAGTGGTCAAAGTCAGCTCCCATATGCCGGAAATCCTGCAATGGCAGCGTGAAAGCTTGGTCGCGAAGCTGGAAGATGCTCAGGTTCAGTTGGAAAACAAACGTCTGGAGCAGGAGTTGGTTCTGCTGGCACAACGAAt                            >  1:194897/1‑251 (MQ=255)
                                                                                                                                                                                                                    aacgtctcGCTGGGCTGGATGGTACGCTGAACGACTTTATGGTCGCGCGCGAAACCGAAGGTCAGGCACTGAAAGCATTTATCGAGCAGCGTCTGGAAGGCGTCACCGCCGAAGTGGTCAAAGTCCGCTCCCATATGCCGGAAATCCTGCAATGGCAGCGTGAACGCTTGGTCGCGAAGCTGGAAGATGCTCAGGTTCAGTTGGAAAACAACCGTCTGGAGCAGGAGTTGGTTCTGCTGGCACAACGAAtt                           <  2:141893/247‑1 (MQ=255)
                                                                                                                                                                                                                    aaaTTCTCGCGGCGCTGGATGGTACGCTGAACGACTTTATTGTCGCGCGCGAAACCGAAGGTCAGGCACTGAACGCATTGATCGAGCAGCGTCTGGAAGGCGTCACCGCCGAAGTGGTCAAAGTCCGCTCCCATATGCCGGAAATCCTGCAATGGCAGCGTGAACGCTTGGTCGCGAAGCTGGAAGATGCTCAGGTTCAGTTGGAAACCAACCGTCTGGAGCAGGAGTTGGTTCTGCTGGCACAACGAAtt                           >  1:389470/1‑251 (MQ=255)
                                                                                                                                                                                                                           cgcgGCGCTGGATGGTACGCTGAACGACTTTATTGTCGCGCGCGAAACCGAAGGTCAGGCACTGAAAGCATTGATCGAGCAGCGTCTGGAAGGCGTCACCGCCGAAGTGGTCAAAGTCCGCTCCCATATGCCGGAAATCCTGCAATGGCAGCGTGAACGCTTGGTCGCGAAGCTGGAAGATGCTCAGGTTCAGTTGGAAAACAACCGTCTGGAGCAGGAGTTGGTTCTGCTGGCACAACGAATTGACGt                      <  1:257778/249‑1 (MQ=255)
                                                                                                                                                                                                                             cggcgCTGGATGGGACGCTGAACGACTTTATTGTCGCGCGCGAAACCGAAGGTCAGGCACTGAAAGCATTGATCGAGCAGCGTCTGGAAGGCGTCACCGCCGAAGTGGTCAAAGTCCGCTCCCATATGCCGGAAATCCTGCAATGGCAGCGTGAACGCTTGGTCGCGAAGCTGGAAGATGCTCAGGTTCAGTTGGAAAACAACCGTCTGGAGCAGTATTTGGTTCTGCTGGCACAACGAATTGACGTTGcc                  <  1:119670/251‑1 (MQ=255)
                                                                                                                                                                                                                              ggcgCTGGATGGTACGCTGAACGACTTTAGTGTCGCGCGCGAAACCGAAGGCCAGGCACTGAAAGCATTGATCGAGCCGCGTCTGGAAGGCGTCACCGCCGAAGTGGTCAAAGTCCGCTCCCATATGCCGGAAATCCTGCAATGGCAGCGTGAACGCTTGGTCGCGAAGCTGGAAGATGCTCAGGTTCAGTTGGAAAACAACCGTCTGGAGCAGGAGTTGGTTCTGCTGGCACAACGAATTGACGTTGCCg                 <  2:358158/251‑1 (MQ=255)
                                                                                                                                                                                                                                      aTGGTACGCTGAACGACTTTATTGTCGCGCGCGAAACCGAAGGTCTGACACTGAAAGCATTGATCGAGCAGAGTCTGGAAGGCGTCACCGCCGACGTGGTCAAAGTCCGCTCCCCTATGCCGGCACTACTGCAATGGAAGCGTGAACGCTTGGTCGCGAAGCTGGAAGATGCTCAGGTTCAGTTGGAAAACAACCGTCTGGAGCAGGAGTTGGTTCTGCTGGCACAACGAATTGACGTTGCAGAAGAACTg         >  1:293051/1‑251 (MQ=255)
                                                                                                                                                                                                                                       ggttACGCTGAACGACTTTATTGTCGCGCGCGAAACCGAAGGTCAGGCACTGAAAGCATTGATCGAGCAGCGTCTGGAAGGCGTCACCGCCGAAGTGGTCAAAGTCCGCTCCCATATGCCGGAAATCCTGCAATGGCAGCGTGAACGCTTGGTCGCGAAGCTGGAAGATGCTCAGGTTCAGTTGGAAAACAACCGTCTGGAGCAGGAGTTGGTTCTGCTGGCACAACGAATTGACGTTGCCGAAGAACTgg        <  1:412192/248‑1 (MQ=255)
                                                                                                                                                                                                                                             gCTGAACGACTGTATTGTCGCGCGCGAAACCGAAGGCCAGGCACTGAAAGCATTGATCGAGCAGCGGCTGGAAGGCGTCACCGCCGAAGTGGTCAAAGTCCGCTCCCATATGCCGGAAATCCTGCAATGGCAGCGTGAACGCTTGGGCGCGAAGCTGGAAGATGCTCAGGTTCAGTTGGAAAACAACCGTCTGGAGCAGGAGTTGGTTCTGCTGGCACAACGAATTGACGTTGCCGAAGAACTGGACCGcc  <  2:361613/251‑1 (MQ=255)
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                       
GCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGCCAGATGTTAGCGCGCAAGGTGAGCTGATCCTCAACGAAAAACTGGCTAAACAGCTGGTAACTGCCGCGAACTGGGTAAAAATGCAGAGTGACGAAGGGGAAATCAACCCGGTTGATATTCTACGCTGGCCGGGCGTGATGGCAGCCCAGGAGCAGGATCTTGACGCCATTGCCGCTGAAATTCTCGCGGCGCTGGATGGTACGCTGGACGACTTTATTGTCGCGCGCGAAACCGAAGGTCAGGCACTGAAAGCATTGATCGAGCAGCGTCTGGAAGGCGTCACCGCCGAAGTGGTCAAAGTCCGCTCCCATATGCCGGAAATCCTGCAATGGCAGCGTGAGCGTCTGGTCGCGAAGCTGGAAGATGCTCAGGTGCAACTGGAAAACAACCGTCTGGAGCAGGAACTGGTTCTGCTGGCACAACGAATTGACGTTGCCGAAGAACTGGATCGCC  >  NC_000913/3816849‑3817336

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