Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,715,260 C→T 5.3% A21T (GCT→ACT) ‡ 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,2600CT 5.3% 328.1 / 4.2 94A21T (GCT→ACT) ‡yfiEputative DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (43/46);  new base T (3/2);  total (46/48)
Fisher's exact test for biased strand distribution p-value = 6.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATAAATGTGCGGCAGCAACTTTTTTCCTTCGCGGGTAAGGCACATCCGTCGGCCAATTTTCTCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATCTAGAATATAAAAACATTCATTTTTTTAAATGTTCCGTGTCGGGTACTG  >  NC_000913/2715114‑2715397
                                                                                                                                                  |                                                                                                                                          
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                                                                                                                             >  4:216788/1‑142 (MQ=255)
                   ttttttCCTTCGCGGGTAAGGCACATCCGTCGGCCAATTTTCTCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAAc                                                                                                                             <  3:216788/142‑1 (MQ=255)
                   ttttttCCTTCGCGGGTAAGGCACATCCGTCGGCCAATTTTCTCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTc                                                                                                                      >  3:83140/1‑149 (MQ=255)
                           ttCGCGGGTAAGGCACATCCGTCGGCCAATTTTCTCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACgg                                                                                                              <  1:1522676/149‑1 (MQ=255)
                           ttCGCGGGTAAGGCACATCCGTCGGCCAATTTTCTCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACgg                                                                                                              >  3:1395702/1‑149 (MQ=255)
                            tCGCGGGTAAGGCACATCCGTCGGCCAATTTTCTCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGt                                                                                                             >  4:1889153/1‑149 (MQ=255)
                                   tAAGGCACATCCGTCGGCCAATTTTCTCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGc                                                                                                      >  2:516403/1‑149 (MQ=255)
                                   tAAGGCACATCCGTCGGCCAATTTTCTCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGc                                                                                                      <  4:83140/149‑1 (MQ=255)
                                   tAAGGCACATCCGTCGGCCAATTTTCTCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGc                                                                                                      <  4:1162177/149‑1 (MQ=255)
                                            tCCGTCGGCCAATTTTCTCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATaaaa                                                                                             <  3:1157338/149‑1 (MQ=255)
                                               gTCGGCCAATTTTCTCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACgg                                                                                          <  3:500899/149‑1 (MQ=255)
                                                    ccAATTTTCTCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGgc                                                                                         >  3:699365/1‑145 (MQ=255)
                                                    ccAATTTTCTCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGgc                                                                                         <  4:699365/145‑1 (MQ=255)
                                                    ccAATTTTCTCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAg                                                                                     >  3:1364915/1‑149 (MQ=255)
                                                     cAATTTTCTCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGaa                                                                                                                              >  4:358962/1‑107 (MQ=255)
                                                     cAATTTTCTCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGaa                                                                                                                              <  3:358962/107‑1 (MQ=255)
                                                     cAATTTTCTCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGaa                                                                                                                              >  4:28179/1‑107 (MQ=255)
                                                     cAATTTTCTCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGaa                                                                                                                              <  3:28179/107‑1 (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                                                                                 <  1:516403/149‑1 (MQ=255)
                                                        ttttCTCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATcc                                                                                 <  4:1395702/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                                                                                 >  3:1109539/1‑130 (MQ=255)
                                                                           tGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATcc                                                                                 <  4:1109539/130‑1 (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                                                                                                                            <  3:1478972/77‑1 (MQ=255)
                                                                                     tGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAAcc                                                                                                                            >  4:1478972/1‑77 (MQ=255)
                                                                                            gctgctGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATCTAGaa                                             <  3:1889153/149‑1 (MQ=255)
                                                                                            gctgctGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAACAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATCTAGaa                                             <  2:1246897/149‑1 (MQ=255)
                                                                                              tgctgAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATCTAGAata                                           >  3:1780642/1‑149 (MQ=255)
                                                                                               gctgAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTaa                                                                                                        <  2:1327848/87‑1 (MQ=255)
                                                                                               gctgAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTaa                                                                                                        >  1:1327848/1‑87 (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                                                                                                        <  2:1151429/83‑1 (MQ=255)
                                                                                                   aaTATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTaa                                                                                                        >  1:1151429/1‑83 (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                                                >  3:853951/1‑134 (MQ=255)
                                                                                                        gAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATCTa                                                <  4:853951/134‑1 (MQ=255)
                                                                                                             gTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTccacca                                                                                                                 >  3:255009/1‑64 (MQ=255)
                                                                                                             gTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTccacca                                                                                                                 <  4:255009/64‑1 (MQ=255)
                                                                                                                aCCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATCTAGAATATAAAAACATTCAtttttt                         <  4:1780642/149‑1 (MQ=255)
                                                                                                                   gTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCg                                                                                                     >  3:1553999/1‑70 (MQ=255)
                                                                                                                   gTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCg                                                                                                     <  4:1553999/70‑1 (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                          <  2:1807078/125‑1 (MQ=255)
                                                                                                                                       tttttGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATCTAGAATATAAAAACATTCAttttt                          >  1:1807078/1‑125 (MQ=255)
                                                                                                                                       tttttGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATCTAGAATATAAAAACATTCAttttt                          >  1:1375995/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                                     <  2:300672/112‑1 (MQ=255)
                                                                                                                                         tttGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATCTAGAATATaaaaa                                     >  1:300672/1‑112 (MQ=255)
                                                                                                                                              gAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATCTAg                                               <  2:1801561/97‑1 (MQ=255)
                                                                                                                                              gAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATCTAg                                               >  1:1801561/1‑97 (MQ=255)
                                                                                                                                                  |                                                                                                                                          
ATAAATGTGCGGCAGCAACTTTTTTCCTTCGCGGGTAAGGCACATCCGTCGGCCAATTTTCTCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATCTAGAATATAAAAACATTCATTTTTTTAAATGTTCCGTGTCGGGTACTG  >  NC_000913/2715114‑2715397

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