Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,148,049 C→G 9.3% 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 9.3% 216.4 / 10.6 64G90A (GGC→GCC) yijOAraC family putative transcriptional activator
Reads supporting (aligned to +/- strand):  ref base C (32/25);  new base G (4/2);  total (36/28)
Fisher's exact test for biased strand distribution p-value = 6.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.29e-01

CTGCGGCTGCATTTGCATTTCATTGAGCGTTTGTAGCAGAAAAGAACCTATGCGTGGGCCGCGTCGGGCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGCACATCGCCGGGAACTAACGTTGCTCCAGCGCCGGTATCAACAAACTCACCTTCCAGCACAATTTCCAGCCGT  >  NC_000913/4147907‑4148195
                                                                                                                                              |                                                                                                                                                  
cTGCGGCTGCATTTGCATTTCATTGAGCGTTTGTAGCAGAAAAGAACCTATGCGTGGGCCGCGTCGGGCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACa                                                                                                                                              >  1:191013/1‑149 (MQ=255)
cTGCGGCTGCATTTGCATTTCATTGAGCGTTTGTAGCAGAAAAGAACCTATGCGTGGGCCGCGTCGGGCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGACAAACa                                                                                                                                              <  2:491070/149‑1 (MQ=255)
     gCTGCATTTGCATTTCATTGAGCGTTTGTAGCAGAAAAGAACCTATGCGTGGGCCGCGTCGGGCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGGCAAACAGCaac                                                                                                                                         <  2:179959/149‑3 (MQ=255)
      cTGCATTTGCATTTCATTGAGCGTTTGTAGCAGAAAAGAACCTATGCGTGGGCCGCGTCGGGCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGc                                                                                                                                        <  4:115179/149‑1 (MQ=255)
                atttCATTGAGCGTTTGTAGCAGAAAAGAACCTATGCGTGGGCCGCGTCGGGCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGGACGCTAAAGgtagt                                                                                                                              <  4:460984/149‑1 (MQ=255)
                atttCATTGAGCGTTTGTAGCAGAAAAGAACCTATGCGTGGGCCGCGTCGGGCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGgtagt                                                                                                                              >  4:21308/1‑149 (MQ=255)
                 tttCATTGAGCGTTTGTAGCAGAAAAGAACCTATGCGTGGGCCGCGTCGGGCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAg                                                                                                                                   >  1:161167/1‑143 (MQ=255)
                 tttCATTGAGCGTTTGTAGCAGAAAAGAACCTATGCGTGGGCCGCGTCGGGCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAg                                                                                                                                   <  2:161167/143‑1 (MQ=255)
                 tttCATTGAGCGTTTGTAGCAGAAAAGAACCTATGCGTGGGCCGCGTCGGGCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGgtagta                                                                                                                             >  1:544894/1‑149 (MQ=255)
                     aTTGAGCGTTTGTAGCAGAAAAGAACCTATGCGTGGGCCGCGTCGGGCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAg                                                                                                                                   <  2:196953/139‑1 (MQ=255)
                     aTTGAGCGTTTGTAGCAGAAAAGAACCTATGCGTGGGCCGCGTCGGGCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAg                                                                                                                                   >  1:196953/1‑139 (MQ=255)
                     aTTGAGCGTTTGTAGCAGAAAAGAACCTATGCGTGGGCCGCGTCGGGCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCgg                                                                                                                         >  1:414443/1‑149 (MQ=255)
                            gTTTGTAGCAGAAAAGAACCTATGCGTGGGCCGCGTCGGGCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTg                                                                                                                  >  1:551471/1‑149 (MQ=255)
                                 tAGCAGAAAAGAACCTATGCGTGGGCCGCGTCGGGCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCAtt                                                                                                             >  2:560053/1‑149 (MQ=255)
                                   gCAGAAAAGAACCTATGCGTGGGCCGCGTCGGGCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATtgt                                                                                                           >  3:304819/1‑149 (MQ=255)
                                   gCAGAAAAGAACCTATGCGTGGGCCGCGTCGGGCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATtgt                                                                                                           <  2:162025/149‑1 (MQ=255)
                                   gCAGAAAAGAACCTATGCGTGGGCCGCGTCGGGCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATtgt                                                                                                           <  4:99687/149‑1 (MQ=255)
                                   gCAGAAAAGAACCTATGCGTGGGCCGCGTCGGGCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATtgt                                                                                                           <  1:608627/149‑1 (MQ=255)
                                      gAAAAGAACCTATGCGTGGGCCGCGTCGGGCGACGTGTTGCTTCGCCACATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGa                                                                                                        >  3:523119/1‑149 (MQ=255)
                                        aaaGAACCTATGCGTGGGCCGCGTCGGGCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGaaa                                                                                                      >  4:45638/1‑149 (MQ=255)
                                         aaGAACCTATGCGTGGGCCGCGTCGGGCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGaaaa                                                                                                     <  1:560053/149‑1 (MQ=255)
                                            aaCCTATGCGTGGGCCGCGTCGGGCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTc                                                                                                  >  1:88982/1‑149 (MQ=255)
                                            aaCCTATGCGTGGGCCGCGTCGGGCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTc                                                                                                  <  3:232589/149‑1 (MQ=255)
                                              ccTATGCGTGGGCCGCGTCGGGCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCa                                                                                                <  1:476091/149‑1 (MQ=255)
                                                   gCGTGGGCCGCGTCGGGCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGccac                                                                                           <  4:665416/149‑1 (MQ=255)
                                                     gTGGGCCGCGTCGGGCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGccacca                                                                                         <  3:321642/149‑1 (MQ=255)
                                                          ccGCGTCGGGCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAgg                                                                                    <  1:647408/149‑1 (MQ=255)
                                                            gcgTCGGGCGACGTGTTGCTTCGCCACATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGaa                                                                                  <  4:523119/149‑1 (MQ=255)
                                                                cGGGCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTa                                                                              <  3:5883/149‑1 (MQ=255)
                                                                  ggCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACa                                                                            <  4:589955/149‑1 (MQ=255)
                                                                  ggCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACa                                                                            >  1:436558/1‑149 (MQ=255)
                                                                     gACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGca                                                                         >  4:479088/1‑149 (MQ=255)
                                                                          gtTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGCACATCg                                                                    >  1:393807/1‑149 (MQ=255)
                                                                                  gCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGCACATCGCCGGGAAc                                                            <  1:105476/149‑1 (MQ=255)
                                                                                  gCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGCACATCGCCGGGAAc                                                            >  4:458947/1‑149 (MQ=255)
                                                                                   ccAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGCACATCGCCGGGAACt                                                           >  3:52745/1‑149 (MQ=255)
                                                                                             atatTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGCACATCGCCGGGAACTAACGTTGCTc                                                 >  4:579755/1‑149 (MQ=255)
                                                                                             atatTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGCACATCGCCGGGAACTAACGTTGCTc                                                 >  1:173644/1‑149 (MQ=255)
                                                                                              tatTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGCACATCGCCGGGAACTAACGTTGCTcc                                                >  4:537247/1‑149 (MQ=255)
                                                                                                ttGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGCACATCGCCGGGAACTAACGTTGCTCCAg                                              >  3:20107/1‑149 (MQ=255)
                                                                                                    ttGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGCACATCGCCGGGAACTAACGTTGCTCCAGCGcc                                          >  4:124337/1‑149 (MQ=255)
                                                                                                         aTCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGCACATCGCCGGGAACTAACGTTGCTCCAGCGCCGGTAt                                     >  3:424022/1‑149 (MQ=255)
                                                                                                           cccATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGCACATCGCCGGGAACTAACGTTGCTCCAGCGCCGGTATca                                   >  3:269875/1‑149 (MQ=255)
                                                                                                              aTTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGCACATCGCCGGGAACTAACGTTGCTCCAGCGCCGGTATcaaca                                >  1:528373/1‑149 (MQ=255)
                                                                                                                tGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGCACATCGCCGGGAACTAACGTTGCTCCAGCGCCGGTATCAACaaa                              <  2:88982/149‑1 (MQ=255)
                                                                                                                tGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGCACATCGCCGGGAACTAACGTTGCTCCAGCGCCGGTATCAACaaa                              <  4:269875/149‑1 (MQ=255)
                                                                                                                     gacgCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGCACATCGCCGGGAACTAACGTTGCTCCAGCGCCGGTATCAACAAACTCAc                         >  3:559158/1‑149 (MQ=255)
                                                                                                                     gacgCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGCACATCGCCGGGAACTAACGTTGCTCCAGCGCCGGTATCAACAAACTCAc                         <  3:21308/149‑1 (MQ=255)
                                                                                                                      acgCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGCACATCGCCGGGAACTAACGTTGCTCCAGCGCCGGTATCAACAAACTCAcc                        <  2:544894/149‑1 (MQ=255)
                                                                                                                          tGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGCACATCGCCGTGAACTAACGTTGCTCCAGCGCCGGTATCAACAAACTCACCTGcc                    <  2:191013/149‑1 (MQ=255)
                                                                                                                          tGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGCACATCGCCGGGAACTAACGTTGCTCCAGCGCCGGCATCAACAAACTCACCTTcc                    >  4:238971/1‑149 (MQ=255)
                                                                                                                           gAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGCACATCGCCGGGAACTAACGTTGCTCCAGCGCCGGTATCAACAAACTCACCTTCCa                   <  2:551471/149‑1 (MQ=255)
                                                                                                                             aGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGCACATCGCCGGGAACTAACGTTGCTCCAGCGCCGGTATCAACAAACTCACCTTCCAGc                 >  4:149800/1‑149 (MQ=255)
                                                                                                                             aGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGCACATCGCCGGGAACTAACGTTGCTCCAGCGCCGGTATCAACAAACTCACCTTCCAGc                 >  1:44528/1‑149 (MQ=255)
                                                                                                                              gCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGCACATCGCCGGGAACTAACGTTGCTCCAGCGCCGGTATCAACAAACTCACCTTCCAGca                >  2:440262/1‑149 (MQ=255)
                                                                                                                                cGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGCACATCGCCGGGAACTAACGTTGCTCCAGCGCCGGTATCAACAAACTCACCTTCCAGcaca              >  4:616016/1‑149 (MQ=255)
                                                                                                                                   gTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGCACATCGCCGGGAACTAACGTTGCTCCAGCGCCGGTATCAACAAACTCACCTTCCAGCACAAtt           >  3:318002/1‑149 (MQ=255)
                                                                                                                                    gtGCTGTTTGGCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGCACATCGCCGGGAACTAACGTTGCTCCAGCGCCGGTATcaac                                 >  3:536766/2‑126 (MQ=255)
                                                                                                                                    gtGCTGTTTGGCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGCACATCGCCGGGAACTAACGTTGCTCCAGCGCCGGTATcaac                                 <  4:536766/125‑1 (MQ=255)
                                                                                                                                    gtGCTGTTTGGCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGCACATCGCCGGGAACTAACGTTGCTCCAGCGCCGGTATCAACAAACTCACCTTCCAGCACAAttt          >  2:625691/2‑149 (MQ=255)
                                                                                                                                     tGCTGTTTGGCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGCACATCGCCGGGAACTAACGTTGCTCCAGCGCCGGTATCAACAAACTCACCTTCCAGCACAATTTc         >  2:47826/1‑149 (MQ=255)
                                                                                                                                     tGCTGTTTGGCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGCACATCGCCGGGAACTAACGTTGCTCCAGCGCCGGTATCAACAAACTCACCTTCCAGCACAATTTc         >  3:575219/1‑149 (MQ=255)
                                                                                                                                       cTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGCACATCGCCGGGAACTAACGTTGCTCCAGCGCCGGTATCAACAAACTCACCTTCCAGCACAATTTCCa       <  2:44528/149‑1 (MQ=255)
                                                                                                                                            tGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGCACATCGCCGGGAACTAACGTTGCTCCAGCGCCGGTATCAACAAACTCACCTTCCAGCACAATTTCCAGCCGt  <  2:393807/149‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                  
CTGCGGCTGCATTTGCATTTCATTGAGCGTTTGTAGCAGAAAAGAACCTATGCGTGGGCCGCGTCGGGCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGCACATCGCCGGGAACTAACGTTGCTCCAGCGCCGGTATCAACAAACTCACCTTCCAGCACAATTTCCAGCCGT  >  NC_000913/4147907‑4148195

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