Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,494,883:1 (A)7→8 100% coding (261/564 nt) idnK → D‑gluconate kinase, thermosensitive

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,494,8761.A93.7% 250.2 / 9.7 79coding (254/564 nt)idnKD‑gluconate kinase, thermosensitive
Reads supporting (aligned to +/- strand):  ref base . (3/2);  new base A (35/39);  total (38/41)
Fisher's exact test for biased strand distribution p-value = 6.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATCTGATGAAGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATTCTCGCGCGAATGC  >  NC_000913/4494778‑4494971
                                                                                                   |                                                                                                
aTCTGATGAAGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTa                                                                                                 <  1:1904193/100‑1 (MQ=255)
 tCTGATGAAGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTaa                                                                                                <  1:968656/100‑1 (MQ=255)
  cTGATGAAGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTaaa                                                                                              >  2:92381/1‑100 (MQ=255)
  cTGATGAAGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTaaa                                                                                              >  1:1693971/1‑100 (MQ=255)
   tgatgaAGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTaaaa                                                                                             >  2:1552722/1‑100 (MQ=255)
    gatgaAGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTaaaaa                                                                                            >  1:804790/1‑100 (MQ=255)
      tgaAGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTaaaaaaa                                                                                          <  1:733079/100‑1 (MQ=255)
      tgaAGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTaaaaaaa                                                                                          >  1:1639346/1‑100 (MQ=255)
        aaGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTaaaaaaaac                                                                                        <  1:1648624/100‑3 (MQ=255)
        aaGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTaaaaaaaac                                                                                        <  1:2109983/100‑3 (MQ=255)
         aGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTaaaaaaaaca                                                                                       <  1:2054855/100‑4 (MQ=255)
          gATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTaaaaaaaacag                                                                                      >  2:1756707/1‑96 (MQ=255)
          gATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTaaaaaaaacag                                                                                      <  2:1954943/100‑5 (MQ=255)
          gATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTaaaaaaaac                                                                                        <  1:702812/98‑3 (MQ=255)
              gACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAACAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAACAAAAACAGTATc                                                                                   >  2:98199/1‑100 (MQ=255)
              gACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATc                                                                                   >  2:1805160/1‑100 (MQ=255)
               aCTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCg                                                                                  <  2:656936/100‑1 (MQ=255)
               aCTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCg                                                                                  >  1:792504/1‑100 (MQ=255)
                cTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGt                                                                                 <  1:406776/100‑1 (MQ=255)
                  tCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTg                                                                                >  1:1907061/1‑99 (MQ=255)
                  tCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGa                                                                               >  1:1751430/1‑100 (MQ=255)
                  tCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGa                                                                               >  1:381529/1‑100 (MQ=255)
                  tCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTATGTTCATCATTAAAAAAAACAGTATAGTGa                                                                               >  1:729259/1‑100 (MQ=255)
                    ccTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGata                                                                             >  2:1776384/1‑100 (MQ=255)
                    ccTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGCGat                                                                              >  2:34701/1‑99 (MQ=255)
                    ccTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAAAAGTATCGTGata                                                                             >  2:535077/1‑100 (MQ=255)
                    ccTGGATGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAATAATGAAACAGTATTTATTGTCTGTTCATAATTAAAAAAAACAGTATCGt                                                                                 <  2:181378/96‑1 (MQ=255)
                        gctggAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATAtttt                                                                         <  1:944460/100‑1 (MQ=255)
                        gattaaaCGCATAAATGATGCTCCATACATTCTTTATAAAAAGAATGAAACAGTATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATAttt                                                                          <  2:1957395/94‑1 (MQ=255)
                          tggAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTAc                                                                       <  1:734019/100‑1 (MQ=255)
                            gAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGc                                                                     <  1:1545666/100‑1 (MQ=255)
                             aaCGCTTAAATGATGCTTCATACAGTCTTTATAACAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCa                                                                    >  1:2190650/1‑100 (MQ=255)
                               cGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAg                                                                  >  1:274174/1‑100 (MQ=255)
                               cGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAACAGTATCGTGATATTTTACGCa                                                                    >  2:1824517/1‑97 (MQ=255)
                                gCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTGATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAg                                                                  >  2:1640521/1‑99 (MQ=255)
                                 cTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAggg                                                                <  1:975956/100‑1 (MQ=255)
                                 cTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAggg                                                                <  2:1751419/100‑1 (MQ=255)
                                 cTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAg                                                                  <  1:535076/98‑1 (MQ=255)
                                  ttAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGTTATGTTACGCAAg                                                                  >  2:1162552/1‑97 (MQ=255)
                                  ttAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGt                                                               <  1:1669637/100‑1 (MQ=255)
                                   tAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTa                                                              >  1:859181/1‑100 (MQ=255)
                                   tAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTa                                                              <  1:555824/100‑1 (MQ=255)
                                      atgatgCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGcc                                                           <  1:644751/100‑1 (MQ=255)
                                        gatgggtCATACAGTCTTTATAAAAAAAATGATACAGGATTTATTGTATGTTTATCATTAAAAAAAACAGTATCTTTATATTTTACGCAAGGGTAGcc                                                           <  2:1884474/92‑1 (MQ=255)
                                         atgCTTCATACAGTTTTTATAAAAAGAATGAAACAAGACTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTAATATTTTATGCAAGGATat                                                             >  2:112893/1‑94 (MQ=255)
                                           gCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCAt                                                      >  2:617419/1‑100 (MQ=255)
                                           gCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCAt                                                      >  1:1735320/1‑100 (MQ=255)
                                               cATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTc                                                  <  2:113466/100‑1 (MQ=255)
                                                ttccgatCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCa                                                 <  1:2099755/94‑1 (MQ=255)
                                                aTACAGTCTTTATAAAAAGAATGAAACAGTATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCa                                                 <  2:720220/100‑1 (MQ=255)
                                                aTACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCa                                                 <  1:146700/100‑1 (MQ=255)
                                                 tACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCAt                                                >  2:1139407/1‑100 (MQ=255)
                                                 tACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCAt                                                >  1:429797/1‑100 (MQ=255)
                                                      tCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTcc                                            <  2:804790/98‑1 (MQ=255)
                                                        tttATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCctct                                         <  1:560146/100‑1 (MQ=255)
                                                         ttATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCaagatcggaa                                        >  2:2099741/1‑91 (MQ=255)
                                                         ttATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTg                                        <  2:313225/100‑1 (MQ=255)
                                                          tataAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTgg                                       >  2:1901726/1‑100 (MQ=255)
                                                          tataAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTgg                                       <  1:400973/100‑1 (MQ=255)
                                                           ataAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGt                                      >  2:1155988/1‑100 (MQ=255)
                                                           agaaaaaGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGAGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGt                                      <  1:655843/98‑1 (MQ=255)
                                                            taAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGtt                                     <  1:961310/100‑1 (MQ=255)
                                                            taAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGAATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGtt                                     >  2:556804/1‑100 (MQ=255)
                                                             aaaaaGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTa                                    <  2:966242/100‑1 (MQ=255)
                                                              aaaaGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAg                                   >  1:1328506/1‑100 (MQ=255)
                                                              aaaaGAATGAAACAGGATTTATTGTATGTTCATCATTAAAAAAAACAGTATCGTGATATTTTAGGAAAGGGTAGCACCCATGTTAATTTACTCTGGTTAg                                   >  2:932634/1‑100 (MQ=255)
                                                               aaaGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGa                                  >  2:676742/1‑100 (MQ=255)
                                                               aaaGAATGAAACAGGATTTATTATCTGTTCATCATTAAAAAAAACAGTATCGTCATATTTTACGCAAGGATAGCCCCCATGTTCATTTCCTCTGGTTaca                                  >  2:2129627/1‑98 (MQ=255)
                                                                 aGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGAt                                 >  2:854473/1‑99 (MQ=255)
                                                                   aaTGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACATTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATgg                               <  2:440775/99‑1 (MQ=255)
                                                                   aaTGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGc                              <  1:722709/100‑1 (MQ=255)
                                                                    aTGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCg                             >  1:1983871/1‑100 (MQ=255)
                                                                    aTGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCg                             >  1:1384100/1‑100 (MQ=255)
                                                                     tGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGa                            >  1:1734639/1‑100 (MQ=255)
                                                                        aaCAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTa                         <  2:12102/100‑1 (MQ=255)
                                                                        aaCAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTa                         <  2:2204298/100‑1 (MQ=255)
                                                                              aTTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAAc                   <  1:866518/100‑1 (MQ=255)
                                                                               tttATTGTCTGTTCATCATTAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAAc                   <  1:214071/98‑1 (MQ=255)
                                                                                    tGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATtct             <  2:949784/100‑1 (MQ=255)
                                                                                     gtatGTTTAACATTAAAAAAAACAGTATCGTGTTATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATtctc            <  2:476169/97‑1 (MQ=255)
                                                                                     gTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGTAACTATtctc            >  2:1481527/1‑100 (MQ=255)
                                                                                       cTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATTCTcgc          <  2:963288/100‑1 (MQ=255)
                                                                                       cTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATTCTcgc          >  1:1282440/1‑100 (MQ=255)
                                                                                       cTGTTCATCATTAAAAAAAACAGTATCATGATATTTTACGCAAGGGTAGCCCCCAAGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATTATcgc          >  2:442480/1‑100 (MQ=255)
                                                                                         gTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATTCTcgcgc        <  1:621470/100‑1 (MQ=255)
                                                                                          ttCATTATTAAAAAAAACAGTATCGTGAAATTTTACGCAAGGGTAGCCACCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATTCTcgcgcg       <  1:143946/100‑1 (MQ=255)
                                                                                          ttCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATTCTcgcgcg       <  2:754997/100‑1 (MQ=255)
                                                                                             atcatTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATTCTCGCGCGAAt    <  2:24457/100‑1 (MQ=255)
                                                                                             atcatTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATTCTCGCGCGAAt    <  1:1246040/100‑1 (MQ=255)
                                                                                             atcatTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATTCTCGCGCGAAt    >  2:12225/1‑100 (MQ=255)
                                                                                             atcatTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATTCTCGCGCGAAt    <  2:18743/100‑1 (MQ=255)
                                                                                              tcatTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATTCTCGCGCGAAt  <  1:912855/100‑2 (MQ=255)
                                                                                              tcatTAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATTCTCGCGCGAATGc  >  2:932213/1‑100 (MQ=255)
                                                                                                   |                                                                                                
ATCTGATGAAGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATTCTCGCGCGAATGC  >  NC_000913/4494778‑4494971

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: