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.A97.5% 430.6 / 2.8 119coding (254/564 nt)idnKD‑gluconate kinase, thermosensitive
Reads supporting (aligned to +/- strand):  ref base . (3/0);  new base A (54/62);  total (57/62)
Fisher's exact test for biased strand distribution p-value = 1.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 137 total aligned reads displayed.

TATCTGATGAAGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑AA‑AAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATTCTCGCGCGAA  >  NC_000913/4494777‑4494968
                                                                                                    |                                                                                              
tATCTGATGAAGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCAtt                                                                                                 <  1:2592263/100‑1 (MQ=255)
 aTCTGATGAAGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑a                                                                                              <  1:2958436/100‑1 (MQ=255)
   cTGATGAAGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑aa‑a                                                                                           >  1:1714230/1‑100 (MQ=255)
   cTGATGAAGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑aa‑a                                                                                           >  2:2564518/1‑100 (MQ=255)
     gatgaAGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑aa‑aaa                                                                                         <  2:347298/100‑1 (MQ=255)
     gatgaAGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑aa‑aaa                                                                                         >  2:860593/1‑100 (MQ=255)
      atgaAGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑aa‑aaaa                                                                                        >  1:1714909/1‑100 (MQ=255)
       tgaAGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTGTTGTCTGTTTATCATT‑‑aa‑aaaga                                                                                       <  2:2075597/100‑3 (MQ=255)
       tgaAGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑aa‑aaaaa                                                                                       <  2:1414128/100‑1 (MQ=255)
         aaGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑aa‑aaaaaac                                                                                     <  2:465734/100‑3 (MQ=255)
          aGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGATTGAAACAGGATTTATTGTCTGTTCATCATT‑‑aa‑aaaaa                                                                                       >  1:748234/1‑97 (MQ=255)
          aGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑aa‑aaaaaaca                                                                                    <  2:1939567/100‑4 (MQ=255)
           gATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑aa‑aaaaaacag                                                                                   >  2:2707706/1‑96 (MQ=255)
             tCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AA‑AAAAACAGTa                                                                                  >  2:2416780/1‑100 (MQ=255)
                  ttCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AA‑AAAAACAGTATCGTg                                                                             <  1:2416769/100‑1 (MQ=255)
                   tCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AA‑AAAAACAGTATCGTGa                                                                            <  1:2378124/100‑1 (MQ=255)
                    cccTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AA‑AAAAACAGTATCGTGat                                                                           >  1:2448278/1‑100 (MQ=255)
                    cccTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AA‑AAAAACAGTATCGTGat                                                                           >  2:735868/1‑100 (MQ=255)
                     ccTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AA‑AAAAACAGTATCGTGata                                                                          >  1:1117141/1‑100 (MQ=255)
                     ccTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AA‑AAAAACAGTATCGTGata                                                                          >  1:231322/1‑100 (MQ=255)
                      ctggctggAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AA‑AAAAACAGTATCGTGatat                                                                         >  1:2341860/1‑100 (MQ=255)
                         gctggAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AA‑AAAAACAGTATCGTGATAtttt                                                                      <  2:712847/100‑1 (MQ=255)
                          ctggAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AA‑AAAAACAGTATCGTGATATTTTa                                                                     <  1:1996219/100‑1 (MQ=255)
                           tggAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AA‑AAAAACAGTATCGTGATATTTTAc                                                                    <  2:1724339/100‑1 (MQ=255)
                                cGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AA‑AAAAACAGTATCGTGATATTTTACGCAAg                                                               >  2:668185/1‑100 (MQ=255)
                                 gCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AA‑AAAAACAGTATCGTGATATTTTACGCAAgg                                                              <  1:2499298/100‑1 (MQ=255)
                                 gCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGACTTATAGTCTGTGCATCATTA‑AA‑AAAAACAGAATCGTGATCTTTTACGCTAg                                                               >  1:809463/1‑99 (MQ=255)
                                  cTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AA‑AAAAACAGTATCGTGATATTTTACGCAAggg                                                             <  1:375891/100‑1 (MQ=255)
                                  cTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AA‑AAAAACAGTATCGTGATATTTTACGCAAggg                                                             >  1:1013286/1‑100 (MQ=255)
                                    tAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AA‑AAAAACAGTATCGTGATATTTTACGCAAGGGTa                                                           >  1:218180/1‑100 (MQ=255)
                                    tAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AA‑AAAAACAGTATCGTGATATTTTACGCAAGGGTa                                                           >  1:374009/1‑100 (MQ=255)
                                     aaaTGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AA‑AAAAACAGTATCGTGATATTTTACGCAAGGGTAg                                                          >  1:2318562/1‑100 (MQ=255)
                                        tgatgCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AA‑AAAAACAGTATCGTGATATTTTACGCAAGGGTAGccc                                                       <  1:2794599/100‑1 (MQ=255)
                                        tgatgCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AA‑AAAAACAGTATCGTGATATTTTACGCAAGGGTAGccc                                                       >  1:1042231/1‑100 (MQ=255)
                                        tgatgCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AA‑AAAAACAGTATCGTGATATTTTACGCAAGGGTAGccc                                                       <  1:2344896/100‑1 (MQ=255)
                                         gatgCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AA‑AAAAACAGTATCGTGATATTTTACGCAAGGGTAGcccc                                                      >  1:3038586/1‑100 (MQ=255)
                                             cTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AA‑AAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATg                                                  <  1:2707690/100‑1 (MQ=255)
                                             cTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AA‑AAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATg                                                  >  1:917137/1‑100 (MQ=255)
                                                gtttcattCTTTTTATAAAGACTGAAACAGGATTTATTGTCTGTTCATCATTA‑AA‑AAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTc                                               >  2:1903840/8‑100 (MQ=255)
                                                cATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AA‑AAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTc                                               <  2:85262/100‑1 (MQ=255)
                                                 aTACAGTTTTTATAAAAAGAATGAAACAGGATTTATTGTCTGAAAAAAAAAA‑AA‑AAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGt                                                 <  1:2522122/97‑1 (MQ=255)
                                                 aTACAGTCTTTATAAAAAGAATGAAACTGAATTTATTGTCTGTTCATCATTA‑AA‑AAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTc                                               <  1:155986/99‑1 (MQ=255)
                                                      gTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AA‑AAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTcc                                         <  1:419871/100‑1 (MQ=255)
                                                       tCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AA‑AAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCct                                        <  2:1013290/100‑1 (MQ=255)
                                                        cTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AA‑AAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCctc                                       <  1:1588699/100‑1 (MQ=255)
                                                         tttATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AA‑AAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCctct                                      <  2:2613244/100‑1 (MQ=255)
                                                         tttATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATAAATA‑AA‑AAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCctct                                      <  1:2164977/100‑1 (MQ=255)
                                                          ttATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AA‑AAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTg                                     >  2:1681509/1‑100 (MQ=255)
                                                           tataAAAAGAATGAAACATGATTTATTGTCTGTTCATCATTA‑AA‑AAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTgg                                    <  2:1007827/100‑1 (MQ=255)
                                                           tataAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AA‑AAAAACAGTATCGTGATATTTTACGCATGGGTAGCCCCCATGTTCATTTCCTCTgg                                    >  1:1330939/1‑100 (MQ=255)
                                                           tataAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AA‑AAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTgg                                    <  1:2734277/100‑1 (MQ=255)
                                                           tataAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AA‑AAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTgg                                    <  2:1714916/100‑1 (MQ=255)
                                                           tataAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AA‑AAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTgg                                    <  2:596397/100‑1 (MQ=255)
                                                            ataAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AA‑AAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGt                                   >  2:2568160/1‑100 (MQ=255)
                                                             taAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AA‑AAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGtt                                  >  1:2240116/1‑100 (MQ=255)
                                                             taAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AA‑AAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCCCTGGtt                                  >  2:2437625/1‑100 (MQ=255)
                                                              aaaaaGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AA‑AAAAACAGTATCGTGATATTTTACGCAAGGTTAGCCCCCATGTTCATTTCCTATGGTTa                                 >  1:1104622/1‑100 (MQ=255)
                                                              aaaaaGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AA‑AAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTa                                 >  1:657288/1‑100 (MQ=255)
                                                              aaaaaGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AA‑AAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCACCATGTTCATTTCCTCTGGTTa                                 >  1:593541/1‑100 (MQ=255)
                                                               aaaaGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AA‑AAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTaa                                >  2:2917935/1‑99 (MQ=255)
                                                               aaaaGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AA‑AAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAg                                >  1:29038/1‑100 (MQ=255)
                                                                 aaGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AA‑AAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGAt                              <  2:898109/100‑1 (MQ=255)
                                                                 aaGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AA‑AAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGAt                              <  2:1424841/100‑1 (MQ=255)
                                                                    aaTGAAACAGGATTTATTGTCTGTTCATCATTAAAA‑AAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATgg                            <  2:2600739/100‑1 (MQ=255)
                                                                    aaTGAAACAGGATTTATTGTCTGTTCATCATTA‑AA‑AAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCTTCTGGTTAGATGGc                           >  1:1656395/1‑100 (MQ=255)
                                                                     aTGAAACAGGATTTATTGTCTGTTCATCATTA‑AA‑AAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCg                          >  1:1656897/1‑100 (MQ=255)
                                                                      tGAAACAGGATTTATTGTCTGTTCATCATTA‑AA‑AAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCg                          <  1:890713/99‑1 (MQ=255)
                                                                      tGAAACAGGATTTATTGTCTGTTCATCATTA‑AA‑AAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGa                         >  1:1688365/1‑100 (MQ=255)
                                                                       gAAACAGGATTTATTGTCTGTTCATCATTA‑AA‑AAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGAc                        <  2:1005326/100‑1 (MQ=255)
                                                                       gAAACAGGATTTATTGTCTGTTCATCATTA‑AA‑AAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGAc                        <  2:1754165/100‑1 (MQ=255)
                                                                       gAAACAGGATTTATTGTCTGTTCATCATTA‑AA‑AAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGAc                        <  1:165578/100‑1 (MQ=255)
                                                                        aaaCAGGATTTATTGTCTGTTCATCATTA‑AA‑AAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACt                       >  1:2314869/1‑100 (MQ=255)
                                                                         aaCAGGATTTATTGTCTGTTCATCATTA‑AA‑AAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTa                      <  2:1714237/100‑1 (MQ=255)
                                                                           cAGGATTTATTGTCTGTTCATCATTA‑AA‑AAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATg                    <  2:120938/100‑1 (MQ=255)
                                                                           cAGGATTTATTGTCTGTTCATCATTA‑AA‑AAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATg                    <  2:374010/100‑1 (MQ=255)
                                                                               aTTTATTGTCTGTTCATCATTA‑AA‑AAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAAc                <  2:792142/100‑1 (MQ=255)
                                                                               aTTTATTGTCTGTTCATCATTA‑AA‑AAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAAc                <  1:2564504/100‑1 (MQ=255)
                                                                               aTTTATTGTCTGTTCATCATTA‑AA‑AAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAAc                >  1:2905430/1‑100 (MQ=255)
                                                                                 ttATTGTCTGTTCATCATTA‑AA‑AAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTa              >  1:2894899/1‑100 (MQ=255)
                                                                                   aTTGTCTGTTCATCATTA‑AA‑AAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTAtt            >  2:2906654/1‑100 (MQ=255)
                                                                                    ttGTCTGTTCATCATTA‑AA‑AAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATtc           >  2:1243273/1‑100 (MQ=255)
                                                                                     tGTCTGTTCATCATTA‑AA‑AAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATtct          >  2:1429899/1‑100 (MQ=255)
                                                                                     tGTCTGTTCATCATTA‑AA‑AAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATtct          >  1:925308/1‑100 (MQ=255)
                                                                                     tGTCTGTTCATCATTA‑AA‑AAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATtct          <  1:3029498/100‑1 (MQ=255)
                                                                                      gTCTGTTCATCATTA‑AA‑AAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATtctc         >  2:2909217/1‑100 (MQ=255)
                                                                                      gTCTGTTCATCATTA‑AA‑AAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATtctc         >  1:1905822/1‑100 (MQ=255)
                                                                                       tCTGTTCATCATTA‑AA‑AAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATTCTcg        <  2:218180/100‑1 (MQ=255)
                                                                                       tCTGTTCATCATTA‑AA‑AAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATTCTcg        <  1:1124768/100‑1 (MQ=255)
                                                                                       tCTGTTCATCATT‑‑AACAAAAACAGTATCGTGATATTTTACGCAATGGAAGACCCCACGCTCATTTACTCTCGTTAGATCGCGCCCATGACACTactatc         >  1:2735053/1‑94 (MQ=255)
                                                                                         tGTTCATCATTA‑AA‑AAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGTGACTATGAAACTATTCTcgcg      >  1:1690349/1‑100 (MQ=255)
                                                                                         tGTTCATCATT‑‑AA‑AAAAACAGTATCGTGATATTTTACGTAATGGTAGCCCCCATGTTCATTTCTTTTGGTTAGATGGCGACTATGAAACTATtc           >  1:638099/1‑94 (MQ=255)
                                                                                         tGTTCATCATT‑‑AA‑AAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATTCTcgcgc     >  1:1220941/1‑100 (MQ=255)
                                                                                          gTTCATCATTA‑AA‑AAAAACAGTCTCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATTCTcgcgc     >  1:207503/1‑100 (MQ=255)
                                                                                          gTTCATCATTA‑AA‑AAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATTCTcgcgc     <  2:930504/100‑1 (MQ=255)
                                                                                           ttCATCATTA‑AA‑AAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATTCTcgcgcg    <  1:3041474/100‑1 (MQ=255)
                                                                                           ttCATCATTA‑AA‑AAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATTCTcgcgcg    <  1:1243269/100‑1 (MQ=255)
                                                                                           ttCATCATTA‑AA‑AAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATTCTcgcgcg    <  2:38804/100‑1 (MQ=255)
                                                                                            tcatcaTTA‑AA‑AAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATTCTCGCGCGa   >  1:1075206/1‑100 (MQ=255)
                                                                                             catcatTA‑AA‑AAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATTCTCGCGCGaa  >  1:2961372/1‑100 (MQ=255)
                                                                                             catcatTA‑AA‑AAAAACAGTATCGTGATATTTTACGCAAGGGCAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATTCTCGCGCGaa  >  1:1682276/1‑100 (MQ=255)
                                                                                                    |                                                                                              
TATCTGATGAAGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑AA‑AAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATTCTCGCGCGAA  >  NC_000913/4494777‑4494968

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