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.A94.2% 460.2 / 19.2 138coding (254/564 nt)idnKD‑gluconate kinase, thermosensitive
Reads supporting (aligned to +/- strand):  ref base . (3/5);  new base A (68/62);  total (71/67)
Fisher's exact test for biased strand distribution p-value = 4.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 165 total aligned reads displayed.

TCTGATGAAGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑AAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATTCTCGCGCGAATG  >  NC_000913/4494779‑4494970
                                                                                                  |                                                                                                
 cTGATGAAGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑aaa                                                                                              >  2:2724620/1‑100 (MQ=255)
 cTGATGAAGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑aaa                                                                                              >  2:247533/1‑100 (MQ=255)
 cTGATGAAGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑aaa                                                                                              <  1:2354640/100‑1 (MQ=255)
  tgatgaAGATCGACTTCCGTGGCTGTAACGCTTAAATGATGCTTCATAAAAAAATAATAAAAAGAATGAAACAGGATTTCTTGTCTGTTCATCATT‑‑aaaa                                                                                             <  2:2340879/100‑1 (MQ=255)
  tgatgaAGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAGAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑aaaa                                                                                             >  2:3421000/1‑100 (MQ=255)
   gatgaAGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑aaaaa                                                                                           >  2:3562678/1‑100 (MQ=255)
     tgaAGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑aaaaaaa                                                                                         <  2:3116205/100‑1 (MQ=255)
     tgaAGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGATTTATTGTCTGTTCAACATT‑‑aaaaaaaa                                                                                        >  2:3631279/1‑99 (MQ=255)
      gaAGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑aaaaaaaa                                                                                        >  1:3411702/1‑99 (MQ=255)
       cgatctgACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑aaaaaaaac                                                                                       <  1:1246961/94‑3 (MQ=255)
       aaGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑aaaaaaaac                                                                                       <  2:2554045/100‑3 (MQ=255)
       aaGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑aaaaaaaac                                                                                       >  2:3822630/1‑98 (MQ=255)
       aaGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑aaaaaaaac                                                                                       <  1:3051988/100‑3 (MQ=255)
        aGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑aaaaaaaaca                                                                                      <  2:3059093/100‑4 (MQ=255)
        aGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑aaaaaaaac                                                                                       <  2:2165027/99‑3 (MQ=255)
         gATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑aaaaaaagcag                                                                                     >  2:286900/1‑96 (MQ=255)
         gATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑aaaaaaaacag                                                                                     >  1:3154165/1‑96 (MQ=255)
            cGACTTCACTGGCTGGAAAGCTTAAATGATACTTCATACAGTCTTTATAAAAAGAATGAAACAGTATTTATTGTATGTTAATCATT‑‑aaaaaaaacagt                                                                                    <  2:3441692/98‑6 (MQ=255)
             gACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATc                                                                                  >  2:1266174/1‑100 (MQ=255)
             gACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTAGTGTCTGTTCATCATTAAAAAAAACAGTaac                                                                                  >  1:1632776/1‑98 (MQ=255)
              aCTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCg                                                                                 <  1:923703/100‑1 (MQ=255)
              aCTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCg                                                                                 <  1:804308/100‑1 (MQ=255)
               cTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGt                                                                                <  1:292306/100‑1 (MQ=255)
                 tCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGa                                                                              <  1:395620/100‑1 (MQ=255)
                   ccTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGata                                                                            >  1:3454565/1‑100 (MQ=255)
                    ctggctggAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGatat                                                                           >  2:2906436/1‑100 (MQ=255)
                    ctggctggAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGata                                                                            >  2:2281814/1‑99 (MQ=255)
                      ggctggAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATAttt                                                                         >  1:1759346/1‑100 (MQ=255)
                      ggctggAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATAttt                                                                         >  1:3413739/1‑100 (MQ=255)
                      ggctggAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATAttt                                                                         <  1:1282773/100‑1 (MQ=255)
                        ctggAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTa                                                                       >  2:3757891/1‑100 (MQ=255)
                         tggAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTAc                                                                      <  1:1759136/100‑1 (MQ=255)
                         tggAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTAAAAAAAAAAAAAAAACAGTATCGTGATATTTTAc                                                                      <  2:2109591/100‑1 (MQ=255)
                          ggAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACg                                                                     >  2:2284511/1‑100 (MQ=255)
                          ggAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACg                                                                     <  1:459709/100‑1 (MQ=255)
                          ggAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACg                                                                     >  2:68362/1‑100 (MQ=255)
                           gAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTAc                                                                      <  2:3476875/98‑1 (MQ=255)
                            aaCGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCa                                                                   >  1:1419871/1‑100 (MQ=255)
                              cGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAg                                                                 <  1:2691977/100‑1 (MQ=255)
                               gCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACg                                                                     >  2:3676272/1‑95 (MQ=255)
                                cTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAggg                                                               >  2:109640/1‑100 (MQ=255)
                                cTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAggg                                                               >  1:1021579/1‑100 (MQ=255)
                                cTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCACCATTAAAAAAAACAGTATCGTGATATTTTACGCAAggg                                                               <  1:2299071/100‑1 (MQ=255)
                                cTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATAAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAggg                                                               >  2:2033809/1‑100 (MQ=255)
                                 ttAAATTATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATGTTTTACGCAAggg                                                               <  2:891910/99‑1 (MQ=255)
                                  tAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAgggn                                                              <  2:1164110/99‑2 (MQ=255)
                                  tAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTa                                                             >  2:3179565/1‑100 (MQ=255)
                                    aaTGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGc                                                           >  2:3347292/1‑100 (MQ=255)
                                     atgatgCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGcc                                                          >  1:2879029/1‑100 (MQ=255)
                                      tgatgCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGccc                                                         <  2:1633942/100‑1 (MQ=255)
                                         tgCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCa                                                      <  2:507155/100‑1 (MQ=255)
                                         tgCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCa                                                      <  1:2654979/100‑1 (MQ=255)
                                         tgCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCa                                                      >  2:1598855/1‑100 (MQ=255)
                                          gCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCAt                                                     >  1:3288991/1‑100 (MQ=255)
                                          gCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑AAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATg                                                    >  1:736862/1‑100 (MQ=255)
                                           cTTCATACATGGTTTATAAAAAGAATGAAACAGTATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCa                                                      <  2:1880970/98‑1 (MQ=255)
                                           cTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTTATCATTAAAAAAAAAAAGTATCGTGATATTTTACGCAAGGGTAGcgccaa                                                      >  2:792217/1‑97 (MQ=255)
                                              cATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTc                                                 >  2:3764604/1‑100 (MQ=255)
                                              cATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTc                                                 >  2:956357/1‑100 (MQ=255)
                                                tACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCAt                                               <  1:2131931/100‑1 (MQ=255)
                                                 acaatatttATAAAAAGAATGAAACAGGACTTATTGTCTGTTCAACATT‑‑AAAACAAACAGTATCGTGATATTTTACTCACGGGTAGCCCCCTTGTTCAtt                                              <  2:1055685/94‑1 (MQ=255)
                                                 aCAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTATCCCCCATGTTCAtt                                              >  2:693735/1‑100 (MQ=255)
                                                   aGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTc                                            <  2:3333108/100‑1 (MQ=255)
                                                      cTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTTATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCctc                                         >  2:1794149/1‑100 (MQ=255)
                                                      cTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCctc                                         <  1:141237/100‑1 (MQ=255)
                                                      cTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCctc                                         <  1:286901/100‑1 (MQ=255)
                                                         tataAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTATGCAAGGGTAGCCCCCATGTTCATTTCCTCTgg                                      <  1:3625559/100‑1 (MQ=255)
                                                           naaaaaGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGtt                                    >  2:418714/2‑100 (MQ=255)
                                                            aaaaaGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTa                                   <  1:2395819/100‑1 (MQ=255)
                                                            aaaaaGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTa                                   >  1:1889572/1‑100 (MQ=255)
                                                             aaaaGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAg                                  >  1:2839479/1‑100 (MQ=255)
                                                             aaaaGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGAGTAGCCCCCATGTTCATTTCCTCTGGTTa                                   >  2:3555419/1‑99 (MQ=255)
                                                             aaaaGAATGAAACAGGATTTATTGTATGTTCATCATT‑‑AAAACAAACAGTATCGTGATATTCTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTa                                   <  2:3810391/99‑1 (MQ=255)
                                                                 gAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATgg                              >  2:3688051/1‑100 (MQ=255)
                                                                  aaTGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGc                             <  1:2640921/100‑1 (MQ=255)
                                                                  aaTGAAACAGAATTTATTTTCTGGTCATAATTAAAAAAAACAGTATATGATATTTTACGCACGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATg                               <  2:1707889/97‑1 (MQ=255)
                                                                   aTGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCg                            <  2:1759336/100‑1 (MQ=255)
                                                                     gAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACCGTATCGTGATATTTTACGCAAGGGTAGCTCCCATGTTCATTTCCTCTGGTTAGATGGCGAc                          <  1:3626106/100‑1 (MQ=255)
                                                                     gAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGAc                          <  1:3283141/100‑1 (MQ=255)
                                                                       aaCAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTa                        <  1:218303/100‑1 (MQ=255)
                                                                          aGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGa                     <  1:789508/100‑1 (MQ=255)
                                                                          aGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGa                     <  2:905138/100‑1 (MQ=255)
                                                                          aGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGa                     <  1:579772/100‑1 (MQ=255)
                                                                          aGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGa                     >  1:3517182/1‑100 (MQ=255)
                                                                           ggATTTATTGTCTGTTCATCATTAAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTGTGa                     >  1:3473690/1‑100 (MQ=255)
                                                                           ggATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGaa                    <  2:3290341/100‑1 (MQ=255)
                                                                            gATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGaaa                   >  1:2826352/1‑100 (MQ=255)
                                                                            gATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGaa                    >  2:2710651/1‑99 (MQ=255)
                                                                             aTTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAAc                  <  2:1930033/100‑1 (MQ=255)
                                                                              tttATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACt                 <  2:2616204/100‑1 (MQ=255)
                                                                                tATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTAt               >  1:1386042/1‑100 (MQ=255)
                                                                                tATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTAt               >  2:3153623/1‑100 (MQ=255)
                                                                                   tGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATtct            >  1:1004516/1‑100 (MQ=255)
                                                                                    gTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATtctc           >  2:1010110/1‑100 (MQ=255)
                                                                                       tGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATTCTcgcg        <  1:247533/100‑1 (MQ=255)
                                                                                       tGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACACAAGGGTAACCCCAATGTTCATTTACGATGGTTAGATGGCGACTATGAAACTATTCTcgc         >  2:2145305/1‑99 (MQ=255)
                                                                                        gTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATTCTcgcgc       <  1:3420336/100‑1 (MQ=255)
                                                                                        gTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATTCTcgcgc       >  1:3329499/1‑100 (MQ=255)
                                                                                        gTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATTCTcgcgc       >  1:3094679/1‑100 (MQ=255)
                                                                                        gTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGACCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATTCTcgcgc       >  2:1545455/1‑100 (MQ=255)
                                                                                                  |                                                                                                
TCTGATGAAGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑AAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATTCTCGCGCGAATG  >  NC_000913/4494779‑4494970

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