Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,182,764 G→A 9.8% H22H (CAC→CAT)  potD ← spermidine/putrescine ABC transporter periplasmic binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,182,7640GA 9.8% 176.3 / 8.7 61H22H (CAC→CAT) potDspermidine/putrescine ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base G (30/25);  new base A (2/4);  total (32/29)
Fisher's exact test for biased strand distribution p-value = 4.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.16e-01

ATGTTTTCAGCTTGGCGTACATGGTTTCGTTCGACTCGTAAGTCGAATAGATAACCTTAATACCGGTTTCTTTGGTGAACTGTTCAAGCAGTCCTGGCGGCACGTACTCGGTCCAGTTGTAGAAATACAGCGTGTTGTTGTCATCGGCGTGAGCGGCGCTCATGCCCAGTGCCAGAGCACCCGCCGCGAGCAGGTGGCGTGACCATTTTTTCATTTAACGTCCCCTGTGTTACCTTTCGTTTTATCACGAGCAATAAGCTGGCTGGCAATAACCATCACCAGCGACAGCACC  >  NC_000913/1182616‑1182907
                                                                                                                                                    |                                                                                                                                               
aTGTTTTCAGCTTGGCGTACATGGTTTCGTTCGACTCGTAAGTCGAATAGATAACCTTAATACCGGTTTCTTTGGTGAACTGTTCAAGCAGTCCTGGCGGCACGTACTCGGTCCAGTTGTAGAAATACAGCGTGTTGTTGTCATCGGCg                                                                                                                                                 >  3:357163/1‑149 (MQ=255)
 tGTTTTCAGCTTGGCGTACATGGTTTCGTTCGACTCGTAAGTCGAATAGATAACCTTAATACCGGTTTCTTTGGTGAACTGTTCAAGCAGTCCTGGCGGCACGTACTCGGTCCAGTTGTAGAAATACAGCGTGTTGTTGTCATCGGCGt                                                                                                                                                >  1:658350/1‑149 (MQ=255)
     ttCAGCTTGGCGTACATGGTTTCGTTCGACTCGTAAGTCGAATAGATAACCTTAATACCGGTTTCTTTGGTGAACTGTTCAAGCAGTCCTGGCGGCACGTACTCGGTCCAGTTGTAGAAATACAGCGTGTTGTTGTCATCGGCGTGAgc                                                                                                                                            <  3:663054/149‑1 (MQ=255)
         gCTTGGCGTACATGGTTTCGTTCGACTCGTAAGTCGAATAGATAACCTTAATACCGGTTTCTTTGGTGAACTGTTCAAGCAGTCCTGGCGGCACGTACTCGGTCCAGTTGTAGAAATACAGCGTGTTGTTGTCATCGGCATGAgcgccg                                                                                                                                        >  3:667628/1‑146 (MQ=255)
              gCGTACATGGTTTCGTTCGACTCGTAAGTCGAATAGATAACCTTAATACCGGTTTCTTTGGTGAACTGTTCAAGCAGTCCTGGCGGCACGTACTCGGTCCAGTTGTAGAAATACAGCGTGTTGTTGTCATCGGCGTGAGCGGCGCTCAt                                                                                                                                   <  1:130889/149‑1 (MQ=255)
              gCGTACATGGTTTCGTTCGACTCGTAAGTCGAATAGATAACCTTAATACCGGTTTCTTTGGTGAACTGTTCAAGCAGTCCTGGCGGCACGTACTCGGTCCAGTTGTAGAAATACAGCGTGTTGTTGTCATCGGCGTGAGCGGCGCTCAt                                                                                                                                   <  2:115093/149‑1 (MQ=255)
                gTACATGGTTTCGTTCGACTCGTAAGTCGAATAGATAACCTTAATACCGGTTTCTTTGGTGAACTGTTCAAGCAGTCCTGGCGGCACGTACTCGGTCCAGTTGTAGAAATACAGCGTGTTGTTGTCATCGGCATGAGCGCCGCTCACGc                                                                                                                                 <  4:667628/149‑1 (MQ=255)
                 tACATGGTTTCGTTCGACTCGTAAGTCGAATAGATAACCTTAATACCGGTTTCTTTGGTGAACTGTTCAAGCAGTCCTGGCGGCACGTACTCGGTCCAGTTGTAGAAATACAGCGTGTTGTTGTCATCGGCATGAGCGCCGCTCACGcc                                                                                                                                <  3:124419/149‑1 (MQ=255)
                     tGGTTTCGTTCGACTCGTAAGTCGAATAGATAACCTTAATACCGGTTTCTTTGGTGAACTGTTCAAGCAGTCCTGGCGGCACGTACTCGGTCCAGTTGTAGAAATACAGCGTGTTGTTGTCATCGGCGTGAGCGGCGCTCATGCCCAGt                                                                                                                            <  2:271321/149‑1 (MQ=255)
                     tGGTTTCGTTCGACTCGTAAGTCGAATAGATAACCTTAATACCGGTTTCTTTGGTGAACTGTTCAAGCAGTCCTGGCGGCACGTACTCGGTCCAGTTGTAGAAATACAGCGTGTTGTTGTCATCGGCGTGAGCGGCGCTCATGCCCAGt                                                                                                                            <  1:142237/149‑1 (MQ=255)
                      ggTTTCGTTCGACTCGTAAGTCGAATAGATAACCTTAATACCGGTTTCTTTGGTGAACTGTTCAAGCAGTCCTGGCGGCACGTACTCGGTCCAGTTGTAGAAATACAGCGTGTTGTTGTCATCGGCGTGAGCGGCGCTCATGCCCAGTg                                                                                                                           <  3:340624/149‑1 (MQ=255)
                      ggTTTCGTTCGACTCGTAAGTCGAATAGATAACCTTAATACCGGTTTCTTTGGTGAACTGTTCAAGCAGTCCTGGCGGCACGTACTCGGTCCAGTTGTAGAAATACAGCGTGTTGTTGTCATCGGCATGAGCGCCGCTCACGccgatga                                                                                                                           <  4:92339/149‑6 (MQ=255)
                        tttCGTTCGACTCGTAAGTCGAATAGATAACCTTAATACCGGTTTCTTTGGTGAACTGTTCAAGCAGTCCTGGCGGCACGTACTCGGTCCAGTTGTAGAAATACAGCGTGTTGTTGTCATCGGCGTGAGCGGCGCTCATGCCCAGTGcc                                                                                                                         >  4:236074/1‑149 (MQ=255)
                        tttCGTTCGACTCGTAAGTCGAATAGATAACCTTAATACCGGTTTCTTTGGTGAACTGTTCAAGCAGTCCTGGCGGCACGTACTCGGTCCAGTTGTAGAAATACAGCGTGTTGTTGTCATCGGCGTGAGCGGCGCTCATGCCCAGTGcc                                                                                                                         <  1:643407/149‑1 (MQ=255)
                         ttcgttcgACTCGTAAGTCGAATAGATAACCTTAATACCGGGTTCTTTGGTGAACTGTTCAAGCAGTCCTGGCGGCACGTACTCGGTCCAGTTGTAGAAATACAGCGTGTTGTTGTCATCGGCGTGAGCGGCGCTCATGCCCAGTGCCa                                                                                                                        <  2:516959/149‑1 (MQ=255)
                            gttcgACTCGTAAGTCGAATAGATAACCTTAATACCGGTTTCTTTGGTGAACTGTTCAAGCAGTCCTGGCGGCACGTACTCGGTCCAGTTGTAGAAATACAGCGTGTTGTTGTCATCGGCGTGAGCGGCGCTCATGCCCAGTGCCagag                                                                                                                     >  4:99858/1‑149 (MQ=255)
                             ttcgACTCGTAAGTCGAATAGATAACCTTAATACCGGTTTCTTTGGTGAACTGTTCAAGCAGTCCTGGCGGCACGTACTCGGTCCAGTTGTAGAAATACAGCGTGTTGTTGTCATCGGCGTGAGCGGCGCTCATGCCCAGTGCCAGAGc                                                                                                                    >  2:68325/1‑149 (MQ=255)
                                gACTCGTAAGTCGAATAGATAACCTTAATACCGGTTTCTTTGGTGAACTGTTCAAGCAGTCCTGGCGGCACGTACTCGGTCCAGTTGTAGAAATACAGCGTGTTGTTGTCATCGGCGTGAGCGGCGCTCATGCCCAGTGCCAGAGCAcc                                                                                                                 >  1:236266/1‑149 (MQ=255)
                                 aCTCGTAAGTCGAATAGATAACCTTAATACCGGTTTCTTTGGTGAACTGTTCAAGCAGTCCTGGCGGCACGTACTCGGTCCAGTTGTAGAAATACAGCGTGTTGTTGTCATCGGCGTGAGCGGCGCTCATGCCCAGTGCCAGAGCAccc                                                                                                                <  1:113236/149‑1 (MQ=255)
                                 aCTCGTAAGTCGAATAGATAACCTTAATACCGGTTTCTTTGGTGAACTGTTCAAGCAGTCCTGGCGGCACGTACTCGGTCCAGTTGTAGAAATACAGCGTGTTGTTGTCATCGGCGTGAGCGGCGCTCATGCCCAGTGCCAGAGCAccc                                                                                                                >  4:440671/1‑149 (MQ=255)
                                   tCGTAAGTCGAATAGATAACCTTAATACCGGTTTCTTTGGTGAACTGTTCAAGCAGTCCTGGCGGCACGTACTCGGTCCAGTTGTAGAAATACAGCGTGTTGTTGTCATCGGCGTGAGCGGCGCTCATGCCCAGTGCCAGAGCACccgc                                                                                                              >  1:230762/1‑149 (MQ=255)
                                          tCGAATAGATAACCTTAATACCGGTTTCTTTGGTGAACTGTTCAAGCAGTCCTGGCGGCACGTACTCGGTCCAGTTGTAGAAATACAGCGTGTTGTTGTCATCGGCGTGAGCGGCGCTCATGCCCAGTGCCAGAGCACCCGCCGCGAGc                                                                                                       <  3:236074/149‑1 (MQ=255)
                                          tCGAATAGATAACCTTAATACCGGTTTCTTTGGTGAACTGTTCAAGCAGTCCTGGCGGCACGTACTCGGTCCAGTTGTAGAAATACAGCGTGTTGTTGTCATCGGCGTGAGCGGCGCTCATGCCCAGTGCCAGAGCACCCGCCGCGAGc                                                                                                       >  1:278773/1‑149 (MQ=255)
                                           cGAATAGATAACCTTAATACCGGTTTCTTTGGTGAACTGTTCAAGCAGTCCTGGCGGCACGTACTCGGTCCAGTTGTAGAAATACAGCGTGTTGTTGTCATCGGCGTGAGCGGCGCTCATGCCCAGTGCCAGAGCACCCGCCGCGAGCa                                                                                                      <  2:627413/149‑1 (MQ=255)
                                                aGATAACCTTAATACCGGTTTCTTTGGTGAACTGTTCAAGCAGTCCTGGCGGCACGTACTCGGTCCAGTTGTAGAAATACAGCGTGTTGTTGTCATCGGCGTGAGCGGCGCTCATGCCCAGTGCCAGAGCACCCGCCGCGAGCAGGTgg                                                                                                 >  1:127518/1‑149 (MQ=255)
                                                      ccTTAATACCGGTTTCTTTGGTGAACTGTTCAAGCAGTCCTGGCGGCACGTACTCGGTCCAGTTGTAGAAATACAGCGTGTTGTTGTCATCGGCGTGAGCGGCGCTCATGCCCAGTGCCAGAGCACCCGCCGCGAGCAGGTGGCGTGAc                                                                                           >  4:257860/1‑149 (MQ=255)
                                                       cgtAATACCGGTTTCTTTGGTGAACTGTTGAAGCAGTCCTGGCGGCACGTACTCGGTCCAGTTGTAGAAATACAGCGTGTTGTTGTCATCGGCGTGAGCGGCGCTCATGCCCAGTGCCAGAGCACCCGCCGCGAGCAGGTGGCGTGAcc                                                                                          >  3:584681/3‑149 (MQ=255)
                                                       cTTAATGCCGGGTTCTTTGGTGAACTGTTCAAGCAGTCCTGGCGGCACGTACTCGGTCCAGTTGTAGAAATACAGCGTGTTGTTGTCATCGGCGTGAGCGGCGCTCATGCCCAGTGCCaga                                                                                                                      >  2:340349/1‑121 (MQ=255)
                                                       cTTAATACCGGTTTCTTTGGTGAACTGTTCAAGCAGTCCTGGCGGCACGTACTCGGTCCAGTTGTAGAAATACAGCGTGTTGTTGTCATCGGCGTGAGCGGCGCTCATGCCCAGTGCCaga                                                                                                                      <  1:340349/121‑1 (MQ=255)
                                                            tACCGGTTTCTTTGGTGAACTGTTCAAGCAGTCCTGGCGGCACGTACTCGGTCCAGTTGTAGAAATACAGCGTGTTGTTGTCATCGGCGTGAGCGGCGCTCATGCCCAGTGCCAGAGCACCCGCCGCGAGCAGGTGGCGTGACCAtttt                                                                                     >  4:650944/1‑149 (MQ=255)
                                                                ggTTTCTTTGGTGAACTGTTCAAGCAGTCCTGGCGGCACGTACTCGGTCCAGTTGTAGAAATACAGCGTGTTGTTGTCATCGGCGTGAGCGGCGCTCATGCCCAGTGCCAGAGCACCCGCCGCGAGCAGGTGGCGTGACCATTTTTTCa                                                                                 >  3:249240/1‑149 (MQ=255)
                                                                      tttGGTGAACTGTTCAAGCAGTCCTGGCGGCACGTACTCGGTCCAGTTGTAGAAATACAGCGTGTTGTTGTCATCGGCGTGAGCGGCGCTCATGCCCAGTGCCAGAGCACCCGCCGCGAGCAGGTGGCGTGACCATTTTTTCATTTAAc                                                                           <  3:257860/149‑1 (MQ=255)
                                                                      tttGGTGAACTGTTCAAGCAGTCCTGGCGGCACGTACTCGGTCCAGTTGTAGAAATACAGCGTGTTGTTGTCATCGGCGTGAGCGGCGCTCATGCCCAGTGCCAGAGCACCCGCCGCGAGCAGGTGGCGTGACCATTTTTTCATTTAAc                                                                           <  3:440671/149‑1 (MQ=255)
                                                                      tttGGTGAACTGTTCAAGCAGTCCTGGCGGCACGTACTCGGTCCAGTTGTAGAAATACAGCGTGTTGTTGTCATCGGCATGAGCGCCGCTCACGccgatga                                                                                                                           <  3:227941/101‑6 (MQ=38)
                                                                      tttGGTGAACTGTTCAAGCAGTCCTGGCGGCACGTACTCGGTCCAGTTGTAGAAATACAGCGTGTTGTTGTCATCGGCATGAGCGCCGCTCACGccgatga                                                                                                                           >  4:227941/1‑96 (MQ=38)
                                                                        tGGTGAACTGTTCAAGCAGTCCTGGCGGCACGTACTCGGTCCAGTTGTAGAAATACAGCGTGTTGTTGTCATCGGCGTGAGCGGCGCTCATGCCCAGTGCCAGAGCACCCGCCGCGAGCAGGTGGCGTGACCATTTTTTCATTTAACGt                                                                         <  3:650944/149‑1 (MQ=255)
                                                                         ggTGAACTGTTCAAGCAGTCCTGGCGGCACGTACTCGGTCCAGTTGTAGAAATACAGCGTGTTGTTGTCATCGGCGTGAGCGGCGCTCATGCCCAGTGCCAGAGCACCCGCCGCGAGCAGGTGGCGTGACCATTTTTTCATTTAACGTc                                                                        <  2:127518/149‑1 (MQ=255)
                                                                           tGAACTGTTCAAGCAGTCCTGGCGGCACGTACTCGGTCCAGTTGTAGAAATACAGCGTGTTGTTGTCATCGGCGTGAGCGGCGCTCATGCCCAGTGCCAGAGCACCCGCCGCGAGCAGGTGGCGTGACCATTTTTTCATTTAACGTccc                                                                      >  3:593635/1‑149 (MQ=255)
                                                                           tGAACTGTTCAAGCAGTCCTGGCGGCACGTACTCGGTCCAGTTGTAGAAATACAGCGTGTTGTTGTCATCGGCGTGAGCGGCGCTCATGCCCAGTGCCAGAGCACCCGCCGCGAGCAGGTGGCGTGACCATTTTTTCATTTAACGTccc                                                                      <  4:249240/149‑1 (MQ=255)
                                                                                   tCAAGCAGTCCTGGCGGCACGTACTCGGTCCAGTTGTAGAAATACAGCGTGTTGTTGTCATCGGCGTGAGCGGCGCTCATGCCCAGTGCCAGAGCACCCGCCGCGAGCAGGTGGCGTGACCATTTTTTCATTTAACGTCCCCTGTGTTa                                                              <  3:635595/149‑1 (MQ=255)
                                                                                      aGCAGTCCTGGCGGCACGTACTCGGTCCAGTTGTAGAAATACAGCGTGTTGTTGTCATCGGCGTGAGCGGCGCTCATGCCCAGTGCCAGAGCACCCGCCGCGAGCAGGTGGCGTGACCATTTTTTCATTTAACGTCCCCTGTGTTACCt                                                           >  4:463616/1‑149 (MQ=255)
                                                                                       gCAGTCCTGGCGGCACGTACTCGGTCCAGTTGTAGAAATACAGCGTGTTGTTGTCATCGGCGTGAGCGGCGCTCATGCCCAGTGCCAGAGCACCCGCCGCGAGCAGGTGGCGTGACCATTTTTTCATTTAACGTCCCCTGTGTTACCtt                                                          >  3:480755/1‑149 (MQ=255)
                                                                                         aGTCCTGGCGGCACGTACTCGGTCCAGTTGTAGAAATACAGCGTGTTGTTGTCATCGGCGTGAGCGGCGCTCATGCCCAGTGCCAGAGCACCCGCCGCGAGCAGGTGGCGTGACCATTTTTTCATTTAACGTCCCCTGTGTTACCTTTc                                                        <  4:476384/149‑1 (MQ=255)
                                                                                          gTCCTGGCGGCACGTACTCGGTCCAGTTGTAGAAATACAGCGTGTTGTTGTCATCGGCGTGAGCGGCGCTCATGCCCAGTGCCAGAGCACCCGCCGCGAGCAGGTGGCGTGACCATTTTTTCATTTAACGTCCCCTGTGTTACCTTTCg                                                       <  2:236266/149‑1 (MQ=255)
                                                                                             cTGGCGGCACGTCCTCGGTCCAGCTGTAGAAATACCGCGTGTTGTTGTCATCGGCGTGGGCGGCGCTCATGCCCCGTGCCAGAGCCCCCGCCGCGAGCCGGTGGCGTGACCATTTTTTCATTTAACGTCCCCTGTGTTACCTTTCGttt                                                    >  4:636392/1‑149 (MQ=255)
                                                                                              tGGCGGCACGTACTCGGTCCAGTTGTAGAAATACAGCGTGTTGTTGTCATCGGCGTGAGCGGCGCTCATGCCCAGTGCCAGAGCACCCGCCGCGAGCAGGTGGCGTGACCATTTTTTCATTTAACGTCCCCTGTGTTACCTTTCGtttt                                                   >  1:148997/1‑149 (MQ=255)
                                                                                              tGGCGGCACGTACTCGGTCCAGTTGTAGAAATACAGCGTGTTGTTGTCATCGGCGTGAGCGGCGCTCATGCCCAGTGCCAGAGCACCCGCCGCGAGCAGGTGGCGTGACCATTTTTTCATTTAACGTCCCCTGTGTTACCTTTCGtttt                                                   <  3:284340/149‑1 (MQ=255)
                                                                                              tGGCGGCACGTACTCGGTCCAGTTGTAGAAATACAGCGTGTTGTTGTCATCGGCGTGAGCGGCGCTCATGCCCAGTGCCAGAGCACCCGCCGCGAGCAGGTGGCGTGACCATTTTTTCATTTAACGTCCCCTGTGTTACCTTTCGttgt                                                   <  3:567707/149‑3 (MQ=255)
                                                                                                   gcACGAACTCGGTCCAGTTGTAGAAATACAGCGAGTTGTTGTCATCGGCGTGAGCGGCGCTCATGCCCAGTGCCAGAGCACCAGCCGCGAGCAGGTGGCGTGACCATTTTTTCATTTAACGTCCCCTGTGTTACCTTTCGTTTTATCAc                                              >  1:200009/1‑149 (MQ=255)
                                                                                                     aCGTACTCGGTCCAGTTGTAGAAATACAGCGTGTTGTTGTCATCGGCGTGAGCGGCGCTCATGCCCAGTGCCAGAGCACCCGCCGCGAGCAGGTGGCGTGACCATTTTTTCATTTAACGTCCCCTGTGTTACCTTTCGTTTTATCACGa                                            >  3:370645/1‑149 (MQ=255)
                                                                                                     aCGTACTCGGTCCAGTTGTAGAAATACAGCGTGTTGTTGTCATCGGCGTGAGCGGCGCTCATGCCCAGTGCCAGAGCACCCGCCGCGAGCAGGTGGCGTGACCATTTTTTCATTTAACGTCCCCTGTGTTACCTTTCGTTTTATCACGa                                            >  2:103019/1‑149 (MQ=255)
                                                                                                     aCGTACTCGGTCCAGTTGTAGAAATACAGCGTGTTGTTGTCATCGGCGTGAGCGGCGCTCATGCCCAGTGCCAGAGCACCCGCCGCGAGCAGGTGGCGTGACCATTTTTTCATTTAACGTCCCCTGTGTTACCTTTCGTTTTATCACGa                                            >  3:255912/1‑149 (MQ=255)
                                                                                                         aCTCGGTCCAGTTGTAGAAATACAGCGTGTTGTTGTCATCGGCGTGAGCGGCGCTCATGCCCAGTGCCAGAGCACCCGCCGCGAGCAGGTGGCGTGACCATTTTTTCATTTAACGTCCCCTGTGTTACCTTTCGTTTTATCACGa                                            <  2:149124/145‑1 (MQ=255)
                                                                                                         aCTCGGTCCAGTTGTAGAAATACAGCGTGTTGTTGTCATCGGCGTGAGCGGCGCTCATGCCCAGTGCCAGAGCACCCGCCGCGAGCAGGTGGCGTGACCATTTTTTCATTTAACGTCCCCTGTGTTACCTTTCGTTTTATCACGa                                            >  1:149124/1‑145 (MQ=255)
                                                                                                         aCTCGGTCCAGTTGTAGAAATACAGCGTGTTGTTGTCATCGGCGTGAGCGGCGCTCATGCCCAGTGCCAGAGCACCCGCCGCGAGCAGGTGGCGTGACCATTTTTTCATTTAACGTCCCCTGTGTTACCTTTCGTTTTATCACGAGCaa                                        >  2:316235/1‑149 (MQ=255)
                                                                                                           tCGGTCCAGTTGTAGAAATACAGCGTGTTGTTGTCATCGGCGTGAGCGGCGCTCATGCCCAGTGCCAGAGCACCCGCCGCGAGCAGGTGGCGTGACCATTTTTTCATTTAACGTCCCCTGTGTTACCTTTCGTTTTATCACGAGc                                          >  3:642325/1‑145 (MQ=255)
                                                                                                           tCGGTCCAGTTGTAGAAATACAGCGTGTTGTTGTCATCGGCGTGAGCGGCGCTCATGCCCAGTGCCAGAGCACCCGCCGCGAGCAGGTGGCGTGACCATTTTTTCATTTAACGTCCCCTGTGTTACCTTTCGTTTTATCACGAGc                                          <  4:642325/145‑1 (MQ=255)
                                                                                                             ggTCCAGTTGTAGAAATACAGCGTGTTGTTGTCATCGGCGTGAGCGGCGCTCATGCCCAGTGCCAGAGCACCCGCCGCGAGCAGGTGGCGTGACCATTTTTTCATTTAACGTCCCCTGTGTTACCTTTCGTTTTATCACGAGCAATAAg                                    <  2:275391/149‑1 (MQ=255)
                                                                                                                              aCAGCGTGTTGTTGTCATCGGCGTGAGCGGCGCTCATGCCCAGTGCCAGAGCACCCGCCGCGAGCAGGTGGCGTGACCATTTTTTCATTTAACGTCCCCTGTGTTACCTTTCGTTTTATCACGAGCAATAAGCTGGCTGGCAATAACCa                   >  1:135863/1‑149 (MQ=255)
                                                                                                                                    tgttgttgtCATCGGCGTGAGCGGCGCTCATGCCCAGTGCCAGAGCACCCGCCGCGAGCAGGTGGCGTGACCATTTTTTCATTTAACGTCCCCTGTGTTACCTTTCGTTTTATCACGAGCAATAAGCTGGCTGGCAATAACCATCACCa             >  2:375297/1‑149 (MQ=255)
                                                                                                                                         ttgtCATCGGCGTGAGCGGCGCTCATGCCCAGTGCCAGAGCACCCGCCGCGAGCAGGTGGCGTGACCATTTTTTCATTTAACGTCCCCTGTGTTACCTTTCGTTTTATCACGAGCAATAAGCTGGCTGGCAATAACCATCACCAGCGAc        >  2:644229/1‑149 (MQ=255)
                                                                                                                                               tcgtcgTCAGCGGCGCTCATGCCCAGTGCCAGAGCACCCGCCGCGAGCAGGTGGCGTGACCATTTTTTCATTTAACGTCCCCTGTGTTACCTTTCGTTTTATCACGAGCAATAAGCTGGCTGGCAATATCCATCACCAGCGACAGCAcc  >  4:485859/1‑149 (MQ=255)
                                                                                                                                               tCGGCGTGAGCGGCGCTCATGCCCAGTGCCAGAGCACCCGCCGCGAGCAGGTGGCGTGACCATTTTTTCATTTAACGTCCCCTGTGTTACCTTTCGTTTTATCACGAGCAATAAGCTGGCTGGCAATAACCATCACCAGCGACAGCAcc  >  1:609532/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                               
ATGTTTTCAGCTTGGCGTACATGGTTTCGTTCGACTCGTAAGTCGAATAGATAACCTTAATACCGGTTTCTTTGGTGAACTGTTCAAGCAGTCCTGGCGGCACGTACTCGGTCCAGTTGTAGAAATACAGCGTGTTGTTGTCATCGGCGTGAGCGGCGCTCATGCCCAGTGCCAGAGCACCCGCCGCGAGCAGGTGGCGTGACCATTTTTTCATTTAACGTCCCCTGTGTTACCTTTCGTTTTATCACGAGCAATAAGCTGGCTGGCAATAACCATCACCAGCGACAGCACC  >  NC_000913/1182616‑1182907

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑