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)
                aCTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AAAAAAACAGTATCg                                                                               >  2:1237884/1‑100 (MQ=255)
                aCTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AAAAAAACAGTATCg                                                                               <  1:2526025/100‑1 (MQ=255)
                 cTTCCCTGGCTGGAACGCTTAAATGATGTTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AAAAAAACAGTATCGt                                                                              <  1:450752/100‑1 (MQ=255)
                  ttCCCTGGCTGGAACGCTTATATTATTCTTTATACAGTCTTTATAAAAATAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑AAAAAAACAATATCGTGa                                                                            <  2:3523966/100‑1 (MQ=255)
                  ttCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AAAAAAACAGTATCg                                                                               <  1:958602/98‑1 (MQ=255)
                    cccTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AAAAAAACAGTATCGTGat                                                                           >  2:1366346/1‑100 (MQ=255)
                      ctggctggAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AAAAAAACAGTATCGTGatat                                                                         >  1:10076/1‑100 (MQ=255)
                      ctggctggAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AAAAAAACAGTATCGTGatat                                                                         >  1:1234074/1‑100 (MQ=255)
                         gctggAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AAAAAAACAGTATCGTGATAtttt                                                                      <  1:618930/100‑1 (MQ=255)
                          ctggAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAAACAGTATCGTGATAtttt                                                                      >  2:2863006/1‑100 (MQ=255)
                           ttgAACGCTTAAATGATGTTTCATACAGTCTTTATAAAAAGAATGAAAAAGGATTTATTGTCTGTTCATCATTA‑AAAAAAACAGTATCGTGATATTTTa                                                                     <  1:1861589/97‑1 (MQ=255)
                           tggAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AAAAAAACAGTATCGTGATATTTTAc                                                                    <  1:1151343/100‑1 (MQ=255)
                           tggAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATAAAACAGGATTTATTGTCTGTTCATCATTA‑AAAAAAACAGTATCGTCATATTTTa                                                                     >  2:833855/1‑99 (MQ=255)
                            ggAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AAAAAAACAGTATCGTGATATTTTACg                                                                   >  1:1767180/1‑100 (MQ=255)
                                cGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTTTTCATCATTA‑AAAAAAACAGTATCGTGATATTTTACGc                                                                  >  2:3087555/1‑97 (MQ=255)
                                  cTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AAAAAAACAGTATCGTGATATTTTACGCAAggg                                                             >  1:2940147/1‑100 (MQ=255)
                                  cTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AAAAAAACAGTATCGTGATATTTTACGCAAggg                                                             <  1:1937752/100‑1 (MQ=255)
                                  cTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AAAAAAACAGTATCGTGATATTTTACGCAAggg                                                             >  2:768726/1‑100 (MQ=255)
                                  cTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AAAAAAACAGTATCGTGATATTTTACGCAAggg                                                             <  1:2452088/100‑1 (MQ=255)
                                  cTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AAAAAAACAGTATCGTGATATTTTACGCAAggg                                                             >  2:3184309/1‑100 (MQ=255)
                                  cTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑AAAAAAACAGTATCGTGATATTTTACGCAAGGGt                                                            <  1:3019158/100‑1 (MQ=255)
                                   ttAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AAAAAAACAGTATCGTGATATTTTACGCAAgg                                                              >  1:3574550/1‑98 (MQ=255)
                                    tAAATGCTGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AAAAAAACAGTATCGTGATATTTTACGCAAGGGTa                                                           >  2:794654/1‑100 (MQ=255)
                                    tAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATTATTA‑AAAAAAACAGTATCGTGATATTTTACGCAAGGGTa                                                           >  1:2167421/1‑100 (MQ=255)
                                    tAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AAAAAAACAGTATCGTGATATTTTACGCAAggga                                                            >  1:1408230/1‑98 (MQ=255)
                                    tAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AAAAAAACAGTATCGTGATATTTTACGCAAGGGTa                                                           <  1:607386/100‑1 (MQ=255)
                                    tAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AAAAAAACAGTATCGTGATATTTTACGCAAGGGTa                                                           >  1:3451852/1‑100 (MQ=255)
                                     aaaTGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AAAAAAACAGTATCGTGATATTTTACGCAAGGGTAg                                                          >  2:679584/1‑100 (MQ=255)
                                     aaaTGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AAAAAAACAGTATCGTGATATTTTACGCAAGGGTAg                                                          >  1:2796983/1‑100 (MQ=255)
                                        tgatgCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGccc                                                       >  1:1996415/1‑100 (MQ=255)
                                        tgatgCTTCAGACAGTTTTTATAAAAAGAATTAAACAGGATTTATTGTCTGTTTATCATTA‑AAAAAAACATTATCGTGATATTTTACGCAAGGGTAg                                                          <  2:1858338/97‑1 (MQ=255)
                                         gatgCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGcccc                                                      >  1:3714191/1‑100 (MQ=255)
                                           tgCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCa                                                    <  1:1702546/100‑1 (MQ=255)
                                            gCTTCATACAGTTTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCa                                                    <  1:2688691/99‑1 (MQ=255)
                                             cTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATg                                                  >  1:1001385/1‑100 (MQ=255)
                                                cATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTc                                               <  2:3659820/100‑1 (MQ=255)
                                                 aTACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCa                                              >  2:147010/1‑100 (MQ=255)
                                                  tACAGTCTTTATAAAAAGAATGAAACAAGATTTATTGTCTGTTCATCATTA‑AAAAAAACAGTATAGTGATATTTTACGCAAGAGTAACCCCAATGTTCa                                              >  2:824717/1‑99 (MQ=255)
                                                   aCAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCAtt                                            <  1:768735/100‑1 (MQ=255)
                                                     aGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTc                                          >  1:287938/1‑100 (MQ=255)
                                                     aGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTc                                          >  2:2764036/1‑100 (MQ=255)
                                                      gTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTcc                                         >  2:625515/1‑100 (MQ=255)
                                                      gTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTcc                                         >  2:2772033/1‑100 (MQ=255)
                                                       tCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑AAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCct                                        >  2:2553878/1‑99 (MQ=255)
                                                        cTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCctc                                       >  2:555640/1‑100 (MQ=255)
                                                        cTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AAAAAAACAGTATCGTGATATTTTACGCAAAGGTAGCCCCCATGTTCATTTCctc                                       >  1:1365315/1‑100 (MQ=255)
                                                         tttATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCctct                                      >  2:342850/1‑100 (MQ=255)
                                                          ttATAAAAAGAATTAAATAGGATTTATTGTCTGTTCATCATTA‑AAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTg                                     <  2:3650124/100‑1 (MQ=255)
                                                           tataAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTgg                                    <  2:398980/100‑1 (MQ=255)
                                                           tataAAAAGAATGAAACAGGATTTATTATCTGTTCATCATTA‑AAAAAAACAGTATCGTGCTATTTTACGTAAGGGGAGCCGCCATGTTCATTTCCTCTg                                     >  2:393371/1‑99 (MQ=255)
                                                            ataAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AAAAAAACAGTATCGTGATATTTTACGCAAGGTTAGCCCCCATGTTCATTTCCTCTGGt                                   <  2:1895970/100‑1 (MQ=255)
                                                            ataAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AAAAAAACAGTATCGTGATATTTTACGCAAAGGTAGCCCCCATGTTCATTTCCTCTg                                     >  1:497739/1‑98 (MQ=255)
                                                             taAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGtt                                  >  1:3392135/1‑100 (MQ=255)
                                                             taAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGtt                                  >  2:193141/1‑100 (MQ=255)
                                                              aaaaaGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGtt                                  <  2:158688/99‑1 (MQ=255)
                                                              aaaaaGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTa                                 <  1:2990775/100‑1 (MQ=255)
                                                              aaaaaGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCACCATGTTCATTTCCTCTGGTTa                                 >  2:1631621/1‑100 (MQ=255)
                                                               aaaaGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AAAAAAACATTATCATGATATTTTACACAAGAGTAGCCCCCATATTCATTTCCTCTGGTTAg                                >  1:194282/1‑100 (MQ=255)
                                                               aaaaGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAg                                >  2:151705/1‑100 (MQ=255)
                                                                       gAAACAGGATTTATTGTCTGTTCATCATTA‑AAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGAc                        >  1:2241397/1‑100 (MQ=255)
                                                                        aaaCAGGATTTATTGTCTGTTCATCATTA‑AAAAAAACAGTATCGTGCTATTTTACGCAAGGGTAACCCCAATGTTCATTTCCTCTGGTTAGATGGCGACt                       >  1:399380/1‑100 (MQ=255)
                                                                        aaaCAGGATTTATTGTCTGTTCATCATTA‑AAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACt                       >  2:554787/1‑100 (MQ=255)
                                                                            aGGATTTATTGTCTGTTCATCATTA‑AAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGa                   >  2:3130283/1‑100 (MQ=255)
                                                                               aTTTATTGTCTGTTCATCATTA‑AAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGa                   <  1:942088/97‑1 (MQ=255)
                                                                               aTTTATTGTCTGTTCATCATTA‑AAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAAc                >  2:1701470/1‑100 (MQ=255)
                                                                               aTTTATTGTCTGTTCATCATTA‑AAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAAc                <  2:1376543/100‑1 (MQ=255)
                                                                               aTTTATTGTCTGTTCATCATT‑‑AAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACt               >  1:2708971/1‑100 (MQ=255)
                                                                               aTTTATTGTCTGTTCATAATTA‑AAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAAc                <  2:285807/100‑1 (MQ=255)
                                                                                 ttATGGTCGGTTCATCATTA‑AAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACt               <  2:2808835/99‑1 (MQ=255)
                                                                                   aTTGTCTGTTCATCATTA‑AAAAAAACAGTATCGTGATATTTTACGCGAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTAtt            >  2:2172582/1‑100 (MQ=255)
                                                                                   aTTGTCTGTTCATCATTA‑AAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTAtt            >  2:3371833/1‑100 (MQ=255)
                                                                                    ttGTGTGTTCATCATTA‑AAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCTTCTGGTTAGATGGCGACTATGAAACTATtc           <  2:2109314/100‑1 (MQ=255)
                                                                                     tGTCTGTTCATCATTA‑AAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATtct          <  2:493633/100‑1 (MQ=255)
                                                                                      gTCTGTTCATCATTA‑AAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATtctc         >  2:1810812/1‑100 (MQ=255)
                                                                                      gTCTGTTCATCATTA‑AAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTATGATTAGATGGCGACTATGAAACTATtctc         >  2:2381121/1‑100 (MQ=255)
                                                                                        cTGTTCATCATTA‑AAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCACCCATGTTCATTTTCTCTGGTTAGATGGCGACTATGAAACTAtttttgc       >  2:3898806/1‑95 (MQ=255)
                                                                                         tGTTCATCATTA‑AAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATTCTcgcg      >  1:1078676/1‑100 (MQ=255)
                                                                                         tGTTCATCATTA‑AAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATTCTcgcg      <  2:3451813/100‑1 (MQ=255)
                                                                                          gTTCATCATTA‑AAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATTCTcgcgc     <  1:625522/100‑1 (MQ=255)
                                                                                            tcatcaTTA‑AAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATTCTCGCGCGa   >  2:89859/1‑100 (MQ=255)
                                                                                            tcatcaTTA‑AAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATTCTAGCGCGa   >  1:2326116/1‑100 (MQ=255)
                                                                                             catcatTA‑AAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATTCTCGCGCGaa  >  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: ‑