Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,683,479 Δ1 bp 14.3% intergenic (‑32/+151) dctA ← / ← yhjK C4‑dicarboxylic acid, orotate and citrate transporter/cyclic‑di‑GMP phosphodiesterase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,683,4790T.14.3% 426.3 / 67.2 147intergenic (‑32/+151)dctA/yhjKC4‑dicarboxylic acid, orotate and citrate transporter/cyclic‑di‑GMP phosphodiesterase
Reads supporting (aligned to +/- strand):  ref base T (64/58);  new base . (12/9);  total (78/69)
Fisher's exact test for biased strand distribution p-value = 8.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.24e-01

Only 100 of 149 total aligned reads displayed.

TCGCCAAGCGGTTTCATTTGCTCGCCTATTTCAGGATAGAAATGGCCAAGGAGAATACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGT  >  NC_000913/3683332‑3683627
                                                                                                                                                   |                                                                                                                                                    
tCGCCAAGCGGTTTCATTTGCTCGCCTATTTCAGGATAGAAATGGCCAAGGAGAATACCAATGGCTATCGCTGTCAGGAACTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGt                                                                                                                                                      >  4:20308/1‑148 (MQ=255)
  gCCAAGCGGTTTCATTTGCTCGCCTATTTCAGGATAGAAATGGCCAAGGAGAATACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGTTCTTGTAAg                                                                                                                                                   >  7:376487/1‑149 (MQ=255)
  gCCAAGCGGTTTCATTTGCTCGCCTATTTCAGGATAGAAATGGCCAAGGAGAATACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTaa                                                                                                                                                    <  6:39716/148‑1 (MQ=255)
      aGCGGTTTCATTTGCTCGCCTATTTCAGGATAGAAATGGCCAAGGAGAATACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTa                                                                                                                                               <  4:65695/149‑1 (MQ=255)
       gCGGTTTCATTTGCTCGCCTATTTCAGGATAGAAATGGCCAAGGAGAATACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTAt                                                                                                                                              <  4:143458/149‑1 (MQ=255)
       gCGGTTTCATTTGCTCGCCTATTTCAGGATAGAAATGGCCAAGGAGAATACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGGGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTAt                                                                                                                                              <  5:232190/149‑1 (MQ=255)
             tCATTTGCTCGCCTATTTCAGGATAGAAATGGCCAAGGAGAATACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATggg                                                                                                                                           <  8:344017/146‑1 (MQ=255)
             tCATTTGCTCGCCTATTTCAGGATAGAAATGGCCAAGGAGAATACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATggg                                                                                                                                           >  7:344017/1‑146 (MQ=255)
               aTTTGCTCGCCTATTTCAGGATAGAAATGGCCAAGGAGAATACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTAcc                                                                                                                                      >  7:361087/1‑149 (MQ=255)
                     tCGCCTATTTCAGGATAGAAATGGCCAAGGAGAATACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTGCGGc                                                                                                                                >  6:99895/1‑149 (MQ=255)
                       gCCTATTTCAGGATAGAAATGGCCAAGGAGAATACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCATAGGTCTTGTAAGGTTATGGGGTACCTGCGGCCt                                                                                                                              >  8:320482/1‑149 (MQ=255)
                       gCCTATTTCAGGATAGAAATGGCCAAGGAGAATACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTGCGGcc                                                                                                                               <  1:403571/148‑1 (MQ=255)
                         cTATTTCAGGATAGAAATGGCCAAGGAGAATACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCATAGGTCTTGTAAGGTTATg                                                                                                                                             >  1:119573/1‑132 (MQ=255)
                         cTATTTCAGGATAGAAATGGCCAAGGAGAATACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCATAGGTCTTGTAAGGTTATg                                                                                                                                             <  2:119573/132‑1 (MQ=255)
                            tttcttgATAGAAATGGCCAAGGAGAATACCAATGGCTATCGCTGTCAGGACCTGATAGTAAAGGCTTTTAACCAGAGAGGTTTTCATGGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTaaaa                                                                                                                         <  2:158043/143‑1 (MQ=255)
                            tttCAGGATAGAAATGGCCAAGGAGAATACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGTTCTTGTAAGGTTATGGGGTACCTGCGGCCTTaa                                                                                                                           <  3:156912/147‑1 (MQ=255)
                            tttCAGGATAGAAATGGCCAAGGAGAATACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTaaaa                                                                                                                         >  1:157860/1‑149 (MQ=255)
                            tttCAGGATAGAAATGGCCAAGGAGAATACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTCTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACATGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTaa                                                                                                                           <  3:18603/147‑1 (MQ=255)
                               cAGGATAGAAATGGCCAAGGAGAATACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGAAGGTTATGGGGTAc                                                                                                                                       >  1:319558/1‑131 (MQ=255)
                                   aTAGAAATGGCCAAGGAGAATACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGTTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACAcc                                                                                                                  <  8:102432/149‑1 (MQ=255)
                                   aTAGAAATGGCCAAGGAGAATACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGTTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACAcc                                                                                                                  >  7:185669/1‑149 (MQ=255)
                                   aTAGAAATGGCCAAGGAGAATACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACAcc                                                                                                                  <  3:34384/149‑1 (MQ=255)
                                   aTAGAAATGGCCAAGGAGAATACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACAcc                                                                                                                  <  3:128160/149‑1 (MQ=255)
                                      gAAATGGCCAAGGAGAATACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGTTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCa                                                                                                                <  5:199629/148‑1 (MQ=255)
                                       aaaTGGCCAAGGAGAATACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTCTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTAACTGCGGCCTTAAAATAACACCCAGa                                                                                                              >  6:327305/1‑149 (MQ=255)
                                            gCCAAGGAGAATACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGCAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGacaac                                                                                                          <  3:283817/148‑1 (MQ=255)
                                            gCCAAGGAGAATACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGAAGGTTATGGGGTACCTGCGTCCTTAAAATAACACCCAGACAACAt                                                                                                        <  8:129417/149‑1 (MQ=255)
                                              cAAGGAGAATACCAATGGCTATCGCTGTCAGGACCTGAAAGTATAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATc                                                                                                       >  6:112902/1‑149 (MQ=255)
                                              cAAGGAGAATACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGGGTCCTTTAGTAAAACCACATGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATc                                                                                                       <  8:268704/149‑1 (MQ=255)
                                                 ggAGAATACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATcaca                                                                                                    <  2:388291/149‑1 (MQ=255)
                                                   agaATACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGa                                                                                                  <  3:25192/149‑1 (MQ=255)
                                                     aaTACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGTTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGaaa                                                                                                <  8:185669/149‑1 (MQ=255)
                                                     aaTACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGaaa                                                                                                <  2:146874/149‑1 (MQ=255)
                                                       tACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCATAGGTCTTGTAAGGTTATggg                                                                                                                                           >  4:333329/1‑104 (MQ=255)
                                                       tACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCATAGGTCTTGTAAGGTTATggg                                                                                                                                           <  3:333329/104‑1 (MQ=255)
                                                           aaTGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGAAGGTTATGGGGTACCTGCGGCCTTaaaa                                                                                                                         >  6:215647/1‑117 (MQ=255)
                                                               gCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACATGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGa                                                                                      >  6:347514/1‑149 (MQ=255)
                                                               gCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGAAGGTTAAGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGa                                                                                      <  1:38927/148‑1 (MQ=255)
                                                                  aTCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCAGAGGGGGTCCTTTAGTAAAACCACAGGTCTGGTAAGGTTATGGGGTACCTGCGGCCTATAATTAACACCCAGACAACATCACAGAAATGTACCTGGATCa                                                                                   >  4:269221/1‑149 (MQ=255)
                                                                      cTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAAtt                                                                              <  8:367772/149‑1 (MQ=255)
                                                                       tGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAAtt                                                                              >  7:109019/1‑149 (MQ=255)
                                                                            ggACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGaaa                                                                          <  8:426599/148‑1 (MQ=255)
                                                                            ggACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGaaa                                                                          >  7:426599/1‑148 (MQ=255)
                                                                                 tGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACATGTCTTGTAAGGTTATGGGTTCCTTCGGCCTTAAAATAAAACCCCGACAAAATAACAGAAATTTTCCTGGATCATAATTTAAACAAAAAg                                                                   >  5:400717/1‑149 (MQ=255)
                                                                                   aaaGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAg                                                                   <  3:69793/148‑1 (MQ=255)
                                                                                     aGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGtt                                                                <  7:6415/149‑1 (MQ=255)
                                                                                      gTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTa                                                               >  1:121830/1‑148 (MQ=255)
                                                                                        aaaGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTAACTGCGGACTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTaa                                                              >  3:119358/1‑148 (MQ=255)
                                                                                          aGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAg                                                          <  1:371065/149‑1 (MQ=255)
                                                                                            gCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTCCCTGGATCATAATTGAAACAAAAAGGGTAAAAAGt                                                         >  2:297269/1‑149 (MQ=255)
                                                                                             cTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGtt                                                        <  8:109019/149‑1 (MQ=255)
                                                                                              ttttAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGttt                                                       >  8:2552/1‑149 (MQ=255)
                                                                                              ttttAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGttt                                                       >  7:142439/1‑148 (MQ=255)
                                                                                              atttAAACAGATAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACATGGATCAGAATAGAAACAAAAAGGTTAAAAAGatt                                                       >  4:290136/2‑148 (MQ=255)
                                                                                                ttAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGttt                                                       <  7:2552/147‑1 (MQ=255)
                                                                                                ttAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGCAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAAt                                                                               <  4:51271/123‑1 (MQ=255)
                                                                                                ttAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGCAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAAt                                                                               >  3:51271/1‑123 (MQ=255)
                                                                                                 tAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAg                                                    <  6:183566/148‑1 (MQ=255)
                                                                                                  aaaCAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACATGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTg                                                      >  5:130661/1‑146 (MQ=255)
                                                                                                  aaaCAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACATGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTg                                                      <  6:130661/146‑1 (MQ=255)
                                                                                                    aCAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACATGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTg                                                 >  4:222353/1‑149 (MQ=255)
                                                                                                      agagagGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAGAAAGTTTGAGCTGGc                                               >  2:60317/1‑149 (MQ=255)
                                                                                                       gagagGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCt                                              <  3:391082/149‑1 (MQ=255)
                                                                                                        agagGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCt                                              >  3:308749/1‑148 (MQ=255)
                                                                                                        agagGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTc                                             >  4:382080/1‑149 (MQ=255)
                                                                                                        agagGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCg                                            >  7:128412/1‑149 (MQ=255)
                                                                                                          agGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGTTCATAATTGAAACAAAAAGGGTAAAAAGTTTGAGCTGGCTcac                                           >  1:351691/1‑147 (MQ=255)
                                                                                                            gTTTTCATAGGGTGTCCTTTAGTAAAACCACAGTTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGca                                          >  7:136707/1‑148 (MQ=255)
                                                                                                            gTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGca                                          <  4:4800/148‑1 (MQ=255)
                                                                                                                 cATAGGGTGTCCTTTAGTAAAACCACAGTTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCAGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAAc                                     <  8:376487/148‑1 (MQ=255)
                                                                                                                 cATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAAc                                     >  6:406773/1‑147 (MQ=255)
                                                                                                                  aTAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAAc                                     >  3:38934/1‑147 (MQ=255)
                                                                                                                  aTAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCACACAAAcc                                    <  8:58797/148‑1 (MQ=255)
                                                                                                                   tAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCACACAAACCAg                                  <  8:219717/149‑1 (MQ=255)
                                                                                                                      ggTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGCTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCACACAAACCAGCAc                               >  2:386163/1‑149 (MQ=255)
                                                                                                                            cTTTAGTAAAACCACAGTTCTTGTAAGGATATGGGGTACCTGCGGACTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACttttt                          >  4:155485/1‑148 (MQ=255)
                                                                                                                             tttAGTAAAACCACATGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTa                         >  5:387957/1‑148 (MQ=255)
                                                                                                                             tttAGTAAAACCACAGTTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTa                         >  7:180655/1‑148 (MQ=255)
                                                                                                                               tAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAGAAAGTTTGAGCTGGCTCACACAAACCAGCACTTTTTaaa                       <  1:60317/148‑1 (MQ=255)
                                                                                                                                aGTAAAACCACAGTTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCACACAAACCAGCACTTTTTaaat                      <  7:297951/148‑2 (MQ=255)
                                                                                                                                 gTAAAACCACAGGTCTTGAAGGTTATGGGGTACCTGCGGCCTTAAAATAAAACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGttt                   >  4:304944/1‑149 (MQ=255)
                                                                                                                                  tAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGCGCTGGCTCAGACAAACCAGCACTTTTTaatgggt                   >  2:237454/1‑144 (MQ=255)
                                                                                                                                  tAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACACAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCACACAAACCAGCACTTTTTAAAGttt                   >  6:342833/1‑149 (MQ=255)
                                                                                                                                   aaaaCCACAGGTCTTGTAAGGTTATGGGGAACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGtttt                  >  6:17532/1‑149 (MQ=255)
                                                                                                                                    aaaCCACATGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTg                 <  3:222353/149‑1 (MQ=255)
                                                                                                                                    aaaCCACATGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCACACAAACCAGCACTTTTTAAAGTTTTg                 <  1:237454/149‑1 (MQ=255)
                                                                                                                                    aaaCCACAGGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTg                 >  3:219262/1‑149 (MQ=255)
                                                                                                                                    aaaCCACAGGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGGTAAAAAGTTTGAGCTGGCTCACACAAACCAGCACTTTTTAAAGtttt                  <  5:189349/148‑1 (MQ=255)
                                                                                                                                    aaaCCACAGGTCTTGTAAGGCTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCACACAAACCAGCACTTTTTAAAGTTTTg                 <  1:386163/149‑1 (MQ=255)
                                                                                                                                       ccACAGGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCACACAAACCAGCACTTTTTAAAGTTTTGTaa              <  2:41668/149‑1 (MQ=255)
                                                                                                                                         acaGGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGtt                    >  1:61250/1‑141 (MQ=255)
                                                                                                                                         acaGGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCACACAAACCAGCACTTTTTAAAGTTTTGTAAt             >  5:193260/1‑148 (MQ=255)
                                                                                                                                         acaGGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACAACCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCACACAAACCAGCACTTTTTAAAGTTTTGTAAt             >  7:84300/1‑148 (MQ=255)
                                                                                                                                         acaGGTCTTGAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCa           <  8:142439/149‑1 (MQ=255)
                                                                                                                                          caGGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAAt             >  7:124373/1‑147 (MQ=255)
                                                                                                                                             gTCTTGAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCAAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCa           <  3:252945/145‑1 (MQ=255)
                                                                                                                                             gTCTTGAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCAAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCa           >  4:252945/1‑145 (MQ=255)
                                                                                                                                               cTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTg      >  5:285724/1‑149 (MQ=255)
                                                                                                                                                 tGTAAGGTTATGGGGTACCTGCGGCCTTAAAATACCCCCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCACACCAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTggg    >  2:277445/1‑149 (MQ=255)
                                                                                                                                                   tAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCACACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGt  <  1:252147/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
TCGCCAAGCGGTTTCATTTGCTCGCCTATTTCAGGATAGAAATGGCCAAGGAGAATACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGT  >  NC_000913/3683332‑3683627

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: