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

ATCTGATGAAGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑AA‑AAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATTCTCGCGCGAATGC  >  NC_000913/4494778‑4494971
                                                                                                   |                                                                                                
aTCTGATGAAGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑a                                                                                                 <  1:1904193/100‑1 (MQ=255)
 tCTGATGAAGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑aa                                                                                                <  1:968656/100‑1 (MQ=255)
  cTGATGAAGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑aa‑a                                                                                              >  2:92381/1‑100 (MQ=255)
  cTGATGAAGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑aa‑a                                                                                              >  1:1693971/1‑100 (MQ=255)
   tgatgaAGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑aa‑aa                                                                                             >  2:1552722/1‑100 (MQ=255)
    gatgaAGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑aa‑aaa                                                                                            >  1:804790/1‑100 (MQ=255)
      tgaAGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑aa‑aaaaa                                                                                          <  1:733079/100‑1 (MQ=255)
      tgaAGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑aa‑aaaaa                                                                                          >  1:1639346/1‑100 (MQ=255)
        aaGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑aa‑aaaaaac                                                                                        <  1:1648624/100‑3 (MQ=255)
        aaGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑aa‑aaaaaac                                                                                        <  1:2109983/100‑3 (MQ=255)
         aGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑aa‑aaaaaaca                                                                                       <  1:2054855/100‑4 (MQ=255)
          gATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑aa‑aaaaaacag                                                                                      >  2:1756707/1‑96 (MQ=255)
          gATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑aa‑aaaaaacag                                                                                      <  2:1954943/100‑5 (MQ=255)
          gATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑aa‑aaaaaac                                                                                        <  1:702812/98‑3 (MQ=255)
              gACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAACAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑AACAAAAACAGTATc                                                                                   >  2:98199/1‑100 (MQ=255)
              gACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAA‑AAAAACAGTATc                                                                                   >  2:1805160/1‑100 (MQ=255)
               aCTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAA‑AAAAACAGTATCg                                                                                  <  2:656936/100‑1 (MQ=255)
               aCTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAA‑AAAAACAGTATCg                                                                                  >  1:792504/1‑100 (MQ=255)
                cTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAA‑AAAAACAGTATCGt                                                                                 <  1:406776/100‑1 (MQ=255)
                  tCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAA‑AAAAACAGTATCGTg                                                                                >  1:1907061/1‑99 (MQ=255)
                  tCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAA‑AAAAACAGTATCGTGa                                                                               >  1:1751430/1‑100 (MQ=255)
                  tCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAA‑AAAAACAGTATCGTGa                                                                               >  1:381529/1‑100 (MQ=255)
                  tCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTATGTTCATCATTAAA‑AAAAACAGTATAGTGa                                                                               >  1:729259/1‑100 (MQ=255)
                    ccTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAA‑AAAAACAGTATCGTGata                                                                             >  2:1776384/1‑100 (MQ=255)
                    ccTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAA‑AAAAACAGTATCGCGat                                                                              >  2:34701/1‑99 (MQ=255)
                    ccTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAA‑AAAAAAAGTATCGTGata                                                                             >  2:535077/1‑100 (MQ=255)
                    ccTGGATGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAATAATGAAACAGTATTTATTGTCTGTTCATAATTAAA‑AAAAACAGTATCGt                                                                                 <  2:181378/96‑1 (MQ=255)
                        gctggAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAA‑AAAAACAGTATCGTGATAtttt                                                                         <  1:944460/100‑1 (MQ=255)
                        gattaaaCGCATAAATGATGCTCCATACATTCTTTATAAAAAGAATGAAACAGTATTTATTGTCTGTTCATCATTAAA‑AAAAACAGTATCGTGATAttt                                                                          <  2:1957395/94‑1 (MQ=255)
                          tggAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAA‑AAAAACAGTATCGTGATATTTTAc                                                                       <  1:734019/100‑1 (MQ=255)
                            gAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAA‑AAAAACAGTATCGTGATATTTTACGc                                                                     <  1:1545666/100‑1 (MQ=255)
                             aaCGCTTAAATGATGCTTCATACAGTCTTTATAACAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAA‑AAAAACAGTATCGTGATATTTTACGCa                                                                    >  1:2190650/1‑100 (MQ=255)
                               cGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAA‑AAAAACAGTATCGTGATATTTTACGCAAg                                                                  >  1:274174/1‑100 (MQ=255)
                               cGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑AA‑AAAAACAGTATCGTGATATTTTACGCa                                                                    >  2:1824517/1‑97 (MQ=255)
                                gCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTGATCATTAAA‑AAAAACAGTATCGTGATATTTTACGCAAg                                                                  >  2:1640521/1‑99 (MQ=255)
                                 cTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAA‑AAAAACAGTATCGTGATATTTTACGCAAggg                                                                <  1:975956/100‑1 (MQ=255)
                                 cTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAA‑AAAAACAGTATCGTGATATTTTACGCAAggg                                                                <  2:1751419/100‑1 (MQ=255)
                                 cTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAA‑AAAAACAGTATCGTGATATTTTACGCAAg                                                                  <  1:535076/98‑1 (MQ=255)
                                  ttAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAA‑AAAAACAGTATCGTGTTATGTTACGCAAg                                                                  >  2:1162552/1‑97 (MQ=255)
                                  ttAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAA‑AAAAACAGTATCGTGATATTTTACGCAAGGGt                                                               <  1:1669637/100‑1 (MQ=255)
                                   tAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAA‑AAAAACAGTATCGTGATATTTTACGCAAGGGTa                                                              >  1:859181/1‑100 (MQ=255)
                                   tAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAA‑AAAAACAGTATCGTGATATTTTACGCAAGGGTa                                                              <  1:555824/100‑1 (MQ=255)
                                      atgatgCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAA‑AAAAACAGTATCGTGATATTTTACGCAAGGGTAGcc                                                           <  1:644751/100‑1 (MQ=255)
                                        gatgggtCATACAGTCTTTATAAAAAAAATGATACAGGATTTATTGTATGTTTATCATTAAA‑AAAAACAGTATCTTTATATTTTACGCAAGGGTAGcc                                                           <  2:1884474/92‑1 (MQ=255)
                                         atgCTTCATACAGTTTTTATAAAAAGAATGAAACAAGACTTATTGTCTGTTCATCATTAAA‑AAAAACAGTATCGTAATATTTTATGCAAGGATat                                                             >  2:112893/1‑94 (MQ=255)
                                           gCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAA‑AAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCAt                                                      >  2:617419/1‑100 (MQ=255)
                                           gCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAA‑AAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCAt                                                      >  1:1735320/1‑100 (MQ=255)
                                               cATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAA‑AAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTc                                                  <  2:113466/100‑1 (MQ=255)
                                                ttccgatCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAA‑AAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCa                                                 <  1:2099755/94‑1 (MQ=255)
                                                aTACAGTCTTTATAAAAAGAATGAAACAGTATTTATTGTCTGTTCATCATTAAA‑AAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCa                                                 <  2:720220/100‑1 (MQ=255)
                                                aTACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAA‑AAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCa                                                 <  1:146700/100‑1 (MQ=255)
                                                 tACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAA‑AAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCAt                                                >  2:1139407/1‑100 (MQ=255)
                                                 tACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAA‑AAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCAt                                                >  1:429797/1‑100 (MQ=255)
                                                      tCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑AA‑AAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTcc                                            <  2:804790/98‑1 (MQ=255)
                                                        tttATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAA‑AAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCctct                                         <  1:560146/100‑1 (MQ=255)
                                                         ttATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAA‑AAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCaagatcggaa                                        >  2:2099741/1‑91 (MQ=255)
                                                         ttATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAA‑AAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTg                                        <  2:313225/100‑1 (MQ=255)
                                                          tataAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAA‑AAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTgg                                       >  2:1901726/1‑100 (MQ=255)
                                                          tataAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAA‑AAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTgg                                       <  1:400973/100‑1 (MQ=255)
                                                           ataAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAA‑AAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGt                                      >  2:1155988/1‑100 (MQ=255)
                                                           agaaaaaGAATGAAACAGGATTTATTGTCTGTTCATCATTAAA‑AAAAACAGTATCGAGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGt                                      <  1:655843/98‑1 (MQ=255)
                                                            taAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAA‑AAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGtt                                     <  1:961310/100‑1 (MQ=255)
                                                            taAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAA‑AAAAACAGAATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGtt                                     >  2:556804/1‑100 (MQ=255)
                                                             aaaaaGAATGAAACAGGATTTATTGTCTGTTCATCATTAAA‑AAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTa                                    <  2:966242/100‑1 (MQ=255)
                                                              aaaaGAATGAAACAGGATTTATTGTCTGTTCATCATTAAA‑AAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAg                                   >  1:1328506/1‑100 (MQ=255)
                                                              aaaaGAATGAAACAGGATTTATTGTATGTTCATCATTAAA‑AAAAACAGTATCGTGATATTTTAGGAAAGGGTAGCACCCATGTTAATTTACTCTGGTTAg                                   >  2:932634/1‑100 (MQ=255)
                                                               aaaGAATGAAACAGGATTTATTGTCTGTTCATCATTAAA‑AAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGa                                  >  2:676742/1‑100 (MQ=255)
                                                               aaaGAATGAAACAGGATTTATTATCTGTTCATCATTAAA‑AAAAACAGTATCGTCATATTTTACGCAAGGATAGCCCCCATGTTCATTTCCTCTGGTTaca                                  >  2:2129627/1‑98 (MQ=255)
                                                                 aGAATGAAACAGGATTTATTGTCTGTTCATCATTAAA‑AAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGAt                                 >  2:854473/1‑99 (MQ=255)
                                                                   aaTGAAACAGGATTTATTGTCTGTTCATCATTAAA‑AAAAACATTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATgg                               <  2:440775/99‑1 (MQ=255)
                                                                   aaTGAAACAGGATTTATTGTCTGTTCATCATTAAA‑AAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGc                              <  1:722709/100‑1 (MQ=255)
                                                                    aTGAAACAGGATTTATTGTCTGTTCATCATTAAA‑AAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCg                             >  1:1983871/1‑100 (MQ=255)
                                                                    aTGAAACAGGATTTATTGTCTGTTCATCATTAAA‑AAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCg                             >  1:1384100/1‑100 (MQ=255)
                                                                     tGAAACAGGATTTATTGTCTGTTCATCATTAAA‑AAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGa                            >  1:1734639/1‑100 (MQ=255)
                                                                        aaCAGGATTTATTGTCTGTTCATCATTAAA‑AAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTa                         <  2:12102/100‑1 (MQ=255)
                                                                        aaCAGGATTTATTGTCTGTTCATCATTAAA‑AAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTa                         <  2:2204298/100‑1 (MQ=255)
                                                                              aTTTATTGTCTGTTCATCATTAAA‑AAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAAc                   <  1:866518/100‑1 (MQ=255)
                                                                               tttATTGTCTGTTCATCATT‑AA‑AAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAAc                   <  1:214071/98‑1 (MQ=255)
                                                                                    tGTCTGTTCATCATTAAA‑AAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATtct             <  2:949784/100‑1 (MQ=255)
                                                                                     gtatGTTTAACATTAAA‑AAAAACAGTATCGTGTTATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATtctc            <  2:476169/97‑1 (MQ=255)
                                                                                     gTCTGTTCATCATTAAA‑AAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGTAACTATtctc            >  2:1481527/1‑100 (MQ=255)
                                                                                       cTGTTCATCATTAAA‑AAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATTCTcgc          <  2:963288/100‑1 (MQ=255)
                                                                                       cTGTTCATCATTAAA‑AAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATTCTcgc          >  1:1282440/1‑100 (MQ=255)
                                                                                       cTGTTCATCATTAAA‑AAAAACAGTATCATGATATTTTACGCAAGGGTAGCCCCCAAGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATTATcgc          >  2:442480/1‑100 (MQ=255)
                                                                                         gTTCATCATTAAA‑AAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATTCTcgcgc        <  1:621470/100‑1 (MQ=255)
                                                                                          ttCATTATTAAA‑AAAAACAGTATCGTGAAATTTTACGCAAGGGTAGCCACCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATTCTcgcgcg       <  1:143946/100‑1 (MQ=255)
                                                                                          ttCATCATTAAA‑AAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATTCTcgcgcg       <  2:754997/100‑1 (MQ=255)
                                                                                             atcatTAAA‑AAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATTCTCGCGCGAAt    <  2:24457/100‑1 (MQ=255)
                                                                                             atcatTAAA‑AAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATTCTCGCGCGAAt    <  1:1246040/100‑1 (MQ=255)
                                                                                             atcatTAAA‑AAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATTCTCGCGCGAAt    >  2:12225/1‑100 (MQ=255)
                                                                                             atcatTAAA‑AAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATTCTCGCGCGAAt    <  2:18743/100‑1 (MQ=255)
                                                                                              tcatTAAA‑AAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATTCTCGCGCGAAtn   <  1:912855/100‑2 (MQ=255)
                                                                                              tcatT‑AA‑AAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATTCTCGCGCGAATGc  >  2:932213/1‑100 (MQ=255)
                                                                                                   |                                                                                                
ATCTGATGAAGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑AA‑AAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATTCTCGCGCGAATGC  >  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: ‑