Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,494,883:1 (A)7→8 100% 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.0% 473.1 / 32.5 150coding (254/564 nt)idnKD‑gluconate kinase, thermosensitive
Reads supporting (aligned to +/- strand):  ref base . (3/9);  new base A (78/60);  total (81/69)
Fisher's exact test for biased strand distribution p-value = 6.62e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 174 total aligned reads displayed.

TATCTGATGAAGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑AAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATTCTCGCGCGAA  >  NC_000913/4494777‑4494968
                                                                                                    |                                                                                             
 aTCTGATGAAGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTTTAAAAAGAATGAAACAAGATATATTGTCTGATCACAATT‑‑a                                                                                             >  1:2595724/1‑100 (MQ=255)
   cTGATGAAGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑aa                                                                                            <  2:2100165/99‑1 (MQ=255)
    tgatgaAGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑aaaa                                                                                          <  1:3129866/100‑1 (MQ=255)
    tgatgaAGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑aaaa                                                                                          >  2:2972599/1‑100 (MQ=255)
    tgatgaAGATCGACTTCCCTGACAGGAACGCTTAAATGATGCTTCATACAGTCTCTCTAAAAAGAATGAAACAGGATTTATTGTCTGTTCCTCATT‑‑a                                                                                             >  2:1288065/1‑97 (MQ=255)
     gatgaAGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑aaaaa                                                                                         >  2:2481638/1‑100 (MQ=255)
       tgaAGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑aaaaaaa                                                                                       >  1:1581420/1‑100 (MQ=255)
         aaGATCGTTTTTCCTGGCTGGAACGCTTAAATGATGCTTCATAAAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑aaaaaaaac                                                                                     <  1:206573/100‑3 (MQ=255)
         aaGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑aaaaaaaac                                                                                     <  1:1793318/100‑3 (MQ=255)
         aaGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑aaaaaaaac                                                                                     >  2:2534036/1‑98 (MQ=255)
          aGATCGACTTCCCTGGCTGGAATGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑aaaaaaaaca                                                                                    <  1:3811230/100‑4 (MQ=255)
          aGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑aaaaaaaaca                                                                                    <  1:1112061/100‑4 (MQ=255)
          aGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTATCTATTCATCATT‑‑aaaaaaaaca                                                                                    >  1:1554151/1‑97 (MQ=255)
           gATCGACTTCCCTGGCTGGAACGTTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑aaaaagaaca                                                                                    <  2:3518117/99‑6 (MQ=255)
           gATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑aaaaaaaacag                                                                                   <  2:2821785/100‑5 (MQ=255)
            aTCGACTTCCCTGGCTGGAACGCTTAGATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑aaaaaaaacag                                                                                   <  2:2658717/99‑5 (MQ=255)
            aTCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑aaaaaaaacagt                                                                                  >  2:3822717/1‑95 (MQ=255)
               gACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑aaaaaaaacagaat                                                                                >  1:904684/1‑92 (MQ=255)
                aCTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCg                                                                               >  2:1237884/1‑100 (MQ=255)
                aCTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCg                                                                               <  1:2526025/100‑1 (MQ=255)
                 cTTCCCTGGCTGGAACGCTTAAATGATGTTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGt                                                                              <  1:450752/100‑1 (MQ=255)
                  ttCCCTGGCTGGAACGCTTATATTATTCTTTATACAGTCTTTATAAAAATAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑AAAAAAACAATATCGTGa                                                                            <  2:3523966/100‑1 (MQ=255)
                  ttCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCg                                                                               <  1:958602/98‑1 (MQ=255)
                    cccTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGat                                                                           >  2:1366346/1‑100 (MQ=255)
                      ctggctggAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGatat                                                                         >  1:10076/1‑100 (MQ=255)
                      ctggctggAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGatat                                                                         >  1:1234074/1‑100 (MQ=255)
                         gctggAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATAtttt                                                                      <  1:618930/100‑1 (MQ=255)
                          ctggAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAAACAGTATCGTGATAtttt                                                                      >  2:2863006/1‑100 (MQ=255)
                           ttgAACGCTTAAATGATGTTTCATACAGTCTTTATAAAAAGAATGAAAAAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTa                                                                     <  1:1861589/97‑1 (MQ=255)
                           tggAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTAc                                                                    <  1:1151343/100‑1 (MQ=255)
                           tggAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATAAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTCATATTTTa                                                                     >  2:833855/1‑99 (MQ=255)
                            ggAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACg                                                                   >  1:1767180/1‑100 (MQ=255)
                                cGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTTTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGc                                                                  >  2:3087555/1‑97 (MQ=255)
                                  cTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAggg                                                             >  1:2940147/1‑100 (MQ=255)
                                  cTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAggg                                                             <  1:1937752/100‑1 (MQ=255)
                                  cTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAggg                                                             >  2:768726/1‑100 (MQ=255)
                                  cTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAggg                                                             <  1:2452088/100‑1 (MQ=255)
                                  cTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAggg                                                             >  2:3184309/1‑100 (MQ=255)
                                  cTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑AAAAAAACAGTATCGTGATATTTTACGCAAGGGt                                                            <  1:3019158/100‑1 (MQ=255)
                                   ttAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAgg                                                              >  1:3574550/1‑98 (MQ=255)
                                    tAAATGCTGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTa                                                           >  2:794654/1‑100 (MQ=255)
                                    tAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATTATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTa                                                           >  1:2167421/1‑100 (MQ=255)
                                    tAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAggga                                                            >  1:1408230/1‑98 (MQ=255)
                                    tAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTa                                                           <  1:607386/100‑1 (MQ=255)
                                    tAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTa                                                           >  1:3451852/1‑100 (MQ=255)
                                     aaaTGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAg                                                          >  2:679584/1‑100 (MQ=255)
                                     aaaTGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAg                                                          >  1:2796983/1‑100 (MQ=255)
                                        tgatgCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGccc                                                       >  1:1996415/1‑100 (MQ=255)
                                        tgatgCTTCAGACAGTTTTTATAAAAAGAATTAAACAGGATTTATTGTCTGTTTATCATTAAAAAAAACATTATCGTGATATTTTACGCAAGGGTAg                                                          <  2:1858338/97‑1 (MQ=255)
                                         gatgCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGcccc                                                      >  1:3714191/1‑100 (MQ=255)
                                           tgCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCa                                                    <  1:1702546/100‑1 (MQ=255)
                                            gCTTCATACAGTTTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCa                                                    <  1:2688691/99‑1 (MQ=255)
                                             cTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATg                                                  >  1:1001385/1‑100 (MQ=255)
                                                cATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTc                                               <  2:3659820/100‑1 (MQ=255)
                                                 aTACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCa                                              >  2:147010/1‑100 (MQ=255)
                                                  tACAGTCTTTATAAAAAGAATGAAACAAGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATAGTGATATTTTACGCAAGAGTAACCCCAATGTTCa                                              >  2:824717/1‑99 (MQ=255)
                                                   aCAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCAtt                                            <  1:768735/100‑1 (MQ=255)
                                                     aGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTc                                          >  1:287938/1‑100 (MQ=255)
                                                     aGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTc                                          >  2:2764036/1‑100 (MQ=255)
                                                      gTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTcc                                         >  2:625515/1‑100 (MQ=255)
                                                      gTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTcc                                         >  2:2772033/1‑100 (MQ=255)
                                                       tCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑AAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCct                                        >  2:2553878/1‑99 (MQ=255)
                                                        cTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCctc                                       >  2:555640/1‑100 (MQ=255)
                                                        cTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAAGGTAGCCCCCATGTTCATTTCctc                                       >  1:1365315/1‑100 (MQ=255)
                                                         tttATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCctct                                      >  2:342850/1‑100 (MQ=255)
                                                          ttATAAAAAGAATTAAATAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTg                                     <  2:3650124/100‑1 (MQ=255)
                                                           tataAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTgg                                    <  2:398980/100‑1 (MQ=255)
                                                           tataAAAAGAATGAAACAGGATTTATTATCTGTTCATCATTAAAAAAAACAGTATCGTGCTATTTTACGTAAGGGGAGCCGCCATGTTCATTTCCTCTg                                     >  2:393371/1‑99 (MQ=255)
                                                            ataAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGTTAGCCCCCATGTTCATTTCCTCTGGt                                   <  2:1895970/100‑1 (MQ=255)
                                                            ataAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAAGGTAGCCCCCATGTTCATTTCCTCTg                                     >  1:497739/1‑98 (MQ=255)
                                                             taAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGtt                                  >  1:3392135/1‑100 (MQ=255)
                                                             taAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGtt                                  >  2:193141/1‑100 (MQ=255)
                                                              aaaaaGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGtt                                  <  2:158688/99‑1 (MQ=255)
                                                              aaaaaGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTa                                 <  1:2990775/100‑1 (MQ=255)
                                                              aaaaaGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCACCATGTTCATTTCCTCTGGTTa                                 >  2:1631621/1‑100 (MQ=255)
                                                               aaaaGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACATTATCATGATATTTTACACAAGAGTAGCCCCCATATTCATTTCCTCTGGTTAg                                >  1:194282/1‑100 (MQ=255)
                                                               aaaaGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAg                                >  2:151705/1‑100 (MQ=255)
                                                                       gAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGAc                        >  1:2241397/1‑100 (MQ=255)
                                                                        aaaCAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGCTATTTTACGCAAGGGTAACCCCAATGTTCATTTCCTCTGGTTAGATGGCGACt                       >  1:399380/1‑100 (MQ=255)
                                                                        aaaCAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACt                       >  2:554787/1‑100 (MQ=255)
                                                                            aGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGa                   >  2:3130283/1‑100 (MQ=255)
                                                                               aTTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGa                   <  1:942088/97‑1 (MQ=255)
                                                                               aTTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAAc                >  2:1701470/1‑100 (MQ=255)
                                                                               aTTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAAc                <  2:1376543/100‑1 (MQ=255)
                                                                               aTTTATTGTCTGTTCATCATT‑‑AAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACt               >  1:2708971/1‑100 (MQ=255)
                                                                               aTTTATTGTCTGTTCATAATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAAc                <  2:285807/100‑1 (MQ=255)
                                                                                 ttATGGTCGGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACt               <  2:2808835/99‑1 (MQ=255)
                                                                                   aTTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCGAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTAtt            >  2:2172582/1‑100 (MQ=255)
                                                                                   aTTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTAtt            >  2:3371833/1‑100 (MQ=255)
                                                                                    ttGTGTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCTTCTGGTTAGATGGCGACTATGAAACTATtc           <  2:2109314/100‑1 (MQ=255)
                                                                                     tGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATtct          <  2:493633/100‑1 (MQ=255)
                                                                                      gTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATtctc         >  2:1810812/1‑100 (MQ=255)
                                                                                      gTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTATGATTAGATGGCGACTATGAAACTATtctc         >  2:2381121/1‑100 (MQ=255)
                                                                                        cTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCACCCATGTTCATTTTCTCTGGTTAGATGGCGACTATGAAACTAtttttgc       >  2:3898806/1‑95 (MQ=255)
                                                                                         tGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATTCTcgcg      >  1:1078676/1‑100 (MQ=255)
                                                                                         tGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATTCTcgcg      <  2:3451813/100‑1 (MQ=255)
                                                                                          gTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATTCTcgcgc     <  1:625522/100‑1 (MQ=255)
                                                                                            tcatcaTTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATTCTCGCGCGa   >  2:89859/1‑100 (MQ=255)
                                                                                            tcatcaTTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATTCTAGCGCGa   >  1:2326116/1‑100 (MQ=255)
                                                                                             catcatTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATTCTCGCGCGaa  >  1:3640520/1‑100 (MQ=255)
                                                                                                    |                                                                                             
TATCTGATGAAGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑AAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATTCTCGCGCGAA  >  NC_000913/4494777‑4494968

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