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.A92.6% 400.4 / 22.5 121coding (254/564 nt)idnKD‑gluconate kinase, thermosensitive
Reads supporting (aligned to +/- strand):  ref base . (5/4);  new base A (57/55);  total (62/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.99e-01

Only 100 of 141 total aligned reads displayed.

TATCTGATGAAGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑‑AAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATTCTCGCGCGAATGCAGC  >  NC_000913/4494777‑4494974
                                                                                                     |                                                                                                       
tATCTGATGAAGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCAtt                                                                                                          <  2:2988741/100‑1 (MQ=255)
   cTGATGAAGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑‑aaa                                                                                                   >  2:1896966/1‑100 (MQ=255)
    tgatgaAGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑‑aaaa                                                                                                  <  1:355971/100‑1 (MQ=255)
    tgatgaAGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑‑aaaa                                                                                                  >  1:25189/1‑100 (MQ=255)
    tgatgaAGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑‑aaaa                                                                                                  >  2:2305449/1‑100 (MQ=255)
       tgaAGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑‑aaaaaaa                                                                                             >  1:897537/1‑100 (MQ=255)
       tgaAGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑‑aaaaaaa                                                                                             >  2:726532/1‑100 (MQ=255)
           gATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑‑aaaaaaaacag                                                                                         <  1:3006593/100‑5 (MQ=255)
           gATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑‑aaaaaaaacag                                                                                         <  2:2212580/100‑5 (MQ=255)
             tCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑‑AAAAAAACAGTa                                                                                        >  1:3315/1‑100 (MQ=255)
                aCTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑‑AAAAAAACAGTATCg                                                                                     >  1:2267761/1‑100 (MQ=255)
                aCTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑‑AAAAAAACAGTATCg                                                                                     >  1:1253649/1‑100 (MQ=255)
                aCTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑‑AAAAAAACAGTATCg                                                                                     <  1:2141442/100‑1 (MQ=255)
                aCTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑‑AAAAAAACAGTATCg                                                                                     <  2:1835921/100‑1 (MQ=255)
                 cTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑‑AAAAAAACAGTATCGt                                                                                    <  2:2479777/100‑1 (MQ=255)
                  ttCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑‑AAAAAAACAGTATCGTg                                                                                   <  2:496125/100‑1 (MQ=255)
                   tCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑‑AAAAAAACAGTATCGTGa                                                                                  <  1:2907946/100‑1 (MQ=255)
                    cccTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑‑AAAAAAACAGTATCGTGat                                                                                 >  2:1399016/1‑100 (MQ=255)
                       tggctggAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑‑AAAAAAACAGTATCGTGATAtt                                                                              <  2:1970328/100‑1 (MQ=255)
                        ggctggAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑‑AAAAAAACAGTATCGTGATAttt                                                                             <  2:1021441/100‑1 (MQ=255)
                        ggctggAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑‑AAAAAAACAGTATCGTGATAttt                                                                             <  2:1307594/100‑1 (MQ=255)
                        ggctggAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑‑AAAAAAACAGTATCGTGATAttt                                                                             >  2:1442297/1‑100 (MQ=255)
                           tggAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑‑AAAAAAACAGTATCGTGATATTTTa                                                                           <  2:1216673/99‑1 (MQ=255)
                           tggAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑‑AAAAAAACAGTATCGTGATATTTTAc                                                                          <  1:1379946/100‑1 (MQ=255)
                            ggAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑‑AAAAAAACAGTATCGTGATATTTTACg                                                                         <  1:2825945/100‑1 (MQ=255)
                                cGCTTAAATGATGGTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑‑AAAAAAACAGTATCGTGATATTTTACGCAAg                                                                     >  2:1877587/1‑100 (MQ=255)
                                cGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑‑AAAAAAACAGTATCGTGATATTTTACGCAAg                                                                     <  1:1352725/100‑1 (MQ=255)
                                 gCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑‑AAAAAAACAGTATCGTGATATTTTACGCAAgg                                                                    >  1:2900262/1‑100 (MQ=255)
                                  cTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑‑AAAAAAACAGTATCGTGATATTTTACGCAAggg                                                                   <  1:30410/100‑1 (MQ=255)
                                  cTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑‑AAAAAAACAGTATCGTGATATTTTACGCAAggg                                                                   >  2:2545434/1‑100 (MQ=255)
                                  cTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑‑AAAAAAACAGTATCGTGATATTTTACGCAAggg                                                                   >  2:1056297/1‑100 (MQ=255)
                                  cTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑‑AAAAAAACAGTATCGTGATATTTTACGCAAg                                                                     >  2:2842446/1‑98 (MQ=255)
                                  cTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑‑AAAAAAACAGTATCGTGATATTTTACGCAAGGGt                                                                  >  1:1539126/1‑100 (MQ=255)
                                   ttAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGTATTTATTGTCTGTTAATCATT‑‑‑AAAAAACACAGTATCCTGATAtttt                                                                            >  2:2733594/1‑90 (MQ=255)
                                    tAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATCA‑‑‑AAAAAAACAGTATCGTGATATTTTACGCAAGGGTa                                                                 >  1:97222/1‑100 (MQ=255)
                                    tAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑‑AAAAAAACAGTATCGTGATATTTTACGCAAg                                                                     >  2:1085217/1‑96 (MQ=255)
                                    tAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑‑AAAAAAACAGTATCGTGATATTTTACGCAAGGGTa                                                                 <  2:2221276/100‑1 (MQ=255)
                                    tAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑‑AAAAAAACAGTATCGTGATATTTTACGCAAGGGTAg                                                                >  2:111162/1‑100 (MQ=255)
                                    tAAATGATCGTTCAGACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑‑AAAAAAACAGTATCGTGATATTTTACGCAAGGGTa                                                                 <  2:1796250/100‑1 (MQ=255)
                                     aaaTGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑‑AAAAAAACAGTATCGTGATATTTTACGCAAGGGTAg                                                                <  1:2968494/100‑1 (MQ=255)
                                         gatgCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑‑AAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGcccc                                                            <  1:1000845/100‑1 (MQ=255)
                                         gatgCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑‑AAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGccc                                                             >  1:123137/1‑99 (MQ=255)
                                         gatgCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑‑AAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGc                                                               >  2:3002645/1‑97 (MQ=255)
                                           tgCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑‑AAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCa                                                          <  1:1442309/100‑1 (MQ=255)
                                           tgCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑‑AAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCa                                                          <  2:2663176/100‑1 (MQ=255)
                                           tgCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑‑AAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCAt                                                         <  2:1539113/100‑1 (MQ=255)
                                           tgCTTCATACAGTCCTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑‑AAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCa                                                          <  1:794125/100‑1 (MQ=255)
                                            gCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑‑AAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCAt                                                         >  2:869840/1‑100 (MQ=255)
                                             cTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑‑AAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATg                                                        >  1:294766/1‑100 (MQ=255)
                                               tCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑‑AAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGtt                                                      <  1:2301003/100‑1 (MQ=255)
                                                  tACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑‑AAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCAt                                                   >  1:2091566/1‑100 (MQ=255)
                                                   aCAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑‑AAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCAtt                                                  >  1:2544121/1‑100 (MQ=255)
                                                   aCAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑‑AAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCAtt                                                  >  2:2959092/1‑100 (MQ=255)
                                                   aCAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑‑AAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCAtt                                                  >  2:610413/1‑100 (MQ=255)
                                                      gTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑‑AAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTcc                                               <  1:2695113/100‑1 (MQ=255)
                                                      gTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑‑AAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTc                                                >  2:724145/1‑99 (MQ=255)
                                                         tttATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑‑AACAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTc                                                >  2:961163/1‑96 (MQ=255)
                                                          ttATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑‑AAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTg                                           <  2:1424832/100‑1 (MQ=255)
                                                           tataAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑‑AAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTgg                                          <  2:1949713/100‑1 (MQ=255)
                                                           tataAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑‑AAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTgg                                          >  2:196582/1‑100 (MQ=255)
                                                             taAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑‑AAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGtt                                        <  2:867458/100‑1 (MQ=255)
                                                             taAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑‑AAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGtt                                        <  1:1071444/100‑1 (MQ=255)
                                                             taAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑‑AAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGtt                                        <  2:944638/100‑1 (MQ=255)
                                                             taAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑‑AAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGtt                                        <  1:2115459/100‑1 (MQ=255)
                                                             taAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑‑AAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGtt                                        >  2:1071434/1‑100 (MQ=255)
                                                             taAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑‑AAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGtt                                        <  1:2017140/100‑1 (MQ=255)
                                                             taAAAAGAATGAAACAGGATATATTGTCTGTTCATCATTA‑‑AAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGtt                                        >  2:1858164/1‑100 (MQ=255)
                                                             taAAAAGAAGGAAACAGGATTTATTGTCTGTTCATAAATA‑‑AAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGtt                                        <  2:683790/100‑1 (MQ=255)
                                                              aaaaaGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑‑AAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGtt                                        >  1:868317/1‑99 (MQ=255)
                                                              aaaaaGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑‑AAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTa                                       >  2:2208261/1‑100 (MQ=255)
                                                              aaaaaGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑‑AAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTa                                       <  1:724152/100‑1 (MQ=255)
                                                              aaaaaGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑‑AAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTa                                       <  2:1385477/100‑1 (MQ=255)
                                                              aaaaaGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑‑AAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTa                                       >  2:511370/1‑100 (MQ=255)
                                                                  aGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑‑AAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATg                                   <  2:451245/100‑1 (MQ=255)
                                                                   gAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAAACAGTATCGTGATATTTTACGCAAGTGTAGCCCCCATGTTCATTTCCTCTGGTTAGATg                                   <  1:671806/100‑1 (MQ=255)
                                                                     aTGAAACAGGATTTATTGTCTGTTCATCATTA‑‑AAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCg                                <  1:1834512/100‑1 (MQ=255)
                                                                     aTGAAACAGGATTTATTGTCTGTTCATCATTA‑‑AAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCg                                <  1:151557/100‑1 (MQ=255)
                                                                      tGAAACAGGATTTATTGTCTGTTCATCATTA‑‑AAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGa                               <  1:883771/100‑1 (MQ=255)
                                                                        aaaCAGGATTTATTGTCTGTTCATCATTA‑‑AAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACt                             >  1:3012421/1‑100 (MQ=255)
                                                                        aaaCAGGATTTATTGTCTGTTCATCATTA‑‑AAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACt                             <  1:44033/100‑1 (MQ=255)
                                                                        aaaCAGGATTTATTGTCTGTTCATCATTA‑‑AAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCGGGTTAGATGGCGACt                             >  1:1039731/1‑100 (MQ=255)
                                                                          aCAGGATTTATTGTCTGTTCATCATTA‑‑AAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTAt                           <  2:2091548/100‑1 (MQ=255)
                                                                          aCAGGATTTATTGTCTGTTCATCATTA‑‑AAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTAt                           <  2:3315/100‑1 (MQ=255)
                                                                          aCAGGATTTATTGTCTGTTCATCATTA‑‑AAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTAt                           <  1:560768/100‑1 (MQ=255)
                                                                           cAGGATTTATTGTCTGTTCATCATTAAAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGa                               >  2:2188519/1‑97 (MQ=255)
                                                                           cAGGATTTATTGTCTGTTCATCATTA‑‑AAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATg                          <  1:1399028/100‑1 (MQ=255)
                                                                            aGGATTTATTGTCTGTTCATCATTA‑‑AAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGa                         <  1:1708259/100‑1 (MQ=255)
                                                                            aGGATTTATTGTCTGTTCATCATT‑‑‑AAAAACAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGa                         <  2:2450121/100‑1 (MQ=255)
                                                                             ggATTTATTGTCTGTTCATCATTA‑‑AAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGaa                        >  2:2461017/1‑100 (MQ=255)
                                                                               aattATTGTCTGTTCATCATTA‑‑AAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAAc                      >  1:2043910/3‑100 (MQ=255)
                                                                               aTTTATTGTCTGTTCATCATTA‑‑AAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAAc                      >  2:1813459/1‑100 (MQ=255)
                                                                               aTTTATTGTCTGTTCATCATTA‑‑AAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAAc                      <  2:2267741/100‑1 (MQ=255)
                                                                                  tATTGTCTGTTCATCATTA‑‑AAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTAt                   >  1:4613/1‑100 (MQ=255)
                                                                                      gTCTGTTCATCATTA‑‑AAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATtctc               >  2:1857245/1‑100 (MQ=255)
                                                                                       tCTGTTCATCATTA‑‑AAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATTCTcg              >  2:1126146/1‑100 (MQ=255)
                                                                                          gTTCATCATTA‑‑AAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATTCTcgcgc           >  2:418165/1‑100 (MQ=255)
                                                                                          gTTCATCATTA‑‑AAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATTCTcgcgc           >  1:332010/1‑100 (MQ=255)
                                                                                          gTTCATCATTA‑‑AAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATTCTcgcgc           >  2:2498979/1‑100 (MQ=255)
                                                                                           ttCATCATTA‑‑AAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATTCTcgcgcg          <  2:1335187/100‑1 (MQ=255)
                                                                                                  tt‑‑‑AAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATTCTCGCGCGAATGCAGc  >  1:2697838/1‑100 (MQ=255)
                                                                                                     |                                                                                                       
TATCTGATGAAGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑‑AAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATTCTCGCGCGAATGCAGC  >  NC_000913/4494777‑4494974

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