Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 553,999 G→A 7.8% H147Y (CAC→TAC)  ppiB ← peptidyl‑prolyl cis‑trans isomerase B (rotamase B)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913553,9990GA 7.8% 233.9 / 7.3 77H147Y (CAC→TAC) ppiBpeptidyl‑prolyl cis‑trans isomerase B (rotamase B)
Reads supporting (aligned to +/- strand):  ref base G (32/39);  new base A (3/3);  total (35/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.85e-01

GGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCATACCGCTACGACCGGTTGCAACACCTTTGATTTTGTCTACCACGTCCATGCCGTCAACCACTTCAGCAAACACGCAGTAGCCCCAACCTTGCAGGCTTTCGCCAGAGAAGTTCAGGAAGTCGTTATCAACCACGTTGATGAAGAA  >  NC_000913/553851‑554146
                                                                                                                                                    |                                                                                                                                                   
ggCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTg                                                                                                                                                     <  5:311191/149‑1 (MQ=255)
 gCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTg                                                                                                                                                     >  6:371161/1‑148 (MQ=255)
 gCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGc                                                                                                                                                    <  2:210299/149‑1 (MQ=255)
  cTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCa                                                                                                                                                   <  8:36503/149‑1 (MQ=255)
  cTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCa                                                                                                                                                   <  8:123659/149‑1 (MQ=255)
     cccccGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCATa                                                                                                                                                 <  1:310963/148‑1 (MQ=255)
     cccccGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTACATAc                                                                                                                                                >  3:370012/1‑149 (MQ=255)
     cccccGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTACATAc                                                                                                                                                >  3:364293/1‑149 (MQ=255)
       cccGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCATACCg                                                                                                                                              >  8:198316/1‑149 (MQ=255)
          gCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCATACCGCTa                                                                                                                                           >  4:244053/1‑149 (MQ=255)
           cTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCATACCGCTAc                                                                                                                                          <  5:336279/149‑1 (MQ=255)
              aaaaaCGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGTTTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCATACCGCTACGAc                                                                                                                                       <  5:324282/149‑1 (MQ=255)
              aaaaaCGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCATACCGCTACGAc                                                                                                                                       <  6:83634/149‑1 (MQ=255)
               aaatCGCAGAACACCGGCTGGGATCGCCGGGTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTCACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCATACCGCTACGAcc                                                                                                                                      <  2:328612/149‑1 (MQ=255)
                  aCGCAGAAAACCGGAGGTGAGAGCCGGTTCTTGCACGCAGAGATGCAGATCTGCAATAAAGAGTTGCGGCCCGATTACTCGCTAACTGTCAGGATTTCAATTATAACGTCTTCTTTTGGCACGTCCTGGTGCATACCGCTACGACCGGt                                                                                                                                   <  3:282846/149‑1 (MQ=255)
                            ccGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCATACCGCTACGACCGGTTGCAACACCt                                                                                                                         >  5:45543/1‑149 (MQ=255)
                              ggcggTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCATACCGCTACGACCGGTTGCAACACCttt                                                                                                                       <  5:226390/149‑1 (MQ=255)
                                 ggTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCATACCGCTACGACCGGTTGCAACACCTTTGAt                                                                                                                    <  5:284137/149‑1 (MQ=255)
                                    gATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCATACCGCTACGACCGGTTGCAACACCTTTGAtttt                                                                                                                 >  2:15664/1‑149 (MQ=255)
                                    gATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTACATACCGCTACGACCGGTTGc                                                                                                                                <  7:361050/134‑1 (MQ=255)
                                    gATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTACATACCGCTACGACCGGTTGc                                                                                                                                >  8:361050/1‑134 (MQ=255)
                                      tCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCATACCGCTACGACCGGTTGCAACACCTTTGATTTTGt                                                                                                               >  2:323574/1‑149 (MQ=255)
                                         ccGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCATACCGCTACGACCGGTTGCAACACCTTTGATTTTGTCTa                                                                                                            <  3:244053/149‑1 (MQ=255)
                                          cGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCATACCGCTACGACCGGTTGCAACACCTTTGATTTTGTCTAc                                                                                                           <  5:226938/149‑1 (MQ=255)
                                            gTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCATACCGCt                                                                                                                                            >  3:180052/1‑114 (MQ=255)
                                            gTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCATACCGCt                                                                                                                                            <  4:180052/114‑1 (MQ=255)
                                                 tcCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATGCTAACGTCTTCTTTTGGCACGTCCTGGTGCATACCGCTACGACCGGTTGCAACACCTTTGATTTTGTCTACCACGTc                                                                                                     >  6:209307/1‑148 (MQ=255)
                                                  ccACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCATACCGCTACGACCGGTTGCAACACCTTTGATTTTGTCTCCCACGTCCa                                                                                                   >  8:298285/1‑149 (MQ=255)
                                                   cACGCAGAGATGATGATCTGCAATAAAGAGTGGCGCCACGTTTACTCGCTTACGTTCACGCTTTCACTGATAACGTCTTCTTTTGGCACGTCCTGGTGCATACCGCTACGACCGGTTGCAACACCTTTGATTTTGTCTACCACGTCCa                                                                                                   <  5:66565/148‑1 (MQ=255)
                                                   cACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCATACCGCTACGACCGGTTGCAACACCTTTGATTTTGTCTACCACGTCCa                                                                                                   <  1:15664/148‑1 (MQ=255)
                                                    actcAGAGATGAAGATCTGAAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCATAGCGCTACGACCGGTTGCAACACCTTTGATTTTGTCTACCACGTCCATg                                                                                                 <  5:15896/146‑1 (MQ=255)
                                                      gCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCATACCGCTACGACCGGTTGCAACACCTTTGATTTTGTCTACCACGTCCATGcc                                                                                               >  1:129447/1‑149 (MQ=255)
                                                           gaTGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCATACCGCTACGACCGGTTGCAACACCTTTGATTTTGTCTACCACGTCCATGCCGTCaa                                                                                          <  3:101800/149‑1 (MQ=255)
                                                            aTGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCATACCGCTACGACCGGTTGCAACACCTTTGATTTTGTCTACCACGTCCATGCCGTCagc                                                                                         >  4:107987/1‑147 (MQ=255)
                                                            aTGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCATACCGCTACGACCGGTTGCAACACCTTTGATTTTGTCTACCACGTCCATGCCGTCAAc                                                                                         >  8:354790/1‑149 (MQ=255)
                                                             tGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCATACCGCTACGACCGGTTGCAACACCTTTGATTTTGTCTACCACGTCCATGCCGTCAAcc                                                                                        >  7:197395/1‑149 (MQ=255)
                                                              gAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCATACCGCTACGACCGGTTGCAACACCTTTGATTTTGTCTACCACGTCCATGCCGTCAACCa                                                                                       >  8:274466/1‑149 (MQ=255)
                                                              gAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCATACCGCTACGACCGGTTGCAACACCTTTGATTTTGTCTACCACGTCCATGCCGTCAACCa                                                                                       >  2:321029/1‑149 (MQ=255)
                                                               aaGATCTGCAATTATGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTGAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCATACCGCTACGACCGGTTGCAACACCTTTGATTTTGTCTACCACGTCCATGCCGTCAACCa                                                                                       <  6:57860/148‑1 (MQ=255)
                                                               aaGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCATACCGCTACGACCGGTTGCAACACCTTTGATTTTGTCTACCACGTCCATGCCGTCAACCAc                                                                                      >  1:266812/1‑149 (MQ=255)
                                                                aGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCATACCGCTACGACCGGTTGCAACACCTTTGATTTTGTCTACCACGTCCATGCCGTCAACCACt                                                                                     >  5:375344/1‑149 (MQ=255)
                                                                 gATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTTACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCATACCGCTACGACCGGTTGCAaca                                                                                                                            >  6:359417/1‑109 (MQ=255)
                                                                 gATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTTACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCATACCGCTACGACCGGTTGCAaca                                                                                                                            <  5:359417/109‑1 (MQ=255)
                                                                     tGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCATACCGCTACGACCGGTTGCAACACCTTTGATTTTGTCTACCACGTCCATGCCGTCAACCACTTCAg                                                                                 >  6:78646/1‑148 (MQ=255)
                                                                      gCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCATACCGCTACGACCGGTTGCAACACCTTTGATTTTGTCTACCACGTCCATGCCGTCAACCACTTCAGc                                                                                >  7:205095/1‑148 (MQ=255)
                                                                       cAATAAAGAGTGTCGCCACGTTTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCATACCGCTACGACCGGTTGCAACACCTTTGATTTTGTCTACCACGTCCATGCCGTCAACCACTTCa                                                                                  <  5:234268/145‑1 (MQ=255)
                                                                       cAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCATACCGCTACGACCGGTTGCAACACCTTTGATTTTGGCTACCACGTCCATGCCGTCAACCACTTCAGCa                                                                               >  7:320791/1‑148 (MQ=255)
                                                                          tAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCAGGTGCATACCGCTACGACCGGTTGCAACACCTTTGATTTTGTCTACCACGTCCATGCCGTCAACCACTTCAGCAAaca                                                                           >  3:23350/1‑149 (MQ=255)
                                                                              gagTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCATACCGCTACGACCGGTTGCAACACCTTTGATTTTGTCTACCACGTCCATGCCGTCAACCACTTCAGCAAACACGCa                                                                       <  5:371161/149‑1 (MQ=255)
                                                                                      ccACGATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCATACCGCTACGACCGGTTGCAACACCTTTGATTTTGTCTACCACGTCCATGCCGTCAACCACTTCAGCAAACACGCAGTAGcccc                                                               <  5:78646/149‑1 (MQ=255)
                                                                                         cGATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCATACCGCTACGACCGGTTGCAACACCTTTGATTTTGTCTACCACGTCCATGCCGTCAACCACTTCAGCAAACACGCAGTAGCCCCAAc                                                            <  5:271079/149‑1 (MQ=255)
                                                                                          gATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCATACCGCTACGACCGGTTGCAACACCTTTGATTTTGTCTACCACGTCCATGCCGTCAACCACTTCAGCAAACACGCAGTAGCCCCAAcc                                                           <  6:45543/149‑1 (MQ=255)
                                                                                             tACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCATACCGCTACGACCGGTTGCAACACCTTTGATTTTGTCTACCACGTCCATGCCGTCAACCACTTCAGCAAACACGCAGTAGCCCCAACCTTg                                                        >  1:321097/1‑149 (MQ=255)
                                                                                               cTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCATACCGCTACGACCGGTTGCAACACCTTTGATTTTGTCTACCACGTCCATGCCGTCAACCACTTCAGCAAACACGCAGTAGCCCCAACCTTg                                                        >  7:114031/1‑147 (MQ=255)
                                                                                                  gCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCATACCGCTACGACCGGTTGCAACACCTTTGATTTTGTCTACCACGTCCATGCCGTCAACCACTTCAGCAAACACGCAGTAGCCCCAACCTTGCAg                                                     <  8:114031/147‑1 (MQ=255)
                                                                                                  gCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCATACCGCTACGACCGGTTGCAACACCTTTGATTTTGTCTACCACGTCCATGCCGTCAACCACTTCAGCAAACACGCAGTAGCCCCAACCTTGCAg                                                     <  2:321097/147‑1 (MQ=255)
                                                                                                  gCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCATACCGCTACGACCGGTTGCAACACCTTTGATTTTGTCTACCACGTCCATGCCGTCAACCACTTCAGCAAACACGCAGTAGCCCCAACCTTGCAg                                                     <  8:320791/147‑1 (MQ=255)
                                                                                                  gCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCATACCGCTACGACCGGGTGCAACACCTTTGATTTTGTCTACCACGTCCATGCCGTCAACCACTTCAGCAAACCCGCAGTAGCCCCAACCTTGCAGGc                                                   >  2:117812/1‑149 (MQ=255)
                                                                                                  gCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCATACCGATACGACCGGTTGCAACACATTTGATTTTGTCTACAACGTCCATGCCGTCAAACACTTCAGCAAACACGCAGTAGCCCAAACCTTGCAGGc                                                   >  8:204948/1‑149 (MQ=255)
                                                                                                      aCGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCATACCGCTACGACCGGTTGCAACACCTTTGATTTTGTCTACCACGTCCATGCCGTCAGCCACTTCAGCAAACACGCAGTAGCCCCAACCTTGCAGGCTTTc                                               <  3:107987/149‑1 (MQ=255)
                                                                                                       cGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTACATACCGCTACGACCGGTTGCAACACCTTTGATTTTGTCTACCACGTCCATGCCGTCAACCACTTCAGCAAACACGCAGTAGCCCCAACCTTGCAGGCTTTCg                                              <  4:364293/149‑1 (MQ=255)
                                                                                                       cGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTACATACCGCTACGACCGGTTGCAACACCTTTGATTTTGTCTACCACGTCCATGCCGTCAACCACTTCAGCAAACACGCAGTAGCCCCAACCTTGCAGGCTTTCg                                              <  4:370012/149‑1 (MQ=255)
                                                                                                        ggTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCATACCGCTACGACCGGTTGCAACACCTTTGATTTTGTCTACCACGTCCATGCCGTCAACCACTTCAGCAAACACGCAGTAGCCCCAACCTTGCAGGCTTTCGc                                             <  8:197395/149‑1 (MQ=255)
                                                                                                           cACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCATACCGCTACGACCGGTTGCAACACCTTTGATTTTGTCTACCACGTCCATGCCGTCAACCACTTCAGCAAACACGCAGTAGCCCCAACCTTGCAGGCTTTCGc                                             >  5:92576/1‑146 (MQ=255)
                                                                                                           cACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCATACCGCTACGACCGGTTGCAACACCTTTGATTTTGTCTACCACGTCCATGCCGTCAACCACTTCAGCAAACACGCAGTAGCCCCAACCTTGCAGGCTTTCGCCa                                           <  6:92576/148‑1 (MQ=255)
                                                                                                            aCGCTTTAAAAGATAACGACTTCTTTTGGCACGTCCTCGTGCATACCCCTACCACCGGATGCAACACCTTTGATTTTGTCTACCACGTCCATGCCGTCAGCCACTTCAGCAAACACGCAGTAGCCCCAACATTGCAAGCTTTCGCCa                                           >  2:224510/1‑147 (MQ=255)
                                                                                                                    aaTGATAACGTCTTCTTTTGGCACGTCCTGGTGCATACCGCTACGACCGGTTGCAACACCTTTGATTTTGTCTACCACGTCCATGCCGTCAACCACTTCAGCAAACACGCAGTAGCCCCAACCTTGCAGGCTTTc                                               <  5:214420/135‑1 (MQ=255)
                                                                                                                    aaTGATAACGTCTTCTTTTGGCACGTCCTGGTGCATACCGCTACGACCGGTTGCAACACCTTTGATTTTGTCTACCACGTCCATGCCGTCAACCACTTCAGCAAACACGCAGTAGCCCCAACCTTGCAGGCTTTc                                               >  6:214420/1‑135 (MQ=255)
                                                                                                                              tcttctTTTGGCACGTCCTGGTGCATACCGCTACGACCGGTTGCAACACCTTTGATTTTGTCTACCACGTCCATGCCGTCAACCACTTCAGCAAACACGCAGTAGCCCCAACCTTGCAGGCTTTCGCCAGAGAAGTTCAGGAAGTCGtt                       <  5:209307/149‑1 (MQ=255)
                                                                                                                              tcttctTTTGGCACGTCCTGGTGCATACCGCTACGACCGGTTGCAACACCTTTGATTTTGTCTACCACGTCCATGCCGTCAACCACTTCAGCAAACACGCAGTAGCCCCAACCTTGCAGGCTTTCGCCAGAGAAGTTCAGGAAGTCGtt                       <  6:380665/149‑1 (MQ=255)
                                                                                                                               cttcttTTGGCACGTCCTGGTGCATACCGCTACGACCGGTTGCAACACCTTTGATTTTGTCTACCACGTCCATGCCGTCAACCACTTCAGCAAACACGCAGTAGCCCCAACCTTGCAGGCTTTCGCCAGAGAAGTTCAGGAAGTCGTTa                      >  5:388277/1‑149 (MQ=255)
                                                                                                                                       ggCACGTCCTGGTGCATACCGCTACGACCGGTTGCAACACCTTTGATTTTGTCTACCAAGTCCATGCCGTCAACCACTTCAGCAAACACGCAGTAGCCCCAACCTTGCAGGCTTTCGCCAGAGAAGTTCAGGAAGTCGTTATCAACCa               <  1:323574/148‑1 (MQ=255)
                                                                                                                                        gCACGTCCTGGTGCATACCGCTACGACCGGTTGCAACACCTTTGATTTTGTCTACCACGTCCATGCCGTCAACCACTTCAGCAAGCACGCAGTAGCCCCAACCTTGCAGGCTTTCGCCAGAGAAGTTCAGGAAGTCGTTATCAACCACg             >  6:24543/1‑149 (MQ=255)
                                                                                                                                        gCACGTCCTGGTGCATACCGCTACGACCGGTTGCAACACCTTTGATTTTGTCTACCACGTCCATGCCGTCAACCACTTCAGCAAACACGCAGTAGCCCCAACCTTGCAGGCTTTCGCCAGAGAAGTTCAGGAAGTCGTTATCAACCACg             >  3:221511/1‑149 (MQ=255)
                                                                                                                                            gTCCTGGTGCATACCGCTACGACCGGTTGCAACACCTTTGATTTTGTCTACCACGTCCATGCCGTCAACCACTTCAGCAAACACGCAGTAGCCCCAACCTTGCAGGCTTTCGCCAGAGAAGTTCAGGAAGTCGTTATCAACCACGTtga         <  8:205095/149‑1 (MQ=255)
                                                                                                                                             tCCTGGTGCATACCGCTACGACCGGTTGCAACACCTTTGATTTTGTCTACCACGTCCATGCCGTCAACCACTTCAGCAAGCACGCAGTAGCCCCAACCTTGCAGGCTTTCGCCAGAGAAGTTCAGGAAGTCGTTATCAACCACGTtgat        <  5:24543/149‑1 (MQ=255)
                                                                                                                                                  gTGCATACCGCTACGACCGGTTGCAACACCTTTGATTTTGTCTACCACGTCCATGCCGTCAACCACTTCAGCAAACACGCAGTAGCCCCAACCTTGCAGGCTTTCGCCaga                                         >  1:30291/1‑111 (MQ=255)
                                                                                                                                                  gTGCATACCGCTACGACCGGTTGCAACACCTTTGATTTTGTCTACCACGTCCATGCCGTCAACCACTTCAGCAAACACGCAGTAGCCCCAACCTTGCAGGCTTTCGCCaga                                         <  2:30291/111‑1 (MQ=255)
                                                                                                                                                   tGCATACCGCTACGACCGGTTGCAACACCTTTGATTTTGTCTACCACGTCCATGCCGTCAACCACTTCAGCAAACACGCAGTAGCCCCAACCTTGCAGGCTTTCGCCAGAGAAGTTCAGGAAGTCGTTATCAACCACGTTGATgaagaa  <  7:198316/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCATACCGCTACGACCGGTTGCAACACCTTTGATTTTGTCTACCACGTCCATGCCGTCAACCACTTCAGCAAACACGCAGTAGCCCCAACCTTGCAGGCTTTCGCCAGAGAAGTTCAGGAAGTCGTTATCAACCACGTTGATGAAGAA  >  NC_000913/553851‑554146

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: ‑