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.A85.0% 285.3 / 43.0 93coding (254/564 nt)idnKD‑gluconate kinase, thermosensitive
Reads supporting (aligned to +/- strand):  ref base . (8/6);  new base A (45/34);  total (53/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 108 total aligned reads displayed.

ATCTGATGAAGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATTCTCGCGCGAATG  >  NC_000913/4494778‑4494970
                                                                                                    |                                                                                              
aTCTGATGAAGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAAGAGGATTTATTGTCTGTTCATCATTa                                                                                               >  2:2801241/1‑100 (MQ=255)
aTCTGATGAAGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCAtt                                                                                                 >  2:2960309/1‑99 (MQ=255)
   tgatgaAGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTaaaa                                                                                            <  1:464916/100‑1 (MQ=255)
   tgatgaAGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTaaaa                                                                                            >  2:297480/1‑100 (MQ=255)
    gatgaAGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTaaaaa                                                                                           <  2:1488649/100‑1 (MQ=255)
     atgaAGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTaaaaaa                                                                                          >  2:3348749/1‑100 (MQ=255)
      tgaAGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTaaaaaaa                                                                                         <  1:1890219/100‑1 (MQ=255)
       gaAGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTaaaaaaaa                                                                                        >  2:2764158/1‑99 (MQ=255)
       gaAGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAc                                                                                        <  1:817073/100‑1 (MQ=255)
        aaGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTaaaaaaaac                                                                                       <  2:1284820/100‑3 (MQ=255)
        aaGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTaaaaaaaac                                                                                       <  1:2670248/100‑3 (MQ=255)
           aTCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTaaaaaaaacagt                                                                                    >  1:2891907/1‑95 (MQ=255)
           aTCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTaaaaaaaaacag                                                                                    <  1:2699070/100‑6 (MQ=255)
           aTCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTAATAaaaaaaaaaaaacag                                                                                    <  2:3118107/100‑6 (MQ=255)
            tCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTa                                                                                    >  2:1178848/1‑100 (MQ=255)
              gACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATc                                                                                  >  1:1101333/1‑100 (MQ=255)
               aCTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCg                                                                                 <  2:417864/100‑1 (MQ=255)
               aCTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCg                                                                                 <  1:44628/100‑1 (MQ=255)
                cTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGt                                                                                >  1:370156/1‑100 (MQ=255)
                   cccTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGat                                                                             <  2:3233535/100‑1 (MQ=255)
                    ccTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAACAGTATCGTGatat                                                                           <  1:1892653/100‑1 (MQ=255)
                        gctggAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATAtttt                                                                        <  1:1685924/100‑1 (MQ=255)
                         ctggAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAACAGTATCGTGATATTTTAc                                                                      >  2:3221332/1‑100 (MQ=255)
                           ggAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTa                                                                       >  2:3003600/1‑98 (MQ=255)
                            taaaatttAAATAATGCTTCATACATTCTTTAAAAAAAGAATGAAATAGTATTTATAGTCTTTTTATAATTAAAAAAAACAGTATCGTGATATTTTACg                                                                     <  2:751793/93‑1 (MQ=255)
                            gAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGc                                                                    >  1:530468/1‑100 (MQ=255)
                               cGCTTAACTGAGGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCGTTAAAAAAACAGTATCGTGATATTTTACGCAAgg                                                                >  2:1767378/1‑100 (MQ=255)
                               cGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAg                                                                 <  1:3268501/100‑1 (MQ=255)
                                gCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAg                                                                 >  2:2181214/1‑99 (MQ=255)
                                 cTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAggg                                                               >  1:2762394/1‑100 (MQ=255)
                                  ttAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGt                                                              <  2:3036053/100‑1 (MQ=255)
                                   tAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAg                                                            >  1:290950/1‑100 (MQ=255)
                                    aaaTGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTTAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAg                                                            <  1:2365047/100‑1 (MQ=255)
                                        gatgCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGcccc                                                        >  1:3257846/1‑100 (MQ=255)
                                         atgCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGccccc                                                       <  2:1399479/100‑1 (MQ=255)
                                          tgCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAACAGGAACGTGATTTTTTTCGCAAGGGTAGCCCCCAt                                                     >  1:2660744/1‑100 (MQ=255)
                                           gCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCAt                                                     >  2:1970869/1‑100 (MQ=255)
                                           gCTTCATACAGTCTTTATAAAAAGAACGAAACAGGATCTATTGTCTGTTCACCATTCAAAAAACAGTATCGTGATATTTTACGAAAGGGTACCCCCCATg                                                    >  2:1119648/1‑100 (MQ=255)
                                            cTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATg                                                    >  1:2615789/1‑100 (MQ=255)
                                            cTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGAGTAGCACCCATg                                                    >  2:2958263/1‑100 (MQ=255)
                                             ttCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATg                                                    <  2:1399094/99‑1 (MQ=255)
                                               cATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTc                                                 <  2:756869/100‑1 (MQ=255)
                                                aTACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCa                                                >  2:1177324/1‑100 (MQ=255)
                                                aTACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCa                                                >  1:2868651/1‑100 (MQ=255)
                                                  aCAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCAtt                                              >  2:18280/1‑100 (MQ=255)
                                                   cAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCAttt                                             <  1:1178852/100‑1 (MQ=255)
                                                    aGTCTTTATAAAAAGAATGAAACAGGCTTTATTGTCTGTTCATAATTAAAAAAAACAGTGTCGTGATATTTTACGGAAGGATAGCCCCCATGTTCAttt                                             >  2:407874/1‑99 (MQ=255)
                                                    aGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTTTTTATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCACCCATGTTCATTTc                                            >  2:1803356/1‑100 (MQ=255)
                                                     gTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTcc                                           <  1:1970877/100‑1 (MQ=255)
                                                     gTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCACCCATGTTCAttt                                             >  2:2815765/1‑98 (MQ=255)
                                                       cTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCctc                                         <  1:387666/100‑1 (MQ=255)
                                                       cTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCctc                                         <  1:1142221/100‑1 (MQ=255)
                                                       cTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCctc                                         >  2:1257339/1‑100 (MQ=255)
                                                       cTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCctct                                        >  1:151616/1‑100 (MQ=255)
                                                        tttATAAAAAGAATGAAACAGGATTTATTGTGTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCctc                                         <  2:2623874/99‑1 (MQ=255)
                                                         ttATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTg                                       <  1:2321295/100‑1 (MQ=255)
                                                          tataAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTgg                                      <  1:818324/100‑1 (MQ=255)
                                                           ataAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTgg                                      >  1:1075551/1‑99 (MQ=255)
                                                           ataAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGt                                     >  2:222474/1‑100 (MQ=255)
                                                            taAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGtt                                    <  1:1385047/100‑1 (MQ=255)
                                                            taAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGtt                                    >  1:1319874/1‑100 (MQ=255)
                                                             aaaaaGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCACCCATGTTCATTTCCTCTGGTTa                                   >  2:482395/1‑100 (MQ=255)
                                                             aaaaaGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAg                                  >  1:12367/1‑100 (MQ=255)
                                                              aaaaGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAg                                  >  2:1078140/1‑100 (MQ=255)
                                                              aaaaGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAg                                  >  2:3167743/1‑100 (MQ=255)
                                                              aaaaGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAg                                  <  1:2947837/100‑1 (MQ=255)
                                                              aaaaGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAg                                  >  1:2440442/1‑100 (MQ=255)
                                                              aaaaGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAg                                  >  2:3297545/1‑100 (MQ=255)
                                                              aaaaGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCACGTTCATTTCCTCTGGTTAg                                  >  1:985427/1‑100 (MQ=255)
                                                               aaaGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGa                                 <  1:286307/100‑1 (MQ=255)
                                                               aaaGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGAt                                <  2:151615/100‑1 (MQ=255)
                                                                 aGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATg                               <  2:2911827/100‑1 (MQ=255)
                                                                 aGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATg                               <  2:859463/100‑1 (MQ=255)
                                                                    aTGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAACAGTATCGTGATATTTTACGCAATGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGa                           >  2:672462/1‑100 (MQ=255)
                                                                     tGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGa                           <  1:2528664/100‑1 (MQ=255)
                                                                      gAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGAc                          <  2:1319870/100‑1 (MQ=255)
                                                                       aaaCAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACt                         >  1:3139275/1‑100 (MQ=255)
                                                                       aaaCAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACt                         >  1:343693/1‑100 (MQ=255)
                                                                        aaCAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTa                        <  1:1700664/100‑1 (MQ=255)
                                                                        aaCAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTa                        >  1:2810008/1‑100 (MQ=255)
                                                                        aaCAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTa                        <  1:2181225/100‑1 (MQ=255)
                                                                          cAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATg                      >  2:2574406/1‑100 (MQ=255)
                                                                           aGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGa                     >  1:3293002/1‑100 (MQ=255)
                                                                            ggATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGaa                    <  2:3257828/100‑1 (MQ=255)
                                                                            ggATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGa                     >  2:2610098/1‑99 (MQ=255)
                                                                             gATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGaaa                   >  2:2890082/1‑100 (MQ=255)
                                                                              aTTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAAc                  <  2:3298583/100‑1 (MQ=255)
                                                                              aTTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCACCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAAc                  >  2:2667319/1‑100 (MQ=255)
                                                                                 tATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTAt               >  2:2984988/1‑100 (MQ=255)
                                                                                  aTTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTAtt              >  2:567378/1‑100 (MQ=255)
                                                                                    tGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATtct            >  2:1782174/1‑100 (MQ=255)
                                                                                     gTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATtctc           >  1:1101884/1‑100 (MQ=255)
                                                                                      tCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATTCTcg          >  1:2718658/1‑100 (MQ=255)
                                                                                        tGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATTCTcgcg        >  1:470962/1‑100 (MQ=255)
                                                                                          ttCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATTCTcgcgcg      <  2:2615777/100‑1 (MQ=255)
                                                                                          ttCATCATTAAAAAAAACAGTATCCTGTTAATTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATTCTcgcg        <  2:746067/98‑1 (MQ=255)
                                                                                            catcatTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATTCTcgcgcg      >  2:3001195/1‑98 (MQ=255)
                                                                                            catcatTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATTCTCGCGCGaa    >  1:1219149/1‑100 (MQ=255)
                                                                                            catcatTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATTCTCGCGCGaa    >  1:2123145/1‑100 (MQ=255)
                                                                                              tcatTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATTCTCGCGCGAATg  <  1:1275935/100‑1 (MQ=255)
                                                                                                    |                                                                                              
ATCTGATGAAGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATTCTCGCGCGAATG  >  NC_000913/4494778‑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: