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.A97.3% 390.9 / 3.2 113coding (254/564 nt)idnKD‑gluconate kinase, thermosensitive
Reads supporting (aligned to +/- strand):  ref base . (1/2);  new base A (54/56);  total (55/58)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.77e-01

Only 100 of 128 total aligned reads displayed.

ATCTGATGAAGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATTCTCGCGCGAAT  >  NC_000913/4494778‑4494969
                                                                                                    |                                                                                              
  cTGATGAAGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTaaa                                                                                             <  2:1328751/100‑1 (MQ=255)
   tgattaAGATCGACTTCCCTGGCTGTAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAATAATGAAACAGGATTTATTGTCTGTTCATCATTaaaa                                                                                            <  2:496990/100‑1 (MQ=255)
   tgatgaAGATCGACTTCCCTGGCTGGAACTCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTaaaa                                                                                            <  2:1589770/100‑1 (MQ=255)
   tgatgaAGATCGACTTCCCTGGCTGGAACGTTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTaaaa                                                                                            <  1:1857935/100‑1 (MQ=255)
   tgatgaAGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTaaaa                                                                                            >  1:1312193/1‑100 (MQ=255)
   tgatgaAGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTaaaa                                                                                            >  2:2162769/1‑100 (MQ=255)
     atgaAGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTaaaaaa                                                                                         >  1:1063918/1‑100 (MQ=255)
      tgaAGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTaaaaa                                                                                           >  1:1995856/1‑98 (MQ=255)
         aGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTaaaaaaaaca                                                                                     >  2:2413143/1‑97 (MQ=255)
         aGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTaaaaaaaaca                                                                                     >  1:868379/1‑97 (MQ=255)
          gATCGACTTCCATGGCTGGAAAGCTTAAATAATGCTTAATACCGTCTTTATAAAAAGAATGAAACAGGAATTATTGTCTGTTCATCATTaaaaaaaa                                                                                       >  2:1476105/1‑96 (MQ=255)
             cGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTAt                                                                                  <  2:790098/100‑1 (MQ=255)
              gACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATc                                                                                 >  1:1584864/1‑100 (MQ=255)
              gACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGGCTGTTCATCATTAAAAAAAACAGTAt                                                                                  >  2:2379125/1‑99 (MQ=255)
                cTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGt                                                                               >  2:711475/1‑100 (MQ=255)
                 ttCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTg                                                                              <  1:1824832/100‑1 (MQ=255)
                 ttCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTg                                                                              <  1:2005255/100‑1 (MQ=255)
                   cccTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTAATGTCTATTCATCATTAAAAAAAAAAGTATAGTGat                                                                            >  1:1514600/1‑100 (MQ=255)
                   cccTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGAATTTATTGTCTGTTCATCATTAAAAAAACAGTATCGTGata                                                                           >  1:740391/1‑100 (MQ=255)
                    ccTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGata                                                                           >  2:637449/1‑100 (MQ=255)
                    ccTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGata                                                                           >  2:750846/1‑100 (MQ=255)
                    ccTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGat                                                                            >  2:455871/1‑99 (MQ=255)
                         ctggAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTa                                                                      >  2:2331099/1‑100 (MQ=255)
                         ctggAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTa                                                                      <  1:1696473/100‑1 (MQ=255)
                          tggAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTAc                                                                     <  1:1031263/100‑1 (MQ=255)
                          tggAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTAc                                                                     >  1:1085079/1‑100 (MQ=255)
                            gAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGAATCGTGATATTTTACGc                                                                   >  1:274100/1‑100 (MQ=255)
                             aaCGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCa                                                                  <  1:1930014/100‑1 (MQ=255)
                             aaCGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCa                                                                  >  1:2633104/1‑100 (MQ=255)
                               cGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAg                                                                >  2:1017283/1‑100 (MQ=255)
                               cGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAg                                                                <  1:479626/100‑1 (MQ=255)
                                gCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAgg                                                               <  2:92957/100‑1 (MQ=255)
                                  ttAAATGATGCTTCATACATTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAggg                                                              <  1:523653/99‑1 (MQ=255)
                                  ttAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAggg                                                              <  1:2495840/99‑1 (MQ=255)
                                   tAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACCGTATCGTGATATTTTACGCAAgg                                                               >  2:744421/1‑97 (MQ=255)
                                    aaaTGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTan                                                           <  1:686250/100‑2 (MQ=255)
                                      atgatgCTTCATACATTCTTTATAAAAAGAATGAAACAGTATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGcc                                                         <  1:1126820/100‑1 (MQ=255)
                                       tgatgCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGccc                                                        <  2:2446649/100‑1 (MQ=255)
                                        gatgCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGcccc                                                       >  2:1645407/1‑100 (MQ=255)
                                        gatgCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGcccc                                                       >  1:962990/1‑100 (MQ=255)
                                        gatgCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGccc                                                        >  2:1902477/1‑99 (MQ=255)
                                           gCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCAt                                                    >  2:547751/1‑100 (MQ=255)
                                           gCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCAt                                                    <  1:1485078/100‑1 (MQ=255)
                                           gCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCAt                                                    >  2:1157956/1‑100 (MQ=255)
                                             ttCATACAGTCTTTATAAAAAGAATGAAACAGTATTTATTGTTTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCa                                                     <  2:1173129/97‑1 (MQ=255)
                                             ttCATACAGTCTTTATAAAAAGAATGAAACAGGATTTGTTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGt                                                  >  2:1412665/1‑100 (MQ=255)
                                             ttCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGt                                                  >  2:310765/1‑100 (MQ=255)
                                                aTACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCa                                               <  2:537752/100‑1 (MQ=255)
                                                aTACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCa                                               <  2:148919/100‑1 (MQ=255)
                                                aTACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCa                                               <  1:2441162/100‑1 (MQ=255)
                                                 tACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCAt                                              >  2:1728039/1‑100 (MQ=255)
                                                  aCAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCAtt                                             <  1:1899628/100‑1 (MQ=255)
                                                  aCAGTATTTATAAAAATAATGAAACAGGATTTATTGTCTGTTCATCAATAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTc                                                <  2:1066781/97‑1 (MQ=255)
                                                     gTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTcc                                          <  2:274098/100‑1 (MQ=255)
                                                       cTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCctc                                        >  1:2292352/1‑100 (MQ=255)
                                                        tttATAAAAAGAATGAAACAGGATTTATTGTCTGTTTATTATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCctct                                       <  1:978747/100‑1 (MQ=255)
                                                        tttATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAACCCCCATATTCATTTCctct                                       >  1:967105/1‑100 (MQ=255)
                                                           ataAAAAGAATGATACAGGATTCACTGTCTGATCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCACCCATGTTCATTTCTTCTGGt                                    >  2:1087503/1‑100 (MQ=255)
                                                           ataAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGt                                    >  2:1351150/1‑100 (MQ=255)
                                                            taAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTTATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGtt                                   <  1:2314783/100‑1 (MQ=255)
                                                            taAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGtt                                   <  2:1950613/100‑1 (MQ=255)
                                                            taAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGtt                                   <  1:327695/100‑1 (MQ=255)
                                                            taAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGtt                                   <  1:357039/100‑1 (MQ=255)
                                                            taAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGtt                                   >  1:650273/1‑100 (MQ=255)
                                                            taAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGAAAGGTTAACCCCCATATTCATATCCTCTGGtt                                   >  1:651627/1‑100 (MQ=255)
                                                             aaaaaGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAACCACCATGTTCATTTCCTCTGGTTa                                  >  1:1194973/1‑100 (MQ=255)
                                                              aaaaGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAg                                 <  2:1765533/100‑1 (MQ=255)
                                                              aaaaGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAg                                 <  2:1208031/100‑1 (MQ=255)
                                                              aaaaGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAg                                 >  2:489501/1‑100 (MQ=255)
                                                               aaaGAATGAAACAGGATTTAGTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGa                                <  2:2204692/100‑1 (MQ=255)
                                                                aaGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGAt                               <  1:250801/100‑1 (MQ=255)
                                                                    aTGAAACAGGATTTATTGTCTGTTCATCATTAAAAATAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCg                           <  1:1820821/100‑1 (MQ=255)
                                                                     tGAAACATGATGTATAGTCTGTTGATCATTAAAAAAAAGAGTATCGTTATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGc                            >  2:386933/1‑98 (MQ=255)
                                                                     tGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGa                          >  2:1497304/1‑100 (MQ=255)
                                                                      gAAACAGGATTTATTGTCTGTAAATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCACATGTTCATTTACTCTGGTTAGATGGCGAc                         >  2:1530863/1‑100 (MQ=255)
                                                                       aaaCAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACt                        >  1:1954746/1‑100 (MQ=255)
                                                                        aaCAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTa                       <  1:2572054/100‑1 (MQ=255)
                                                                            ggATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGaa                   >  2:1877893/1‑100 (MQ=255)
                                                                            ggATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGaa                   <  1:2022922/100‑1 (MQ=255)
                                                                             gATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGaa                   <  1:1503777/99‑1 (MQ=255)
                                                                              aTTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAAc                 <  2:1085069/100‑1 (MQ=255)
                                                                              aTTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAAc                 <  2:1916/100‑1 (MQ=255)
                                                                              aTTTATTCTATGTTCATAATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGaaa                  <  2:115779/99‑1 (MQ=255)
                                                                                ttATTGTCTGTTAATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAAc                 <  2:467971/98‑1 (MQ=255)
                                                                                ttATTGTCCTTGCATCATTAAAAAAAACAGTATCGTGATATTTTATGTAAGTGTAGCCCCCACGTTCATTTCCTCTGTTTAGATGGCGACTATGAAAc                 <  2:2579361/98‑1 (MQ=255)
                                                                                    tGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATtct           <  1:1681534/100‑1 (MQ=255)
                                                                                     gTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATtctc          >  2:1272683/1‑100 (MQ=255)
                                                                                      tCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATTCTcg         <  1:507064/100‑1 (MQ=255)
                                                                                      tCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGGTTATTTCTTCTGGTTTGATGGGGACTATGAAACTATtct           >  2:1581281/1‑98 (MQ=255)
                                                                                       cTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATTCTcg         <  1:783577/99‑1 (MQ=255)
                                                                                        tGTTCATCATTAAAAAAAAAAGTATCGTGATATTTTACTCAAGTGTATCCACCATGTTCATTTACTCTGTTTAGATTTCGACTATGAAACTATTCTcg         >  2:81054/1‑98 (MQ=255)
                                                                                        ngTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCAAGTTCATTTCCTCTGGATAGATGGCGACTATGAAACTACtctc          >  2:2570092/2‑97 (MQ=255)
                                                                                         gTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATTCTcgcgc      >  2:40696/1‑100 (MQ=255)
                                                                                         gTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATTCTcgcgc      >  1:880522/1‑100 (MQ=255)
                                                                                         gTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATTCTcgcgc      >  2:1747859/1‑100 (MQ=255)
                                                                                         gTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATTCTcgcgc      >  2:75777/1‑100 (MQ=255)
                                                                                          ttCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATTCTcgcgcg     >  1:1401491/1‑100 (MQ=255)
                                                                                             atcatTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATTCTCGCGCGAAt  <  2:2066860/100‑1 (MQ=255)
                                                                                             atcatTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATTCTCGCGCGAAt  <  2:2588850/100‑1 (MQ=255)
                                                                                             atcatTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATTCTCGCGCGAAt  <  1:1877910/100‑1 (MQ=255)
                                                                                                    |                                                                                              
ATCTGATGAAGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATTCTCGCGCGAAT  >  NC_000913/4494778‑4494969

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: