Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,213,179 Δ1 bp 5.9% intergenic (+20/‑55) rrfE → / → yjaA 5S ribosomal RNA of rrnE operon/stress‑induced protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,213,1790T. 5.9% 659.7 / 31.0 187intergenic (+20/‑55)rrfE/yjaA5S ribosomal RNA of rrnE operon/stress‑induced protein
Reads supporting (aligned to +/- strand):  ref base T (91/85);  new base . (8/3);  total (99/88)
Fisher's exact test for biased strand distribution p-value = 2.22e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.47e-01

Only 100 of 189 total aligned reads displayed.

ACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAA‑GGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGCAAACGTGAAGTGTCAGGAGACGCTGCCTTCAGT  >  NC_000913/4213032‑4213326
                                                                                                                                                   |                                                                                                                                                     
aCAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAACAACCCAGGTc                                                                                                                                                      >  1:681078/1‑149 (MQ=255)
aCAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTc                                                                                                                                                      >  1:1062481/1‑149 (MQ=255)
aCAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTc                                                                                                                                                      >  4:1134968/1‑149 (MQ=255)
     aTTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATaa                                                                                                                                                 <  2:1628615/149‑1 (MQ=255)
     aTTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGCGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAAGCATCAAATTAGAAAAACCCCGGTCCATaa                                                                                                                                                 >  2:393094/1‑149 (MQ=255)
      tttGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAg                                                                                                                                                >  3:1162957/1‑149 (MQ=255)
      tttGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAg                                                                                                                                                >  3:467741/1‑149 (MQ=255)
       ttGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAgg                                                                                                                                               <  4:1162957/149‑1 (MQ=255)
           cTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCgg                                                                                                                                          >  1:1109866/1‑149 (MQ=255)
            tGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCggg                                                                                                                                         >  1:1189888/1‑149 (MQ=255)
            tGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCggg                                                                                                                                         <  2:336125/149‑1 (MQ=255)
                  cAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGttttt                                                                                                                                   <  3:427100/149‑1 (MQ=255)
                   aGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGtttttt                                                                                                                                  <  2:1189888/149‑1 (MQ=255)
                   aGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGtttttt                                                                                                                                  <  4:1391600/149‑1 (MQ=255)
                         gcgGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCatat                                                                                                                            >  1:1722473/1‑149 (MQ=255)
                          cgGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATc                                                                                                                           <  3:870340/149‑1 (MQ=255)
                            gtggtCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCaa                                                                                                                         <  4:467741/149‑1 (MQ=255)
                                  ccACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTAttt                                                                                                                   >  4:371787/1‑149 (MQ=255)
                                  ccACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTAttt                                                                                                                   >  2:591492/1‑149 (MQ=255)
                                  ccACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTAttt                                                                                                                   <  3:1102822/149‑1 (MQ=255)
                                      cTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCat                                                                                                               >  2:1145446/1‑149 (MQ=255)
                                       tGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGc                                                                                                                                <  1:1590942/131‑1 (MQ=255)
                                        gACCCAATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGAAGGTAGCGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAACAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCatga                                                                                                             >  4:1588050/1‑149 (MQ=255)
                                               tGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAA‑ggg                                                                                                     >  3:1404424/1‑149 (MQ=255)
                                                gCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAA‑GGGa                                                                                                    <  3:1266941/149‑1 (MQ=255)
                                                  cGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAA‑GGGAAt                                                                                                  >  4:927844/1‑149 (MQ=255)
                                                      cTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAA‑GGGAATCTCa                                                                                              <  4:301848/149‑1 (MQ=255)
                                                       tCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAA‑GGGAATCTCAt                                                                                             <  1:1627901/149‑1 (MQ=255)
                                                        cAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCg                                                                                                                                           <  1:1385991/103‑1 (MQ=255)
                                                         aGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGATCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCat                                                                                                               <  1:1271563/130‑1 (MQ=255)
                                                           aaGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAA‑GGGAATCTCATGTCa                                                                                         <  2:1588619/149‑1 (MQ=255)
                                                              tGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAA‑GGGAATCTCATGTCAGtt                                                                                      <  3:927844/149‑1 (MQ=255)
                                                               gAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAA‑GGGAATCTCATGTCAGTTc                                                                                     <  2:322001/149‑1 (MQ=255)
                                                                aaaCGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAA‑GGGAATCTCATg                                                                                            <  3:250647/141‑1 (MQ=255)
                                                                aaaCGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAA‑GGGAATCTCATGTCAGTTCt                                                                                    <  4:1081305/149‑1 (MQ=255)
                                                                aaaCGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAA‑GGGAATCTCATGTCAGTTCt                                                                                    <  4:1404424/149‑1 (MQ=255)
                                                                 aaCGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAA‑GGGAATCTCATGTCAGTTCTg                                                                                   >  1:942100/1‑149 (MQ=255)
                                                                 aaCGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAA‑GGGAATCTCATGTCAGTTCTg                                                                                   >  4:1081884/1‑149 (MQ=255)
                                                                      cGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAA‑GGGAATCTCATGTCAGTTCTGTatat                                                                              <  3:1317974/149‑1 (MQ=255)
                                                                          gcgcCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAA‑GGGAATCTCATGTCAGTTCTGTATATCCaa                                                                          >  1:1562146/1‑149 (MQ=255)
                                                                            gcCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAA‑GGGAATCTCATGTCAGTTCTGTATATCCAAAt                                                                        >  1:1467434/1‑149 (MQ=255)
                                                                             ccGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAA‑GGGAATCTCATGTCAGTTCTGTATATCCAAAtt                                                                       >  1:409770/1‑149 (MQ=255)
                                                                                aTGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAA‑GGGAATCTCATGTCAGTTCTGTATATCCAAATtcgt                                                                    <  2:409770/149‑1 (MQ=255)
                                                                                 tGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAA‑GGGAATCTCATGTCAGTTCTGTATATCCAAATtcgtc                                                                   <  3:371787/149‑1 (MQ=255)
                                                                                 tGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAA‑GGGAATCTCATGTCAGTTCTGTATATCCAAATtcgtc                                                                   <  4:458043/149‑1 (MQ=255)
                                                                                  ggTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAA‑GGGAATCTCATGTCAGTTCTGTATATCCAAATtcgtcg                                                                  >  4:599915/1‑149 (MQ=255)
                                                                                         tgGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAA‑GGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCaa                                                           >  3:633511/1‑149 (MQ=255)
                                                                                           gggTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTg                                                                                                                                 <  2:801642/78‑1 (MQ=255)
                                                                                           gggTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTg                                                                                                                                 >  1:801642/1‑78 (MQ=255)
                                                                                                  cccATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAA‑GGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGtt                                                  >  1:384536/1‑149 (MQ=255)
                                                                                                  cccATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGG‑CCTTATGGACCGGGGtttttc                                                                                                                                  >  2:694925/1‑69 (MQ=17)
                                                                                                   ccATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAA‑GGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTc                                                 >  2:125282/1‑149 (MQ=255)
                                                                                                     aTGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAA‑GGGAATCTCATGTCAGTCCTGTATATCCAATTTCGTCGTAATCAAAGTACTTTTcgc                                               <  2:947400/149‑1 (MQ=255)
                                                                                                      tGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAA‑g                                                                                                       >  2:302858/1‑92 (MQ=255)
                                                                                                              tAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAA‑GGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGaa                                      <  2:1314109/149‑1 (MQ=255)
                                                                                                               aGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAA‑GGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGaaa                                     <  1:1021850/149‑1 (MQ=255)
                                                                                                                gggAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAA‑GGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGaaa                                     >  3:1685853/1‑148 (MQ=255)
                                                                                                                gggAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAA‑GGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAg                                    >  3:771192/1‑149 (MQ=255)
                                                                                                                 ggAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAA‑GGGAAtc                                                                                                 >  3:990894/1‑87 (MQ=255)
                                                                                                                 ggAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAA‑GGGAAtc                                                                                                 <  4:990894/87‑1 (MQ=255)
                                                                                                                    aCTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAA‑GGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGCaaa                                >  2:1067030/1‑149 (MQ=255)
                                                                                                                     cTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAA‑GGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGCAAAc                               <  3:1464668/149‑1 (MQ=255)
                                                                                                                         cAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAC‑GGGAATCGCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGGAAACGTGa                           >  1:1336924/1‑149 (MQ=255)
                                                                                                                         cAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAA‑GGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGCAAACGTGa                           >  4:487703/1‑149 (MQ=255)
                                                                                                                         cAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAA‑GGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGCAAACGTGa                           >  3:691280/1‑149 (MQ=255)
                                                                                                                          aGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAA‑GGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGCAAACGTGaa                          >  1:761944/1‑149 (MQ=255)
                                                                                                                          aGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAA‑GGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGCAAACGTGaa                          <  2:1467434/149‑1 (MQ=255)
                                                                                                                          aGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAA‑GGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGCAAACGTGaa                          >  3:908494/1‑149 (MQ=255)
                                                                                                                           ggCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAA‑GGGAATCTCAt                                                                                             <  3:739918/81‑1 (MQ=255)
                                                                                                                           ggCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAA‑GGGAATCTCAt                                                                                             >  4:739918/1‑81 (MQ=255)
                                                                                                                            gCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAA‑GGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGCAAACGTGAAgt                        <  1:355997/149‑1 (MQ=255)
                                                                                                                              aTCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAA‑GGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGCAAACGTGAAgtgt                      >  1:1702701/1‑149 (MQ=255)
                                                                                                                               tCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAA‑GGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGCAAACGTGAAGTGTc                     >  3:346894/1‑149 (MQ=255)
                                                                                                                                 aaaTTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAA‑GGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACt                                                     <  1:1241047/115‑1 (MQ=255)
                                                                                                                                 aaaTTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAA‑GGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGCAAACGTGAAGTGTCAg                   <  2:761944/149‑1 (MQ=255)
                                                                                                                                 aaaTTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAA‑GGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGCAAACGTGAAGTGTCAg                   <  3:484851/149‑1 (MQ=255)
                                                                                                                                   aTTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAA‑GGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGCAAACGTGAAGTGTCAGga                 <  4:346894/149‑1 (MQ=255)
                                                                                                                                      aGAAAAACCCCGGTCCATTAAG‑CCGGGGTTTTTTGCAAATCAATTATTTTCAATATGAAGGGGAATGT‑TTGTCTGTTTTTTATATCCAAAATCCTCGTAATCAAATTACTGTTCGCGATCTTGAAAGCAAACGTGAAGTGTCAGGAGAc              <  2:1062481/149‑1 (MQ=255)
                                                                                                                                      aGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAA‑GGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGCAAACGTGAAGTGTCAGGAGAc              <  1:431896/149‑1 (MQ=255)
                                                                                                                                       gAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAA‑GGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGCAAACGTGAAGTGTCAGGAGACg             >  1:481876/1‑149 (MQ=255)
                                                                                                                                       aaaaaaCCCCGG‑CCTTATGGACCGGGGTTTTTTGCATATCAATTATTTGCATGATGAA‑GGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGCAAACGTGAAGTGTCAGGAGACg             >  2:421373/2‑149 (MQ=255)
                                                                                                                                       aaaaaaCCCCGG‑CCTTATGGACCGGGGTTTTTTGCATATCAATTATTTGCATGATGAA‑GGGAATCTCATGTCAGTTCTCTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGCAAACGTGAAGTGTCAGGAGACg             >  4:1213116/2‑149 (MQ=255)
                                                                                                                                        aaaaaCCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAA‑GGGAATCTCATGTCAGTTCTGTATATCCAGATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGCAAACGTGAAGTGTCAGGAGACGc            >  3:1222688/1‑149 (MQ=255)
                                                                                                                                        aaaaaCCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAA‑GGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGCAAACGTGAAGTGTCAGGAGACGc            >  3:651480/1‑149 (MQ=255)
                                                                                                                                        aaaaaCCCCGG‑CCTTATGGACCGGGGTTTTTTGCATATCAATTATTTGCATGATGAA‑GGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGCAAACGTGAAGTGTCa                    >  4:75787/1‑141 (MQ=255)
                                                                                                                                        aaaaaCCCCGG‑CCTTATGGACCGGGGTTTTTTGCATATCAATTATTTGCATGATGAA‑GGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGCAAACGTGAAGTGTCAGGAGACGc            >  3:1007446/1‑149 (MQ=255)
                                                                                                                                           aaCCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAA‑GGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGt                                                   >  3:877826/1‑107 (MQ=255)
                                                                                                                                            aCCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAA‑GGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGCAAACGTGAAGTGTCAGGAGACGCTGcc        >  1:304278/1‑149 (MQ=255)
                                                                                                                                            aCCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAA‑GGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGCAAACGTGAAGTGTCAGGAGACGCTGcc        <  1:44490/149‑1 (MQ=255)
                                                                                                                                              cccGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAA‑GGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTc                                                 >  2:670706/1‑106 (MQ=255)
                                                                                                                                              cccGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAA‑GGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTc                                                 <  1:670706/106‑1 (MQ=255)
                                                                                                                                              cccGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAA‑GGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTTAAAGCAAACGTGAAGTGTCAGGAGACGCTGCCtt      >  4:22840/1‑149 (MQ=255)
                                                                                                                                              cccGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAA‑GGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGCAAACGTGaa                          >  1:835607/1‑129 (MQ=255)
                                                                                                                                              cccGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAA‑GGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGCAAACGTGaa                          <  2:835607/129‑1 (MQ=255)
                                                                                                                                               ccGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAA‑GGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGCAAACGTGAAGTGTCAGGAGACGCTGCCTTc     >  2:500804/1‑149 (MQ=255)
                                                                                                                                                cGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAA‑GGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGCAAACGTGAAGTGTCAGGAGACGCTGCCTTCa    >  4:263408/1‑149 (MQ=255)
                                                                                                                                                 ggTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAA‑GGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGCAAACGTGAAGTGTCAGGAGACGCTGCCTTCAg   <  3:263408/149‑1 (MQ=255)
                                                                                                                                                  gTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAA‑GGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGCAAACGTGAAGTGTCAGGAGACGCTGCCTTCAGt  >  4:424714/1‑149 (MQ=255)
                                                                                                                                                   tCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAA‑GGGAATCTCATGTCAGTTCTGTATATCCAAATtcgt                                                                    <  4:1289783/82‑1 (MQ=255)
                                                                                                                                                   tCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAA‑GGGAATCTCATGTCAGTTCTGTATATCCAAATtcgt                                                                    >  3:1289783/1‑82 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
ACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAA‑GGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGCAAACGTGAAGTGTCAGGAGACGCTGCCTTCAGT  >  NC_000913/4213032‑4213326

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