Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,148,049 C→G 5.1% G90A (GGC→GCC)  yijO ← AraC family putative transcriptional activator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,148,0490CG 5.1% 326.3 / 6.7 113G90A (GGC→GCC) yijOAraC family putative transcriptional activator
Reads supporting (aligned to +/- strand):  ref base C (43/63);  new base G (3/3);  total (47/66)
Fisher's exact test for biased strand distribution p-value = 6.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 116 total aligned reads displayed.

TGCTCCTGCGGCTGCATTTGCATTTCATTGAGCGTTTGTAGCAGAAAAGAACCTATGCGTGGGCCGCGTCGGGCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGCACATCGCCGGGAACTAACGTTGCTCCAGCGCCGGTATCAACAAACTCACCTTCCAGCACAATTTCCAGCCGTG  >  NC_000913/4147902‑4148196
                                                                                                                                                   |                                                                                                                                                   
tgctCCTGCGGCTGCATTTGCATTTCATTGAGCGTTTGTAGCAGAAAAGAACCTATGCGTGGGCCGCGTCGGGCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTTCTGTTTGca                                                                                                                                                    >  2:1510650/1‑148 (MQ=255)
tgctCCTGCGGCTGCATTTGCATTTCATTGAGCGTTTGTAGCAGAAAAGAACCTATGCGTGGGCCGCGTCGGGCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGATGTTTGcc                                                                                                                                                    >  1:1023236/1‑149 (MQ=255)
  ctCCTGCGGCTGCATTTGCATTTCTTTGAGCGTTTGTAGCAGAAAAGATCCTATGAGTGGGCCGCGTCGGGCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCaa                                                                                                                                                  >  2:1039225/1‑149 (MQ=255)
   tCCTGCGGCTGCATTTGCATTTCATTGAGCGTTTGTAGCAGAAAAGAACCTATGCGTGGGCCGCGTCGGGCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGGCaaa                                                                                                                                                 <  1:814915/149‑1 (MQ=255)
    ccTGCGGCCGCATTTGGATTTCCTTGTGTGTTTGTAGCTAAAAATAACCTATGCGGTGGGAGCGTTTGGCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAAc                                                                                                                                                <  2:1388866/149‑1 (MQ=255)
     cTGCGGCTGCATTTGCATTTCATTGAGCGTTTGTAGCAGAAAAGAACCTTTGCGTGGGCCGCGTCTGGCGACGTGTTGCTTCGCAAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACa                                                                                                                                               <  1:1442921/149‑1 (MQ=255)
     cTGCGGCTGCATTTGCATTTCATTGAGCGTTTGTAGCAGAAAAGAACCTATGCGTGGGCCGCGTCGGGCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACa                                                                                                                                               >  2:369243/1‑149 (MQ=255)
        cGGCTGCATTTGCATTTCATTGAGCGTTTGTAGCAGAAAAGAACCTATGCGTGGGCCGCGTCGGGCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCa                                                                                                                                            <  2:969272/149‑1 (MQ=255)
        cGGCTGCATTTGCATTTCATTGAGCGTTTGTAGCAGAAAAGAACCTATGCGTGGGCCGCGTCGGGCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCa                                                                                                                                            <  2:389304/149‑1 (MQ=255)
         ggCTGCATTTGCATTTCATTGAGCGTTTGTAGCAGAAAAGAACCTATGCGTGGGCCGCGTCGGGCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCAc                                                                                                                                           <  1:1527474/149‑1 (MQ=255)
          gCTGCATTTTGATTTCATTGAGCGTTTGTAGCAGAAAAGAACCTATGCGTGGGCCGCGTCGGGCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGGCAAACAGCaac                                                                                                                                          <  2:148275/149‑3 (MQ=255)
          gCTGCATTTGCATTTCATTGAGCGTTTGTAGCAGAAAAGAACCTATGCGTGGGCCGCGTCGGGCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACg                                                                                                                                          <  1:544852/149‑1 (MQ=255)
            tgcatttgcattTCATTGAGCGTTTGTAGCAGAAAAGAACCTATGCGTGGGCCGCGTCGGGCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGCCGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCt                                                                                                                                        >  1:811721/1‑149 (MQ=255)
               atttgcatttCATTGAGCGTTTGTAGCAGAAAAGAACCTATGCGTGGGCCGCGTCGGGCGACGTGTTGCTTCGCCAGATTTTGATATTGTATGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTaaa                                                                                                                                     <  1:944047/149‑1 (MQ=255)
               atttgcatttCATTGAGCGTATGTAGCAGAAAAGAACCTATGCGTGGGCCGCGTCGGGCGACGTGTTGCTTCGCCAGTTTTTGATATTGTTTGCAATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGGCAAACAGCACGCTaaa                                                                                                                                     >  1:6233/1‑149 (MQ=255)
                tttgcatttCATTGAGCGTTTGTAGCAGAAAAGAACCTATGCGTGGGCCGCGTCGGGCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAg                                                                                                                                    <  1:708972/149‑1 (MQ=255)
                 ttgcatttCATTGAGCGTTTGTAGCAGAAAAGAACCTATGCGTGGGCCGCGTCGGGCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAgg                                                                                                                                   >  1:463490/1‑149 (MQ=255)
                  tgcatttCATTGAGCGTTTGTAGCAGAAAAGAACCTATGCGTGGGCCGCGTCGGGCGCCGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGgt                                                                                                                                  >  2:77061/1‑149 (MQ=255)
                      tttCATTGAGCGTTTGTAGCAGAAAAGAACCTATGCGTGGGCCGCGTCGGGCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGgtagta                                                                                                                              >  1:1347581/1‑149 (MQ=255)
                      tttCATTGAGCGTTTGTAGCAGAAAAGAACCTATGCGTGGGCCGCGTCGGGCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGgtagta                                                                                                                              >  1:1371464/1‑149 (MQ=255)
                       ttCATTGAGCGTATGTAGCAGAAAAGAACCAATGCGTGGGCCGCGTCGGGCGACGTTTTGCTTCGCCAGATTTTGATATAGTTTGCCATCCAATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCAAAAGGtagtag                                                                                                                             <  2:895118/149‑1 (MQ=255)
                          aTTGAGCGTTTGTAGCAGAAAAGAACCTATGCGTGGGCCGCGTCGGGCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCgg                                                                                                                          <  1:410393/149‑1 (MQ=255)
                          aTTGAGCGTTTGTAGCAGAAAAGAACCTATGCGTGGGCCGCGTCGGGCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCgg                                                                                                                          >  2:1178886/1‑149 (MQ=255)
                           ttGAGCGTTTGTAGCAGAAAAGAACCTATGCGTGGGCCGCGTCGGGCGACGTGTTGCTTCGCCAGTTTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCggg                                                                                                                         <  1:12129/149‑1 (MQ=255)
                           ttGAGCGTTTGTAGCAGAAAAGAACCTATGCGTGGGCCGCGTCGGGCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCggg                                                                                                                         >  2:370368/1‑149 (MQ=255)
                           ttGAGCGTGTGTAGCACAAAAGAACCTATGCCGGGGCCTCGTCGGGCGACGGGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCggg                                                                                                                         <  2:1023236/149‑1 (MQ=255)
                            tGAGCGTTTGTAGCAGAAAAGAACCTATGCGTGGGCCGCGTCGGGCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCgggg                                                                                                                        <  1:1538964/149‑1 (MQ=255)
                             gAGCGTTTGTAGCAGAAAAGAACCTATGCGTGGGCCGCGTCGGGCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCggggg                                                                                                                       <  1:652037/149‑1 (MQ=255)
                                cGTTTGTAGCAGAAAAGAACCTATGCGTGGGCCGCGTCGGGCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCtt                                                                                                                    <  1:162838/149‑1 (MQ=255)
                                cGTTTGTAGCAGAAAAGAACCTATGCGTGGGCCGCGTCGGGCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCtt                                                                                                                    <  1:989677/149‑1 (MQ=255)
                                  tttGTAGCAGAAAAGAACCTATGCGTGGGCCGCGTCGGGCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGc                                                                                                                  <  1:35365/149‑1 (MQ=255)
                                   ttGTAGCAGAAAAGAACCTATGCGTGGGCCGCGTCGGGCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGcc                                                                                                                 <  1:278998/149‑1 (MQ=255)
                                     gTAGCAGAAAATAACCTATGCGTGGGCCGCGTCGGGCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGCCGCTGAACCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCAt                                                                                                               >  2:1394139/1‑149 (MQ=255)
                                      tggCAGAAAAGAACCTATGCGTGGGCCGCGTCGGGCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCAtt                                                                                                              <  1:541865/147‑1 (MQ=255)
                                      tcgCAGAAACGAACCGATGCGTGGGCCGCGTCGGGCGACGGGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCAtt                                                                                                              <  1:281014/147‑1 (MQ=255)
                                      tAGCAGAAAAGAACCTATGCGTGGGCCGCGTCGGGCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCAtt                                                                                                              >  2:1533750/1‑149 (MQ=255)
                                       aGCAGAAAAGAACCTATGCGTGGGCCGCGTCGGGCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGGAGTAGCGGGGGCTTGCCATtg                                                                                                             >  2:72986/1‑149 (MQ=255)
                                        gCAGAAAAGAACCTATGCGTGGGCCGCGTCGGGCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATtgt                                                                                                            <  1:345401/149‑1 (MQ=255)
                                        gCAGAAAAGAACCTATGCGTGGGCCGCGTCGGGCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCGGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATtgt                                                                                                            <  2:1529733/149‑1 (MQ=255)
                                         cAGAAAAGAACCTATGCGTGGGCCGCGTCGGGCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATtgtg                                                                                                           <  1:598833/149‑1 (MQ=255)
                                           gAAAAGAACCTATGCGTGGGCCGCGTCGGGCGACGTGTTGCTTCGCCAGATTTTGATATTGTTAGCCATCCCATTGTACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGa                                                                                                         >  2:1179815/1‑149 (MQ=255)
                                              aaGACCCTATGCGGGGGCCGCGTCGGGCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGaaaa                                                                                                      <  1:153182/149‑1 (MQ=255)
                                              aaGAACCTATGCGTGGGCCGCGTCGGGCGACGTGTTGCTACGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGaaaa                                                                                                      >  2:3800/1‑149 (MQ=255)
                                                     tATGCGTGGGCCGCGTCGGGCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCGATTGGGGAAAATTCCAGc                                                                                               >  1:229802/1‑149 (MQ=255)
                                                       tGCGTGGGCCGCGTCGGGCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCAGTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACACCCCGCTAAAGGTAGTAGCGGGGGCGTGCCATTGTGGAAAATTCCAGcgc                                                                                             >  2:1445320/1‑147 (MQ=255)
                                                        gCGTGGGCCGCGTCGGGCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGccac                                                                                            <  1:988193/149‑1 (MQ=255)
                                                          gTGGGCCGCGTCGGGCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGccacca                                                                                          <  1:28375/149‑1 (MQ=255)
                                                           tGGGCCGCGTCGGGCGACGTGGTGCGTCGCCAGATTTTGATATTGTTTGCCATCCCAATGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACcag                                                                                         <  1:150718/149‑1 (MQ=255)
                                                            gggCCGCGTCGGGCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACcagc                                                                                        <  1:369243/149‑1 (MQ=255)
                                                             ggCCGCGTCGGGCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACcagca                                                                                       <  2:666320/149‑1 (MQ=255)
                                                             ggCCGCGTCGGGCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACcagca                                                                                       <  1:370368/149‑1 (MQ=255)
                                                             ggCCGCGTCGGGCGACGGGTAGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACcagca                                                                                       <  1:883937/149‑1 (MQ=255)
                                                                      gggCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTAc                                                                              <  1:996082/149‑1 (MQ=255)
                                                                      gggCGACGTGTTGCTTCGCCAGATTTTGATATTGTTAGCCATCCCATTGCACGACGCAGAACCCGAGTTGCTGTTTGCCAAACAGCCCGCTAAAGGGAGTAGCGGGGGCTTGCCATTCTGGAAAATTCCACCCACCAGCAGGAACGTAc                                                                              >  1:1262709/1‑149 (MQ=255)
                                                                         cGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAATTTCCAGCCACCAGCAGGCACGTACAGc                                                                           >  2:196842/1‑149 (MQ=255)
                                                                         cGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGc                                                                           <  2:521494/149‑1 (MQ=255)
                                                                          gACGTGTTGCTTCGCCATATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGca                                                                          <  1:391370/149‑1 (MQ=255)
                                                                            cGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGcaca                                                                        <  1:1510650/149‑1 (MQ=255)
                                                                             gtgtTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGCACAt                                                                       >  1:253628/1‑149 (MQ=255)
                                                                               gtTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGAAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCCGGAACGTCCCGCACATCg                                                                     >  2:753150/1‑149 (MQ=255)
                                                                               gtTGCTTCGCCAGATTTTGATATTGTTTGCCATCACATTGCACGACGCTGAAGACGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGCACATCg                                                                     >  2:159054/1‑149 (MQ=255)
                                                                                ttGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGCACATCGc                                                                    <  2:776270/149‑1 (MQ=255)
                                                                                ttGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGCACATCGc                                                                    >  1:554836/1‑149 (MQ=255)
                                                                                  gCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGCACATCGCCg                                                                  <  2:440132/149‑1 (MQ=255)
                                                                                  gCTGCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTTGAAAATTCCAGCCACCAGCAGGAACGTACAGCACATCGCCg                                                                  <  1:512453/149‑1 (MQ=255)
                                                                                   cTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGCACATCGCCgg                                                                 <  2:1311978/149‑1 (MQ=255)
                                                                                      cGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGCACATCGCCGGGaa                                                              <  1:1394139/149‑1 (MQ=255)
                                                                                       gCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCCCGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGCACGTACAGCACCTCGCCGGGAAc                                                             >  2:876635/1‑149 (MQ=255)
                                                                                        ccAGATTTTGATATTGTTTGCCATCACATTGCAAGACGCTGAAGCCGAGTTGCTGTTTGACAAACAGCACGCTAAATGTAGTAGCTGGGTCTTGCCATTGTGGAAAATTCCAGCAACCAGCAGGAACGTAGAGCAAATCGCCGGGAACt                                                            >  2:1152933/1‑149 (MQ=255)
                                                                                         cAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGCACATCGCCGGGAACTa                                                           >  1:995627/1‑149 (MQ=255)
                                                                                           gATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGCACATCGCCGGGAACTAAc                                                         <  1:684319/149‑1 (MQ=255)
                                                                                            aTTTTGATATTGTTTGCCATCCCAATGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGCACATCGCCGGGAACTAACg                                                        >  2:164368/1‑149 (MQ=255)
                                                                                             ttttGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGCACATCGCCGGGAACTAACGt                                                       >  1:888303/1‑149 (MQ=255)
                                                                                                  atatTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGCACATCGCCGGGAACTAACGTTGCTc                                                  <  1:1353441/149‑1 (MQ=255)
                                                                                                     ttGTTTGCCATCCCATTGCACGACGCTGAAGCCGCGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGACCGTACAGCACCTCGGCGGGAACTTACGTTGATCCAg                                               >  2:1054790/1‑149 (MQ=255)
                                                                                                        tttGCCATCCCATTGCACGCCGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGCACATCGCCGGGAACTAACGTTGCTCCAgcgc                                            >  1:554539/1‑149 (MQ=255)
                                                                                                        tttGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGCACATCGCCGGGAACTAACGTTGCTCCAgcgc                                            <  2:253628/149‑1 (MQ=255)
                                                                                                        tttGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGCACATCGCCGGGAACTAACGTTGCTCCAgcgc                                            <  1:1178886/149‑1 (MQ=255)
                                                                                                          tGCCATCCGATTGCACGACGCTGTAGCCTAGTTGATGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGCACATCGCCGGGAACTAACGTTGCTCCAGCGCCg                                          <  2:6233/149‑1 (MQ=255)
                                                                                                          tGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGCACATCGCCGGGAACTAACGTTGCTCCAGCGCCg                                          <  2:1198640/149‑1 (MQ=255)
                                                                                                            ccATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTCCAGCACATCGCCGGGAACTAACGTTGCTCCAGCGCCGGt                                        >  2:982863/1‑149 (MQ=255)
                                                                                                             cATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCCAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGCACATCGCCGGGAACTAACGTTGCTCCAGCGCCGGTa                                       <  2:26554/149‑1 (MQ=255)
                                                                                                             cATCCCATTGCACGACGATGAAGCCCCGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTTTGGAAAATTCCCGCCCCCAGCAGGAACGTACAGCACATCGCCGGGCACAAACGTTGCTCCCGCGCCggaa                                       >  1:844630/1‑147 (MQ=255)
                                                                                                              aTCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGCACATCGCCGGGAACTAACGTTGCTCCAGCGCCGGTAt                                      <  2:440202/149‑1 (MQ=255)
                                                                                                                 ccATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGCACATCGCCGGGAACTAACGTTGCTCCAGCGCCGGTATcaa                                   <  1:421524/149‑1 (MQ=255)
                                                                                                                  cagtGCCCGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAATGGTAGTAGCGGGGGCTTGCCATTGTGGAACATTCCAGCCCCAAGCAGGAACGTACAGCACATCGCCGGGAGCTAACGTTGCTCCAGCGCCGGTTGcaac                                  >  2:857933/4‑149 (MQ=255)
                                                                                                                  cATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGCTCATCGCCGGGAACTAACGTTGCTCCAGCGCCGGTATcaac                                  <  1:1384618/149‑1 (MQ=255)
                                                                                                                     tGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGCACATCGCCGGGAACTAACGTTGCTCCAGCGCCGGTATCAACaaa                               <  2:951560/149‑1 (MQ=255)
                                                                                                                      gCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTACAGGTAGTCGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGCACATCGCCGGGAACTAACGTTGCTCCAGCGCCGGTATCAACAAAc                              <  1:285600/149‑1 (MQ=255)
                                                                                                                            cgCTGATGGCGAGTTGCTGTTTGCCTAACAGCACGCTATAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGCACATCGACGGGAACTAACGTTGCTCCAGCGCCGGTATCAACAAACTCACCt                        <  2:1371464/149‑1 (MQ=255)
                                                                                                                            cgCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCCGCAGGAACGTACAGCACATCGCCGGGAACTAACGTTGCTCCAGCGCCGGTATCAACAAACTCACCt                        <  2:1347581/149‑1 (MQ=255)
                                                                                                                            cgCTGAAGCCGAGTTGCTGCTTGCCAAACAGCACGCTAAATGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGCACATCGCCGGGAACTAACGTTGCTCCAGCGCCGGAATCAACAAACTCACCt                        <  2:12066/149‑1 (MQ=255)
                                                                                                                                 aaGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGCACATCGCCGGGAACTAACGTTGCTCCAGCGCCGGTATCAACAAACTCACCTTCCAg                   <  1:1232205/149‑1 (MQ=255)
                                                                                                                                    ccGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGCCCATCGCCGGGAACTAACGTTGCTCCAGCGCCGGTATCAACAAACTCACCTTCCAGcac                >  2:134851/1‑149 (MQ=255)
                                                                                                                                         ttGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTAGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGCACATCGCCGGCAACTAACGTTGCTACAGCGCCTGTATCAACACACTAACCTTCCAGCACAAttt           >  2:741154/1‑149 (MQ=255)
                                                                                                                                         gtGCTGTTTGGCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGCACGTACAGCACATCGCCGGGAACTAACGTTGCTCCAGCGCCGGTATCAACAAACTCACCTTCCAGCACAAttt           >  1:1005889/2‑149 (MQ=255)
                                                                                                                                         gtGCTGTTTGGCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGCACATCGCCGGGAACTAACGTTGCTCCAGCGCCGGTATCAACAAACTCACCTTCCAGCGCAAttt           >  1:615704/2‑149 (MQ=255)
                                                                                                                                            cTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCCCCAGCAGGAACGTACAGCACATCGCCGGGAACTAACGTTGCTCCAGCGCCGGTATCAACAATCTCACCTTCCAGCACAATTTccc        >  2:769777/1‑148 (MQ=255)
                                                                                                                                              gTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTGCAGCACATCGCCGGGAACTAACGTTGCTCCAGCGCCGGTATCAACAAACTCACCTTCCAGCACAATTTCCAGc      >  2:492323/1‑149 (MQ=255)
                                                                                                                                                  gCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGCACATCGCCGGGAACTAACGTTGCTCCAGCGCCGGTATCAACAAACTCACATTCCAGCACAATTTCCAGCCGTg  >  1:10701/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
TGCTCCTGCGGCTGCATTTGCATTTCATTGAGCGTTTGTAGCAGAAAAGAACCTATGCGTGGGCCGCGTCGGGCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGCACATCGCCGGGAACTAACGTTGCTCCAGCGCCGGTATCAACAAACTCACCTTCCAGCACAATTTCCAGCCGTG  >  NC_000913/4147902‑4148196

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