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.8% 422.5 / 5.4 124coding (254/564 nt)idnKD‑gluconate kinase, thermosensitive
Reads supporting (aligned to +/- strand):  ref base . (2/2);  new base A (63/57);  total (65/59)
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.26e-01

Only 100 of 145 total aligned reads displayed.

TATCTGATGAAGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑AAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATTCTCGCGCGAATG  >  NC_000913/4494777‑4494970
                                                                                                    |                                                                                                
tATCTGATGAAGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCAtt                                                                                                   <  2:1144915/100‑1 (MQ=255)
   cTGATGAAGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑aaa                                                                                             <  1:1456205/100‑1 (MQ=255)
   cTGATGAAGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCGGTTCATCATT‑‑aaa                                                                                             >  2:26073/1‑100 (MQ=255)
    tgatgaAGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑aaaa                                                                                            >  2:2150022/1‑100 (MQ=255)
    tgatgaAGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑aaaa                                                                                            >  2:2558144/1‑100 (MQ=255)
     gatgaAGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑aaaaa                                                                                           >  2:834602/1‑100 (MQ=255)
      atgaAGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑aaaaaa                                                                                          <  2:594072/100‑1 (MQ=255)
       tgaAGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑aaaaaaa                                                                                         >  2:1648955/1‑100 (MQ=255)
       tgaAGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑aaaaaaa                                                                                         >  1:602681/1‑100 (MQ=255)
         aaGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑aaaaaaaac                                                                                       <  1:1659013/100‑3 (MQ=255)
         aaGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑aaaaaaaac                                                                                       <  2:1427626/100‑3 (MQ=255)
         aaGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑aaaaaaaac                                                                                       <  2:624959/100‑3 (MQ=255)
         aaGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑aaaaaaaac                                                                                       <  2:2472463/100‑3 (MQ=255)
         aaGATCGACTTCCCTAGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCAGCATT‑‑aaaaaaaac                                                                                       <  1:1797790/100‑3 (MQ=255)
          aGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑aaaaaaaaac                                                                                      <  2:1273932/100‑4 (MQ=255)
           gATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑aaaaaaaacag                                                                                     <  1:1587379/100‑5 (MQ=255)
              cGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTAt                                                                                   <  1:1623789/100‑1 (MQ=255)
                aCTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCg                                                                                 <  1:2554698/100‑1 (MQ=255)
                aCTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCg                                                                                 <  2:2056774/100‑1 (MQ=255)
                 cTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGCAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCg                                                                                 >  2:932942/1‑99 (MQ=255)
                 cTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGt                                                                                >  1:2251337/1‑100 (MQ=255)
                     ccTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGata                                                                            >  2:402388/1‑100 (MQ=255)
                     ccTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGata                                                                            >  1:2102461/1‑100 (MQ=255)
                      ctggctggAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGatat                                                                           >  1:2789161/1‑100 (MQ=255)
                      ctggctggAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGatat                                                                           >  2:1803281/1‑100 (MQ=255)
                       tggctggAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATAtt                                                                          <  1:2859878/100‑1 (MQ=255)
                       tggctggAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATAtt                                                                          <  2:259596/100‑1 (MQ=255)
                       tggctggAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATAtt                                                                          >  1:1401093/1‑100 (MQ=255)
                          ctggAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTa                                                                       >  2:2796197/1‑100 (MQ=255)
                           tggAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTAc                                                                      <  2:1229427/100‑1 (MQ=255)
                            ggAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACg                                                                     <  2:2106888/100‑1 (MQ=255)
                                cGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAg                                                                 >  2:906728/1‑100 (MQ=255)
                                 gCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAgg                                                                >  2:4381/1‑100 (MQ=255)
                                  cTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAggg                                                               <  1:2557053/100‑1 (MQ=255)
                                  cTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAggg                                                               >  1:315984/1‑100 (MQ=255)
                                   ttAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGt                                                              <  2:602687/100‑1 (MQ=255)
                                    tAAATGATGCTTCATACAGTCTTTCTAAAAAGAATGAAACAGGTTTTATTGTCTGTTCATCATTAAAAAAAACAGTCTTGTGATATTTTACGCCAAGGt                                                              >  1:924489/1‑99 (MQ=255)
                                    tAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTa                                                             <  2:2637595/100‑1 (MQ=255)
                                    tAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTa                                                             >  2:2151614/1‑100 (MQ=255)
                                     aaaTGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAg                                                            >  1:656789/1‑100 (MQ=255)
                                     aaaTGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAg                                                            >  1:2237532/1‑100 (MQ=255)
                                     aaaTGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAg                                                            >  2:806866/1‑100 (MQ=255)
                                       atgatgCTTCATACAGTCTTTATAAAAAGAATGACACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGAAAGGGTAg                                                            >  2:1045967/1‑98 (MQ=255)
                                        tgatgCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGccc                                                         <  1:858828/100‑1 (MQ=255)
                                            gCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCa                                                      <  2:636917/100‑1 (MQ=255)
                                            gCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCAt                                                     >  1:1599598/1‑100 (MQ=255)
                                             cTTCATACAGTCTTTATAAAAAGAATGAAATAGGATTTATTGTCTGTTAATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATg                                                    <  1:489523/100‑1 (MQ=255)
                                               tCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGtt                                                  <  2:1074673/100‑1 (MQ=255)
                                                  tACAGTCTTTGTAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACTGTATCGTGATATTTTACGCAAGGGAAGCACCCAGGTTCAt                                               >  1:820235/1‑100 (MQ=255)
                                                  tACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCa                                                >  1:49769/1‑99 (MQ=255)
                                                  tACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCAt                                               >  2:2214883/1‑100 (MQ=255)
                                                  tACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCAt                                               >  1:1498835/1‑100 (MQ=255)
                                                  tACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCAt                                               >  1:981688/1‑100 (MQ=255)
                                                  tACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCAt                                               >  1:2618283/1‑100 (MQ=255)
                                                    cAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCAGCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCAttt                                             <  1:2771420/100‑1 (MQ=255)
                                                     aGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTAAttt                                             >  1:1263108/1‑99 (MQ=255)
                                                      gTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTcc                                           >  2:695260/1‑100 (MQ=255)
                                                      gTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTcc                                           >  1:877007/1‑100 (MQ=255)
                                                      gTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTc                                            >  1:877313/1‑99 (MQ=255)
                                                       tCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTcc                                           <  1:1467398/100‑1 (MQ=255)
                                                       tCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCct                                          <  1:2537339/100‑1 (MQ=255)
                                                         tttATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCctct                                        >  1:2087464/1‑100 (MQ=255)
                                                          ttATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCATGTCTACACCCCATGTTCATTTCCTCTg                                       >  1:1379661/1‑100 (MQ=255)
                                                           tataAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTgg                                      >  1:949079/1‑100 (MQ=255)
                                                           tataAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTgg                                      <  2:688076/100‑1 (MQ=255)
                                                           tataAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑AAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGt                                     <  2:2571706/100‑1 (MQ=255)
                                                            ataAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGt                                     <  2:2545544/100‑1 (MQ=255)
                                                            ataAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGt                                     >  1:2585939/1‑100 (MQ=255)
                                                             taAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGtt                                    >  1:69409/1‑100 (MQ=255)
                                                             taAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGtt                                    <  1:1194706/100‑1 (MQ=255)
                                                             taAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGtt                                    >  1:1050051/1‑100 (MQ=255)
                                                              aaaaaGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTa                                   >  1:1406782/1‑100 (MQ=255)
                                                              aaaaaGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTa                                   >  2:1267855/1‑100 (MQ=255)
                                                               aaaaGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCACCCATGTTCATTTCCTCTGGTTAg                                  >  1:433565/1‑100 (MQ=255)
                                                                aaaGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGa                                 >  2:439117/1‑100 (MQ=255)
                                                                aaaGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGa                                 <  2:1037511/100‑1 (MQ=255)
                                                                aaaGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑ACAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGa                                 >  1:1556342/1‑100 (MQ=255)
                                                                   gAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATgg                              <  1:1689546/100‑1 (MQ=255)
                                                                    aaTGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGc                             <  1:142514/100‑1 (MQ=255)
                                                                    aaTGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGc                             <  1:898460/100‑1 (MQ=255)
                                                                    aaTGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCACCCATGTTCATTTCCTCTGGTTAGATGGc                             >  1:1837519/1‑100 (MQ=255)
                                                                     aTGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCg                            >  1:325010/1‑100 (MQ=255)
                                                                     aTGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGCAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCg                            <  1:1046472/100‑1 (MQ=255)
                                                                       gAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGAc                          <  2:2605808/100‑1 (MQ=255)
                                                                       gAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGAc                          <  2:1950602/100‑1 (MQ=255)
                                                                        aaaCAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACt                         >  2:1092225/1‑100 (MQ=255)
                                                                         aaCAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGAcca                        >  1:1760518/1‑98 (MQ=255)
                                                                             ggATTTATTGTCTGTTCATCATTAAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGa                     <  2:335894/100‑1 (MQ=255)
                                                                             ggATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGaa                    <  2:656795/100‑1 (MQ=255)
                                                                             ggATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGaa                    <  1:740256/100‑1 (MQ=255)
                                                                                tttATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACt                 <  1:2851611/100‑1 (MQ=255)
                                                                                      gTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATtctc           <  1:662570/100‑1 (MQ=255)
                                                                                      gTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATtctc           >  1:1268904/1‑100 (MQ=255)
                                                                                       tCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATTCTcg          <  1:106511/100‑1 (MQ=255)
                                                                                        cTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATTCTcgc         >  1:1145810/1‑100 (MQ=255)
                                                                                          gTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGAAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATTCTcgcgc       >  1:1795353/1‑100 (MQ=255)
                                                                                            tcatcaTTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATTCTCGCGCGa     >  1:1895974/1‑100 (MQ=255)
                                                                                            tcatcaTTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATTCTCGCGCGa     >  2:1783798/1‑100 (MQ=255)
                                                                                             catcatTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATTCTCGCGCGaa    <  2:152858/100‑1 (MQ=255)
                                                                                               tcatTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATTCTCGCGCGAATg  <  2:383815/100‑1 (MQ=255)
                                                                                                    |                                                                                                
TATCTGATGAAGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑AAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATTCTCGCGCGAATG  >  NC_000913/4494777‑4494970

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