Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,117,004 T→A 100% D73V (GAT→GTT)  glpK ← glycerol kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,117,0040TA100.0% 192.3 / ‑6.7 62D73V (GAT→GTT) glpKglycerol kinase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base A (33/28);  minor base G (0/1);  total (33/29)
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

CTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCA  >  NC_000913/4116858‑4117147
                                                                                                                                                  |                                                                                                                                               
ccaaaCCCGCACGTTTTAAATGCACGCAGATTTCTGCGGAACCACGGCACTGCCAAACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACg                                                                                                                                               <  2:463481/147‑1 (MQ=255)
 tAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACgg                                                                                                                                              >  1:861776/1‑149 (MQ=255)
         cACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACGTCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTTCCAGACAATAGTGGTTTCACGCTGGTTCGTAATGCCGATAGCAGCAATTTGAACGGAACTGata                                                                                                                                      <  1:600973/149‑1 (MQ=255)
          aCGTTTTAAATGCTCGCAGATTTCTGGGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGatat                                                                                                                                     <  2:71391/149‑1 (MQ=255)
            gTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGGTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACGGATAGCTGCAATTTGAACGGAACTGATATCg                                                                                                                                   >  2:441693/1‑149 (MQ=255)
              tttAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGCACTGCTATCGGc                                                                                                                                 >  2:384898/1‑149 (MQ=255)
                  aaTGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTc                                                                                                                             <  1:1086070/149‑1 (MQ=255)
                   aTGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGAAAAATTGAAAGCTGCAATTTGAACGGAACTGATATCGGCTTTCg                                                                                                                            <  1:1138934/149‑1 (MQ=255)
                      cTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCa                                                                                                                         <  1:524062/149‑1 (MQ=255)
                       tCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAg                                                                                                                        <  1:136052/149‑1 (MQ=255)
                       tCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAg                                                                                                                        <  1:90510/149‑1 (MQ=255)
                        cGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGCTAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGc                                                                                                                       >  1:903673/1‑149 (MQ=255)
                          cAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCAc                                                                                                                     >  2:623136/1‑149 (MQ=255)
                            gATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACtt                                                                                                                   <  2:35547/149‑1 (MQ=255)
                                tCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAATTGATATCGGCTTTCGCCAGCACTTCTAc                                                                                                               <  1:692254/149‑1 (MQ=255)
                                tCTGCGGTACCGCGTCACGGCCATACAATGGCGTTATAGATAGGCTTGACGTTTTTTTTGGACCATACAATGTTGGTTACACGTTGGTTCGTAATACCGGTAGCTGCAATTTGAGCGGCACTGATAGCGGCTTTCGCCAGCACTTCTAc                                                                                                               <  2:557983/149‑1 (MQ=255)
                                  tGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCCTCTCTACCa                                                                                                             >  2:1037342/1‑149 (MQ=255)
                                       tACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTg                                                                                                        <  1:673866/149‑1 (MQ=255)
                                        acgacgGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCCATTCTACCAGCGTgg                                                                                                       >  1:896685/1‑149 (MQ=255)
                                         cgacgGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTCCCAGCGTGGa                                                                                                      >  1:511104/1‑149 (MQ=255)
                                             ggCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCGGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCtt                                                                                                  >  2:353188/1‑149 (MQ=255)
                                                    cccgACAATTGCGTTATCGACAATCCCGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTgg                                                                                           >  2:255567/4‑149 (MQ=255)
                                                    ccAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTgg                                                                                           >  1:955076/1‑149 (MQ=255)
                                                      agtcgaTGGTGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGcc                                                                                         <  1:431011/144‑1 (MQ=255)
                                                      aggcAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTCTCCAGCCAATAGTGGTTTCACTCTGGTTCGTAATACCGATAGCTGCAACTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGcc                                                                                         <  2:291373/146‑1 (MQ=255)
                                                       gACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCAAGCTCGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGccc                                                                                        <  1:554071/149‑1 (MQ=255)
                                                        aCAATGTCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCa                                                                                       <  1:649600/149‑1 (MQ=255)
                                                        aCAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATCCCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCGCTTCTACCAGCGTGGAGATTCGGGTGGcccg                                                                                       >  2:608839/1‑148 (MQ=255)
                                                         cAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAg                                                                                      <  1:734821/149‑1 (MQ=255)
                                                         cAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGATAGCGGCTTTCGCCAGCACTTTTACCAGCGTGGAGCTTTGGGTGGCCCAg                                                                                      >  2:189965/1‑149 (MQ=255)
                                                           aTGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTTGTTCGTACTACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGAt                                                                                    <  2:862637/149‑1 (MQ=255)
                                                               cGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATCGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTcc                                                                                <  2:907237/149‑1 (MQ=255)
                                                                    cagatagGCTTGCCGGTTTCTTATTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAGTACCGATAGCTGCAATTTTGACGTTACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTCGAGCTTTGGGTGACCCAGATTTCCATTgg                                                                           <  1:728142/148‑1 (MQ=255)
                                                                       atagGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTc                                                                        <  1:756365/149‑1 (MQ=255)
                                                                         agGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCCATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCgt                                                                      >  1:605902/1‑149 (MQ=255)
                                                                            cTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGtt                                                                   <  1:796666/149‑1 (MQ=255)
                                                                             ttGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCGTTGGGGTGTGTTc                                                                  >  1:218189/1‑149 (MQ=255)
                                                                                 cGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTAcc                                                              >  2:1048207/1‑149 (MQ=255)
                                                                                      gcTTTTTCCCGGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTTCCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGCCGTGTTCTACCCAAcc                                                         <  2:353783/148‑1 (MQ=255)
                                                                                       cTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCCGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCATTGGGTGGCCCAGATTTCCCCTGGGTCGTGTTCTCCCCCAccc                                                        >  2:269841/1‑148 (MQ=255)
                                                                                        tttttCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTg                                                       >  2:679591/1‑149 (MQ=255)
                                                                                        tttttCCCAGACAATAGTGGTTTCACGCTGGTTCATAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGGTCTACCCAACCTg                                                       >  1:583021/1‑149 (MQ=255)
                                                                                         ttttCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGCTATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTgg                                                      >  1:44736/1‑149 (MQ=255)
                                                                                              ccAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTg                                                 >  2:368381/1‑149 (MQ=255)
                                                                                               cAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCTCCAGCACTTGTACCAGCGTGGAGCTTTGGGTGGCCCCGATTTCCAGTGGCTCGTGCTGTTACCAGCCTGGttttag                                                >  2:766323/1‑147 (MQ=255)
                                                                                               cAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCCGGTTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTgg                                                >  2:236433/1‑149 (MQ=255)
                                                                                               cAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAAGCTGGTTTggg                                                >  2:665176/1‑149 (MQ=255)
                                                                                                      tAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCGGATTTCCATTGGGTGGTGTTCTACCCAACCAGGTCGTGGGTAGAtt                                         >  2:178543/1‑149 (MQ=255)
                                                                                                        gTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTg                                       >  1:1017128/1‑149 (MQ=255)
                                                                                                         tGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGc                                      >  1:1064035/1‑149 (MQ=255)
                                                                                                            tttCACGCTGGTTCGTAATACCGATAGATGCAATTTGAACGGAACTGATATCGGCTTTCGACAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGGTTTCCATTGGGTCGTGTTCTACCCAACCTGGTCTTGGGTAGATTTGCTCa                                   >  2:821300/1‑149 (MQ=255)
                                                                                                             ttCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCaa                                  >  1:909852/1‑149 (MQ=255)
                                                                                                                        ttgTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGc                       <  2:1064673/147‑1 (MQ=255)
                                                                                                                            aaTACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGaca                   <  1:679591/149‑1 (MQ=255)
                                                                                                                             gtACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCAGTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGAGTTGCTCAAATTCGCGCTGCGacac                  <  1:1141909/148‑1 (MQ=255)
                                                                                                                             aTACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGacac                  <  1:1141227/149‑1 (MQ=255)
                                                                                                                                 cGAGAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTa              >  2:1060152/1‑149 (MQ=255)
                                                                                                                                       cTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGata        <  2:44736/149‑1 (MQ=255)
                                                                                                                                        tGCAATTTGAACGGAACTGATATCTGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGatat       >  1:801844/1‑149 (MQ=255)
                                                                                                                                         gCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATAtt      <  1:353188/149‑1 (MQ=255)
                                                                                                                                         gCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATAtt      >  2:431940/1‑149 (MQ=255)
                                                                                                                                             tttGAACGGAACTGATATCGGCTTTCGCCAGCAATACTACCAGCGTGAAACTTTGGATGGCCCAGATTCCCATAGGGTCGTGTTCTACCCAACCTGGATTCGAGCAGACATAATCAAAATCACCCTCCGAAACACAAATGATATTGGCa  >  2:321890/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                               
CTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCA  >  NC_000913/4116858‑4117147

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: ‑