Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,148,049 C→G 5.7% 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.7% 822.4 / 30.9 230G90A (GGC→GCC) yijOAraC family putative transcriptional activator
Reads supporting (aligned to +/- strand):  ref base C (110/107);  new base G (7/6);  total (117/113)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.22e-01

Only 100 of 240 total aligned reads displayed.

CTGCTCCTGCGGCTGCATTTGCATTTCATTGAGCGTTTGTAGCAGAAAAGAACCTATGCGTGGGCCGCGTCGGGCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGCACATCGCCGGGAACTAACGTTGCTCCAGCGCCGGTATCAACAAACTCACCTTCCAGCACAATTTCCAGCCGT  >  NC_000913/4147901‑4148195
                                                                                                                                                    |                                                                                                                                                  
ctgctCCTGCGGCTGCGTTTGCATTTCATTGAGCGTTTGTTGCAGAAAAGAACCTATGCGTGGGCCGCGTCGGGCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGc                                                                                                                                                    <  2:629105/149‑1 (MQ=255)
ctgctCCTGCGGCTGCATTTGCATTTCATTGAGCGTTTGTAGCCGAAAAGAACCTATGCGTGGGCCGCGTCGGGCGACGTGTTGCTTCGCCAGAGTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGc                                                                                                                                                    >  1:541425/1‑149 (MQ=255)
 tgctCCTGCGGCTGCATTTGCATTTCATTGAGCGTTTGTAGCAGAAAAGAACCTATGCGTGGGCCGCGTCGGGCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGcc                                                                                                                                                   >  2:1612372/1‑149 (MQ=255)
 tgctCCTGCGGCTGCATTTGCATTTCATTGAGCGTTTGTAGCAGAAAAGAACCTATGCGTGGGCCGCGTCGGGCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGcc                                                                                                                                                   >  2:1699020/1‑149 (MQ=255)
    tCCTGCGGCTGCATTTGCATTTCATTGAGCGTTTGTAGCAGAAAAGAACCTATGCGTGGGCCGCGTCGGGCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCaaa                                                                                                                                                <  1:1604549/149‑1 (MQ=255)
      cTGCGGCTGCATTGGCATTTCATTGAGCGTTTGTAGCAGAAAAGAACATATGCGTGGGCCGCGTCGGGCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACa                                                                                                                                              <  3:145106/149‑1 (MQ=255)
       tGCGGCTGCATTTGCATTTCATTGAGCGTTTGTAGCAGAAAAGAACCTATGCGTGGGCCGCGTCGGGCGACGTGTTGCTCCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAg                                                                                                                                             <  2:1461453/149‑1 (MQ=255)
                 tttgcatttCATTGAGCGTTTGTAGCAGAAAAGAACCTATGCGTGGGCCGCGTCGGGCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAg                                                                                                                                   <  4:442470/149‑1 (MQ=255)
                 tttgcatttCATTGAGCGTTTGTAGCAGAAAAGAACCTATGCGTGGGCCGCGTCGGGCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAg                                                                                                                                   <  4:1309168/149‑1 (MQ=255)
                   tgcatttCATTGAGCGTTTGTAGCAGAAAAGAACCTATGCGTGGGCCGTGTCGGGCGACGTATTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGgt                                                                                                                                 >  4:1391890/1‑149 (MQ=255)
                    gcatttCATTGAGCGTTTGTAGCAGAAAAGAACCTATGCGTGGGCCGCGTCGGGCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGgta                                                                                                                                >  2:1784790/1‑149 (MQ=255)
                    gcatttCATTGAGCGTTTGTAGCAGAAAAGAACCTATGCGTGGGCCGCGTCGGGCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGgta                                                                                                                                <  2:1494561/149‑1 (MQ=255)
                           aTTGAGCGTTTGTAGCAGAAAAGAACCTATGCGTGGGCCGCGTCGGGCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAAc                                                                                                                                               <  2:282739/127‑1 (MQ=255)
                              gAGCGTTTGTAGCAGAAAAGAACCTATGCGTGGGCCGCGTCGGGCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCggggg                                                                                                                      <  1:1699020/149‑1 (MQ=255)
                                gCGTTTGTAGCAGAAAAGAACCTATGCGTGGGCCGCGTCGGGCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCt                                                                                                                    >  2:507291/1‑149 (MQ=255)
                                       tAGCAGAAAAGAACCTATGCGTGGGCCGCGTCGGGCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCAtt                                                                                                             <  1:639765/149‑1 (MQ=255)
                                       tAGCAGAAAAGAACCTATGCGTGGGCCGCGTCGGGCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCAtt                                                                                                             <  4:1638184/149‑1 (MQ=255)
                                        aGCAGAAAAGAACCTATGCGTGGGCCGCGTCGGGCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATtg                                                                                                            <  2:541425/149‑1 (MQ=255)
                                         gCAGAAAAGAACCTATGCGTGGGCCGCGTCGGGCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATtgt                                                                                                           <  2:359755/149‑1 (MQ=255)
                                         gCAGAAAAGAACCTATGCGTGGGCCGCGTCGGGCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATtgt                                                                                                           <  4:160169/149‑1 (MQ=255)
                                         gCAGAAAAGAACCTATGCGTGGGCCGCGTCGGGCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATtgt                                                                                                           <  4:1547466/149‑1 (MQ=255)
                                            gAAAAGAACCTATGCGTGGGCCGCGTCGGGCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATtgt                                                                                                           >  4:486031/1‑146 (MQ=255)
                                            gAAAAGAACCTATGCGTGGGCCGCGTCGGGCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATtgt                                                                                                           <  3:1339915/146‑1 (MQ=255)
                                            gAAAAGAACCTATGCGTGGGCCGCGTCGGGCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATtgt                                                                                                           <  3:486031/146‑1 (MQ=255)
                                            gAAAAGAACCTATGCGTGGGCCGCGTCGGGCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGa                                                                                                        >  1:1487538/1‑149 (MQ=255)
                                              aaaGAACCTATGCGTGGGCCGCGTCGGGCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGaaa                                                                                                      >  4:1482433/1‑149 (MQ=255)
                                              aaaGAACCTATGCGTGGGCCGCGTCGGGCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGaaa                                                                                                      >  4:1235934/1‑149 (MQ=255)
                                                                  gcgTCGGGCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGaa                                                                                  <  1:1466557/149‑1 (MQ=255)
                                                                    gTCGGGCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACcagcag                                                                                     <  3:1041253/144‑1 (MQ=255)
                                                                    gTCGGGCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACcagcag                                                                                     >  4:1041253/1‑144 (MQ=255)
                                                                         gCGACGTGTTGCTTCGCCAGATTTTGCTATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAg                                                                           >  2:857060/1‑149 (MQ=255)
                                                                         gCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGtagtag                                                                                                                            <  3:672578/100‑1 (MQ=255)
                                                                         gCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGtagtag                                                                                                                            >  4:947827/1‑100 (MQ=255)
                                                                         gCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGtagtag                                                                                                                            <  3:947827/100‑1 (MQ=255)
                                                                         gCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGtagtag                                                                                                                            >  4:672578/1‑100 (MQ=255)
                                                                         gCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAg                                                                           <  3:1235934/149‑1 (MQ=255)
                                                                            aCGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGcac                                                                        >  4:594824/1‑149 (MQ=255)
                                                                               tgtTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGgta                                                                                                                                >  2:1747475/1‑90 (MQ=255)
                                                                                gtTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATtgtg                                                                                                          <  4:1874710/111‑1 (MQ=255)
                                                                                   gCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGc                                                                                              >  1:1004839/1‑120 (MQ=255)
                                                                                   gCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGc                                                                                              <  2:1004839/120‑1 (MQ=255)
                                                                                     ttCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGCACATCGCCggg                                                               >  1:1353751/1‑149 (MQ=255)
                                                                                       cGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGCACATCGCCGGGaa                                                             <  1:211561/149‑1 (MQ=255)
                                                                                       cGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGCACATCGCCGGGaa                                                             <  3:1391890/149‑1 (MQ=255)
                                                                                        gCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGgta                                                                                                                                >  2:760287/1‑81 (MQ=255)
                                                                                        gCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGCACATCGCCGGGAAc                                                            <  1:857060/149‑1 (MQ=255)
                                                                                         ccAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTaaa                                                                                                                                    >  1:359238/1‑76 (MQ=255)
                                                                                         ccAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTaaa                                                                                                                                    <  2:1054741/76‑1 (MQ=255)
                                                                                         ccAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTaaa                                                                                                                                    >  1:1054741/1‑76 (MQ=255)
                                                                                         ccAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCa                                                                                                <  3:1505947/112‑1 (MQ=255)
                                                                                         ccAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACa                                                                            >  3:597707/1‑132 (MQ=255)
                                                                                         ccAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGCACATCGCCGGGAACt                                                           >  3:1264079/1‑149 (MQ=255)
                                                                                         ccAGATTTTGATATTGTTTGCCATCCCATTGCACGACGATGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACa                                                                            <  4:597707/132‑1 (MQ=255)
                                                                                                ttGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGCACATCGCCGGGAACTAACGTTg                                                    >  3:72246/1‑149 (MQ=255)
                                                                                                 tGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGCACATCGCCGGGAACTAACGTTGc                                                   <  3:1482433/149‑1 (MQ=255)
                                                                                                 tGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGCACATCGCCGGGAACTAACGTTGc                                                   <  4:1264079/149‑1 (MQ=255)
                                                                                                 tGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGCACATCGCCGGGAACTAACGTTGc                                                   <  2:1487538/149‑1 (MQ=255)
                                                                                                  gATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGCACATCGCCGGGAACTAACGTTGCt                                                  >  1:372272/1‑149 (MQ=255)
                                                                                                   atatTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGCACATCGCCGGGAACTAAc                                                        >  3:422543/1‑142 (MQ=255)
                                                                                                    tatTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGCACATCGCCGGGAACTAACGTTGCt                                                  >  1:972035/1‑147 (MQ=255)
                                                                                                         tttGCCCTCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGt                                                                               <  4:215353/113‑1 (MQ=255)
                                                                                                         tttGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGt                                                                               >  3:215353/1‑113 (MQ=255)
                                                                                                              cATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCgg                                                                                                                         >  2:1643646/1‑66 (MQ=255)
                                                                                                              cATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCa                                                                                                <  1:1582438/91‑1 (MQ=255)
                                                                                                              cATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCa                                                                                                >  2:1582438/1‑91 (MQ=255)
                                                                                                              cATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGCACATCGCCGGGAACTAACGTTGCTCCAGCGCCGGTa                                      >  1:1647206/1‑149 (MQ=255)
                                                                                                              cATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGCACATCGCCGGGAACTAACGTTGCTCCAGCGCCGGTa                                      >  2:1342838/1‑149 (MQ=255)
                                                                                                               aTCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGCACATCGCCGGGAACTAACGTTGCTcc                                                <  3:765332/138‑1 (MQ=255)
                                                                                                                tCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGCACATCGCCGGGAACTAACGTTGCTCCAGCGCCGGTATc                                    >  3:799917/1‑149 (MQ=255)
                                                                                                                tCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGCACATCGCCGGGAACTAACGTTGCTCCAGCGCCGGTATc                                    >  2:207032/1‑149 (MQ=255)
                                                                                                                   cATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGc                                                                                                                 <  2:1212079/69‑1 (MQ=255)
                                                                                                                   cATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAg                                                                           >  2:1836069/1‑107 (MQ=255)
                                                                                                                   cATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAg                                                                           <  1:1836069/107‑1 (MQ=255)
                                                                                                                   cATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGCACATCGCCGGGAACTAACGTTGCTCCAGCGCCGGTAGcaac                                 >  2:373635/1‑149 (MQ=255)
                                                                                                                        cACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGCACATCGCCGGGAACTAACGTTGCTCCAGCGCCGGTATCAACAAACt                            >  2:964252/1‑149 (MQ=255)
                                                                                                                         acgacgCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGCACATCGCCGGGAACTAACg                                                       >  3:1858746/1‑121 (MQ=255)
                                                                                                                          cgacgCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGTGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGCACATCGCCGGGAACTAACGTTGCTCCAGCGCCGGTATCAACAAACTCa                          <  3:1218323/149‑1 (MQ=255)
                                                                                                                           gacgCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGCACATCGCCGGGAACTAACGTTGCTCCAGCGCCGGTa                                      >  4:1427216/1‑136 (MQ=255)
                                                                                                                             cgCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATtgtg                                                                                                          <  1:1258318/66‑1 (MQ=255)
                                                                                                                                tGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAGCGTACAGCACATCGCCGGGAACTAACGTTGCTCCAGCGcc                                          >  4:751977/1‑127 (MQ=255)
                                                                                                                                tGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGCACATCGCCGGGAACTAACGTTGCTCCAGCGCCGGTATCAACAAACTCACCTTcc                    >  4:813971/1‑149 (MQ=255)
                                                                                                                                tGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGCACATCGCCGGGAACTAACGTTGCTCCAGCGCCGGTATCAACAAACTCACCTTcc                    <  4:799796/149‑1 (MQ=255)
                                                                                                                                 gAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGCACATCGCCGGGAACTAACGTTGCTCCAGCGCCGGTATCAACAAACTCACCTTCCa                   <  4:1502874/149‑1 (MQ=255)
                                                                                                                                   aGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACcagca                                                                                      >  1:1732429/1‑80 (MQ=255)
                                                                                                                                   aGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGCACATCGCCGGGAACTAACGTTGCTCCAGCGCCGGTATCAACAAACTCACCTTCCAGc                 >  3:52705/1‑149 (MQ=255)
                                                                                                                                    gCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGCACATCGCCGGGAACTAACGTTGCTCCa                                               <  3:1723916/118‑1 (MQ=255)
                                                                                                                                     ccGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGCACATCGCCGGGAACTAACGTTGCTCCAGCGCCGGTATCAACAAACTCACCTTCCAGcac               >  1:1023466/1‑149 (MQ=255)
                                                                                                                                      cGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGCACATCGCCGGGAACTAACGTTGCTCCAGCGCCGGTATCAACAAACTCACCTTCCAGcaca              >  1:1798102/1‑149 (MQ=255)
                                                                                                                                       gAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGCACATCGCCGGGAACTAACGTTGCTCCAGCGCCGGTATCAACAAACTCACCTTCCAGCACaa             >  4:1200042/1‑149 (MQ=255)
                                                                                                                                       gAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGCACATCGCCGGGAACTAACGTTGCTCCAGCGCCGGTATCAACAAACTCACCTTCCAGCACaa             >  4:26860/1‑149 (MQ=255)
                                                                                                                                         gTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGCACATCGCCGGGAACTAACGTTGCTCCAGCGCCGGTATCAACAAACTCACCTTCCAGCACAAtt           >  4:1836765/1‑149 (MQ=255)
                                                                                                                                          ttGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGCACATCGCCGGGAACTAACGTTGCTCCAGCGCCGGTATCAACAAACTCACCTTCCa                   <  4:1169704/140‑1 (MQ=255)
                                                                                                                                          ttGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGCACATCGCCGGGAACTAACGTTGCTCCAGCGCCGGTATCAACAAACTCACCTTCCa                   >  3:1169704/1‑140 (MQ=255)
                                                                                                                                           tGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGCACATCGCCGGGAACTAACGTTGCTCCAGCGCCGGTATCAACAAACTc                           <  4:461442/131‑1 (MQ=255)
                                                                                                                                             cTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGCACATCGCCGGGAACTAACGTTGCTCCAGCGCCGGTATCAACAAACTCACCTTCCAGCACAATTTCCa       <  3:26860/149‑1 (MQ=255)
                                                                                                                                              tGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGCACATCGCCGGGAACTAACGTTGCTCCAGCGCCGGTATCAACAAACTCACCTTCCAGCACAATTTCCAg      >  4:204164/1‑149 (MQ=255)
                                                                                                                                                tttGCCAAACAGCACGCTAAAGGTAGTTGCGGGGGCTTGCCATTGTGCAAAATTCAAGCCACCGGCAGGAGCGGGCAGGTCATCGCCGGGAACTAACGTTGCTCCAGCGCGGGTATCCACGAACTCACCGCCAAGCACAATTTCCAgac    >  2:1651993/1‑147 (MQ=255)
                                                                                                                                                  tGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGCACATCGCCGGGAACTAACGTTGCTCCAGCGCCGGTATCAACAAACTCACCTTCCAGCACAATTTCCAGCCGt  <  3:97583/149‑1 (MQ=255)
                                                                                                                                                    ccAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGCACATCGCCGGGAACTAACGTTGc                                                   >  2:1301179/1‑98 (MQ=255)
                                                                                                                                                    ccAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGCACATCGCCGGGAACTAACGTTGc                                                   <  1:1301179/98‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
CTGCTCCTGCGGCTGCATTTGCATTTCATTGAGCGTTTGTAGCAGAAAAGAACCTATGCGTGGGCCGCGTCGGGCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGCACATCGCCGGGAACTAACGTTGCTCCAGCGCCGGTATCAACAAACTCACCTTCCAGCACAATTTCCAGCCGT  >  NC_000913/4147901‑4148195

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