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‑‑‑AAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATTCTCGCGCGAATG  >  NC_000913/4494782‑4494970
                                                                                               |                                                                                                 
gatgaAGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATGCAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑‑aaaaa                                                                                           >  1:1008431/1‑100 (MQ=255)
gatgaAGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑‑aaaaa                                                                                           >  2:808158/1‑100 (MQ=255)
gatgaAGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑‑aaaaa                                                                                           >  1:340694/1‑100 (MQ=255)
gatgaAGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑‑aaaaa                                                                                           >  1:340704/1‑100 (MQ=255)
gatgaAGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATAATT‑‑‑aaaaa                                                                                           >  1:1115249/1‑100 (MQ=255)
 atgaAGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑‑aaaaaa                                                                                          <  1:533182/100‑1 (MQ=255)
    aaGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑‑aaaaaaaac                                                                                       >  1:24414/1‑98 (MQ=255)
    aaGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑‑aaaaaaaac                                                                                       >  1:212599/1‑98 (MQ=255)
    aaGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑‑aaaaaaaac                                                                                       <  1:1782524/100‑3 (MQ=255)
     aGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑‑aaaaaaaaca                                                                                      <  1:322703/100‑4 (MQ=255)
     aGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑‑aaaaaaaaca                                                                                      >  1:1394943/1‑97 (MQ=255)
      gATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAATAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑‑aaaaaaaacag                                                                                     >  1:726466/1‑96 (MQ=255)
       aTCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑‑aaaaaaaacagt                                                                                    >  2:1559373/1‑95 (MQ=255)
        tCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑‑AAAAAAACAGTa                                                                                    >  2:1042570/1‑100 (MQ=255)
           aCTTCCCTGGCTGGAACGCTTAAATGATGCTTCATTCAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑‑AAAAAAACAGTATCg                                                                                 >  1:1422507/1‑100 (MQ=255)
             ttCCCTGGCTGGATCGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑‑AAAAAAACAGTATCGt                                                                                <  2:1091776/99‑1 (MQ=255)
             ttCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑‑AAAAAAACAGTATCGTGa                                                                              <  2:1551347/100‑1 (MQ=255)
               cccTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑‑AAAAAAACAGTATCGTGat                                                                             <  1:337383/100‑1 (MQ=255)
                 ctggctggAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑‑AAAAAAACAGTATCGTGatat                                                                           <  2:575386/100‑1 (MQ=255)
                   ggctggAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑‑AAAAAAACAGTATCGTGATAttt                                                                         >  2:613272/1‑100 (MQ=255)
                     ctggAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑‑AAAAAAACAGTATCGTGATATTTTa                                                                       >  2:775219/1‑100 (MQ=255)
                             cTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCTTCATT‑‑‑AAAAAAACAGTATCGTGATATTTTACGCAAggg                                                               <  1:1115891/99‑1 (MQ=255)
                             cTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑‑AAAAAAACAGTATCGTGATATTTTACGCAAgtg                                                               >  1:1181760/1‑98 (MQ=255)
                             cTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑‑AAAAAAACAGTATCGTGATATTTTACGCAAggg                                                               >  1:137317/1‑100 (MQ=255)
                             cTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑‑AAAAAAACAGTATCGTGATATTTTACGCAAggg                                                               <  2:1619981/100‑1 (MQ=255)
                             cTTAAATGATGCTTCATACAGTCTTACTAAAAAGCATGAAACAGGCTTTATTGTCTGTTCATCATTA‑‑AAAAAAACAATATCGTGATATTTTACGCAAg                                                                 >  2:106027/1‑98 (MQ=255)
                               tAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑‑AAAAAAACAGTATCGTGATATTTTACGCAAGGGTa                                                             >  1:1786684/1‑100 (MQ=255)
                                aaaTGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑‑AAAAAAACAGTATCGTGATATTTTACGCAAGGGTAg                                                            <  2:1364391/100‑1 (MQ=255)
                                   tgatgCTTCATACAGTCGTTATAAAAAGAATGAAACAGTATTTATTGTCTGTTCATCATTA‑‑AAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGccc                                                         <  1:1413269/100‑1 (MQ=255)
                                    gatgCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑‑AAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGcccc                                                        >  2:1455787/1‑100 (MQ=255)
                                    gatgCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑‑AAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGcccc                                                        >  2:694816/1‑100 (MQ=255)
                                     atgCTTCATATAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑‑AAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGccccc                                                       <  1:696002/100‑1 (MQ=255)
                                      tgCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑‑AAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCa                                                      <  2:1188133/100‑1 (MQ=255)
                                       gCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑‑AAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCAt                                                     >  1:1361431/1‑100 (MQ=255)
                                        cTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑‑AAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATg                                                    <  2:1008440/100‑1 (MQ=255)
                                        cTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑‑AAAAAAACAGTATCGAGATATTTTACGCAAGGGTAGCCCCCATg                                                    >  2:482961/1‑100 (MQ=255)
                                         ttCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑‑AAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGt                                                   >  1:1553374/1‑100 (MQ=255)
                                            aTACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑‑AAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCa                                                >  1:1070621/1‑100 (MQ=255)
                                             tACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑‑AAAAAAACAGTATTGTGATATTTTACGCAATTGTAGCCCCCATGTTCAt                                               >  1:1178225/1‑100 (MQ=255)
                                             tACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑‑AAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCAt                                               >  1:172461/1‑100 (MQ=255)
                                             tACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑‑AAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCAt                                               >  1:1233938/1‑100 (MQ=255)
                                                aGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAAAACAGTATCGTGAAACTTTACGCAAGGGAAACCCCCATGTTCAtt                                              >  2:1769321/1‑100 (MQ=255)
                                                   cTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑‑AAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCctc                                         <  2:1216861/100‑1 (MQ=255)
                                                   cTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑‑AAACAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCctc                                         >  1:1434929/1‑100 (MQ=255)
                                                    tttATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATAATTA‑‑AAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCACCATGTTCATTTCctct                                        >  1:777748/1‑100 (MQ=255)
                                                      tataAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑‑AAAAAAACAGTATCGTGATATTTTACGCATGGGTAGCCCCCATGTTCATTTCctct                                        >  1:399116/1‑98 (MQ=255)
                                                      tataAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑‑AAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCTCATGTTCATTTCCTCTgg                                      <  1:884745/100‑1 (MQ=255)
                                                      tataAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑‑AAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTgg                                      <  2:1442258/100‑1 (MQ=255)
                                                      tataAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑‑AAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTgg                                      <  2:629226/100‑1 (MQ=255)
                                                      tataAAAAGAATGAAACAGGATTTATTGTCTGTTCATAATTA‑‑AAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTgg                                      <  1:871046/100‑1 (MQ=255)
                                                       ataAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑‑AAAAAAACAGTATCGTGATATTTTACGCAAGGGTAACCCCCATAATCAttaacac                                         >  2:65910/1‑91 (MQ=255)
                                                        taAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑‑AAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGtt                                    <  1:613268/100‑1 (MQ=255)
                                                        taAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑‑AAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGtt                                    >  2:131385/1‑100 (MQ=255)
                                                          aaaaGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑‑AAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGtt                                    >  1:157708/1‑98 (MQ=255)
                                                          aaaaGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑‑AAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAg                                  <  2:1720815/100‑1 (MQ=255)
                                                          aaaaGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑‑AAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAg                                  <  2:75142/100‑1 (MQ=255)
                                                          aaaaGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑‑AAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCACCATGTTCATTTCCTCTGGTTAg                                  >  1:1372808/1‑100 (MQ=255)
                                                           aaaGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAAACAGTATCGTGATATTTTACGCAAGTGTAGCACCCATGCTCATGTCCTCTg                                       >  1:139677/1‑95 (MQ=255)
                                                           aaaGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑‑AAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTGTGGTTAGa                                 >  1:1281832/1‑100 (MQ=255)
                                                           aaaGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑‑AAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGa                                 >  2:1011680/1‑100 (MQ=255)
                                                           aaaGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑‑AAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGa                                 <  2:607048/100‑1 (MQ=255)
                                                           aaaGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑‑AAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGa                                 >  1:1351017/1‑100 (MQ=255)
                                                             aGAATGAAACAGGATTTTTTGTCTGTTCATCTTTA‑‑AAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATg                               <  2:55189/100‑1 (MQ=255)
                                                             aGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑‑AAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATg                               <  1:569171/100‑1 (MQ=255)
                                                              gAATGAAACAGGATTTATTGTCTGTTCATCATTA‑‑AAAAAAACAGTATCGGGATATTTTACGCAAGGGTAGCCCCCATTTTTATTTCCTCTGGTTAGATgg                              <  1:745601/100‑1 (MQ=255)
                                                               aaTGAAACAGGATTTATTGTCTGTTCATCATTA‑‑AAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATgg                              <  2:1548077/99‑1 (MQ=255)
                                                               aaTGAAACAGGATTTATTGTCTGTTCATCATTA‑‑AAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGc                             >  1:977142/1‑100 (MQ=255)
                                                               aaTGAAACAGGATTTATTGTCTGTTCATCATTA‑‑AAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGc                             <  2:212600/100‑1 (MQ=255)
                                                                cttccgatctaTTTATTGTCTGTTCATCATTA‑‑AAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCg                            <  2:1115192/90‑1 (MQ=255)
                                                                aTGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCAAGTTCATTTCCTCTGGTTGGATgg                              >  1:227611/1‑99 (MQ=255)
                                                                aTGAAACAGGATTTATTGTCTGTCCATCATTA‑‑AAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGACCCCATGTTCATTTCCTCTGGTTAGATGGCg                            >  1:902339/1‑100 (MQ=255)
                                                                  gAAACAGGATTTATTGTCTGTTCATCATTA‑‑AAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGAc                          <  1:262703/100‑1 (MQ=255)
                                                                  gAAACAGGATTTATTGTCTGTTCATCATTA‑‑AAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGAc                          >  2:1300141/1‑100 (MQ=255)
                                                                  gAAACAGGATTTATTATCTATTCATCATTA‑‑AAAAAAACAGTATCGTTATATTTTACACAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGa                           >  1:230043/1‑99 (MQ=255)
                                                                   aaaCAGGATTTATTGTCTGTTCATCATTA‑‑AAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACt                         <  1:131385/100‑1 (MQ=255)
                                                                     aCAGGATTTATTGTCTGTTCATCATTA‑‑AAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTAt                       <  2:1153208/100‑1 (MQ=255)
                                                                      cAGGATTTATTGTCTTTTCATTATTA‑‑AAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATg                      <  2:884175/100‑1 (MQ=255)
                                                                       aGGATTTATTGTCTGTTCATCATTA‑‑AAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGa                     >  2:1084574/1‑100 (MQ=255)
                                                                        ggATTTATTGTCTGTTCATCATTA‑‑AAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGaa                    <  2:1354907/100‑1 (MQ=255)
                                                                        ggATTTATTGTCTGTTCATCATTA‑‑AAAAAAACAGTATCGTGATATTTTACGCAAGGGTAACCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGaa                    >  1:1812590/1‑100 (MQ=255)
                                                                          aTTTATTGTCTGTTCATCATTA‑‑AAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGagatcggaag                  >  1:1115182/1‑91 (MQ=255)
                                                                          aTTTATTGTCTGTTCATCATTA‑‑AAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAAc                  <  1:1559965/100‑1 (MQ=255)
                                                                            ttATTGTCTGTTCATCATTA‑‑AAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCATCTGGTTAGATGGCGACTATGAAACTa                >  2:572714/1‑100 (MQ=255)
                                                                              aTTGTCTGTTCATCATTA‑‑AAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTAtt              >  1:1040163/1‑100 (MQ=255)
                                                                              aTTGTCTGTTCATCATTA‑‑AAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTAtt              >  1:1406803/1‑100 (MQ=255)
                                                                               ttGTCTGTTCATCATTA‑‑AAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATtc             <  1:1042560/100‑1 (MQ=255)
                                                                                tGTCTGTTCATCATTA‑‑AAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATtct            >  2:19848/1‑100 (MQ=255)
                                                                                  tCTGTTCATCATTA‑‑AAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATTCTcg          <  2:1618589/100‑1 (MQ=255)
                                                                                   cTGTTCATCATTA‑‑AAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATTCTcgc         <  1:189501/100‑1 (MQ=255)
                                                                                    tGTTCATCATTA‑‑AAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCGCCATGTTCATTTCCTCTGGTTAGATGGGGACTATGAAACTATTCTcgcg        >  1:908616/1‑100 (MQ=255)
                                                                                    tGTTCATCATTA‑‑AAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATTCTcgcg        <  1:17224/100‑1 (MQ=255)
                                                                                      ttCATCATTA‑‑AAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATTCTcgcgcg      <  2:621351/100‑1 (MQ=255)
                                                                                       tcatcaTTA‑‑AAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATTCTCGCGCGa     >  1:1671109/1‑100 (MQ=255)
                                                                                        catcatTA‑‑AAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCTTCTGGTTAGATGGCGACTATGAAACTATTCTCGCGCGaa    >  1:1230769/1‑100 (MQ=255)
                                                                                        catcatTA‑‑AAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATTCTCGCGCGaa    >  2:553369/1‑100 (MQ=255)
                                                                                          tcatTA‑‑AAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATTCTCGCGCGAATg  >  1:1809403/1‑100 (MQ=255)
                                                                                          tcatTA‑‑AAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATTCTCGCGCGAATg  <  2:1180814/100‑1 (MQ=255)
                                                                                          tcatTA‑‑AAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATTCTCGCGCGAATg  <  2:137317/100‑1 (MQ=255)
                                                                                          tcatTA‑‑AAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCACCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATTCTCGCGCGAATg  >  1:966160/1‑100 (MQ=255)
                                                                                               |                                                                                                 
GATGAAGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑‑AAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATTCTCGCGCGAATG  >  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: