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

GCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGCCAGATGTTAGCGCGCAAGGTGAGCTGATCCTCAACGAAAAACTGGCTAAACAGCTGGTAACTGCCGCGAACTGGGTAAAAATGCAGAGTGACGAAGGGGAAATCAACCCGGTTGATATTCTACGCTGGCCGGGCGTGATGGCAGCCCAGGAGCAGGATCTTGACGCCATTGCCGCTGAAATTCTCGCGGCGCTGGATGGTACGCTGGA‑CGACTTTATTGTCGCGCGCGAAACCGAAGGTCAGGCACTGAAAGCATTGATCGAGCAGCGTCTGGAAGGCGTCACCGCCGAAGTGGTCAAAGTCCGCTCCCATATGCCGGAAATCCTGCAATGGCAGCGTGAGCGTCTGGTCGCGAAGCTGGAAGATGCTCAGGTGCAACTGGAAAACAACCGTCTGGAGCAGGAACTGGTTCTGCTGGCACAACGAATTGACGTTGCCGAAGAACTGGATCGCC  >  NC_000913/3816849‑3817336
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                       
gcgcgGTAAAGTGGAATGTACCCTGCGCTATGAGCCAGATGTTAGCGCGCAAGGTGAGCTGATCCTCAACGAAAAACTGGCTAAACAGCTGGTAACTGCCGCGAACTGGGTAAAAATGCAGAGTGACGAAGGGGAAATCAACCCGGTTGATATTCTACGCTGGCCGGGCGTGATGGCAGCCCAGGAGCAGGATCTTGGCGCCATTGCCGCTGAAATTCTCGCGGCGCTGGATGGTACGCTGAA‑CGACTTTa                                                                                                                                                                                                                                               <  2:429259/251‑1 (MQ=255)
                 gTACCCTGCGCTATGAGCCAGATGTTAGCGCGCAAGGTGAGCTGATCCTCAACGAAAAACTGGCTAAACAGCTGGAAACTGCCGCGAACTGGGTAAAAATGCAGAGTGACGAAGGGGAAATCAACCCGGTTGATATTCTACGCTGGCCGGGCGTGATGGCAGCCCAGGAGCAGGATCTTGACGCCATTGCCGCTGAAATTCTCGCGGCGCTGGATGGTACGCTGAA‑CGACTTTATTGTCGCGCGCGAAAcc                                                                                                                                                                                                                              <  2:375188/251‑1 (MQ=255)
                    cccTGCGCTATGAGCCAGATGTTAGCGCGCAAGGTGAGCTGATCCTCAACGAAAAACTGGCTAAACAGCTGGTAACTGCCGCGAACTGGGTAAAAAGGCAGAGTGACGAAGGGGAAATCAACCCGGTTGATATTCTACGCTGGCCGGGCGTGATGGCAGCCCAGGAGCAGGATCTTGACGTCATTGCCGCTGAAATTCTCTCGGCGCTGGATGGTACGCTGAA‑ATTCTTTATTGTCGCGTGCGAAACCGaa                                                                                                                                                                                                                           <  2:57293/251‑1 (MQ=255)
                        gcgcTATGAGCCAGATGTTAGCGCGCAAGGTGAGCTGATCCTCAACGAAAAACTGGCTAAACAGCTGGTAACTGCCGCGAACTGGGTAAAAATGCAGAGTGACGAAGGGGAAATCAACCCGGTTGATATTCTACGCTGGCCGGGCGTGATGGCAGCCCAGGAGCAGGATCTTGACGCCATTGCCGCTGAAATTCTCGCGGCGCTGGATGGTACGCTGAA‑CGACTTTATTGTCGCGCGCGAAACCGAAGGt                                                                                                                                                                                                                        <  1:4846/250‑1 (MQ=255)
                                      aTGTTAGCGCGCAAGGTGAGCTGATCCTCAACGAAAAACTGGCTAAACAGCTGGTAACTGCCGCGAACTGGGTAAAAATGCAGAGTGACGAAGGGGAAATCAACCCGGTTGATATTCTACGCTGGCCGGGCGTGATGGCAGCCCAGGAGCAGGATCTTGACGCCATTGCCGCTGAAATTCTCGCGGCGCTGGATGGTACGCTGAA‑CGACTTTATTGTCGCGCGCGAAACCGAAGGTCAGGCACTGAAAGc                                                                                                                                                                                                          >  1:108216/1‑250 (MQ=255)
                                       tGTTAGCGCGCAAGGTGAGCTGATCCTCAACGAAAAACTGGCTAAACAGCTGGTAACTGCCGCGAACTGGGTAAAAATGCAGAGTGACGAAGGGGAAATCAACCCGGTTGATATTCTACGCTGGCCGGGCGTGATGGCAGCCCAGGAGCAGGATCTTGACGCCATTGCCGCTGAAATTCTCGCGGCGCTGGATGGTACGCTGAA‑CGACTTTATTGTCGCGCGCGAAACCGAAGGTCAGGCACTGAAAGCAt                                                                                                                                                                                                        >  1:11617/1‑251 (MQ=255)
                                             caagcacggTGAGCTGATCCTCAACGAAAAACTGGCTAAACAGCTGGTAACTGCCGCGAACTGGGTAAAAATGAAGAGAGACGAAGGGGAAATCAACCCGGTTGATATTCTACGCTGGCCGGGCGTGATGGCAGCCCAGGAGCAGGATCTTGACGCCATTGCCGCTGAAATTCTCGCGGCGCTGGATGGTACGCTGAA‑CGACTTTATTGTCGCGCGCGAAACCGACGGTCAGGCACTGACAGCATTGATCg                                                                                                                                                                                                  >  2:56867/8‑251 (MQ=255)
                                                       gAGCTGATCCTCAACGAAAAACTGGCTAAACAGCTGGTAACTGCCGCGAACTGGGTAAAAATGCAGAGTGACGAAGGGGAAATCAACCCGGTTGATATTCTACGCTGGCCGGGCGTGATGGCAGCCCAGGAGCAGGATCTTGACGCCATTGCCGCTGAAATTCTCGCGGCGCTGGATGGTACGCTGAA‑CGACTTTATTGTCGCGCGCGAAACCGAAGGTCAGGCACTGAAAGCATTGATCGAGCAGCGTCt                                                                                                                                                                                        <  1:239521/251‑1 (MQ=255)
                                                            gATCCGCAACGAAAAACCGGCTAAACAGCTGGTAACTGCCGCGAACTGTGTAAAAATGCAGAGTGACGAAGGGGAAATCAACCCGGTTGATATTCTACGCTGGCCGGGCGTGATGGCAGCCCAGGAGCAGGATCTTGACGCCATTGCCGCTGAAATTCTCGCGGCGCTGGATGGTACGCTGAA‑CGACTTTATTGTCGCGCGCGAAACCGAAGGTCAGGCACTGAAAGCATTGATCGAGCAGCGTCTGGaa                                                                                                                                                                                    <  2:129834/250‑1 (MQ=255)
                                                                    aCGAAAAACTGGCTAAACAGCTGGTAACTGCCGCGAACTGGGTAAAAATGCAGAGTGACGAAGGGGAAATCAACCCGGTTGATATTCTACGCTGGCCGGGCGTGATGGCAGCCCAGGAGCAGGATCTTGACGCCATTGCCGCTGAAATTCTCGCGGCGCTGGATGGTACGCTGAA‑CGACTTTATTGTCGCGCGCGAAACCGAAGGTCAGGCACTGAAAGCATTGATCGAGCAGCGTCTTGAAGGCGTCAcc                                                                                                                                                                           <  1:88211/251‑1 (MQ=255)
                                                                         aaaCTGGCTAAACAGCTGGTAACTGCCGCGAACTGGGTAAAAATGCAGAGTGACGAAGGGGAAATCAACCCGGTTGATATTCTACGCTGGCCGGGCGTGATGGCAGCCCAGGAGCAGGATCTTGACGCCATTGCCGCTGAAATTCTCGCGGCGCTGGATGGTACGCTGAA‑CGACTTTATTGTCGCGCGCGAAACCGAAGGTCAGGCACTGAAAGCATTGATCGAGCAGCGTCTGGAAGGCGTCAccgccg                                                                                                                                                                       >  1:330193/1‑250 (MQ=255)
                                                                            cTGGCTAAACAGCTGGTAACTGCCGCGAACTGGGTAAACATGCAGAGTGACGAAGGGGAAATCAACCCGGTTGATATTCTACGCTGGCCGGGCGTGATGGCAGCCCAGGAGCAGGATCTTGACGCCATTGCCGCTGAAATTCTCGCGGCGCTGGATGGTACGCTGAA‑CGCCTTTATTGTCGCGCGCGAAACCGAAGGTCAGGCACTGAAAGCATTGCTCGAGCAGCGTCTGGCAGGCGTCACCGCCGAAgg                                                                                                                                                                   >  2:119670/1‑250 (MQ=255)
                                                                                  aaaCAGCTGGTACCTGCCGCGAACTGGGTAAAAATGCAGAGTGGCGAAGGGGAAATCAACCCGGTTGATATTCTACGCTGGCCGGGCGTGATGGCAGCCCAGGAGCAGGATCTGGACGCCATTGCCGCTGAAATTCTCGCGGCGCTGGAGGGTACGCTGAA‑CGACTTTATTGTCGCGCGCGAAACCGAAGGTCAGGCACTGAAAGCATTGATCGAGCAGCGTCTGGAAGGCGTCACCGCCGAAGTGGTCa                                                                                                                                                              <  2:11617/250‑1 (MQ=255)
                                                                                    aCAGCTGGTAACTGCCGCGAACTGGGTAAAAATGCAGAGTGACGAAGGGGAAATCAACCCGGTTGATATTCTACGCTGGCCGGGCGTGATGGCAGCCCAGGAGCAGGATCTTGACGCCATTGCCGCTGAAATTCTCGCGGCGCTGGATGGTACGCTGAA‑CGACTTTATTGTCGCGCGCGAAACCGAAGGTCAGGCACTGAAAGCATTGATCGAGCAGCGTCTGGAAGGCGTCACCGCCGAAGTGGTCAAAg                                                                                                                                                           <  1:413145/251‑1 (MQ=255)
                                                                                     cAGCTGGTAACTGCCGCGAACTGGGTAAAAATGCAGAGTGACGAAGGGGAAATCAACCCGGTTGATATTCTACGCTGGCCGGGCGTGATGGCAGCCCAGGAGCAGGATCTTGACGCCATTGCCGCTGAAATTCTCGCGGCGCTGGATGGTACGCTGAA‑CGACTTTATTGTCGCGCGCGAAACCGAAGGTCAGGCACTGAAAGCATTGATCGAGCAGCg                                                                                                                                                                                           >  1:407509/1‑218 (MQ=255)
                                                                                     cAGCTGGTAACTGCCGCGAACTGGGTAAAAATGCAGAGTGACGAAGGGGAAATCAACCCGGTTGATATTCTACGCTGGCCGGGCGTGATGGCAGCCCAGGAGCAGGATCTTGACGCCATTGCCGCTGAAATTCTCGCGGCGCTGGATGGTACGCTGAA‑CGACTTTATTGTCGCGCGCGAAACCGAAGGTCAGGCACTGAAAGCATTGATCGAGCAGCg                                                                                                                                                                                           <  2:407509/218‑1 (MQ=255)
                                                                                           gTAACTGCCGCGAACTGGGTAAAAATGCAGAGTGACGAAGGGGAAATCAACCCGGTTGATATTCTACGCTGGCCGGGCGTGATGGCAGCCCAGGAGCAGGATCTTGACGCCATTGCCGCTGAAATTCTCGCGGCGCTGGATGGTACGCTGAA‑CGACTTTATTGTCGCGCGCGAAACCGAAGGTCAGGCACTGAAAGCATTGATCGAGCAGCGTCTGGAAGGCGTCACCGCCGAAGTGGTCAAAGTCCGCTc                                                                                                                                                    >  2:314783/1‑251 (MQ=255)
                                                                                              aCTGCCGCGAACTGGGTAAAAATGCAGAGTGACGAAGGGGAACTCAACCCGGTTGATATTCTACGCTGGCCGGGCGTGATGGCAGCCCAGGAGCAGGATCTTGACGCCATTGCCGCTGAAATTCTCGCGGCGCTGGATGGTACGCTGAA‑CGACTTTATTGTCGCGCGCGAAACCGAAGGTCAGGCACTGAAAGCATTGATCGAGCAGCGTCTGGAAGGCGTCACCGCCGCAGTGGTCAAAGGCCGCTccc                                                                                                                                                  >  1:355443/1‑250 (MQ=255)
                                                                                                    gcgAACTGGGTAAAAATGCAGAGTGACGACGGGGAAATCAACACGGTTGATATTCTACGCTGGCCGGGCGCGATGGCCGCCAAGGAGCAGGATCTTGACGCCATCGCCGCTGAAATTATCGAGGCGCTGGATGGTACGCTGAA‑CGACTTTATTGTCGCGCGCGAAACCGACGGTCAGGCACTGAAAGCATTGATCGAGCAACGTCTGGAAGGCGTCACCGCCGAGGTGGTCAAAGTCCGCTCCCATATGcc                                                                                                                                           >  1:206677/1‑251 (MQ=255)
                                                                                                         cTGGGTAAAAATGCAGAGTGACGAAAGGGAAATCAACCCGGTTGATATTCTACGCTGGCAGGGCGTGATGGCCGCCCAGGCGCAGGATCTTGACGCCATTGCCGCTGAAATTCTCGCGGCGCTGGATGGTACGCTGAA‑CGACTTTATTGTCGCGCGCGAAACCGACGGTCAGGCACTGAAAGCATTGATCGAGCAGCGTCTGGAAGGCGTCACCGCCGAAGTGGTCAAAGTCCGCTCCCATCTGCCGGaa                                                                                                                                       >  2:368457/1‑250 (MQ=255)
                                                                                                             gTAAAAATGCAGAGTGACGAAGGGGAAATCAACCCGGTTGATATTCTACGCTGGCCGGGCGTGATGGCAGCCCAGGAGCAGGATCTTGACGCCATTGCCGCTGAAATTCTCGCGGCGCTGGATGGTACGCTGAA‑CGACTTTATTGTCGCGCGCGAAACCGAAGGTCAGGCACTGAAAGCATTGATCGAGCAGCGTCTGGAAGGCGTCACCGCCGAAGTGGTCAAAGTCCGCTCCCATATGCCGGAAATCCt                                                                                                                                  <  2:330828/251‑1 (MQ=255)
                                                                                                                aaaaTGCAGAGTGACGAAGGGGAAATCAACCCGGTTGATATTCTACGCTGGCCGGGCGTGATGGCAGCCCAGGAGCAGGATCTTGACGCCATTGCCGCTGAAATTCTCGCGGCGCTGGATGGTACGCTGAA‑CGACTTTATTGTCGCGCGCGAAACCGAAGGTCAGGCACTGAAAGCATTGATCGAGCAGCGTCTGGAAGGCGTCACCGCCGAAGTGGTCAAAGTCCGCTCCCATATGCCGGAAATCCTGCa                                                                                                                               <  1:185230/251‑1 (MQ=255)
                                                                                                                 aaaTGCAGAGGGACGACGGGGAAATCAACCCGGTTGATATTCTACGCTGGCCGGGCGTGATGGCAGCCCAGGAGCAGGATCTTGACGCCATTGCCGCTGAAATTCTCGCGGCGCTGGATGGTACGCTGAA‑CGACTTTATTGTCGCGCGCGAAACCGAAGGTCAGGCACTGAAAGCATTGATCGAGCAGCGTCTGGAAGGCGTCACCGCCGAAGTGGTCAAAGTCCGCTCCCATATGCCGGAAATCCTGCa                                                                                                                               >  2:418963/1‑250 (MQ=255)
                                                                                                                       agagTGACGAAGGGGAAATCAACCCGGTTGATATTCTACGCTGGCCGGGCGTGATGGCAGCCCAGGAGCAGGATCTTGACGCCATTGCCGCTGAAATTCTCGCGGCGCTGGATGGTACGCTGAA‑CGACTTTATTGTCGCGCGCGAAACCGAAGGTCAGGCACTGAAAGCATTGATCGAGCAGCGTCTGGAAGGCGTCACCGCCGAAGTGGTCAAAGTCCGCTCCCATATGCCGGAAATCCTGCACTGGc                                                                                                                          >  1:173208/1‑249 (MQ=255)
                                                                                                                         agTGACGAAGGGGAAATCAACCCGGTTGATATTCTACGCTGGCCGGGCGTGATGGCAGCCCAGGAGCAGGATCTTGACGCCATTGCCGCTGAAATTCTCGCGGCGCTGGATGGTACGCTGAA‑CGACTTTATTGTCGCGCGCGAAACCGAAGGTCAGGCACTGAAAGCATTGATCGAGCAGCGTCTGGAAGGCGTCACCGCCGAAGTGGTCAAAGTCCGCTCCCATATGCCGGAAATCCTGCAATGGCAGCg                                                                                                                      >  1:228324/1‑251 (MQ=255)
                                                                                                                            gACGAAGGGGAAACCAACCCGGTTGATATTCTACGCGGGCCGGGCGTGATGGCAGCCCAGGAGCCGGATCTTGACGCCATTGCCGCTGAAAGTCTCGCGGCGCTGGATGGTACGCTGAA‑CGACTTTATTGTCGCGCGCGAAACCGAAGGTCAGGCACTGAAAGCATTGATCGAGCAGCGTCTGGAAGGCGTCACCGCCGAAGTGGTCAAAGTCCGCTCCCATATGCCGGAAATCCTGCAATGGCAGCGTGa                                                                                                                   <  2:335984/251‑1 (MQ=255)
                                                                                                                              cGAAGGGGAAATCAACCCGGTTGATATTCTACGCTGGCCGGGCGTGATGGCAGCCCAGGAGCAGGATCTTGACGCCATTGCCGCTGAAATTCTCGCGGCGCTGGATGGTACGCTGAA‑CGACTTTATTGTCGCGCGCGAAACCGAAGGTCAGGCACTGAAAGCATTGATCGAGCAGCGTCTGGAAGGCGTCACCGCCGAAGTGGTCAAAGTCCGCTCCCATATGCCGGAAATCCTGCAATGGCAGCGTGaac                                                                                                                 <  2:330193/251‑3 (MQ=255)
                                                                                                                              cGAAGGGGAAATCAACCCGGTTGATATTCTACGCTGGCCGGGCGTGATGGCAGCCCAGGAGCAGGATCTTGACGCCATTGCCGCTGAAATTCTCGCGGCGCTGGATGGTACGCTGAA‑CGACTTTATTGTCGCGCGCGAAACCGAAGGTCAGGCACTGAAAGCATTGATCGAGCAGCGTCTGGAAGGCGTCACCGCCGAAGTGGTCAAAGTCCGCTCCCATATGCCGGAAATCCTGCAATGGCAGCGTGaac                                                                                                                 >  1:358158/1‑249 (MQ=255)
                                                                                                                                 aGGGGAAATCAACCCGGTTGATATTCTACGCTGGCCGGGCGTTATGGCAGCCCAGGAGCAGGATCTTGACGCCATTGCCGCTGAAATTCTCGCGGCGCTGGATGGTACGCTGAA‑CGACTTTATTGTCGCGCGCGAAACCGAAGGTCAGGCACTGAAAGCATTGATCGAGCAGCGGCTGGAAGGCGTCACCGCCGAAGTGGTCAAAGTCCGCTCCCATATGCCGGAAATCCTGCAATGGCAGCGTGaacgct                                                                                                              <  2:228324/251‑6 (MQ=255)
                                                                                                                                     gAAATCAACCCGGTTGATATTCTACGCTGGCCGGGCGTGATGGCAGCCCAGGAGCAGGATCTTGACGCCATTGCCGCTGAAATTCTCGCGGCGCTGGATGGTACGCTGAA‑CGACTTTATTGTCGCGCGCGAAACCGAAGGTCAGGCACTGAAAGCATTGATCGAGCAGCGTCTGGAAGGCGTCACCGCCGAAGTGGTCAAAGTCCGCTCCCATATGCCGGAAATCCTGCAATGGCAGCGTGaacgcttggt                                                                                                          >  1:147233/1‑242 (MQ=255)
                                                                                                                                           agcccGGTTGATATTCTACGCGGGCCGGGCGTGATGGCAGCCCAGGAGCAGGATCTTGACGCCATTGCCGCTGAAATTCTCGCGGCGCTGGATGGTACGCTGAA‑CGACTTTATTGTCGCGCGCGAAACCGAAGGTCAGGCACTGAAAGCATTGATCGAGCAGCGCCTGGAAGGCGTCACCGCCGAAGTGGTCAAAGTCCGCTCCCATATGCCGGAAATCCTGCAATGGCAGCGTGAACGCTTGGTCGCGaa                                                                                                    <  2:256546/249‑1 (MQ=255)
                                                                                                                                              ccGGTTGATATTCTACGCTGGCCGGGCGTGATGGCAGCCCAGGAGCAGGATCTTGACGCCATTGCCGCTGAAATTCTCGCGGCGCTGGATGGTACGCTGAA‑CGACTTTATTGTCGCGCGCGAAACCGAAGGTCAGGCACTGAAAGCATTGATCGAGCAGCGTCTGGAAGGCGTCACCGCCGAAGTGGTCAAAGTCCGCTCCCATATGCCGGAAATCCTGCAATGGCAGCGTGAACGCTTGGTCGCGAAGCt                                                                                                 >  1:292365/1‑251 (MQ=255)
                                                                                                                                               cGGTTGATATTCTACGCTGGCCGGGCGTGATGGCAGCCCAGGAGCAGGATCTTGACGCCCTTGCCGCTGAAATTCTCGAGGCGCTGGATGGTCCGCTGCA‑CGACTTTATTGTCGCGCGCGAAACCGAAGGTCAGGCACTGAAAGCATTGATCGAGCAGCGTCTGGAAGGCGTCACCGCCGAAGTGGTCAAAGTCCGCTCCCATATGCCGGAAATCCTGCAATGGCAGCgggaacgcttggacgcgtagctg                                                                                                >  2:159501/1‑229 (MQ=255)
                                                                                                                                                   agATATTCTACGCTGGCCGGGCGTGCTGGCAGCCCAGGAGCAGGATCTTGAAGCCATTGCCGCTGAAATTCTCGCGGCGCTGGATGGTACGCTGAA‑CGACTTTATTGTCGCGCGCGAAACCGAAGGTCAGGCACTGAAAGCATTGATCGAGCAGCGTCTGGAAGGCGTCACCGCCGAAGTGGTCAAAGTCCGCTCCCATATGCCGGAAATACTGCAATGGCAGCGTGAACGCTTGGTAGCGAAGCTGGacg                                                                                            >  2:58230/2‑249 (MQ=255)
                                                                                                                                                    gATATTCTACGCTGGCCGGGCGTGATGGCAGCCCAGGAGCAGGATCTTGACGCCATTGCCGCTGAAATTCTCGCGGCGCTGGATGGTACGCTGAA‑CGACTTTATTGTCGCGCGCGAAACCGAAGGTCAGGCACTGAAAGCATTGATCGAGCAGCGTCTGGAAGGCGTCACCGCCGAAGTGGTCAAAGTCCGCTCCCATATGCCGGAAATCCTGCAATGGCAGCGTGAACGCTTGGTCGCGAAGCTGGAAGa                                                                                           <  2:83501/251‑1 (MQ=255)
                                                                                                                                                    gATATTCTACGCGGGCCGGGCGGGATGGCAGCCCAGGAGCAGGATCTGGACGCCAGTGCCGCTGAAATTCTCGCGGCGCTGGATGGTACGCTGAA‑CGACTTTAGTGGCGCGCGCGAAACCGAAGGTCAGGCACTGAAAGCATTGATCGAGCAGCGTCTGGAAGGCGTCACCGCCGAAGTGGTCAAAGTCCGCTCCCATATGCCGGAAATCCTGCAATGGCAGCGTGAACGCTTGGTAGCGAAGCTGGAAg                                                                                            <  2:173208/250‑1 (MQ=255)
                                                                                                                                                      tatTCTACGCTGGCCGGGCGTGATGGCAGCCCAGGAGCAGGATCTTGACGCCATTGCCGCTGAAATTCTCGCGGCGCTGGATGGTACGCTGAA‑CGACTTTATTGTCGCGCGCGAAACCGAAGGTCAGGCACTGAAAGCATTGATCGAGCAGCGTCTGGAAGGCGTCACCGCCGAAGTGGTCAAAGTCCGCTCCCATATGCCGGAAATCCTGCAATGGCAGCGTGAACGCTTGGTCGCGAAGCTGGAAGATg                                                                                         >  2:68698/1‑251 (MQ=255)
                                                                                                                                                      tatTCTACGCTGGCCGGGCGTGATGGCAGCCCAGGAGCAGGATCTTGACGCCATTGCCGCTGAAATTCTCGCGGCGCTGGATGGTACGCTGAA‑CGACTTTATTGTCGCGCGCGAAACCGAAGGTCAGGCACTGAAAGCATTGATCGAGCAGCGTCTGGAAGGCGTCACCGCCGAAGTGGTAAAAGTCCGCTCCCATATGCCGGAAATCCTGCAATGGCAGCGTGAACGCTTGGTCGCGAAGCTGGACGATg                                                                                         >  2:107278/1‑251 (MQ=255)
                                                                                                                                                            aCGCTGGCCGGGCGTGATGGCAGCCCAGGAGCAGGATCTTGACGCCATTGCCGCTGAAATTCTCGCGGCGCTGGATGGTACGCTGAA‑CGACTTTATTGTCGCGCGCGAAACCGAAGGTCAGGCACTGAAAGCATTGATCGAGCAGCGTCTGGAAGGCGTCACCGCCGAAGTGGTCAAAGTCCGCTCCCATATGCCGGAAATCCTGCAATGGCAGCGTGACCGCTTGGTCGCGAAGCTGGAAGATGCTCa                                                                                     <  2:126746/249‑1 (MQ=255)
                                                                                                                                                            aCGCTGGCCGGGCGTGATGGCAGCCCAGGAGCAGGATCTTGACGCCATTGCCGCTGAAATGCTCGCGGCGCTGGATGGTACGCTGAA‑CGACTTTATTGTCGCGCGCGAAACCGAAGGTCAGGCACTGAAAGCATTGATCGAGCAGCGTCTGGAAGGCGTCACCGCCGAAGTGGTCAAAGTCCGCTCCCATATGCCGGAAATCCTGCAATGGCAGCGTGAACGCTTGGTCGCGAAGCTGGAAGATGCTCa                                                                                     <  2:376341/249‑1 (MQ=255)
                                                                                                                                                                 ggCCGGGCGTGATGGCAGCCCAGGAGCAGGATCTTGACGCCATTGCCGCTGGAATTCTCGCGGCGCTGGATGGTACGCTGAA‑CGACTTTATTGTCGCGCGCGAAACCGAAGGTCAGGCACTGAAAGCATTGATCGAGCAGCGTCTGGAAGGCGTCACCGCCGAAGTGGTCAAAGTCCGCTCCCATATGCCGGAAATCCTGCAATGGCAGCGTGAACGCTTGGTCGCGAAGCTGGAAGATGCTCAGGttcag                                                                              <  1:147541/251‑5 (MQ=255)
                                                                                                                                                                   ccGGGCGTGATGGCAGCCCAGGAGCAGGATCTTGACGCCATTGCCGCTGAAATTCTCGCGGCGCTGGATGGTACGCTGAA‑CGTCTTTATTGTCGCGCGCGAAACCGAAGGTCAGGCACTGAAAGCATTGATCGAGCAGCGTCTGGAAGGCTTCACCGCCGAAGTGGTCAAAGTCCGCTCCCATATGCCGGAAATCCTGCAATGGCAGCGTGAACGCTTGGTCGCGAAGCTGGAAGATGCTCAGGttcagtt                                                                            <  1:29951/251‑7 (MQ=255)
                                                                                                                                                                     gggCGTGATGGCAGCCCAGGAGCAGGATCTTGACGCCATTGCCGCTGAAATTCTCGCGGCGCTGGATGGTACGCTGAA‑CGACTTTATTGTCGCGCGCGAAACCGAAGGTCAGGCACTGAAAGCATTGATCGAGCAGCGTCTGGAAGGCGTCACCGCCGAAGTGGTCAAAGTCCGCTCCCATATGCCGGAAATCCTGCAATGGCAGCGTGAACGCTTGGTCGCGAAGCTGGAAGATGCTCAgggtcagttgg                                                                          >  2:101866/1‑242 (MQ=255)
                                                                                                                                                                      ggCGTGATGGCAGCCCAGGAGCAGGATCTTGACGCCATTGCCGCTGAAATTCTCGCGGCGCTGGATGGTACGCTGAA‑CGACTTTATTGTCGCGCGCGAAACCGAAGGTCAGGCACTGAAAGCATTGATCGAGCAGCGTCTGGAAGGCGTCACCGCCGAAGTGGTCAAAGTCCGCTCCCATATGCCGGAAATCCTGCAATGGCAGCGTGAACGCTTGGTCGCGAAGCTGGAAGATGCTCAGGttcagttgga                                                                         <  2:147233/251‑10 (MQ=255)
                                                                                                                                                                      ggCGTGATGGCAGCCCAGGAGCAGGATCTTGACGCCATTGCCGCTGAAATTCTCGCGGCGCTGGATGGTACGCTGAA‑CGACTTTATTGTCGCGCGCGAAACCGAAGGTCAGGCACTGAAAGCATTGATCGAGCAGCGTCTGGAAGGCGTCACCGCCGAAGTGGTCAAAGTCCGCTCCCATATGCCGGAAATCCTGCAATGGCAGCGTGAACGCTTGGTCGCGAAGCTGGAAGATGCTCAGGttcagttgg                                                                          >  2:21702/1‑242 (MQ=255)
                                                                                                                                                                        cgagATGGCCGCCCAGGAGCAGGATCTTGACGCCATTGCCGCTGAAATTCTCGCGGCGCTGGATGGTGCGCTGAA‑CGACTTTATTGTCGCGCGCGAAACCGAAGGTCAGGCACTGCAAGCATTGATCGAGCAGCGTCTGGAAGGCGTCACCGCCGAAGTGGTCAAAGTCCGCTCCCATATGCCGGAAATCCTGCAATGGCAGCGTGAACGCTTGGTCGCGGAGCTGGACGATGCTCAGGttcagttggaa                                                                        >  2:54922/4‑240 (MQ=255)
                                                                                                                                                                                     cAGGAGCAGGATCTTGACGCCATTGCCGCTGAAATTCTCGCGGCGCTGGATGGTACGCTGAA‑CGACTTTATTGTCGCGCGCGAAACCGAAGGTCAGGCACTGAAAGCATTGATCGAGCAGCGTCTGGAAGGCGTCACCGCCGAAGTGGTCAAAGTCCGCTCCCATATGCCGGAAATCCTGCAATGGCAGCGTGAACGCTTGGTCGCGAAGCTGGAAGATGCTCAGGTTCAGTTGGAAAACAACCGTCTgg                                                           >  1:309861/1‑250 (MQ=255)
                                                                                                                                                                                      aGGAGCAGGATCTTGACGCCATTGCCGCTGAAATTCTCGCGGCGCTGGATGGTACGCTGAA‑CGACTTTATTGTCGCGCGCGAAACCGAAGGTCAGGCACTGAAAGCATTGATCGAGCAGCGTCTGGAAGGCGTCACCGCCGAAGTGGTCAAAGTCCGCTCCCATATGCCGGAAATCCTGCAATGGCAGCGTGAACGCTTGGTCGCGAAGCTGGAAGATGCTCAGGTTCAGTTGGAAAACAACCGTCTgg                                                           >  1:52493/1‑249 (MQ=255)
                                                                                                                                                                                       ggAGCAGGATCTTGACGCCATTGCCGCTGAAATTCTCGCGGCGCTGGATGGTACGCTGAA‑TGACTTTATTGTCGCGCGCGAAACCGAAGGTCAGGCACTGAAAGCATTGATCGAGCAGCGTCTGGAAGGCGTCACCGCCGAAGTGGTCAAAGTCCGCTCCCATATGCCGGAAATCCTGCAATGGCAGCGTGAACGCTTGGTCGCGAAGCTGGAAGATGCTCAGGTTCAGTTGGAAAACAACCGTCTGGAGc                                                        <  1:38142/251‑1 (MQ=255)
                                                                                                                                                                                          gCAGGATCTTGACGCCATTGCCGCTGAAATTCTCGCGGCGCTGGATGGTACGCTGAA‑CGACTTTATTGTCGCGCGCGAAACCGAAGGTCAGGCACTGAAAGCATTGATCGAGCAGCGTCTGGAAGGCGTCACCGCCGAAGTGGTCAAAGTCCGCTCCCATATGCCGGAAATCCTGCAATGGCAGCGTGAACGCTTGGTCGCGAAGCTg                                                                                                >  1:99614/1‑208 (MQ=255)
                                                                                                                                                                                          gCAGGATCTTGACGCCATGGCCGCTGAAATTCTCGCGGCGCTGGATGGTACGCTGAA‑CGACTTTATTGTCGCGCGCGAAACCGAAGGTCAGGCACTGAAAGCATTGATCGAGCAGCGTCTGGAAGGCGTCACCGCCGAAGTGGTCAAAGTCCGCTCCCATATGCCGGAAATCCTGCAATGGCAGCGTGAACGCTTGGTCGCGAAGCTg                                                                                                <  2:99614/208‑1 (MQ=255)
                                                                                                                                                                                             ggATCTTGACGCCATTGCCGCTGAAATTCTCGCGGCGCTGGATGGTACGCTGAA‑CGACTTTATTGTCGCGCGCGAAACCGAAGGTCAGGCACTGAAAGCATTGATCGAGCAGCGTCGGGAAGGCGTCACCGCCGAAGTGGTCAAAGTCCGCTCCCATATGCCGGAAATCCTGCAATGGCAGCGTGAACGCTTGGTCGCGAAGCTGGAAGATGCTCAGGTTCAGTTGGAAAACAACCGTCTGGAGCAGGagt                                                  <  2:108216/251‑3 (MQ=255)
                                                                                                                                                                                                   ggACGCCAGTGCCGCTGAAATGCTCGCGGCGCTGGATGGTACGCTGAA‑CGACTTTATTGTCGCGCGCGAAACCGAAGGTCAGGCACTGAAAGCATTGATCGAGCAGCGTCTGGAAGGCGTCACCGCCGAAGTGGTCAAAGTCCGCTCCCATATGCCGGAAATCCTGCAAGGGCAGCGTGAACGCTTGATCGCGAAGCTGGAAGATGCTCAGGTTCAGTTGGAAAACAACCGTCTGGAGCAGGagttggtt                                             <  2:125693/249‑8 (MQ=255)
                                                                                                                                                                                                             gCCGCTGAAATGCTCGCGGCGCTGGATGGTACGC‑GGAGCTCCTTTATTGTCGCGCGCGAAACCGAAGGTCAGGCACTGAAAGCAGTGTTCGCGCAGCGTCTGGAAGGCGTCACCGCCGAAGTGGTCAAAGTCCGCGCCCATATGCCGGAAATCCTGCAATGGCAGCGTGAACGCTTGGTCGCGAAGCTGGAAGATGCTCAGGTTCAGTTGGAAAACAACCGTCTGGAGCAGGAGTTGGTTCTGCTGGcaca                                  <  2:99004/251‑1 (MQ=255)
                                                                                                                                                                                                              ccGCTGAAATTCTCGCGGCGCTGGATGGTACGCTGAA‑CGACTTTATTGTCGCGCGCGAAACCGAAGGTCAGGCACTGAAAGCATTGATCGAGCAGCGTCTGGAAGGCGTCACCGCCGAAGTGGTCAAAGTCCGCTCCCATATGCCGGAAATCCTGCAATGGCAGCGTGAACGCTTGGTCGCGAAGCTGGAAGATGCTCAGGTTCAGTTGGAAAACAACCGTCTGGAGCAGGagtt                                                 <  2:368359/235‑4 (MQ=255)
                                                                                                                                                                                                              ccGCTGAAATTCTCGCGGCGCTGGATGGTACGCTGAA‑CGACTTTATTGTCGCGCGCGAAACCGAAGGTCAGGCACTGAAAGCATTGATCGAGCAGCGTCTGGAAGGCGTCACCGCCGAAGTGGTCAAAGTCCGCTCCCATATGCCGGAAATCCTGCAATGGCAGCGTGAACGCTTGGTCGCGAAGCTGGAAGATGCTCAGGTTCAGTTGGAAAACAACCGTCTGGAGCAGGagtt                                                 >  1:368359/1‑232 (MQ=255)
                                                                                                                                                                                                                   gAAATTCTCGCGGCGCTGGATGGTACGCTGAA‑CGACTTTATTGTCGCGCGCGAAACCGAAGGTCAGGCACTGAAAGCATTGATCGAGCCGCGTCTGGAAGGCGTAACCGCCGAAGTGGTCAAAGTCAGCTCCCATATGCCGGAAATCCTGCAATGGCAGCGTGAAAGCTTGGTCGCGAAGCTGGAAGATGCTCAGGTTCAGTTGGAAAACAAACGTCTGGAGCAGGAGTTGGTTCTGCTGGCACAACGAAt                            >  1:194897/1‑251 (MQ=255)
                                                                                                                                                                                                                    aacgtctcGCTGGGCTGGATGGTACGCTGAA‑CGACTTTATGGTCGCGCGCGAAACCGAAGGTCAGGCACTGAAAGCATTTATCGAGCAGCGTCTGGAAGGCGTCACCGCCGAAGTGGTCAAAGTCCGCTCCCATATGCCGGAAATCCTGCAATGGCAGCGTGAACGCTTGGTCGCGAAGCTGGAAGATGCTCAGGTTCAGTTGGAAAACAACCGTCTGGAGCAGGAGTTGGTTCTGCTGGCACAACGAAtt                           <  2:141893/247‑1 (MQ=255)
                                                                                                                                                                                                                    aaaTTCTCGCGGCGCTGGATGGTACGCTGAA‑CGACTTTATTGTCGCGCGCGAAACCGAAGGTCAGGCACTGAACGCATTGATCGAGCAGCGTCTGGAAGGCGTCACCGCCGAAGTGGTCAAAGTCCGCTCCCATATGCCGGAAATCCTGCAATGGCAGCGTGAACGCTTGGTCGCGAAGCTGGAAGATGCTCAGGTTCAGTTGGAAACCAACCGTCTGGAGCAGGAGTTGGTTCTGCTGGCACAACGAAtt                           >  1:389470/1‑251 (MQ=255)
                                                                                                                                                                                                                           cgcgGCGCTGGATGGTACGCTGAA‑CGACTTTATTGTCGCGCGCGAAACCGAAGGTCAGGCACTGAAAGCATTGATCGAGCAGCGTCTGGAAGGCGTCACCGCCGAAGTGGTCAAAGTCCGCTCCCATATGCCGGAAATCCTGCAATGGCAGCGTGAACGCTTGGTCGCGAAGCTGGAAGATGCTCAGGTTCAGTTGGAAAACAACCGTCTGGAGCAGGAGTTGGTTCTGCTGGCACAACGAATTGACGt                      <  1:257778/249‑1 (MQ=255)
                                                                                                                                                                                                                             cggcgCTGGATGGGACGCTGAA‑CGACTTTATTGTCGCGCGCGAAACCGAAGGTCAGGCACTGAAAGCATTGATCGAGCAGCGTCTGGAAGGCGTCACCGCCGAAGTGGTCAAAGTCCGCTCCCATATGCCGGAAATCCTGCAATGGCAGCGTGAACGCTTGGTCGCGAAGCTGGAAGATGCTCAGGTTCAGTTGGAAAACAACCGTCTGGAGCAGTATTTGGTTCTGCTGGCACAACGAATTGACGTTGcc                  <  1:119670/251‑1 (MQ=255)
                                                                                                                                                                                                                              ggcgCTGGATGGTACGCTGAA‑CGACTTTAGTGTCGCGCGCGAAACCGAAGGCCAGGCACTGAAAGCATTGATCGAGCCGCGTCTGGAAGGCGTCACCGCCGAAGTGGTCAAAGTCCGCTCCCATATGCCGGAAATCCTGCAATGGCAGCGTGAACGCTTGGTCGCGAAGCTGGAAGATGCTCAGGTTCAGTTGGAAAACAACCGTCTGGAGCAGGAGTTGGTTCTGCTGGCACAACGAATTGACGTTGCCg                 <  2:358158/251‑1 (MQ=255)
                                                                                                                                                                                                                                      aTGGTACGCTGAA‑CGACTTTATTGTCGCGCGCGAAACCGAAGGTCTGACACTGAAAGCATTGATCGAGCAGAGTCTGGAAGGCGTCACCGCCGACGTGGTCAAAGTCCGCTCCCCTATGCCGGCACTACTGCAATGGAAGCGTGAACGCTTGGTCGCGAAGCTGGAAGATGCTCAGGTTCAGTTGGAAAACAACCGTCTGGAGCAGGAGTTGGTTCTGCTGGCACAACGAATTGACGTTGCAGAAGAACTg         >  1:293051/1‑251 (MQ=255)
                                                                                                                                                                                                                                       ggttACGCTGAA‑CGACTTTATTGTCGCGCGCGAAACCGAAGGTCAGGCACTGAAAGCATTGATCGAGCAGCGTCTGGAAGGCGTCACCGCCGAAGTGGTCAAAGTCCGCTCCCATATGCCGGAAATCCTGCAATGGCAGCGTGAACGCTTGGTCGCGAAGCTGGAAGATGCTCAGGTTCAGTTGGAAAACAACCGTCTGGAGCAGGAGTTGGTTCTGCTGGCACAACGAATTGACGTTGCCGAAGAACTgg        <  1:412192/248‑1 (MQ=255)
                                                                                                                                                                                                                                             gCTGAA‑CGACTGTATTGTCGCGCGCGAAACCGAAGGCCAGGCACTGAAAGCATTGATCGAGCAGCGGCTGGAAGGCGTCACCGCCGAAGTGGTCAAAGTCCGCTCCCATATGCCGGAAATCCTGCAATGGCAGCGTGAACGCTTGGGCGCGAAGCTGGAAGATGCTCAGGTTCAGTTGGAAAACAACCGTCTGGAGCAGGAGTTGGTTCTGCTGGCACAACGAATTGACGTTGCCGAAGAACTGGACCGcc  <  2:361613/251‑1 (MQ=255)
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                       
GCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGCCAGATGTTAGCGCGCAAGGTGAGCTGATCCTCAACGAAAAACTGGCTAAACAGCTGGTAACTGCCGCGAACTGGGTAAAAATGCAGAGTGACGAAGGGGAAATCAACCCGGTTGATATTCTACGCTGGCCGGGCGTGATGGCAGCCCAGGAGCAGGATCTTGACGCCATTGCCGCTGAAATTCTCGCGGCGCTGGATGGTACGCTGGA‑CGACTTTATTGTCGCGCGCGAAACCGAAGGTCAGGCACTGAAAGCATTGATCGAGCAGCGTCTGGAAGGCGTCACCGCCGAAGTGGTCAAAGTCCGCTCCCATATGCCGGAAATCCTGCAATGGCAGCGTGAGCGTCTGGTCGCGAAGCTGGAAGATGCTCAGGTGCAACTGGAAAACAACCGTCTGGAGCAGGAACTGGTTCTGCTGGCACAACGAATTGACGTTGCCGAAGAACTGGATCGCC  >  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: ‑