Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,147,517 A→G 6.8% Y267Y (TAT→TAC)  yijO ← AraC family putative transcriptional activator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,147,5170AG 6.8% 217.4 / 3.1 87Y267Y (TAT→TAC) yijOAraC family putative transcriptional activator
Reads supporting (aligned to +/- strand):  ref base A (41/40);  new base G (3/3);  total (44/43)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.41e-01

CTATGTGCAATGCAGCATCTACGCATTCCTGCGTGAGCCGCAGCGGGTAATGAGCGCGGTGCGCAAAGTGCTTTCTGCGCCGCAGCAAACCCATCTTATTCTGGAGTAGTCGGTTTTTCTGTCAGTTGGCTGTGGTACTGCCGGCGATATTCCGACGGCGAGCGTTCGGTGTTTTTACGAAACAGACGGCAGAAGTAGTTGCTGTCGACAAAGCCGCAGGCGTGCGCCACTTCTTTTACCTTCAGGTCATAGCCTTTCAATAACGTCTTAGCGTGCTCCAGTCGCGTGTGATT  >  NC_000913/4147371‑4147663
                                                                                                                                                  |                                                                                                                                                  
cTATGTGCAATGCAGCATCTACGCATTCCTGCGTGAGCCGCAGCGGGTAATGAGCGCGGTGCGCAAAGTGCTTTCTGCGCCGCAGCAAACCCATCTTATTCTGGAGTAGTCGGTTTTTCTGTCAGTTGGCTGTGGTACTGCCGGCGata                                                                                                                                                  <  4:271616/149‑1 (MQ=255)
 tATGTGCAATGCAGCATCTACGCATTCCTGCGTGAGCCGCAGCGGGTAATGAGCGCGGTGCGCAAAATGCTTTCTGCGCCGCAGCAAACCCATCTTATTCTGGAGTAGTCGGTTTTTCTGTCAGTTGGCTGTGGTACTGCCGGCGatat                                                                                                                                                 >  8:156050/1‑149 (MQ=255)
     tgCAATGCAGCATCTACGCATTCCTGCGTGAGCCGCAGCGGGTAATGAGCGCGGTGCGCAAAGTGCTTTCTGCGCCGCAGCAAACCCATCTTATTCTGGAGTAGTCGGTTTTTCTGTCAGTTGGCTGTGGTACTGCCGGCGATATTcc                                                                                                                                              >  1:301576/1‑148 (MQ=255)
          tGCAGCATCTACGCATTCCTGCGTGAGCCGCAGCGGGTAATGAGCGCGGTGCGCAAAGTGCTTTCTGCGCCGCAGCACACCCATCTTATTCTGGAGTAGTCGGTTTTTCTGTCAGTTGGCTGTGGTACTGCCGGCGATATTCCGACGGc                                                                                                                                        >  1:29760/1‑149 (MQ=255)
          tGCAGCATCTACGCATTCCTGCGTGAGCCGCAGCGGGTAATGAGCGCGGTGCGCAAAGTGCTTTCTGCGCCGCAGCAAACCCATCTTATTCTGGAGTAGTCGGTTTTTCTGTCAGTTGGCTGTGGTACTGCCGGCGATATTCCGACGGc                                                                                                                                        >  7:211494/1‑149 (MQ=255)
                 tCTACGCATTCCTGCGTGAGCCGCAGCGGGTAATGAGCGCGGTGCGCAAAGTGCTTTCTGCGCCGCAGCAAACCCATCTTATTCTGGAGTAGTCGGTTTTTCTGTCAGTTGGCTGTGGTACTGCCGGCGATATTCCGACGGCGAGCGtt                                                                                                                                 >  4:77825/1‑149 (MQ=255)
                  cTACGCATTCCTGCGTGAGCCGCAGCGGGTAATGAGCGCGGTGCGCAAAGTGCTTTCTGCGCCGCAGCAAACCCATCTTATTCTGGAGTAGTCGGTTTTTCTGTCAGTTGGCTGTGGTACTGCCGGCGATATTCCGACGGCGAGCGTTc                                                                                                                                >  3:44451/1‑149 (MQ=255)
                  cTACGCATTCCTGCGTGAGCCGCAGCGGGTAATGAGCGCGGTGCGCAAAGTGCTTTCTGCGCCGCAGCAAACCCATCTTATTCTGGAGTAGTCGGTTTTTCTGTCAGTTGGCTGTGGTACTGCCGGCGATATTCCGACGGCGAGCGTTc                                                                                                                                >  1:201308/1‑149 (MQ=255)
                   tACGCATTCCTGCGTGAGCCGCAGCGGGTAATGAGCGCGGTGCGCAAAGTGCTTTCTGCGCCGCAGCAAACCCATCTTATTCTGGAGTAGTCGGTTTTTCTGTCAGTTGGCTGTGGTACTGCCGGCGATATTCCGACGGCGAGCGTTCg                                                                                                                               >  3:188012/1‑149 (MQ=255)
                   tACGCATTCCTGCGTGAGCCGCAGCGGGTAATGAGCGCGGTGCGCAAAGTGCTTTCTGCGCCGCAGCAAACCCATCTTATTCTGGAGTAGTCGGTTTTTCTGTCAGTTGGCTGTGGTACTGCCGGCGATATTCCGACGGCGAGCGTTCg                                                                                                                               >  8:30030/1‑149 (MQ=255)
                    aCGCATTCCTGCGTGAGCCGCAGCGGGTAATGAGCGCGGTGCGCAAAGTGCTTTCTGCGCCGCAGCAAACCCATCTTATTCTGGAGTAGTCGGTTTTTCTGTCAGTTGGCTGTGGTACTGCCGGAGATATTCCGCCGGCGAGCGGTCgg                                                                                                                              >  8:367320/1‑149 (MQ=255)
                     cGCATTCCTGCGTGAGCCGCAGCGGGTAATGAGCGCGGTGCGCAAAGTGCTTTCTGCGCCGCAGCAAACCCATCTTATTCTGGAGTAGTCGGTTTTTCTGTCAGTTGGCTGTGGTACTGCCGGCGGTATTCCGACGGCGAGCGTTCGgt                                                                                                                             >  6:330324/1‑149 (MQ=255)
                     cGCATTCCTGCGTGAGCCGCAGCGGGTAATGAGCGCGGTGCGCAAAGTGCTTTCTGCGCCGCAGCAAACCCATCTTATTCTGGAGTAGTCGGTTTTTCTGTCAGTTGGCTGTGGTACTGCCGGCGGTATTCCGACGGCGAGCGTTCGgt                                                                                                                             >  6:324430/1‑149 (MQ=255)
                     cGCATTCCTGCGTGAGCCGCAGCGGGTAATGAGCGCGGTGCGCAAAGTGCTTTCTGCGCCGCAGCAAACCCATCTTATTCTGGAGTAGTCGGTTTTTCTGTCAGTTGGCTGTGGTACTGCCGGCGGTATTCCGACGGCGAGCGTTCGgt                                                                                                                             >  6:324429/1‑149 (MQ=255)
                     cGCATTCCTGCGTGAGCCGCAGCGGGTAATGAGCGCGGTGCGCAAAGTGCTTTCTGCGCCGCAGCAAACCCATCTTATTCTGGAGTAGTCGGTTTTTCTGTCAGTTGGCTGTGGTACTGCCGGCGATATTCCGACGGCGAGCGTTCg                                                                                                                               <  4:8663/147‑1 (MQ=255)
                     cGCATTCCTGCGTGAGCCGCAGCGGGTAATGAGCGCGGTGCGCAAAGTGCTTTCTGCGCCGCAGCAAACCCATCTTATTCTGGAGTAGTCGGTTTTTCTGTCAGTTGGCTGTGGTACTGCCGGCGATATTCCGACGGCGAGCGTTCGgt                                                                                                                             >  1:18082/1‑149 (MQ=255)
                     cGCATTCCTGCGTGAGCCGCAGCGGGTAATGAGCGCGGTGCGCAAAGTGCTTTCTGCGCCGCAGCAAACCCATCTTATTCTGGAGTAGTCGGTGTTTCTGTCAGGTGGCTGTGGTACTGCCGGCGATATTCCGACGGCGAGCGtt                                                                                                                                 <  3:218637/145‑1 (MQ=255)
                     cGCATTCCTGCGTGAGCCGCAGCGGGTAATGAGCGCGGTGCGCAAAGTGCTTTCTGCGCCGCAGCAAACCCATCTTATTCTGGAGTAGTCGGGTTTTCTGTCAGTTGGCTGTGGTACTGCCGGCGATATTCCGACGGCGAGCGTTCGgt                                                                                                                             >  4:61931/1‑149 (MQ=255)
                      gCATTCCTGCGTGAGCCGCAGCGGGTAATGAGCGCGGTGCGCAAAGTGCTTTCTGCGCCGCAGCAAACCCATCTTATTCTGGAGTAGACGGTTTTTCTGTCAGTTGGCTGTGGTACTGCCGGCGATATTCCGACGGCGAGCGTTCGgtg                                                                                                                            >  4:40293/1‑149 (MQ=255)
                         ttCCTGCGTGAGCCGCAGCGGGTAATGAGCGCGGTGCGCAAAGTGCTTTCTGCGCCGCAGCAAACCCATCTTATTCTGGAGTAGTCGGTTTTTCTGTCAGTTGGCTGTGGTACTGCCGGCGATATTCCGACGGCGAGCGTTCGGTGttt                                                                                                                         >  2:105923/1‑149 (MQ=255)
                            cTGCGTGAGCCGCAGCGGGTAATGAGCGCGGTGCGCAAAGTGCTTTCTGCGCCGCAGCAAACCCATCTTATTCTGGAGTAGTCGGTTTTTCTGTCAGTTGGCTGTGGTACTGCCGGCGATATTCCGACGGCGAGCGTTCGGTGttttt                                                                                                                       <  6:147310/148‑1 (MQ=255)
                            cTGCGTGAGCCGCAGCGGGTAATGAGCGCGGTGCGCAAAGTGCTTTCTGCGCCGCAGCAAACCCATCTTATTCTGGAGTAGTCGGTTTTTCTGTCAGTTGGCTGTGGTACTGCCGGCGATATTCCGACGGCGAGCGTTCGGTGTTTTTa                                                                                                                      >  2:176571/1‑149 (MQ=255)
                            cTGCGTGAGCCGCAGCGCGTAATGAGCGCGGTGCGCAAAGTGCTTTCTGCGCCGCAGCAAACCCATCTTATTCAGGAGTAGTCGGTTTTTCTGTCAGTTCGCTGTGGTACTGCCGGCGAAATTCCGACGGCGAGCGTTCGGTGTTTTTa                                                                                                                      >  2:74979/1‑149 (MQ=255)
                               cGTGAGCCGCAGCGGGTAATGAGCGCGGTGCGCAAAGTGCTTTCTGCGCCGCAGCAAACCCATCTTATTCTGGAGTAGTCGGTTTTTCTGTCAGTTGGCTGTGGTACGGCCGGCGATATTCCGACGGCGAGCGTTCGGTGTTTTTAc                                                                                                                     <  1:114354/147‑1 (MQ=255)
                                   aGCCGCAGCGGGTAATGAGCGCGGTGCGCAAAATGCTTTCTGCGCCGCAGCAAACCCATCTTATTCTGGAGTAGTCGGTTTTTCTGTCAGTTGGCTGTGGTACTGCCGGCGATATTCCGACGGCGAGCGTTCGGTGTTTTTACGAAACa                                                                                                               <  7:156050/149‑1 (MQ=255)
                                      cGCAGCGGGTAATGAGCGCGGTGCGCAAAGTGCTTTCTGCGCCGCAGCAAACTCATCTTATTCTGGAGTAGTCGGTTTTTCTGTCAGTTGGCTGTGGTACTGCCGGCGGTATTCCGACGGCGAGCGTTCGGTGTTTTTACGAAACATAc                                                                                                            <  5:324430/149‑1 (MQ=255)
                                      cGCAGCGGGTAATGAGCGCGGTGCGCAAAGTGCTTTCTGCGCCGCAGCAAACCCATCTTATTCTGGAGTAGTCGGTTTTTCTGTCAGTTGGCTGTGGTACTGCCGGCGGTATTCCGACGGCGAGCGTTCGGTGTTTTTACGAAACAGAc                                                                                                            <  5:330324/149‑1 (MQ=255)
                                      cGCAGCGGGTAATGAGCGCGGTGCGCAAAGTGCTTTCTGCGCCGCAGCAAACCCATCTTATTCTGGAGTAGTCGGTTTTTCTGTCAGTTGGCTGTGGTACTGCCGGCGGTATTCCGACGGCGAGCGTTCGGTGTTTTTACGAAACAGAc                                                                                                            <  5:324429/149‑1 (MQ=255)
                                       gCAGCGGGTAATGAGCGCGGTGCGCAAAGTGCTTTCTGCGCCGCAGCAAACCCATCTTATTCTGGAGTAGTCGGTTTTTCTGTCAGTTGGCTGTGGTACTGCCGGCGATATTCCGACGGCGAGCGTTCGGTGTTTTTACGAAACAGACg                                                                                                           >  6:287405/1‑149 (MQ=255)
                                                 aTGAGCGCGGTGCGCAAAGTGCTTTCTGCGCCGCAGCAAACCCATCTTATTCTGGAGTAGTCGGTTTTTCTGTCAGTTGGCTGTGGTACTGCCGGCGATATTCCGACGGCGAGCGTTCGGTGTTTTTACGAAACAGACGGCAGAagtag                                                                                                 <  4:44451/149‑1 (MQ=255)
                                                 aTGAGCGCGGTGCGCAAAGTGCTTTCTGCGCCGCAGCAAACCCATCTTATTCTGGAGTAGTCGGTTTTTCTGTCAGTTGGCTGTGGTACTGCCGGCGATATTCCGACGGCGAGCGTTCGGTGTTTTTACGAAACAGACGGCAGAagtag                                                                                                 <  3:77825/149‑1 (MQ=255)
                                                 aTGAGCGCGGTGCGCAAAGTGCTTTCTGCGCCGCAGCAAACCCATCTTATTCTGGAGTAGTCGGTTTTTCTGTCAGTTGGCTGTGGTACTGCCGGCGATATTCCGACGGCGAGCGTTCGGTGTTTTTACGAAACAGACGGCAGAagta                                                                                                  <  5:300706/148‑1 (MQ=255)
                                                    aGCGCGGTGCGCAAAGTGCTTTCTGCGCCGCAGCAAAGCCATCTTATTCTGGAGTAGTCGGGTTTTGTGTGAGTTGGCTGTGGTACTGCCGGCGATATTCCGAGGGCGAGGGTTCGGTGTTTTTACGAAACAGACGGCAGAAGTAGTTg                                                                                              <  1:74979/149‑1 (MQ=255)
                                                    aGCGCGGTGCGCAAAGTGCTTTCTGCGCCGCAGCAAACCCATCTTATTCTGGAGTAGTCGGTTTTTCTGTCAGTTGGCTGTGGTACTGCCGGCGATATTCCGACGGCGAGCGTTCGGTGTTTTTACGAAACAGACGGCAGAAGTAGtt                                                                                               <  1:105923/148‑1 (MQ=255)
                                                    aGCGCGGTGCGCAAAGTGCTTTCTGCGCCGCAGCAAACCCATCTTATTCTGGAGTAGTCGGTTTTTCTGTCAGTTGGCTGTGGTACTGCCGGCGATATTCCGACGGCGAGCGTTCGGTGTTTTTACGAAACAGACGGCAGAAGTAGTTg                                                                                              <  1:176571/149‑1 (MQ=255)
                                                        cgGTGCGCAAAGTGCTTTCTGCGCCGCAGCAAACCCATCTTATTCTGGAGTAGTCGGTTTTTCTGTCAGTTGGCTGTGGTACTGCCGGCGATATTCCGACGGCGAGCGTTCGGTGTTTTTACGAAACAGACGGCAGAAGTAGTTGCTGt                                                                                          <  5:247938/149‑1 (MQ=255)
                                                         ggTGCGCAAAGTGCTTTCTGCGCCGCAGCAAACCCATCTTATTCTGGAGTAGTCGGTTTTTCTGTCAGTTGGCTGTGGTACTGCCGGCGATATTCCGACGGCGAGCGTTCGGTGTTTTTACGAAACAGACGGCAGAAGTAGTTGCTGTc                                                                                         >  2:93107/1‑149 (MQ=255)
                                                         ggTGCGAAAAGTGGTTTCTGGGCCGGAGGAAAACCAACTTTTTCTGGAGTACGCGGTTTTTGTGTCAGTTTTTAGGGTTCCTGCCGGGGATATTGATTCGGGGTTCGGTCGGTGTTTTTAAGAAACAGACGGCAGAAGTAGTTGCTGTc                                                                                         <  4:63197/149‑1 (MQ=255)
                                                               cAAAGTGCTTTCTGCGCCGCAGCAAACCCATCTTATTCTGGAGTAGTCGGTTTTTCTGTCAGTTGGCTGTGGTACTGCCGGCGATATTCCGACGGCGAGCGTTCGGTGTTTTTACGAAACAGACGGCAGAAGTAGTTGCTGTCGACa                                                                                     >  5:13182/1‑147 (MQ=255)
                                                                aaaGTGCTTTCTGCGCCGCAGCAAACCCATCTTATTCTGGAGTAGTCGGTTTTTCTGTCAGTTGGCTGTGGTACTGCCGGCGATATTCCGACGGCGAGCGTTCGGTGTTTTTACGAAACAGACGGCAGAAGTAGTTGCTGTCGACAAAg                                                                                  >  1:237490/1‑149 (MQ=255)
                                                                 aaGTGCTTTCTGCGCCGCAGCAAACCCATCTTATTCTGGAGTAGTCGGTTTTTCTGTCAGTTGGCTGTGGTACTGCCGGCGATATTCCGACGGCGAGCGTTCGGTGTTTTTACGAAACAGACGGCAGAAGTAGTTGCTGTCGACAAAGc                                                                                 >  5:84269/1‑149 (MQ=255)
                                                                      cTTTCTGCGCCGCAGCAAACCCATCTTATTCTGGAGTAGTCGGTTTTTCTGTCAGTTGGCTGTGGTACTGCCGGCGATATTCCGACGGCGAGCGTTCGGTGTTTTTACGAAACAGACGGCAGAAGTAGTTGCTGTCGACAAAGCCGCAg                                                                            >  1:153223/1‑149 (MQ=255)
                                                                             cgccgcAGCAAACCCATCTTATTGTGGAGTAGTCGGTTTTTCTGTCAGTTGGCTGTGGTACTGCCGGCGATATTCCGACGGCGAGCGTTCGGTGTTTTTACGAAACAGACGGCAGAAGTAGTTGCTGTCGACAAAGCCGCAGGCGTgc                                                                      <  2:301576/148‑1 (MQ=255)
                                                                              gccgcAGCAAACCCATCTTATTCTGGAGTAGTCGGTTTTTCTGTCAGTTGGCTGTGGTACTGCCGGCGATATTCCGACGGCGAGCGTTCGGTGTTTTTACGAAACAGACGGCAGAAGTAGTTGCTGTCGa                                                                                       >  3:361660/1‑130 (MQ=255)
                                                                              gccgcAGCAAACCCATCTTATTCTGGAGTAGTCGGTTTTTCTGTCAGTTGGCTGTGGTACTGCCGGCGATATTCCGACGGCGAGCGTTCGGTGTTTTTACGAAACAGACGGCAGAAGTAGTTGCTGTCGa                                                                                       <  4:361660/130‑1 (MQ=255)
                                                                              gccgcAGCAAACCCATCTTATTCTGGAGTAGTCGGTTTTTCTGTCAGTTGGCTGTGGTACTGCCGGCGATATTCCGACGGCGAGCGTTCGGTGTTTTTACGAAACAGACGGCAGAAGTAGTTGCTGTCGACAAAGCCGCAGGCGTgcgc                                                                    <  3:61931/149‑1 (MQ=255)
                                                                                  cagcaAACCCATCTTATTCTGGAGTAGTCGGTTTTTCTGTCAGTTGGCTGTGGTACTGCCGGCGATATTCCGACGGCGAGCGTTCGGTGTTTTTACGAAACAGACGGCAGAAGTAGTTGCTGTCGACAAAGCCGCAGGCGTGCGCCAct                                                                <  1:89794/149‑1 (MQ=255)
                                                                                       aaCCCATCTTATTCTGGAGTAGTCGGTTTTTCTGTCAGTTGGCTGTGGTACTGCCGGCGATATTCCGACGGCGAGCGTTCGGTGTTTTTACGAAACAGACGGCAGAAGTAGTTGCTGTCGACAAAGCCGCAGGCGTGCGCCACTTCttt                                                           <  2:153223/149‑1 (MQ=255)
                                                                                        aCCCATCTTATTCTGGAGTAGTCGGTTTTTCTGTCAGTTGGCTGTGGTACTGCCGGCGATATTCCGCCGGCGAGCGTTCGGTGTTTTTACGAAACAGACGGCAGAAGTAGTTGCTGTCGACAAAGCCGCAGGCGTGCGCCACTTCtttt                                                          <  2:211514/149‑1 (MQ=255)
                                                                                        aCCCATCTTATTCTGGAGTAGTCGGTTTTTCTGTCAGTTGGCTGTGGTACTGCCGGCGATATTCCGACGGCGAGCGTTCGGTGTTTTTACGAAACAGACGGCAGAAGTAGTTGCTGTCGACAAAGCCGCAGGCGTGCGCCACTTCtttt                                                          <  2:169550/149‑1 (MQ=255)
                                                                                          ccATCTTATTCTGGAGTAGTCGGTTTTTCTGTCAGTTGGCTGTGGTACTGCCGGCGATATTCCGACGGCGAGCGTTCGGTGTTTTTACGAAACAGACGGCAGAAGTAGTTGCTGTCGACaaa                                                                                   >  5:47687/1‑122 (MQ=255)
                                                                                          ccATCTTATTCTGGAGTAGTCGGTTTTTCTGTCAGTTGGCTGTGGTACTGCCGGCGATATTCCGACGGCGAGCGTTCGGTGTTTTTACGAAACAGACGGCAGAAGTAGTTGCTGTCGACaaa                                                                                   <  6:47687/122‑1 (MQ=255)
                                                                                           cATCTTATTCTGGAGTAGTCGGTTTTTCTGTCAGTTGGCTGTGGTACTGCCGGCGATATTCCGACGGCGAGCGTTCGGTGTTTTTACGAAACAGACGGCAGAAGTAGTTGCTGTCGACAAAGCCGCAGGCGTGCGCCACTTCTTTTAcc                                                       <  5:381202/149‑1 (MQ=255)
                                                                                             tCTTATTCTGGAGTAGTCGGTTTTTCTGTCAGTTGGCTGTGGTACTGCCGGCGATATTCCGACGGCGAGCGTTCGGTGTTTTTACGAAACAGACGGCAGAAGTAGTTGCTGTCGACAAAGCCGCAGGCGTGCGCCACTTCTTTTACCtt                                                     >  2:93327/1‑149 (MQ=255)
                                                                                                 aTTCTGGAGTAGTCGGTTTTTCTGTCAGTTGGCTGTGGTACTGCCGGCGATATTCCGACGGCGAGCGTTCGGTGTTTTTACGAAACAGACGGCAGAAGTAGTTGCTGTCGACAAAGCCGCAGGCGTGCGCCACTTCTTTTACCTTc                                                    <  6:302542/146‑1 (MQ=255)
                                                                                                 aTTCTGGAGTAGTCGGTTTTTCTGTCAGTTGGCTGTGGTACTGCCGGCGATATTCCGACGGCGAGCGTTCGGTGTTTTTACGAAACAGACGGCAGAAGTAGTTGCTGTCGACAAAGCCGCAGGCGTGCGCCACTTCTTTTACCTTCAg                                                  <  2:201308/148‑1 (MQ=255)
                                                                                                  ttCTGGAGTAGTCGGTTTTTCTGTCAGTTGGCTGTGGTACTGCCGGCGATATTCCGACGGCGAGCGTTCGGTGTTTTTACGAAACAGACGGCAGAAGTAGTTGCTGTCGAc                                                                                      <  6:310187/111‑1 (MQ=255)
                                                                                                  ttCTGGAGTAGTCGGTTTTTCTGTCAGTTGGCTGTGGTACTGCCGGCGATATTCCGACGGCGAGCGTTCGGTGTTTTTACGAAACAGACGGCAGAAGTAGTTGCTGTCGAc                                                                                      >  5:310187/1‑111 (MQ=255)
                                                                                                  ttCTGGAGTAGTCGGTTTTTCTGTCAGTTGGCTGTGGTACTGCCGGCGATATTCCGACGGCGAGCGTTCGGTGTTTTTACGAAACAGACGGCAGAAGTAGTTGCTGTCGACAAAGCCGCAGGCGTGCGCCACTACTTTTACCTTCAGGt                                                >  2:39172/1‑149 (MQ=255)
                                                                                                      ggAGTAGTCGGTTTTTCTGTCAGTTGGCTGTGGTACTGCCGGCGATATTCCGACGGCGAGCGTTCGGTTTTTTTACGAAACAGACGGCAGAAGTAGTTGCTGTCGACAAAGCCGCAGTCGTGCCCCACATCTATTACCTTCAGGTCa                                              >  1:177965/1‑147 (MQ=255)
                                                                                                       gAGTAGTCGGTTTTTCTGTCAGTTGGCTGTGGTACTGCCGGCGATATTCCGACGGCGAGCGTTCGGTGTTTTTACGAAACAGACGGCAGAAGTAGTTGCTGTCGACAAAGCCGCAGGCGTGCGCCACTTCTTTTACCTTCAGGTCATAg                                           >  8:128212/1‑149 (MQ=255)
                                                                                                       gAGTAGTCGGTTTTTCTGTCAGTTGGCTGTGGTACTGCCGGCGATATTCCGACGGCGAGCGTTCGGTGTTTTTACGAAACAGACGGCAGAAGTAGTTGCTGTCGACAAAGCCGCAGGCGTGCGCCACTACTTTTACCTTCAGGTCa                                              <  1:39172/146‑1 (MQ=255)
                                                                                                         gtagtCGGTTTTTCTGTCAGTTGGCTGTGGTACTGCCGGCGATATTCCGACGGCGAGCGTTCGGTGTTTTTACGAAACAGACGGCAGAAGTAGTTGCTGTCGACAAAGCCGCAGGCGTGCGCCACTTCTTTTACCTTCAGGTCATAGcc                                         <  6:32021/149‑1 (MQ=255)
                                                                                                            gtCGGTTTTTCTGTCAGTTGGCTGTGGTACTGCCGGCGATATTCCGACGGCGAGCGTTCGGTGTTTTTACGAAACAGACGGCAGAAGTAGTTGTTGTCGACAAAGCCGCAGGCGTGCGCCACTTCTTTTACCTTCAGGTCATAGCCttt                                      <  2:287615/149‑1 (MQ=255)
                                                                                                            gtCGGTTTTTCTGTCAGTTGGCTGTGGTACTGCCGGCGATATTCCGACGGCGAGCGTTCGGTGTTTTTACGAAACAGACGGCAGAAGTAGTTGCTGTCGACAAAGCCGCAGGCGTGCGCCACTTCTTTTACCTTCAGGTCATAGCCttt                                      <  4:188012/149‑1 (MQ=255)
                                                                                                             tCGGTTTTTCTGTCAGTTGGCTGTGGTACTGCCGGCGATATTCCGACGGCGAGCGTTCGGTGTTTTTACGAAACAGACGGCAGAAGTAGTTGCTGTCGACAAAGCCGCAGGCGTGCGCCACTTCTTTTACCTTCAGGTCATAGCCttt                                      >  2:381678/1‑148 (MQ=255)
                                                                                                             tCGGTTTTTCTGTCAGTTGGCTGTGGTACTGCCGGCGATATTCCGACGGCGAGCGTTCGGTGTTTTTACGAAACAGACGGCAGAAGTAGTTGCTGTCGACAAAGCCGCAGGCGTGCGCCACTTCTTTTACCTTCAGGTCATAGCCttt                                      >  1:287531/1‑148 (MQ=255)
                                                                                                             tCGGTTTTTCTGTCAGTTGGCTGTGGTACTGCCGGCGATATTCCGACGGCGAGCGTTCGGTGTTTTTACGAAACAGACGGCAGAAGTAGTTGCTGTCGACAAAGCCGCAGGCGTGCGCCACTTCTTTTACCTTCAGGTCATAGCCttt                                      >  7:139847/1‑148 (MQ=255)
                                                                                                             tCGGTTTTTCTGTCAGTTGGCTGTGGTACTGCCGGCGATATTCCGACGGCGAGCGTTCGGTGTTTTTACGAAACAGACGGCAGAAGTAGTTGCTGTCGACAAAGCCGCAGGCGTGCGCCACTTCTTTTACCTTCAGGTCATAGCCttt                                      >  7:326126/1‑148 (MQ=255)
                                                                                                             tCGGTTTTTCTGTCAGTTGGCTGTGGTACTGCCGGCGATATTCCGACGGCGAGCGTTCGGTGTTTTTACGAAACAGACGGCAGAAGTAGTTGCTGTCGACAAAGCCGCAGGCGTGCGCCACTTCTTTTACCTTCAGGTCATAGCCTTTc                                     <  6:13182/149‑1 (MQ=255)
                                                                                                             tCGGTTTTTCTGTCAGTTGGCTGTGGTACTGCCGGCGATATTCCGACGGCGAGCGTTCGGTGTTTTTACGAAACAGACGGCAGAAGTAGTTGCTGTCGACAAAGCCGCAGGCGTGCGCCACTTCTTTTACCTTCAGGTCATAGCCTTTc                                     <  5:287405/149‑1 (MQ=255)
                                                                                                              cGGTTTTTCTGTCAGTTGGCTGTGGTACTGCCGGCGATATTCCGACGGCGAGCGTTCGGTGTTTTTACGAAACAGACGGCAGAAGTAGTTGCTGTCGACAAAGCCGCAGGCGTGCGCCACTTCTTTTACCTTCAGGTCATAGCCTTTc                                     >  7:172079/1‑148 (MQ=255)
                                                                                                                   tttCTGTCAGTTGGCTGTGGTACTGCCGGCGATATTCCGACGGCGAGCGTTCGGTGTTTTTACGAAACAGACGGCAGAAGTAGTTGCTGTCGACAAAGCCGCAGGCGTGCGCCACTTCTTTTACCTTCAGGTCATAGCCTTTCAATAAc                               >  4:178150/1‑149 (MQ=255)
                                                                                                                    ttCTGTCAGTTGGCTGTGGTACTGCCGGCGATATTCCGACGGCGAGCGTTCGGTGTTTTTACGAAACAGACGGCAGAAGTAGTTGCTGTCGACAAAGCCGCAGGCGTGCGCCACTTCTTTTACCTTCAGGTCATAGc                                          >  2:319103/1‑137 (MQ=255)
                                                                                                                    ttCTGTCAGTTGGCTGTGGTACTGCCGGCGATATTCCGACGGCGAGCGTTCGGTGTTTTTACGAAACAGACGGCAGAAGTAGTTGCTGTCGACAAAGCCGCAGGCGTGCGCCACTTCTTTTACCTTCAGGTCATAGc                                          <  1:319103/137‑1 (MQ=255)
                                                                                                                          cAGTTGGCTGTGGTACTGCCGGCGATATTCCGACGGCGAGCGTTCGGTGTTTTTACGAAACAGACGGCAGAAGTAGTTGCTGTCGACAAAGCCGCAGGCGTGCGCCACTTCTTTTACCTTCAGGTCATAGCCTTTCAATAACGTCTTAg                        >  3:208231/1‑149 (MQ=255)
                                                                                                                            gTTGGCTGTGGTACTGCCGGCGATATTCCGACGGCGAGCGTTCGGTGTTTTTACGAAACAGACGGCAGAAGTAGTTGCTGTCGACAAAGCCGCAGGCGTGCGCCACTTCTTTTACCTTCAGGTCATAGCCTTTCAATAACGTCTTAGCg                      >  2:181669/1‑149 (MQ=255)
                                                                                                                                 cTGTGGTACTGCCGGCGATATTCCGACGGCGAGCGTTCGGTGTTTTTACGAAACAGACGGCAGAagtagt                                                                                                <  3:391/70‑1 (MQ=255)
                                                                                                                                 cTGTGGTACTGCCGGCGATATTCCGACGGCGAGCGTTCGGTGTTTTTACGAAACAGACGGCAGAagtagt                                                                                                >  4:391/1‑70 (MQ=255)
                                                                                                                                  tgtgGTACTGCCGGCGATATTCCGACGGCGAGCGTTCGGTGTTTTTACGAAACAGACGGCAGAAGTAGTTGCTGTCGACAAAGCCGCAGGCGTGCGCCACTTCTTTTACCTTCAGGTCATAGCCTTTCAATAACGTCTTAGCGTGCTc                 >  7:11476/1‑148 (MQ=255)
                                                                                                                                    tgGTACTGCCGGCGATATTCCGACGGCGAGCGTTCGGTGTTTTTACGAAACAGACGGCAGAAGTAGTTGCTGTCGACAAAGCCGCAGGCGTGCGCCACTTCTTTTACCTTCAGGTCATAGCCTTTCAATAACGTCTTAGCGTGCTCCAg              >  2:47062/1‑149 (MQ=255)
                                                                                                                                    tgGTACTGCCGGCGATATTCCGACGGCGAGCGTTCGGTGTTTTTACGAAACAGACGGCAGAAGTAGTTGCTGTCGACAAAGCCGCAGGCGTGCGCCACTTCTTTTACCTTCAGGTCATAGCCTTTCAATAACGTCTTAGCGTGCTCCAg              <  3:292961/149‑1 (MQ=255)
                                                                                                                                         cTGCCGGCGATATTCCGACGGCGAGCGTTCGGTGTTTTTACGAAACAGACGGCAGAAGTAGTTGCTGTCGACAAAGCCGCAGGCGTGCGCCACTTCTTTTACCTTCAGGTCATAGCCTTTCAATAACGTCTTAGCGTGCTCCAGTcgc          <  3:40293/148‑1 (MQ=255)
                                                                                                                                         cTGCCGGCGATATTCCGACGGCGAGCGTTCGGTGTTTTTACGAAACAGACGGCAGAAGTAGTTGCTGTCGACAAAGCCGCAGGCGTGCGCCACTTCTTTTACCTTCAGGTCATAGCCTTTCAATAACGTCTTAGCGTGCTCCAGTcgc          <  2:29760/148‑1 (MQ=255)
                                                                                                                                         cTGCCGGCGATATTCCGACGGCGAGCGTTCGGTGTTTTTACGAAACAGACGGCAGAAGTAGTTGCTGGCGACAAAGCCGCAGGCGTGCGCCACTTCTTTTACCTTCAGGTCATAGCCTTTCAATAACGTCTTAGCGTGCTCCAGTcgc          <  2:18082/148‑1 (MQ=255)
                                                                                                                                            ccGGCGATATTCCGACGGCGAGCGTTCGGTGTTTTTACGAAACAGACGGCAGAAGTAGTTGCTGTCGACAAAGCCGCAGGCGTGCGCCACTTCTTTTACCTTCAGGTCATAGCCTTTCAATAACGTCTTAGCGTGCTCCAGTCGCgtg       <  7:30030/148‑1 (MQ=255)
                                                                                                                                            ccGGCGATATTCCGACGGCGAGCGTTCGGTGTTTTTACGAAACAGACGGCAGAAGTAGTTGCTGTCGACAAAGCCGCAGGCGTGCGCCACTTCTTTTACCTTCAGGTCATAGCCTTTCAATAACGTCTTAGCGTGCTCCAGTCGCgt        <  6:84269/147‑1 (MQ=255)
                                                                                                                                             cGGCTATATTCCGACGGCGAGGGTTCGGTGTTTTTAAGAATCAGACGGCAGAAGTAGTTGCTGTCGACAAAGCCGCAGGCGTGCGCCACTTCTTTTACCTTCAGGTCATAGCCTTTCAATAACGTCTTAGCGTGCTCCAGTCGCGtgtg     <  2:177965/149‑1 (MQ=255)
                                                                                                                                                 gATATTCCGACGGCGAGCGTTCGGTGTTTTTACGAAACAGACGGCAGAAGTAGTTGCTGTCGACAAAGCCGCAGGCGTGCGCCACTTCTTTTACCTTCAGGTCATAGCCTTTCAATAACGTCTTAGCGTGCTCCAGTCGCGTGTGAtt  >  2:47994/1‑148 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
CTATGTGCAATGCAGCATCTACGCATTCCTGCGTGAGCCGCAGCGGGTAATGAGCGCGGTGCGCAAAGTGCTTTCTGCGCCGCAGCAAACCCATCTTATTCTGGAGTAGTCGGTTTTTCTGTCAGTTGGCTGTGGTACTGCCGGCGATATTCCGACGGCGAGCGTTCGGTGTTTTTACGAAACAGACGGCAGAAGTAGTTGCTGTCGACAAAGCCGCAGGCGTGCGCCACTTCTTTTACCTTCAGGTCATAGCCTTTCAATAACGTCTTAGCGTGCTCCAGTCGCGTGTGATT  >  NC_000913/4147371‑4147663

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