Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,494,883 (A)7→8 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.A96.5% 297.5 / 2.7 86coding (254/564 nt)idnKD‑gluconate kinase, thermosensitive
Reads supporting (aligned to +/- strand):  ref base . (1/2);  new base A (46/37);  total (47/39)
Fisher's exact test for biased strand distribution p-value = 5.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GATGAAGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑‑AAA‑AAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATTCTCGCGCGAATG  >  NC_000913/4494782‑4494970
                                                                                               |                                                                                                 
gatgaAGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATGCAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑‑aaa‑aa                                                                                           >  1:1008431/1‑100 (MQ=255)
gatgaAGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑‑aaa‑aa                                                                                           >  2:808158/1‑100 (MQ=255)
gatgaAGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑‑aaa‑aa                                                                                           >  1:340694/1‑100 (MQ=255)
gatgaAGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑‑aaa‑aa                                                                                           >  1:340704/1‑100 (MQ=255)
gatgaAGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATAATT‑‑‑aaa‑aa                                                                                           >  1:1115249/1‑100 (MQ=255)
 atgaAGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑‑aaa‑aaa                                                                                          <  1:533182/100‑1 (MQ=255)
    aaGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑‑aaa‑aaaaac                                                                                       >  1:24414/1‑98 (MQ=255)
    aaGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑‑aaa‑aaaaac                                                                                       >  1:212599/1‑98 (MQ=255)
    aaGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑‑aaa‑aaaaac                                                                                       <  1:1782524/100‑3 (MQ=255)
     aGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑‑aaa‑aaaaaca                                                                                      <  1:322703/100‑4 (MQ=255)
     aGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑‑aaa‑aaaaaca                                                                                      >  1:1394943/1‑97 (MQ=255)
      gATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAATAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑‑aaa‑aaaaacag                                                                                     >  1:726466/1‑96 (MQ=255)
       aTCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑‑aaa‑aaaaacagt                                                                                    >  2:1559373/1‑95 (MQ=255)
        tCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑‑AAA‑AAAACAGTa                                                                                    >  2:1042570/1‑100 (MQ=255)
           aCTTCCCTGGCTGGAACGCTTAAATGATGCTTCATTCAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑‑AAA‑AAAACAGTATCg                                                                                 >  1:1422507/1‑100 (MQ=255)
             ttCCCTGGCTGGATCGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑‑AAA‑AAAACAGTATCGt                                                                                <  2:1091776/99‑1 (MQ=255)
             ttCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑‑AAA‑AAAACAGTATCGTGa                                                                              <  2:1551347/100‑1 (MQ=255)
               cccTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑‑AAA‑AAAACAGTATCGTGat                                                                             <  1:337383/100‑1 (MQ=255)
                 ctggctggAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑‑AAA‑AAAACAGTATCGTGatat                                                                           <  2:575386/100‑1 (MQ=255)
                   ggctggAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑‑AAA‑AAAACAGTATCGTGATAttt                                                                         >  2:613272/1‑100 (MQ=255)
                     ctggAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑‑AAA‑AAAACAGTATCGTGATATTTTa                                                                       >  2:775219/1‑100 (MQ=255)
                             cTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCTTCATT‑‑‑AAA‑AAAACAGTATCGTGATATTTTACGCAAggg                                                               <  1:1115891/99‑1 (MQ=255)
                             cTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑‑AAA‑AAAACAGTATCGTGATATTTTACGCAAgtg                                                               >  1:1181760/1‑98 (MQ=255)
                             cTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑‑AAA‑AAAACAGTATCGTGATATTTTACGCAAggg                                                               >  1:137317/1‑100 (MQ=255)
                             cTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑‑AAA‑AAAACAGTATCGTGATATTTTACGCAAggg                                                               <  2:1619981/100‑1 (MQ=255)
                             cTTAAATGATGCTTCATACAGTCTTACTAAAAAGCATGAAACAGGCTTTATTGTCTGTTCATCATTA‑‑AAA‑AAAACAATATCGTGATATTTTACGCAAg                                                                 >  2:106027/1‑98 (MQ=255)
                               tAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑‑AAA‑AAAACAGTATCGTGATATTTTACGCAAGGGTa                                                             >  1:1786684/1‑100 (MQ=255)
                                aaaTGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑‑AAA‑AAAACAGTATCGTGATATTTTACGCAAGGGTAg                                                            <  2:1364391/100‑1 (MQ=255)
                                   tgatgCTTCATACAGTCGTTATAAAAAGAATGAAACAGTATTTATTGTCTGTTCATCATTA‑‑AAA‑AAAACAGTATCGTGATATTTTACGCAAGGGTAGccc                                                         <  1:1413269/100‑1 (MQ=255)
                                    gatgCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑‑AAA‑AAAACAGTATCGTGATATTTTACGCAAGGGTAGcccc                                                        >  2:1455787/1‑100 (MQ=255)
                                    gatgCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑‑AAA‑AAAACAGTATCGTGATATTTTACGCAAGGGTAGcccc                                                        >  2:694816/1‑100 (MQ=255)
                                     atgCTTCATATAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑‑AAA‑AAAACAGTATCGTGATATTTTACGCAAGGGTAGccccc                                                       <  1:696002/100‑1 (MQ=255)
                                      tgCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑‑AAA‑AAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCa                                                      <  2:1188133/100‑1 (MQ=255)
                                       gCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑‑AAA‑AAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCAt                                                     >  1:1361431/1‑100 (MQ=255)
                                        cTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑‑AAA‑AAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATg                                                    <  2:1008440/100‑1 (MQ=255)
                                        cTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑‑AAA‑AAAACAGTATCGAGATATTTTACGCAAGGGTAGCCCCCATg                                                    >  2:482961/1‑100 (MQ=255)
                                         ttCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑‑AAA‑AAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGt                                                   >  1:1553374/1‑100 (MQ=255)
                                            aTACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑‑AAA‑AAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCa                                                >  1:1070621/1‑100 (MQ=255)
                                             tACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑‑AAA‑AAAACAGTATTGTGATATTTTACGCAATTGTAGCCCCCATGTTCAt                                               >  1:1178225/1‑100 (MQ=255)
                                             tACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑‑AAA‑AAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCAt                                               >  1:172461/1‑100 (MQ=255)
                                             tACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑‑AAA‑AAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCAt                                               >  1:1233938/1‑100 (MQ=255)
                                                aGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAA‑AAAACAGTATCGTGAAACTTTACGCAAGGGAAACCCCCATGTTCAtt                                              >  2:1769321/1‑100 (MQ=255)
                                                   cTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑‑AAA‑AAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCctc                                         <  2:1216861/100‑1 (MQ=255)
                                                   cTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑‑AAACAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCctc                                         >  1:1434929/1‑100 (MQ=255)
                                                    tttATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATAATTA‑‑AAA‑AAAACAGTATCGTGATATTTTACGCAAGGGTAGCCACCATGTTCATTTCctct                                        >  1:777748/1‑100 (MQ=255)
                                                      tataAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑‑AAA‑AAAACAGTATCGTGATATTTTACGCATGGGTAGCCCCCATGTTCATTTCctct                                        >  1:399116/1‑98 (MQ=255)
                                                      tataAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑‑AAA‑AAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCTCATGTTCATTTCCTCTgg                                      <  1:884745/100‑1 (MQ=255)
                                                      tataAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑‑AAA‑AAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTgg                                      <  2:1442258/100‑1 (MQ=255)
                                                      tataAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑‑AAA‑AAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTgg                                      <  2:629226/100‑1 (MQ=255)
                                                      tataAAAAGAATGAAACAGGATTTATTGTCTGTTCATAATTA‑‑AAA‑AAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTgg                                      <  1:871046/100‑1 (MQ=255)
                                                       ataAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑‑AAA‑AAAACAGTATCGTGATATTTTACGCAAGGGTAACCCCCATAATCAttaacac                                         >  2:65910/1‑91 (MQ=255)
                                                        taAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑‑AAA‑AAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGtt                                    <  1:613268/100‑1 (MQ=255)
                                                        taAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑‑AAA‑AAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGtt                                    >  2:131385/1‑100 (MQ=255)
                                                          aaaaGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑‑AAA‑AAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGtt                                    >  1:157708/1‑98 (MQ=255)
                                                          aaaaGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑‑AAA‑AAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAg                                  <  2:1720815/100‑1 (MQ=255)
                                                          aaaaGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑‑AAA‑AAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAg                                  <  2:75142/100‑1 (MQ=255)
                                                          aaaaGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑‑AAA‑AAAACAGTATCGTGATATTTTACGCAAGGGTAGCCACCATGTTCATTTCCTCTGGTTAg                                  >  1:1372808/1‑100 (MQ=255)
                                                           aaaGAATGAAACAGGATTTATTGTCTGTTCATCATTAA‑AAA‑AAAACAGTATCGTGATATTTTACGCAAGTGTAGCACCCATGCTCATGTCCTCTg                                       >  1:139677/1‑95 (MQ=255)
                                                           aaaGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑‑AAA‑AAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTGTGGTTAGa                                 >  1:1281832/1‑100 (MQ=255)
                                                           aaaGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑‑AAA‑AAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGa                                 >  2:1011680/1‑100 (MQ=255)
                                                           aaaGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑‑AAA‑AAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGa                                 <  2:607048/100‑1 (MQ=255)
                                                           aaaGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑‑AAA‑AAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGa                                 >  1:1351017/1‑100 (MQ=255)
                                                             aGAATGAAACAGGATTTTTTGTCTGTTCATCTTTA‑‑AAA‑AAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATg                               <  2:55189/100‑1 (MQ=255)
                                                             aGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑‑AAA‑AAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATg                               <  1:569171/100‑1 (MQ=255)
                                                              gAATGAAACAGGATTTATTGTCTGTTCATCATTA‑‑AAA‑AAAACAGTATCGGGATATTTTACGCAAGGGTAGCCCCCATTTTTATTTCCTCTGGTTAGATgg                              <  1:745601/100‑1 (MQ=255)
                                                               aaTGAAACAGGATTTATTGTCTGTTCATCATTA‑‑AAA‑AAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATgg                              <  2:1548077/99‑1 (MQ=255)
                                                               aaTGAAACAGGATTTATTGTCTGTTCATCATTA‑‑AAA‑AAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGc                             >  1:977142/1‑100 (MQ=255)
                                                               aaTGAAACAGGATTTATTGTCTGTTCATCATTA‑‑AAA‑AAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGc                             <  2:212600/100‑1 (MQ=255)
                                                                cttccgatctaTTTATTGTCTGTTCATCATTA‑‑AAA‑AAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCg                            <  2:1115192/90‑1 (MQ=255)
                                                                aTGAAACAGGATTTATTGTCTGTTCATCATTAA‑AAA‑AAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCAAGTTCATTTCCTCTGGTTGGATgg                              >  1:227611/1‑99 (MQ=255)
                                                                aTGAAACAGGATTTATTGTCTGTCCATCATTA‑‑AAA‑AAAACAGTATCGTGATATTTTACGCAAGGGTAGACCCCATGTTCATTTCCTCTGGTTAGATGGCg                            >  1:902339/1‑100 (MQ=255)
                                                                  gAAACAGGATTTATTGTCTGTTCATCATTA‑‑AAA‑AAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGAc                          <  1:262703/100‑1 (MQ=255)
                                                                  gAAACAGGATTTATTGTCTGTTCATCATTA‑‑AAA‑AAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGAc                          >  2:1300141/1‑100 (MQ=255)
                                                                  gAAACAGGATTTATTATCTATTCATCATTA‑‑AAA‑AAAACAGTATCGTTATATTTTACACAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGa                           >  1:230043/1‑99 (MQ=255)
                                                                   aaaCAGGATTTATTGTCTGTTCATCATTA‑‑AAA‑AAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACt                         <  1:131385/100‑1 (MQ=255)
                                                                     aCAGGATTTATTGTCTGTTCATCATTA‑‑AAA‑AAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTAt                       <  2:1153208/100‑1 (MQ=255)
                                                                      cAGGATTTATTGTCTTTTCATTATTA‑‑AAA‑AAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATg                      <  2:884175/100‑1 (MQ=255)
                                                                       aGGATTTATTGTCTGTTCATCATTA‑‑AAA‑AAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGa                     >  2:1084574/1‑100 (MQ=255)
                                                                        ggATTTATTGTCTGTTCATCATTA‑‑AAA‑AAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGaa                    <  2:1354907/100‑1 (MQ=255)
                                                                        ggATTTATTGTCTGTTCATCATTA‑‑AAA‑AAAACAGTATCGTGATATTTTACGCAAGGGTAACCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGaa                    >  1:1812590/1‑100 (MQ=255)
                                                                          aTTTATTGTCTGTTCATCATTA‑‑AAA‑AAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGagatcggaag                  >  1:1115182/1‑91 (MQ=255)
                                                                          aTTTATTGTCTGTTCATCATTA‑‑AAA‑AAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAAc                  <  1:1559965/100‑1 (MQ=255)
                                                                            ttATTGTCTGTTCATCATTA‑‑AAA‑AAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCATCTGGTTAGATGGCGACTATGAAACTa                >  2:572714/1‑100 (MQ=255)
                                                                              aTTGTCTGTTCATCATTA‑‑AAA‑AAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTAtt              >  1:1040163/1‑100 (MQ=255)
                                                                              aTTGTCTGTTCATCATTA‑‑AAA‑AAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTAtt              >  1:1406803/1‑100 (MQ=255)
                                                                               ttGTCTGTTCATCATTA‑‑AAA‑AAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATtc             <  1:1042560/100‑1 (MQ=255)
                                                                                tGTCTGTTCATCATTA‑‑AAA‑AAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATtct            >  2:19848/1‑100 (MQ=255)
                                                                                  tCTGTTCATCATTA‑‑AAA‑AAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATTCTcg          <  2:1618589/100‑1 (MQ=255)
                                                                                   cTGTTCATCATTA‑‑AAA‑AAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATTCTcgc         <  1:189501/100‑1 (MQ=255)
                                                                                    tGTTCATCATTA‑‑AAA‑AAAACAGTATCGTGATATTTTACGCAAGGGTAGCCGCCATGTTCATTTCCTCTGGTTAGATGGGGACTATGAAACTATTCTcgcg        >  1:908616/1‑100 (MQ=255)
                                                                                    tGTTCATCATTA‑‑AAA‑AAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATTCTcgcg        <  1:17224/100‑1 (MQ=255)
                                                                                      ttCATCATTA‑‑AAA‑AAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATTCTcgcgcg      <  2:621351/100‑1 (MQ=255)
                                                                                       tcatcaTTA‑‑AAA‑AAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATTCTCGCGCGa     >  1:1671109/1‑100 (MQ=255)
                                                                                        catcatTA‑‑AAA‑AAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCTTCTGGTTAGATGGCGACTATGAAACTATTCTCGCGCGaa    >  1:1230769/1‑100 (MQ=255)
                                                                                        catcatTA‑‑AAA‑AAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATTCTCGCGCGaa    >  2:553369/1‑100 (MQ=255)
                                                                                          tcatTA‑‑AAA‑AAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATTCTCGCGCGAATg  >  1:1809403/1‑100 (MQ=255)
                                                                                          tcatTA‑‑AAA‑AAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATTCTCGCGCGAATg  <  2:1180814/100‑1 (MQ=255)
                                                                                          tcatTA‑‑AAA‑AAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATTCTCGCGCGAATg  <  2:137317/100‑1 (MQ=255)
                                                                                          tcatTA‑‑AAA‑AAAACAGTATCGTGATATTTTACGCAAGGGTAGCCACCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATTCTCGCGCGAATg  >  1:966160/1‑100 (MQ=255)
                                                                                               |                                                                                                 
GATGAAGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑‑AAA‑AAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATTCTCGCGCGAATG  >  NC_000913/4494782‑4494970

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