Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,715,258 A→G 5.3% A21A (GCT→GCC) ‡ yfiE ← putative DNA‑binding transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,715,2580AG 5.3% 332.3 / 5.9 95A21A (GCT→GCC) ‡yfiEputative DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (44/46);  new base G (3/2);  total (47/48)
Fisher's exact test for biased strand distribution p-value = 6.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.10e-01

TTCATAAATGTGCGGCAGCAACTTTTTTCCTTCGCGGGTAAGGCACATCCGTCGGCCAATTTTCTCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATCTAGAATATAAAAACATTCATTTTTTTAAATGTTCCGTGTCGGGTACTG  >  NC_000913/2715111‑2715397
                                                                                                                                                   |                                                                                                                                            
ttCATAAATGTGCGGCAGCAACTTTTTTCCTTCGCGGGTAAGCCACATCCGTCGGCCAATTTTCTCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAg                                                                                                                                             >  3:1162177/1‑149 (MQ=255)
   aTAAATGTGCGGCAGCAACTTTTTTCCTTCGCGGGTAAGGCACATCCGTCGGCCAATTTTCTCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTc                                                                                                                                         >  4:500899/1‑149 (MQ=255)
     aaaTGTGCGGCAGCAACTTTTTTCCTTCGCGGGTAAGGCACATCCGTCGGCCAATTTTCTCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGc                                                                                                                                       <  1:1266170/149‑1 (MQ=255)
          tgCGGCAGCAACTTTTTTCCTTCGCGGGTAAGGCACATCCGTCGGCCAATTTTCTCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCaaaaa                                                                                                                                  <  2:87697/149‑1 (MQ=255)
              gcagcaACTTTTTTCCTTCGCGGGTAAGGCACATCCGTCGGCCAATTTTCTCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGaa                                                                                                                              >  1:1154778/1‑149 (MQ=255)
                   aaCTTTTTTCCTTCGCGGGTAAGGCACATCCGTCGGCCAATTTTCTCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCctc                                                                                                                         >  2:1522676/1‑149 (MQ=255)
                     cTTTTTTCCTTCGCGGGTAAGGCACATCCGTCGGCCAATTTTCTCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCtt                                                                                                                       <  2:118665/149‑1 (MQ=255)
                      ttttttCCTTCGCGGGTAAGGCACATCCGTCGGCCAATTTTCTCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTAGAGCATAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTc                                                                                                                      >  4:1231778/1‑149 (MQ=255)
                      ttttttCCTTCGCGGGTAAGGCACATCCGTCGGCCAATTTTCTCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAAc                                                                                                                             <  3:216788/142‑1 (MQ=255)
                      ttttttCCTTCGCGGGTAAGGCACATCCGTCGGCCAATTTTCTCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAAc                                                                                                                             >  4:216788/1‑142 (MQ=255)
                      ttttttCCTTCGCGGGTAAGGCACATCCGTCGGCCAATTTTCTCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTc                                                                                                                      >  3:83140/1‑149 (MQ=255)
                           tCCTTCGCGGGTAAGGCACATCCGTCGGCCAATTTTCTCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGaa                                                                                                                                              <  2:1082534/121‑1 (MQ=255)
                           tCCTTCGCGGGTAAGGCACATCCGTCGGCCAATTTTCTCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGaa                                                                                                                                              >  1:1082534/1‑121 (MQ=255)
                              ttCGCGGGTAAGGCACATCCGTCGGCCAATTTTCTCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACgg                                                                                                              <  1:1522676/149‑1 (MQ=255)
                              ttCGCGGGTAAGGCACATCCGTCGGCCAATTTTCTCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACgg                                                                                                              >  3:1395702/1‑149 (MQ=255)
                               tCGCGGGTAAGGCACATCCGTCGGCCAATTTTCTCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGt                                                                                                             >  4:1889153/1‑149 (MQ=255)
                                      tAAGGCACATCCGTCGGCCAATTTTCTCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGc                                                                                                      <  4:1162177/149‑1 (MQ=255)
                                      tAAGGCACATCCGTCGGCCAATTTTCTCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGc                                                                                                      <  4:83140/149‑1 (MQ=255)
                                      tAAGGCACATCCGTCGGCCAATTTTCTCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGc                                                                                                      >  2:516403/1‑149 (MQ=255)
                                               tCCGTCGGCCAATTTTCTCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATaaaa                                                                                             <  3:1157338/149‑1 (MQ=255)
                                                  gTCGGCCAATTTTCTCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACgg                                                                                          <  3:500899/149‑1 (MQ=255)
                                                       ccAATTTTCTCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGgc                                                                                         <  4:699365/145‑1 (MQ=255)
                                                       ccAATTTTCTCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGgc                                                                                         >  3:699365/1‑145 (MQ=255)
                                                       ccAATTTTCTCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAg                                                                                     >  3:1364915/1‑149 (MQ=255)
                                                        cAATTTTCTCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGaa                                                                                                                              <  3:358962/107‑1 (MQ=255)
                                                        cAATTTTCTCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGaa                                                                                                                              >  4:28179/1‑107 (MQ=255)
                                                        cAATTTTCTCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGaa                                                                                                                              <  3:28179/107‑1 (MQ=255)
                                                        cAATTTTCTCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGaa                                                                                                                              >  4:358962/1‑107 (MQ=255)
                                                        cAATTTTCTCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGa                                                                                    >  2:595038/1‑149 (MQ=255)
                                                         aaTTTTCTCAAATAACTGGACTGAAAACTCCTGCTCATGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACGCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGAt                                                                                   <  3:1231778/149‑1 (MQ=255)
                                                         aaTTTTCTCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGAt                                                                                   <  2:1154778/149‑1 (MQ=255)
                                                          aTTTTCTCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATc                                                                                  <  4:1364915/149‑1 (MQ=255)
                                                           ttttCTCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATcc                                                                                 <  4:1395702/149‑1 (MQ=255)
                                                           ttttCTCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATcc                                                                                 <  1:516403/149‑1 (MQ=255)
                                                            tttCTCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATaaaa                                                                                             >  1:81131/1‑136 (MQ=255)
                                                            tttCTCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATaaaa                                                                                             <  2:81131/136‑1 (MQ=255)
                                                               ctcAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATaa                                                                             <  1:595038/149‑1 (MQ=255)
                                                                  aaaTAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATc                                                                                  <  1:1850983/141‑1 (MQ=255)
                                                                  aaaTAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATc                                                                                  >  2:1850983/1‑141 (MQ=255)
                                                                  aaaTAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAAccc                                                                          >  1:887059/1‑149 (MQ=255)
                                                                         tGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCGCAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCt                                                                   >  1:1247359/1‑149 (MQ=255)
                                                                         tGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCt                                                                   >  1:1776446/1‑149 (MQ=255)
                                                                         tGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGCTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAAc                                                                     >  1:1272915/1‑147 (MQ=255)
                                                                         tGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCCCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAAc                                                                     <  2:1272915/147‑1 (MQ=255)
                                                                            aCTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAACAAAACGGCGCAGATCCATAACCCCAAACCTATc                                                                >  1:1246897/1‑149 (MQ=255)
                                                                              tGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAAc                                                                                                                             <  2:13429/86‑1 (MQ=255)
                                                                              tGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAAc                                                                                                                             >  1:13429/1‑86 (MQ=255)
                                                                              tGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGtaa                                                                                                  >  1:594832/1‑113 (MQ=255)
                                                                              tGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGtaa                                                                                                  <  2:594832/113‑1 (MQ=255)
                                                                              tGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATcc                                                                                 <  4:1109539/130‑1 (MQ=255)
                                                                              tGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATcc                                                                                 >  3:1109539/1‑130 (MQ=255)
                                                                               gAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTaa                                                                                                        <  2:1300810/106‑1 (MQ=255)
                                                                               gAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTaa                                                                                                        >  1:1300810/1‑106 (MQ=255)
                                                                                 aaaCTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCGCAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGaaaa                                                           <  2:1247359/149‑1 (MQ=255)
                                                                                  aaCTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAat                                                          >  3:1638396/1‑149 (MQ=255)
                                                                                     tCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATc                                                       <  4:1565054/149‑1 (MQ=255)
                                                                                     tCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATc                                                       >  3:1565054/1‑149 (MQ=255)
                                                                                        tGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAAcc                                                                                                                            >  4:1478972/1‑77 (MQ=255)
                                                                                        tGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAAcc                                                                                                                            <  3:1478972/77‑1 (MQ=255)
                                                                                               gctgctGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATCTAGaa                                             <  3:1889153/149‑1 (MQ=255)
                                                                                               gctgctGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAACAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATCTAGaa                                             <  2:1246897/149‑1 (MQ=255)
                                                                                                 tgctgAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATCTAGAata                                           >  3:1780642/1‑149 (MQ=255)
                                                                                                  gctgAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTaa                                                                                                        >  1:1327848/1‑87 (MQ=255)
                                                                                                  gctgAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTaa                                                                                                        <  2:1327848/87‑1 (MQ=255)
                                                                                                   ctgAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATc                                                       <  1:1854712/135‑1 (MQ=255)
                                                                                                   ctgAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATc                                                       >  2:1854712/1‑135 (MQ=255)
                                                                                                     gAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATc                                                                <  2:1536868/124‑1 (MQ=255)
                                                                                                     gAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATCTAGAATATaaa                                       >  3:226263/1‑149 (MQ=255)
                                                                                                     gAATATGAAAAGTCACCGTCGATTGTGTACAGCACAACTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCGAAACCTATc                                                                >  1:1536868/1‑124 (MQ=255)
                                                                                                      aaTATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTaa                                                                                                        >  1:1151429/1‑83 (MQ=255)
                                                                                                      aaTATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTaa                                                                                                        <  2:1151429/83‑1 (MQ=255)
                                                                                                          tGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAGCTTCGCaaaaat                                                                                                                                 >  2:506780/1‑53 (MQ=25)
                                                                                                          tGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAGCTTCGCaaaaat                                                                                                                                 <  1:506780/54‑2 (MQ=25)
                                                                                                           gAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAGCTTCGCaaaaatt                                                                                                                                >  1:236404/1‑52 (MQ=25)
                                                                                                           gAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAGCTTCGCaaaaatt                                                                                                                                <  2:236404/54‑3 (MQ=25)
                                                                                                           gAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATCTa                                                <  4:853951/134‑1 (MQ=255)
                                                                                                           gAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATCTa                                                >  3:853951/1‑134 (MQ=255)
                                                                                                                gTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTccacca                                                                                                                 >  3:255009/1‑64 (MQ=255)
                                                                                                                gTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTccacca                                                                                                                 <  4:255009/64‑1 (MQ=255)
                                                                                                                   aCCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATCTAGAATATAAAAACATTCAtttttt                         <  4:1780642/149‑1 (MQ=255)
                                                                                                                      gTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCg                                                                                                     <  4:1553999/70‑1 (MQ=255)
                                                                                                                      gTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCg                                                                                                     >  3:1553999/1‑70 (MQ=255)
                                                                                                                             gtgtACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATCTAGAATATAAAAACATTCATTTTTTTAAATGTTcc               <  4:1638396/149‑1 (MQ=255)
                                                                                                                                   aGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATCTAGAATATAAAAACATTCATTTTTTTAAATGTTCCGTGTCg         >  3:1208628/1‑149 (MQ=255)
                                                                                                                                    gCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATCTAGAATATAAAAACATTCATTTTTTTAAATGTTCCGTGTCgg        <  2:1776446/149‑1 (MQ=255)
                                                                                                                                          tttttGCGAGCTTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATCTAGAATATAAAAACATTCATTTTTTTAAATGTTCCGTGTCGGGTACTg  >  3:315186/1‑149 (MQ=255)
                                                                                                                                          tttttGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATCTAGAATATAAAAACATTCAttttt                          <  2:1375995/125‑1 (MQ=255)
                                                                                                                                          tttttGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATCTAGAATATAAAAACATTCAttttt                          >  1:1375995/1‑125 (MQ=255)
                                                                                                                                          tttttGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATCTAGAATATAAAAACATTCAttttt                          <  2:1807078/125‑1 (MQ=255)
                                                                                                                                          tttttGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATCTAGAATATAAAAACATTCAttttt                          >  1:1807078/1‑125 (MQ=255)
                                                                                                                                          tttttGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATCTAGAATATAAAAACATTCATTTTTTTAAATGTTCCGTGTCGGGTACTg  <  4:226263/149‑1 (MQ=255)
                                                                                                                                           ttttGCGAGCTTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATCTAGaa                                             <  2:1706514/105‑1 (MQ=255)
                                                                                                                                           ttttGCGAGCTTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATCTAGaa                                             >  1:1706514/1‑105 (MQ=255)
                                                                                                                                            tttGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATCTAGAATATaaaaa                                     >  1:300672/1‑112 (MQ=255)
                                                                                                                                            tttGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATCTAGAATATaaaaa                                     <  2:300672/112‑1 (MQ=255)
                                                                                                                                                 gAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATCTAg                                               >  1:1801561/1‑97 (MQ=255)
                                                                                                                                                 gAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATCTAg                                               <  2:1801561/97‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                            
TTCATAAATGTGCGGCAGCAACTTTTTTCCTTCGCGGGTAAGGCACATCCGTCGGCCAATTTTCTCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATCTAGAATATAAAAACATTCATTTTTTTAAATGTTCCGTGTCGGGTACTG  >  NC_000913/2715111‑2715397

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: