Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 455,632 A→G 5.5% D167G (GAC→GGC)  tig → peptidyl‑prolyl cis/trans isomerase (trigger factor)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913455,6320AG 5.5% 253.5 / 3.7 73D167G (GAC→GGC) tigpeptidyl‑prolyl cis/trans isomerase (trigger factor)
Reads supporting (aligned to +/- strand):  ref base A (35/34);  new base G (2/2);  total (37/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.96e-01

GGTCTGGAAGCGATCGAAGTTGAAAAACCGATCGTTGAAGTGACCGACGCTGACGTTGACGGCATGCTGGATACTCTGCGTAAACAGCAGGCGACCTGGAAAGAAAAAGACGGCGCTGTTGAAGCAGAAGACCGCGTAACCATCGACTTCACCGGTTCTGTAGACGGCGAAGAGTTCGAAGGCGGTAAAGCGTCTGATTTCGTACTGGCGATGGGCCAGGGTCGTATGATCCCGGGCTTTGAAGACGGTATCAAAGGCCACAAAGCTGGCGAAGAGTTCACCATCGACGTG  >  NC_000913/455487‑455777
                                                                                                                                                 |                                                                                                                                                 
ggTCTGGAAGCGATCGAAGTTGAAAAACCGATCGTTGAAGTGACCGACGCTGACGTTGACGGCATGCTGGATACTCTGCGTAAACAGCAGGCGACCTGGAAAGAAAAAGACGGCGCTGTTGAAGCAGAAGACCGCGTAACCATCGACt                                                                                                                                                 >  4:775802/1‑148 (MQ=255)
    tGGAAGCGATCGAAGTTGAAAAACCGATCGTTGAAGTGACCGACGCTGACGTTGACGGCATGCTGGATACTCTGCGTAAACAGCAGGCGACCTGGAAAGAAAAAGACGGCGCTGTTGAAGCAGAAGACCGCGTAACCATCGACTTCAc                                                                                                                                             <  3:387014/148‑1 (MQ=255)
           gATCGAAGTTGAAAAACCGATCGTTGAAGTGACCGACGCTGACGTTGACGGCATGCTGGATACTCTGCGTAAACAGCAGGCGACCTGGAAAGAAAAAGACGGCGCTGTTGAAGCAGAAGACCGCGTAACCATCGACTTCACCGGTTCt                                                                                                                                      <  3:775802/148‑1 (MQ=255)
             tCGAAGTTGAAAAACCGATCGTTGAAGTGACCGACGCTGACGTTGACGGCATGCTGGATACTCTGCGTAAACAGCAGGCGACCTGGAAAGAAAAAGACGGCGCTGTTGAAGCAGAAGACCGCGTAACCATCGACTTCACCGGTTCTGt                                                                                                                                    <  1:332794/148‑1 (MQ=255)
                        aaaCCGATCGTTGAAGTGACCGACGCTGACGTTGACGGCATGCTGGATACTCTGCGTAAACAGCAGGCGACCTGGAAAGAAAAAGACGGCGCTGTTGAAGCAGAAGACCGCGTAACCATCGACTTCACCGGTTCTGTAGACGGCGAag                                                                                                                         >  3:391456/1‑148 (MQ=255)
                                  ttGAAGTGACCGACGCTGACGTTGACGGCATGCTGGATACTCTGCGTAAACAGCAGGCGACCTGGAAAGAAAAAGACGGCGCTGTTGAAGCAGAAGACCGCGTAACCATCGACTTCACCGGTTCTGTAGACGGCGAAGAGTTCGAAgg                                                                                                               >  1:1105266/1‑148 (MQ=255)
                                      aGTGACCGACGCTGACGTTGACGGCATGCTGGATACTCTGCGTAAACAGCAGGCGACCTGGAAAGAAAAAGACGGCGCTGTTGAAGCAGAAGACCGCGTAACCATCGACTTCACCGGTTCTGTAGACGGCGAag                                                                                                                         <  4:947877/134‑1 (MQ=255)
                                      aGTGACCGACGCTGACGTTGACGGCATGCTGGATACTCTGCGTAAACAGCAGGCGACCTGGAAAGAAAAAGACGGCGCTGTTGAAGCAGAAGACCGCGTAACCATCGACTTCACCGGTTCTGTAGACGGCGAag                                                                                                                         >  3:947877/1‑134 (MQ=255)
                                        tGACCGACGCTGACGTTGACGGCATGCTGGATACTCTGCGTAAACAGCAGGCGACCTGGAAAGAAAAAGACGGCGCTGTTGAAGCAGAAGACCGCGTAACCATCGACTTCACCGGTTCTGTAGACGGCGAAGAGTTCGAAGGCGGTaa                                                                                                         <  4:391456/148‑1 (MQ=255)
                                              aCGCTGACGTTGACGGCATGCTGGATACTCTGCGTAAACAGCAGGCGACCTGGAAAGAAAAAGACGGCGCTGTTGAAGCAGAAGACCGCGTAACCATCGACTTCACCGGTTCTGt                                                                                                                                    >  2:383088/1‑115 (MQ=255)
                                              aCGCTGACGTTGACGGCATGCTGGATACTCTGCGTAAACAGCAGGCGACCTGGAAAGAAAAAGACGGCGCTGTTGAAGCAGAAGACCGCGTAACCATCGACTTCACCGGTTCTGt                                                                                                                                    <  1:383088/115‑1 (MQ=255)
                                              aCGCTGACGTTGACGGCATGCTGGATACTCTGCGTAAACAGCAGGCGACCTGGAAAGAAAAAGACGGCGCTGTTGAAGCAGAAGACCGCGTAACCATCGACTTCACCGGTTCTGTAGACGGCGAAGAGTTCGAAGGCGGTAAAGCGTc                                                                                                   >  4:138847/1‑148 (MQ=255)
                                                 cTGACGTTGACGGCATGCTGGATACTCTGCGTAAACAGCAGGCGACCTGGAAAGAAAAAGACGGCGCTGTTGAAGCAGAAGACCGCGTAACCATCGACTTCACCGGTTCTGTAGACGGCGAAGAGTTCGAAGGCGGTAAAGCGTCTGa                                                                                                <  3:138847/148‑1 (MQ=255)
                                                       ttGACGGCATGCTGGATACTCTGCGTAAACAGCAGGCGACCTGGAAAGAAAAAGACGGCGCTGTTGAAGCAGAAGACCGCGTAACCATCGACTTCACCGGTTCTGTAGACGGCGAAGAGTTCGAAGGCGGTAAAGCGTCTGATTTCGt                                                                                          <  2:1105266/148‑1 (MQ=255)
                                                        tGACGGCATGCTGGATACTCTGCGTAAACAGCAGGCGACCTGGAAAGAAAAAGACGGCGCTGTTGAAGCAGAAGACCGCGTAACCATCGACTTCACCGGTTCTGTAGACGGCGAAGAGTTCGAAGGCGGTAAAGCGTCTGATTTCGTa                                                                                         >  2:1089580/1‑148 (MQ=255)
                                                                     gATACTCTGCGTAAACAGCAGGCGACCTGGAAAGAAAAAGACGGCGCTGTTGAAGCAGAAGACCGCGTAACCATCGACTTCACCGGTTCTGTAGACGGCGAAGAGTTCGAAGGCGGTAAAGCGTCTGATTTCGTACTGGCGATGGGcc                                                                            >  2:757540/1‑148 (MQ=255)
                                                                        aCTCTGCGTAAACAGCAGGCGACCTGGAAAGAAAAAGACGGCGCTGTTGAAGCAGAAGACCGCGTAACCATCGACTTCACCGGTTCTGTAGACGGCGAAGAGTTCGAAGGCGGTAAAGCGTCTGATTTCGTACTGGCGATGGGCCAg                                                                          <  2:558532/147‑1 (MQ=255)
                                                                        aCTCTGCGTAAACAGCAGGCGACCTGGAAAGAAAAAGACGGCGCTGTTGAAGCAGAAGACCGCGTAACCATCGACTTCACCGGTTCTGTAGACGGCGAAGAGTTCGAAGGCGGTAAAGCGTCTGATTTCGTACTGGCGATGGGCCAg                                                                          >  1:558532/1‑147 (MQ=255)
                                                                          tctGCGTAAACAGCAGGCGACCTGGAAAGAAAAAGACGGCGCTGTTGAAGCAGAAGACCGCGTAACCATCGACTTCACCGGTTCTGTAGACGGCGAag                                                                                                                         <  1:198729/98‑1 (MQ=255)
                                                                          tctGCGTAAACAGCAGGCGACCTGGAAAGAAAAAGACGGCGCTGTTGAAGCAGAAGACCGCGTAACCATCGACTTCACCGGTTCTGTAGACGGCGAag                                                                                                                         >  2:198729/1‑98 (MQ=255)
                                                                           ctGCGTAAACAGCAGGCGACCTGGAAAGAAAAAGACGGCGCTGTTGAAGCAGAAGACCGTGTAACCATCGACTTCACCGGTTCTGTAGACGGCGAAGAGTTCGAAGGCGGTAAAGCGTCTGATTTCGTACTGGCg                                                                                   >  4:235263/1‑135 (MQ=255)
                                                                           ctGCGTAAACAGCAGGCGACCTGGAAAGAAAAAGACGGCGCTGTTGAAGCAGAAGACCGTGTAACCATCGACTTCACCGGTTCTGTAGACGGCGAAGAGTTCGAAGGCGGTAAAGCGTCTGATTTCGTACTGGCg                                                                                   <  3:235263/135‑1 (MQ=255)
                                                                             gCGTAAACAGCAGGCGACCTGGAAAGAAAAAGACGGCGCTGTTGAAGCAGAAGACCGCGTAACCATCGACTTCACCGGTTCTGTAGACGGCGAAGAGTTCGAAGGCGGTAAAGCGTCTGATTTCGTACTGGCGATGGGCCAGGGTCGt                                                                    >  4:630481/1‑148 (MQ=255)
                                                                             gCGTAAACAGCAGGCGACCTGGAAAGAAAAAGACGGCGCTGTTGAAGCAGAAGACCGCGTAACCATCGACTTCACCGGTTCTGTAGACGGCGAAGAGTTCGAAGGCGGTAAAGCGTCTGATTTCGTACTGGCGATGGGCCAGGGTCGt                                                                    <  3:630481/148‑1 (MQ=255)
                                                                                  aaCAGCAGGCGACCTGGAAAGAAAAAGACGGCGCTGTTGAAGCAGAAGACCGCGTAACCATCGACTTCACCGGTTCTGTAGACGGCGAAGAGTTCGAAGGCGGTAAAGCGTCTGATTTCGTACTGGCGATGGGCCAGGGTCGTATGAt                                                               >  4:1086994/1‑148 (MQ=255)
                                                                                      gcagGCGACCTGGAAAGAAAAAGACGGCGCTGTTGAAGCAGAAGACCGCGTAACCATCGACTTCACCGGTTCTGTAGACGGCGAAGAGTTCGAAGGCGGTAAAGCGTCTGATTTCGTACTGGCGATGGGCCAGGGTCg                                                                     <  3:27151/138‑1 (MQ=255)
                                                                                      gcagGCGACCTGGAAAGAAAAAGACGGCGCTGTTGAAGCAGAAGACCGCGTAACCATCGACTTCACCGGTTCTGTAGACGGCGAAGAGTTCGAAGGCGGTAAAGCGTCTGATTTCGTACTGGCGATGGGCCAGGGTCg                                                                     <  3:599303/138‑1 (MQ=255)
                                                                                      gcagGCGACCTGGAAAGAAAAAGACGGCGCTGTTGAAGCAGAAGACCGCGTAACCATCGACTTCACCGGTTCTGTAGACGGCGAAGAGTTCGAAGGCGGTAAAGCGTCTGATTTCGTACTGGCGATGGGCCAGGGTCg                                                                     >  4:27151/1‑138 (MQ=255)
                                                                                      gcagGCGACCTGGAAAGAAAAAGACGGCGCTGTTGAAGCAGAAGACCGCGTAACCATCGACTTCACCGGTTCTGTAGACGGCGAAGAGTTCGAAGGCGGTAAAGCGTCTGATTTCGTACTGGCGATGGGCCAGGGTCg                                                                     >  4:599303/1‑138 (MQ=255)
                                                                                         ggCGACCTGGAAAGAAAAAGACGGCGCTGTTGAAGCAGAAGACCGCGTAACCATCGACTTCACCGGTTCTGTAGACGGCGAAGAGTTCGAAGGCGGTAAAGCGTCTGATTTCGTACTGGCGATGGGCCAGGGTCGt                                                                    <  1:301099/136‑1 (MQ=255)
                                                                                         ggCGACCTGGAAAGAAAAAGACGGCGCTGTTGAAGCAGAAGACCGCGTAACCATCGACTTCACCGGTTCTGTAGACGGCGAAGAGTTCGAAGGCGGTAAAGCGTCTGATTTCGTACTGGCGATGGGCCAGGGTCGt                                                                    >  2:301099/1‑136 (MQ=255)
                                                                                          gCGACCTGGAAAGAAAAAGACGGCGCTGTTGAAGCAGAAGACCGCGTAACCATCGACTTCACCGGTTCTGTAGACGGCGAAGAGTTCGAAGGCGGTAAAGCGTCTGAt                                                                                               >  2:32011/1‑108 (MQ=255)
                                                                                          gCGACCTGGAAAGAAAAAGACGGCGCTGTTGAAGCAGAAGACCGCGTAACCATCGACTTCACCGGTTCTGTAGACGGCGAAGAGTTCGAAGGCGGTAAAGCGTCTGAt                                                                                               <  1:32011/108‑1 (MQ=255)
                                                                                           cGACCTGGAAAGAAAAAGACGGCGCTGTTGAAGCAGAAGACCGCGTAACCATCGACTTCACCGGTTCt                                                                                                                                      >  1:420347/1‑68 (MQ=255)
                                                                                           cGACCTGGAAAGAAAAAGACGGCGCTGTTGAAGCAGAAGACCGCGTAACCATCGACTTCACCGGTTCt                                                                                                                                      <  2:420347/68‑1 (MQ=255)
                                                                                                  gaaagaaaAAGACGGCGCTGTTGAAGCAGAAGACCGCGTAACCATCGACTTCACCGGTTCt                                                                                                                                      >  4:1078115/1‑61 (MQ=255)
                                                                                                  gaaagaaaAAGACGGCGCTGTTGAAGCAGAAGACCGCGTAACCATCGACTTCACCGGTTCt                                                                                                                                      >  1:1018077/1‑61 (MQ=255)
                                                                                                  gaaagaaaAAGACGGCGCTGTTGAAGCAGAAGACCGCGTAACCATCGACTTCACCGGTTCt                                                                                                                                      <  3:1078115/61‑1 (MQ=255)
                                                                                                  gaaagaaaAAGACGGCGCTGTTGAAGCAGAAGACCGCGTAACCATCGACTTCACCGGTTCt                                                                                                                                      <  2:1018077/61‑1 (MQ=255)
                                                                                                  gaaagaaaAAGACGGCGCTGTTGAAGCAGAAGACCGCGTAACCATCGACTTCACCGGTTCTGTAGACGGCGAAGAGTTCGAAGGCGGTAAAGCGTCTGAt                                                                                               >  3:258451/1‑100 (MQ=255)
                                                                                                  gaaagaaaAAGACGGCGCTGTTGAAGCAGAAGACCGCGTAACCATCGACTTCACCGGTTCTGTAGACGGCGAAGAGTTCGAAGGCGGTAAAGCGTCTGAt                                                                                               <  4:258451/100‑1 (MQ=255)
                                                                                                   aaagaaaAAGACGGCGCTGTTGAAGCAGAAGACCGCGTAACCATCGGCTTCACCGGTTCTGTAGACGGCGAagag                                                                                                                       >  4:597718/1‑75 (MQ=255)
                                                                                                   aaagaaaAAGACGGCGCTGTTGAAGCAGAAGACCGCGTAACCATCGGCTTCACCGGTTCTGTAGACGGCGAagag                                                                                                                       >  4:632532/1‑75 (MQ=255)
                                                                                                   aaagaaaAAGACGGCGCTGTTGAAGCAGAAGACCGCGTAACCATCGGCTTCACCGGTTCTGTAGACGGCGAagag                                                                                                                       <  3:597718/75‑1 (MQ=255)
                                                                                                   aaagaaaAAGACGGCGCTGTGGAAGCAGAAGACCGCGTAACCATCGGCTTCACCGGTTCTGTAGACGGCGAagag                                                                                                                       <  3:632532/75‑1 (MQ=255)
                                                                                                        aaaaGCCGGCGCTGTTGAAGCAGAAGACCGCGTAACCATCGACTTCACCGGTTCTGTAGACGGCGAAGAGTTCGAAGGCGGTAAAGCGTCTGATTTCGTACTGGCGATGGGCCAGGGTCGt                                                                    <  3:699269/121‑1 (MQ=255)
                                                                                                        aaaaGACGGCGCTGTTGAAGCAGAAGACCGCGTAACCATCGACTTCACCGGTTCTGTAGACGGCGAAGAGTTCGAAGGCGGTAAAGCGTCTGATGTCGTACTGGCGATGGGCCAGGGTCGt                                                                    >  4:699269/1‑121 (MQ=255)
                                                                                                          aaGACGGCGCTGTTGAAGCAGAAGACCGCGTAACCATCGACTTCACCGGTTCTGTAGACGGCGAAGAGTTCGAAGGCGGTAAAGCGTCTGATTTCGTACTGGCGATGGGCCAGGGTCGTATGATCCCGGGCTTTGAAGACGGTATCaa                                       <  1:1089580/148‑1 (MQ=255)
                                                                                                           aGACGGCGCTGTTGAAGCAGAAGACCGCGTAACCATCGACTTCACCGGTTCTGTAGACGGCGAAGAGTTCGAAGGCGGTAAAGCGTCTGATTTCGTACTGGCGATGGGCCAGGGTCGTATGATCCCGGGCTTTGAAGACGGTATCaaa                                      <  3:1086994/148‑1 (MQ=255)
                                                                                                             aCGGCGCTGTTGAAGCAGAAGACCGCGTAACCATCGACTTCACCGGTTCTGTAGACGGCGAAGAGTTCGAAGGCGGTAAAGCGTCTGATTTCGTACTGGCg                                                                                   <  1:606019/101‑1 (MQ=255)
                                                                                                             aCGGCGCTGTTGAAGCAGAAGACCGCGTAACCATCGACTTCACCGGTTCTGTAGACGGCGAAGAGTTCGAAGGCGGTAAAGCGTCTGATTTCGTACTGGCg                                                                                   >  2:606019/1‑101 (MQ=255)
                                                                                                               ggCGCTGTTGAAGCAGAAGACCGCGTAACCATCGACTTCACCGGTTCTGTAGACGGCGAAGAGTTCGAAGGCGGTAAAGCGTCTGATTTCGTACTGGCg                                                                                   <  4:90792/99‑1 (MQ=255)
                                                                                                               ggCGCTGTTGAAGCAGAAGACCGCGTAACCATCGACTTCACCGGTTCTGTAGACGGCGAAGAGTTCGAAGGCGGTAAAGCGTCTGATTTCGTACTGGCg                                                                                   >  3:90792/1‑99 (MQ=255)
                                                                                                                     gTTGAAGCAGAAGACCGCGTAACCATCGACTTCACCGGTTCTGTAGACGGCGAAGAGTTCGAAGGCGGTAAAGCGTCTGATTTCGTACTGGCGATGGGCCAGGGTCGTATGATCCCGGGCTTTGAAGACGGTATCAAAGGCCACAAAg                            <  1:757540/148‑1 (MQ=255)
                                                                                                                       tGAAGCAGAAGACCGCGTAACCATCGACTTCACCGGTTCTGTAGACGGCGAAGAGTTCGAAGGCGGTAAAGCGTCTGATTTCGTACTGGCGATGGGCCAGGGTCGTATGATCCCGGGCTTTGAAGACg                                              >  1:415642/1‑128 (MQ=255)
                                                                                                                       tGAAGCAGAAGACCGCGTAACCATCGACTTCACCGGTTCTGTAGACGGCGAAGAGTTCGAAGGCGGTAAAGCGTCTGATTTCGTACTGGCGATGGGCCAGGGTCGTATGATCCCGGGCTTTGAAGACg                                              <  2:415642/128‑1 (MQ=255)
                                                                                                                          aGCAGAAGACCGCGTAACCATCGACTTCACCGGTTCTGTAGACGGCGAAGAGTTCGAAGGCGGTAAAGCGTCTGATTTCGTACTGGCGATGGGCCAGGGTCGTATGATCCCGGGCTTTGAAGACGGTATCAAAGGCCACAAAGCTGGc                       >  2:664812/1‑148 (MQ=255)
                                                                                                                               aagaCCGCGTAACCATCGACTTCACCGGTTCTGTAGACGGCGAAGAGTTCGAAGGCGGTAAAGCGTCTGATTTCGTACTGGCGATGGGCCAg                                                                          <  3:1121194/92‑1 (MQ=255)
                                                                                                                               aagaCCGCGTAACCATCGACTTCACCGGTTCTGTAGACGGCGAAGAGTTCGAAGGCGGTAAAGCGTCTGATTTCGTACTGGCGATGGGCCAg                                                                          >  4:1121194/1‑92 (MQ=255)
                                                                                                                               aagaCCGCGTAACCATCGACTTCACCGGTTCTGTAGACGGCGAAGAGTTCGAAGGCGGTAAAGCGTCTGATTTCGTACTGGCGATGGGCCAGGGTCGTATGATCCCGGGCTTTGAAGACGGTATCAAAGGCCACAAAGCTGGCGAaga                  >  4:114520/1‑148 (MQ=255)
                                                                                                                               aagaCCGCGTAACCATCGACTTCACCGGTTCTGTAGACGGCGAAGAGTTCGAAGGCGGTAAAGCGTCTGATTTCGTACTGGCGATGGGCCAGGGTCGTATGATCCCGGGCTTTGAAGACGGTATCAAAGGCCACAAAGCTGGCGAaga                  >  2:383611/1‑148 (MQ=255)
                                                                                                                                agaCCGCGTAACCATCGACTTCACCGGTTCTGTAGACGGCGAAGAGTTCGAAGGCGGTAAAGCGTCTGATTTCGTACTGGCGATggg                                                                              <  3:1002200/87‑1 (MQ=255)
                                                                                                                                agaCCGCGTAACCATCGACTTCACCGGTTCTGTAGACGGCGAAGAGTTCGAAGGCGGTAAAGCGTCTGATTTCGTACTGGCGATggg                                                                              >  4:1002200/1‑87 (MQ=255)
                                                                                                                                 gaCCGCGTAACCATCGACTTCACCGGTTCTGTAGACGGCGAAGAGTTCGAAGGCGGTAAAGCGTCTGATTTCGTACTGGCGATg                                                                                <  2:137759/84‑1 (MQ=255)
                                                                                                                                 gaCCGCGTAACCATCGACTTCACCGGTTCTGTAGACGGCGAAGAGTTCGAAGGCGGTAAAGCGTCTGATTTCGTACTGGCGATg                                                                                >  1:137759/1‑84 (MQ=255)
                                                                                                                                          aCCATCGACTTCACCGGTTCTGTAGACGGCGAAGAGTTCGAAGGCGGTAAAGCGTCTGATTTCGTACTGGCGATGGGCCAGGGTCGTATGATCCCGGGCTTTGAAGACGGTATCAAAGGCCACAAAg                            <  1:1035913/127‑1 (MQ=255)
                                                                                                                                          aCCATCGACTTCACCGGTTCTGTAGACGGCGAAGAGTTCGAAGGCGGTAAAGCGTCTGATTTCGTACTGGCGATGGGCCAGGGTCGTATGATCCCGGGCTTTGAAGACGGTATCAAAGGCCACAAAg                            >  2:1035913/1‑127 (MQ=255)
                                                                                                                                          aCCATCGACTTCACCGGTTCTGTAGACGGCGAAGAGTTCGAAGGCGGTAAAGCGTCTGATTTCGTACTGGCGATGGGCCAGGGTCGTATGATCCCGGGCTTTGAAGACGGTATCAAAGGCCACAAAGCTGGCGAAGAGTTCACCATCg       <  1:383611/148‑1 (MQ=255)
                                                                                                                                             aTCGACTTCACCGGTTCTGTAGACGGCGAAGAGTTCGAAGGCGGTAAAGCGTCTGATTTCGTACTGGCGATGGGCCAGGGTCGTATGATCCCggg                                                         >  2:998679/1‑95 (MQ=255)
                                                                                                                                             aTCGACTTCACCGGTTCTGTAGACGGCGAAGAGTTCGAAGGCGGTAAAGCGTCTGATTTCGTACTGGCGATGGGCCAGGGTCGTATGATCCCggg                                                         <  1:998679/95‑1 (MQ=255)
                                                                                                                                              tCGACTTCACCGGTTCTGTAGACGGCGAAGAGTTCGAAGGCGGTAAAGCGTCTGATTTCGTACTGGCGATGGGCCAGGGt                                                                       <  2:372497/80‑1 (MQ=255)
                                                                                                                                              tCGACTTCACCGGTTCTGTAGACGGCGAAGAGTTCGAAGGCGGTAAAGCGTCTGATTTCGTACTGGCGATGGGCCAGGGt                                                                       >  1:372497/1‑80 (MQ=255)
                                                                                                                                               cGACTTCACCGGTTCTGTAGACGGCGAAGAGTTCGAAGGCGGTAAAGCGTCTGATTTCGTACTGGCGATGGGCCAGGGTCGTATGATCCCGGGCTTTGAAGACGGTATCAAAGGCCACAAAGCTGGCGAAGAGTTCACCATCGACGTg  >  4:1079722/1‑148 (MQ=255)
                                                                                                                                                 aCTTCACCGGTTCTGTAGACGGCGAAGAGTTCGAAGGCGGTAAAGCGTCTGATTTCGTACTGGCGAt                                                                                 >  2:110577/1‑67 (MQ=255)
                                                                                                                                                 aCTTCACCGGTTCTGTAGACGGCGAAGAGTTCGAAGGCGGTAAAGCGTCTGATTTCGTACTGGCGAt                                                                                 <  1:110577/67‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                 
GGTCTGGAAGCGATCGAAGTTGAAAAACCGATCGTTGAAGTGACCGACGCTGACGTTGACGGCATGCTGGATACTCTGCGTAAACAGCAGGCGACCTGGAAAGAAAAAGACGGCGCTGTTGAAGCAGAAGACCGCGTAACCATCGACTTCACCGGTTCTGTAGACGGCGAAGAGTTCGAAGGCGGTAAAGCGTCTGATTTCGTACTGGCGATGGGCCAGGGTCGTATGATCCCGGGCTTTGAAGACGGTATCAAAGGCCACAAAGCTGGCGAAGAGTTCACCATCGACGTG  >  NC_000913/455487‑455777

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