Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,183,882 A→T 5.4% I188N (ATC→AAC)  potB ← spermidine/putrescine ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,183,8820AT 5.4% 235.2 / 2.9 74I188N (ATC→AAC) potBspermidine/putrescine ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base A (38/32);  new base T (2/2);  total (40/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.06e-01

ATATTAAGGAACTGGACCTTGATGACGTTACCGATCAGCAGGTTTTTCGCACCGCCCATCAGGTCGGATACATAGAACAGGCCCATCGCTGGCAGCATCACCAGCAGACATCCGGCAATAATTCCCGGCATCGTCAGTGGAATAATGATACGGATAAAAGTCTGTAACTTGCTGGCACCGAGATCGCGCGCCGCCTCAAGCAACGGCTTATCCAGTTTTTCGATACTGGAGTACAGCGGCATCACCATAAACGGCAGCAGAATGTAAACCAGACCGATAATCACCGCACTGGGCGT  >  NC_000913/1183735‑1184030
                                                                                                                                                   |                                                                                                                                                    
attatatggAACTGGACCTTGATGACGTTACGGATCAGCAGGTTGTTCGCACCGCCCATCAGGTCGGATACATAGTACAGGCCCATCGCTGGCAGCATCACCAGCAGACATCCGGCAATAATTCCCGGCATCGTCAGTGGAATAATGAt                                                                                                                                                     <  8:158906/142‑1 (MQ=255)
atatTAAGGAACTGGACCTTGATGACGTTACCGATCAGCAGGTTTTTCGCACCGCCCATCAGGTCGGATACATAGAACAGGCCCATCGCTGGCAGCATCACCAGCAGACATCCGGCAATAATTCCCGGCATCGTCAGTGGAATAATGAt                                                                                                                                                     <  2:247418/149‑1 (MQ=255)
  atTAAGGAACTGGACCTTGATGACGTTACCGATCAGCAGGTTTTTCGCACCGCCCATCAGGTCGGATACATAGAACAGGCCCATCGCTGGCAGCATCACCAGCAGACATCCGGCAATAATTCCCGGCATCGTCAGTGGAATAATGATAc                                                                                                                                                   >  1:187079/1‑149 (MQ=255)
   ttAAGGAACTGGACCTTGATGACGTTACCGATCAGCAGGTTTTTCGCACCGCCCATCAGGTCGGATACATAGAACAGGCCCATCGCTGGCAGCATCACCAGCAGACATCCGGCAATAATTCCCGGCATCGTCAGTGGAATAATGATACg                                                                                                                                                  <  5:52843/149‑1 (MQ=255)
    tAAGGAACTGGACCTTGATGACGTTACCGATCAGCAGGTTTTTCGCACCGCCCATCAGGTCGGATACATAGAACAGGCCCATCGCTGGCAGCATCACCAGCAGACATCCGGCAATAATTCCCGGCATCGTCAGTGGAATAATGATACg                                                                                                                                                  >  1:388333/1‑148 (MQ=255)
          aCTGGACCTTGATGACGTTACCGATCAGCAGGTTTTTCGCACCGCCCATCAGGTCGGATACATAGAACAGGCCCATCGCTGGCAGCATCACCAGCAGACATCCGGCAATAATTCCCGGCATCGTCAGTGGAATAATGATACGGATaaaa                                                                                                                                           <  1:300251/149‑1 (MQ=255)
              gACCTTGATGACGTTACCGATCAGCAGGTTTTTCGCACCGACCATCAGGTCGGATACATAGAACAGGACCATCGCTGGCAGCATCACCAGCAGACATCCGGCAATAATTCCCGGCATCGTCAGTGGAATAATGATACGGATAAAAGTCt                                                                                                                                       >  7:44900/1‑149 (MQ=255)
               aCCTTGATGACGTTACCGATCAGCAGGTTTTTCGCACCGCCCATCAGGTCGGATACATAGAACAGGCCCATCGCTGGCAGCATCACCAGCAGACATCCGGCAATAATTCCCGGCATCGTCAGTGGAATAATGATACGGATAAAAGTCTg                                                                                                                                      <  2:187079/149‑1 (MQ=255)
                 cTTGATGACGTTACCGATCAGCAGGTTTTTCGCACCGCCCATCAGGTCGGATACATAGAACAGGCCCATCGCTGGCAGCATCACCAGCAGACATCCGGCAATAATTCCCGGCATCGTCAGTGGAATAATGATACGGATAAAAGTCTGTa                                                                                                                                    <  8:319832/149‑1 (MQ=255)
                 cTTGATGACGTTACCGATCAGCAGGTTTTTCGCACCGCCCATCAGGTCGGATACATAGAAAAGGCCCATCGCTGGCAGCATCACCAGCAGACATCCGGCAATAATTCCCGGCATAGTCAGTGGAATAATGTTACGGATAAACGTCTGt                                                                                                                                     >  6:86645/1‑148 (MQ=255)
                            tACCGATCAGCAGGTTTTTCGCACCGCCCATCAGGTCGGATACATAGAACAGGCCCATCGCTGGCAGCATCACCAGCAGACATCCGGCAATAATTCCCGGCATCGTCAGTGGAATAATGATACGGATAAAAGTCTGTAACTTGCTGGCa                                                                                                                         >  8:273599/1‑149 (MQ=255)
                               cGATCAGCAGGTTTTTCGCACCGCCCATCAGGTCGGATACATAGAACAGGCCCATCGCTGGCAGCATCACCAGCAGACATCCGGCAATAATTCCCGGCATCGTCAGTGGAATAATGATACGGATAAAAGTCTGTAACTTGCTGGCAcc                                                                                                                       >  1:188700/1‑148 (MQ=255)
                               cGATCAGCAGGTTTTTCGCACCGCCCATCAGGTCGGATACATAGAACAGGCCCATCGCTGGAAGCATCACCAGCAGACATCCGGCAATAATTCCCGGCATCGTCAGGGGAATAATGATACGGATAAAAGTCTGTAACTTGCTGGcccc                                                                                                                       >  1:280107/1‑148 (MQ=255)
                                 aTCAGCAGGTTTTTCGCACCGCCCATCAGGTCGGATACATAGAACAGGCCCATCGCTGGCAGCATCACCAGCAGACATCCGGCAATAATTCCCGGCATCGTCAGTGGAATAATGATACGGATAAAAGTCTGTAACTTGCTGGCAcc                                                                                                                       >  5:357199/1‑146 (MQ=255)
                                     gcagGTTTTTCGCACCGCCCATCAGGTCGGATACATAGAACAGGCCCATCGCTGGCAGCATCACCAGCAGACATCCGGCAATAATTCCCGGCATCGTCAGTGGAATAATGATACGGATAAAAGTCTGTACCTTGCTGGCACCGAGATcg                                                                                                                >  4:12659/1‑149 (MQ=255)
                                     gcagGTTTTTCGCACCGCCCATCAGGTCGGATACATAGAACAGGCCCATCGCTGGCAGCATCACCAGCAGACATCCGGCAATAATTCCCGGCATCGTCAGTGGAATAATGATACGGATAAAAGTCTGTAACTTGCTGGCACCGAGATcg                                                                                                                >  4:351846/1‑149 (MQ=255)
                                      cagGTTTTTCGCACCGCCCATCAGGTCGGATACATAGAACAGGCCCATCGCTGGCAGCATCACCAGCAGACATCCGGCAATAATTCCCGGCATCGTCAGTGGAATAATGATACGGATAAAAGTCTGTAACTTGCTGGCAcc                                                                                                                       >  1:281770/1‑141 (MQ=255)
                                      cagGTTTTTCGCACCGCCCATCAGGTCGGATACATAGAACAGGCCCATCGCTGGCAGCATCACCAGCAGACATCCGGCAATAATTCCCGGCATCGTCAGTGGAATAATGATACGGATAAAAGTCTGTAACTTGCTGGCAcc                                                                                                                       <  2:281770/141‑1 (MQ=255)
                                      cagGTTTTTCGCACCGCCCATCAGGTCGGATACATAGAACAGGCCCATCGCTGGCAGCATCACCAGCAGACATCCGGCAATAATTCCCGGCATCGTCAGTGGAATAATGATACGGATAAAAGTCTGTAACTTGCTGGCACCGAGATcgc                                                                                                               <  2:388333/149‑1 (MQ=255)
                                            tttCGCACCGCCCATCAGGTCGGATACATAGAACAGGCCCATCGCTGGCAGCATCACCAGCAGACATCCGGCAATAATTCCCGGCATCGTCAGTGGAATAATGATACGGATAAAAGTCTGTAACTTGCTGGCACCGAGATCGCGcgcc                                                                                                          <  5:268370/148‑1 (MQ=255)
                                                gCACCGCCCATCAGGTCGGATACATAGAACAGGCCCATCGCTGGCAGCATCACCAGCAGACATCCGGCAATAATTCCCGGCATCGTCAGTGGAATAATGATACGGATAAAAGTCTGTAACTTGCTGGCACCGAGATCGCGCGCCGCCt                                                                                                      >  5:252330/1‑148 (MQ=255)
                                                  aCCGCCCATCAGGTCGGATACATAGAACAGGCCCATCGCGGGCAGCATCACCAGCAGACATCCGGCAATAATTCCCGGCAGCGACAGTGGAATAATGATACGGATAAAAGTCTGTAACTTGCTGGCACCGAGATCGCGCGCCGCCTc                                                                                                     >  5:352659/1‑147 (MQ=255)
                                                      cccATCAGGTCGGATACATAGAACAGGCCCATCGCTGGCAGCATCACCAGCAGACATCCGGCAATAATTCCCGGCATCGTCAGTGGAATAATGATACGGATAAAAGTCTGTAACTTGCTGGCACCGAGATCGCGCGCCGCCTCAAGCaa                                                                                               <  5:345542/149‑1 (MQ=255)
                                                         aTCAGGTCGGATACATAGAACAGGCCCATCGCTGGCAGCATCACCAGCAGACATCCGGCAATAATTCCCGGCATCGTCAGTGGAATAATGATACGGATAAAAGTCTGTAACTTGCTGGCACCGAGATCGCGCGCCGCCTCAAGCAACgg                                                                                            >  6:201802/1‑149 (MQ=255)
                                                          tCAGGTCGGATACATAGAACAGGCCCATCGCTGGCAGCATCACCAGCAGACATCCGGCAATAATTCCCGGCATCGTCAGTGGAATAATGATACGGATAAAAGTCTGTAACTTGCTGGCACCGAGATCGCGCGCCGCCTCAAGCAACGGc                                                                                           >  5:331884/1‑149 (MQ=255)
                                                          tCAGGTCGGATACATAGAACAGGCCCATCGCTGGCAGCATCACCAGCAGACATCCGGCAATAATTCCCGGCATCGTCAGTGGAATAATGATACGGATAAAAGTCTGTAACTTGCTGGCACCGAGATCGCGCGCCGCCTCAAGCAACGGc                                                                                           >  3:275243/1‑149 (MQ=255)
                                                                 ggATACATAGAACAGGCCCATCGCTTGCAGCATCACCAGCAGACATCCGGCAATAATTCCCGGCATCGTCAGTGGAATAATGATACGGATAAAAGTCTGTAACTTGCTGGCACCGAGATCGCGCGCCGCCTCAAGCAACGGCTTATcc                                                                                     <  6:252330/148‑1 (MQ=255)
                                                                       aTAGAACAGGCCCATCGCTGGCAGCATCACCAGCAGACATCCGGCAATAATTCCCGGCATCGTCAGTGGAATAATGATACGGATAAAAGTCTGTAACTTGCTGGCACCGAGATCGCGCGCCGCCTCAAGCAACGGCTTATCCAGttttt                                                                              >  3:1104/1‑149 (MQ=255)
                                                                       aTAGAACAGGCCCATCGCTGGCAGCATCACCAGCAGACATCCGGCAATAATTCCCGGCATCGTCAGTGGAATAATGATACGGATAAAAGTCTGTAACTTGCTGGCACCGAGATCGCGCGCCGCCTCAAGCAACGGCTTATCCAGttttt                                                                              >  8:44307/1‑149 (MQ=255)
                                                                        tggAACAGGCCCATCGCTGGCAGCATCACCAGCAGACATCCGGCAATAATTCCCGGCATCGTCAGTGGAATAATGATACGGATAAAAGTCTGTAACTTGCTGGCACCGAGATCGCGCGCCGCCTCAAGCAACGGCTTATCCAGTTTTTc                                                                             >  8:287663/3‑149 (MQ=255)
                                                                        tAGAACAGTCCCATCGCTGGCAGCATCAACAGCAGACATCCGGCAATAATTACCGGCATCGTCAGTGGAATAATGATACGGATAAAAGTCTGTAACTTGCTGGCACCGAGATCGCGCGCCGCCTCAAGCAACGGCTTATCCAGTTTTTc                                                                             >  8:336919/1‑149 (MQ=255)
                                                                        tAGAACAGGCCCATCGCTGGCAGCATCACCAGCAGACATCCGGCAATAATTCCCGGCATCGTCAGTGGAATAATGATACGGATAAAAGTCTGTAACTTGCTGGCACCGAGATCGCGCGCCGCCTCAAGCAACGGCTTATCCAGTTTTTc                                                                             >  3:192838/1‑149 (MQ=255)
                                                                         aGAACAGGCCCATCGCTGGCAGCATCACCAGCAGACATCCGGCAATAATTCCCGGCATCGTCAGTGGAATAATGATACGGATAAAAGTCTGTAACTTg                                                                                                                               >  6:109820/1‑98 (MQ=255)
                                                                         aGAACAGGCCCATCGCTGGCAGCATCACCAGCAGACATCCGGCAATAATTCCCGGCATCGTCAGTGGAATAATGATACGGATAAAAGTCTGTAACTTg                                                                                                                               <  5:109820/98‑1 (MQ=255)
                                                                            aCAGGCCCATCGCTGGCAGCATCACCAGCAGACATCCGGCAATAATTCCCGGCATCGTCAGTTGAATAATGTTACGGATAAAAGTCTGTAACTTGCTGGCACCGAGATCGCTCGCCGCCTCAAGCAACGGCTTATCCAGTTTTTCGATa                                                                         <  5:86645/149‑1 (MQ=255)
                                                                               ggCCCATCGCTGGCAGCATCACCAGCAGACATCCGGCAATAATTCCCGGCATCGTCAGTGGAATAATGATACGGATAAAAGTCTGTAACTTGCTGGCACCGAGATCGCGCGCCGCCTCAAGCAACGGCTTATCCAGTTTTTCGATACTg                                                                      <  6:331884/149‑1 (MQ=255)
                                                                                gCCCATCGCTGGCAGCATCACCAGCAGACATCCGGCAATAATTCCCGGCATCGTCAGTGGAATAATGATACGGATAAAAGTCTGTAACTTGCTGGCACCGAGATCGCGCGCCGCCTCAAGCAACGGCTTATCCAGTTTTTCGATACTgg                                                                     <  6:334926/149‑1 (MQ=255)
                                                                                 cccATCGCTGGCAGCATCACCAGCAGACATCCGGCAATAATTCCCGGCATCGTCAGTGGAATAATGATACGGATAAAAGTCTGTAACTTGCTGGCACCGAGATCGCGCGCCGCCTCAAGCAACGGCTTATCCAGTTTTTCGATACTGGa                                                                    <  6:357199/149‑1 (MQ=255)
                                                                                  ccATCGCTGGCAGCATCACCAGCAGACATCCGGCAATAATTCCCGGCATCGTCAGTGGAATAATGATACGGATAAAAGTCTGTAACTTGCTGGCACCGAGATCGCGCGCCGCCTCAAGCAACGGCTTATCCAGTTTTTCGATAATGGAg                                                                   >  6:142922/1‑149 (MQ=255)
                                                                                  ccATCGCTGGCAGAATCACCAGCAGACATCCGGCAATAAATACCGGCATCGTCAGTGGAATAATGATACGGATAAAACTCTGTTACTTGCTGGCACCGAGATCGCGCGCCGCCTCAAGCAACGGCTTATAGAGTTTTTCGATACTGGAg                                                                   >  8:371018/1‑149 (MQ=255)
                                                                                    aTCGCTGGCAGCATCACCAGCAGACATCCGGCAATAATTCCCGGCATCGTCAGTGGAATAATGATACGGATAAAAGTCTGTAACTTGCTGGCACCGAGATCGCGCGCCGCCTCAAGCAACGGCTTATCCAGTTTTTCGATACTGGAGTa                                                                 <  5:73514/149‑1 (MQ=255)
                                                                                      cGCTGGCAGCATCACCAGCAGACATCCGGCAATAATTCCCGGCATCGTCAGTGGAATAATGATACGGATAAAAGTCTGTAACTTGCTGGCACCGAGATCGCGCGCCGCCTCAAGCAACGGCTTATCCAGTTTTTCGa                                                                           >  2:350213/1‑137 (MQ=255)
                                                                                      cGCTGGCAGCATCACCAGCAGACATCCGGCAATAATTCCCGGCATCGTCAGTGGAATAATGATACGGATAAAAGTCTGTAACTTGCTGGCACCGAGATCGCGCGCCGCCTCAAGCAACGGCTTATCCAGTTTTTCGa                                                                           <  1:350213/137‑1 (MQ=255)
                                                                                       gCTGGCAGCATCACCAGCAGACATCCGGCAATAATTCCCGGCATCGTCAGTGGAATAATGATACGGATAAAAGTCTGTAACTTGCTGGCACCGAGATCGCGCGCCGCCTCAAGCAACGGCTTATCCAGTTTTTCGATACTGGAGTACa                                                               <  7:273599/148‑1 (MQ=255)
                                                                                           gcagcaTCACCAGCAGACATCCGGCAATAATTCCCGGCATCGTCAGTGGAATAATGATACGGATAAAAGTCTGTAACTTGCTGGCACCGAGATCGCGCGCCGCCTCAAGCAACGGCTTATCCAGTTTTTCGATACTGGAGTACAGCgg                                                           >  7:38768/1‑148 (MQ=255)
                                                                                           gcagcaTCACCAGCAGACATCCGGCAATAATTCCCGGCATCGTCAGTGGAATAATGATACGGATAAAAGTCTGTAACTTGCTGGCACCGAGATCGCGCGCCGCCTCAAGCAACGGCTTATCCAGTTTTTCGATACTGGAGTACAGCGGc                                                          <  4:1104/149‑1 (MQ=255)
                                                                                                   aCCAGCAGACATCCGGCAATAATTCCCGGCATCGTCAGTGGAATAATGATACGGATAAAAGTCTGTAACTTGCTGGCACCGAGATCGCGCGCCGCCTCAAGCAACGGCTTATCCAGTTTTTCGATa                                                                         <  3:231541/126‑1 (MQ=255)
                                                                                                   aCCAGCAGACATCCGGCAATAATTCCCGGCATCGTCAGTGGAATAATGATACGGATAAAAGTCTGTAACTTGCTGGCACCGAGATCGCGCGCCGCCTCAAGCAACGGCTTATCCAGTTTTTCGATa                                                                         >  4:231541/1‑126 (MQ=255)
                                                                                                   aCCAGCAGACATCCGGCAATAATTCCCGGCATCGTCAGTGGAATAATGATACGGATAAAAGTCTGTAACTTGCTGGCACCGAGATCGCGCGCCGCCTCAAGCAACGGCTTATCCAGTTTTTCGATACTGGAGTACAGCGGCATCACCa                                                   >  7:126308/1‑148 (MQ=255)
                                                                                                   aCCAGCAGACATCCGGCAATAATTCCCGGCATCGTCAGTGGAATAATGATACGGATAAAAGTCTGTAACTTGCTGGCACCGAGATCGCGCGCCGCCTCAAGCAACGGCTTATCCAGTTTTTCGATACTGGAGTACAGCGGCATCACCAt                                                  <  3:351846/149‑1 (MQ=255)
                                                                                                       gcagACATCCGGCAATAATTCCCGGCATCGTCAGTGGAATAATGATACGGATAAAAGTCTGTAACTTGCTGGCACCGAGATAGCGCGCCGCCTCAAGCAACGGCTTATCAAGTTTTTCGATACTGGAGTACAGCGGCATCACCATAAAc                                              >  1:278907/1‑149 (MQ=255)
                                                                                                              tCCGGCAATAATTCCCGGCATCGTCAGTGGAATAATGATACGGATAAAAGTCTGTAACTTGCTGGCACCGAGATCGCGCGCCGCCTCAAGCAACGGCTTATCCAGTTTTTCGATACTGGAGTACAGCGGCATCACCATAAACGgtagca                                       <  2:313062/149‑1 (MQ=255)
                                                                                                                  gCAATAATTCCCGGCATCGTCAGTGGAATAATGATACGGATAAAAGTCTGTAACTTGCTGGCACCGAGATCGCGCGCCGCCTCAAGCAACGGCTTATCCAGTTTTTCGATACTGGAGTACAGCGGCATCACCATAAACGGCAGCAGAAt                                   <  4:275243/149‑1 (MQ=255)
                                                                                                                   cAATAATTCCCGGCATCGTCAGTGGAATAATGATACGGATAAAAGTCTGTAACTTGCTGGCACCGAGATCGCGCGCCGCCTCAAGCAACGGCTTATCCAGTTTTTCGATACTGGAGTACAGCGGCATCACCATAAACGGCAGCAGAATg                                  <  3:12659/149‑1 (MQ=255)
                                                                                                                   cAATAATTCCCGGCATCGTCAGTGGAATAATGATACGGATAAAAGTCTGTAACTTGCTGGCACCGAGATCGCGCGCCGCCTCAAGCAACGGCTTATCCAGTTTTTCGATACTGGAGTACAGCGGCATCACCATAAACGGCAGCAGAATg                                  <  5:201802/149‑1 (MQ=255)
                                                                                                                    aataatTCCCGGCATCGTCAGTGGAATAATGATACGGATAAAAGTCTGTAACTTGCTGGCACCGAGATCGCGCGCCGCCTCAAGCAACGGCTTATCCAGTTTTTCGATACTGGAGTACAGCGGCATCACCATAAACGGCAGCAGAATg                                  >  1:403562/1‑148 (MQ=255)
                                                                                                                    aataatTCCCGGCATCGTCAGTGGAATAATGATACGGATAAAAGTCTGTAACTTGCTGGCACCGAGATCGCGCGCCGCCTCAAGCAACGGCTTATCCAGTTTTTCGATACTGGAGTACAGCGGCATCACCATAAACGGCAGCAGAATg                                  >  3:107935/1‑148 (MQ=255)
                                                                                                                     ataatTCCCGGCATCGTCAGTGGAATAATGATACGGATAAAAGTCTGTAACTTGCTGGCACCGAGATCGCTCGACGCCTCAAGCACCGGCTTATCCAGTGTTTCGATACTGGAGTACAGCGGCATCACCATAAACGGCAGCAGAATGt                                 >  4:38467/1‑148 (MQ=255)
                                                                                                                        atTCCCGGCATCGTCAGTGGAATAATGATACGGATAAAAGTCTGTAACTTGCTGGCACCGAGATCGCGCGCCGCCTCAAGCAACGGCTTATCCAGTTTTTCGATACTGGAGTACAGCGGCATCACCATAAACGGCAGCAGAATGTAAAc                             <  8:126308/149‑1 (MQ=255)
                                                                                                                            ccGGCATCGTCAGTGGAATAATGTTACGGATAAAAGTCTGTAACTTGCTGGCACCGAGATCGCGCGCCGCCTCAAGCAACGGCTTATCCAGTTTTTCGATACTGGAGTACAGCGGCATCACCATAAACGGCAGCAGAATGTaa                               >  3:390164/1‑143 (MQ=255)
                                                                                                                            ccGGCATCGTCAGTGGAATAATGTTACGGATAAAAGTCTGTAACTTGCTGGCACCGAGATCGCGCGCCGCCTCAAGCAACGGCTTATCCAGTTTTTCGATACTGGAGTACAGCGGCATCACCATAAACGGCAGCAGAATGTaa                               <  4:390164/143‑1 (MQ=255)
                                                                                                                                cATCGTCAGTGGAATAATGATACGGATAAAAGTCTGTAACTTGCTGGCACCGAGATCGCGCGCCGCCTCAAGCAACGGCTTATCCAGTTTTTCGATACTGGAGTACAGCGGCATCACCATAAACGGCAGCAGAATGTAAACCAGACCGa                     <  7:287663/149‑1 (MQ=255)
                                                                                                                                  tCGTCAGTGGAATAATGATACGGATAAAAGTCTGTAACTTGCTGGCACCGAGATCGCGCGCCGCCTCAAGCAACGGCTTATCCAGTTTTTCGATACTGGAGTACAGCGGCATCACCATAAACGGCAGCAGAATGTAAACCAGACCGAt                    >  7:4379/1‑148 (MQ=255)
                                                                                                                                        gTGGAATAATGATACGGATAAAAGTCTGTAACTTGCTGGCACCGAGATCGCGCGCCGCCTCAAGCAACGGCTTATCCAGTTTTTCGATACTGGAGTACAGCGGCATCACCATAAACGGCAGCAGAATGTAAACCAGACCGATAATCAcc             >  7:2701/1‑149 (MQ=255)
                                                                                                                                          ggAATAATGATACGGATAAAAGTCTGTAACTTGCTGGCACCGAGATCGCGCGCCGCCTCAAGCAACGGCTTATCCAGTTTTTCGATACTGGAGTACAGCGGCATCACCATAAACGgcagca                                       >  2:26854/1‑121 (MQ=255)
                                                                                                                                          ggAATAATGATACGGATAAAAGTCTGTAACTTGCTGGCACCGAGATCGCGCGCCGCCTCAAGCAACGGCTTATCCAGTTTTTCGATACTGGAGTACAGCGGCATCACCATAAACGgcagca                                       <  1:26854/121‑1 (MQ=255)
                                                                                                                                           gAATAATGATACGGATAAAAGTCTGTAACTTGCTGGCACCGAGATCGCGCGCCGCCTCAAGCAACGGCTTATCCAGTTTTTCGATACTGGAGTACAGCg                                                            >  2:7942/1‑99 (MQ=255)
                                                                                                                                           gAATAATGATACGGATAAAAGTCTGTAACTTGCTGGCACCGAGATCGCGCGCCGCCTCAAGCAACGGCTTATCCAGTTTTTCGATACTGGAGTACAGCg                                                            <  1:7942/99‑1 (MQ=255)
                                                                                                                                           gAATAATGATACGGATAAAAGTCTGTAACTTGCTGGCACCGAGATCGCGCGCCGCCTCAAGCAACGGCTTATCCAGTTTTTCGATACTGGAGTACAGCGGCATCACCATAAACGGCAGCAGAATGTAAACCAGACCGATAATCACCGCa          >  7:141247/1‑149 (MQ=255)
                                                                                                                                           gAATAATGATACGGATAAAAGTCTGTAACTTGCTGGCACCGAGATCGCGCGCCGCCTCAAGCAACGGCTTATCCAGTTTTTCGATACTGGAGTACAGCGGCATCACCATAAACGGCAGCAGAATGTAAACCAGACCGATAATCACCGCa          >  7:135071/1‑149 (MQ=255)
                                                                                                                                                 tGATACGGATAAAAGTCTGTAACTTGCTGGCACCGAGATCGCGCGCCGCCTCAAGCAACGGCTTATCCAGTTTTTCGATACTGGAGTACAGCGGCATCACCATAAACGGCAGCAGAATGTAAACCAGACCGATAATCACCGCACTGGGc    <  2:403562/149‑1 (MQ=255)
                                                                                                                                                 tGATACGGATAAAAGTCTGTAACTTGCTGGCACCGAGATCGCGCGCCGCCTCAAGCAACGGCTTATCCAGTTTTTCGATACTGGAGTACAGCGGCATCACCATAAACGGCAGCAGAATGTAAACCAGACCGATAATCACCGCACTGGGc    >  3:301066/1‑149 (MQ=255)
                                                                                                                                                 tGATACGGATAAAAGTCTGTAACTTGCTGGCACCGAGATCGCGCGCCGCCTCAAGCAACGGCTTATCCAGTTTTTCGATACTGGAGTACAGCGGCATCACCATAAACGGCAGCAGAATGTAAACCAGACCGATAATCACCGCACTGGGc    <  7:336919/149‑1 (MQ=255)
                                                                                                                                                 tGATACGGATAAAAGTCTGTAACTTGCTGGCACCGAGATCGCGCGCCGCCTCAAGCAACGGCTTATCCAGTTTTTCGATACTGGAGGACAGCGGAATCACCATAAACGGCAGCAGAATGTAAACCTGACCGATAATCACCGCACTGGGc    <  6:352659/149‑1 (MQ=255)
                                                                                                                                                   aTACGGATAAAAGTCTGTAACTTGCTGGCACCGAGATCGCGCGCCGCCTCAAGCAACGGCTTATCCAGTTTTTCGATACTGGAGTACAGCGGCATCACCATAAACGGCAGCAGAATGTAAACCAGACCGATAATCACCGCACTGGGCGt  >  6:366862/1‑149 (MQ=255)
                                                                                                                                                   aTACGGATAAAAGTCTGTAACTTGCTGGCACCGAGATCGCGCGCCGCCTCAAGCAACGGCTTATCCAGTTTTTCGATACTGGAGTACAGCGGCATCACCATAAACGGCAGCAGAATGTAAACCAGACCGATAATCACCGCACTGGGCGt  <  8:38768/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
ATATTAAGGAACTGGACCTTGATGACGTTACCGATCAGCAGGTTTTTCGCACCGCCCATCAGGTCGGATACATAGAACAGGCCCATCGCTGGCAGCATCACCAGCAGACATCCGGCAATAATTCCCGGCATCGTCAGTGGAATAATGATACGGATAAAAGTCTGTAACTTGCTGGCACCGAGATCGCGCGCCGCCTCAAGCAACGGCTTATCCAGTTTTTCGATACTGGAGTACAGCGGCATCACCATAAACGGCAGCAGAATGTAAACCAGACCGATAATCACCGCACTGGGCGT  >  NC_000913/1183735‑1184030

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