Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,116,703 T→G 100% K173N (AAA→AAC glpK ← glycerol kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,116,7030TG100.0% 279.8 / ‑6.7 86K173N (AAA→AACglpKglycerol kinase
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base G (47/38);  total (47/39)
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.

GAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATTGCT  >  NC_000913/4116557‑4116843
                                                                                                                                                  |                                                                                                                                            
gaagaaCGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCAtgtt                                                                                                                                            >  1:1179936/1‑149 (MQ=255)
 aagaaCGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATGTTc                                                                                                                                           <  1:1242033/149‑1 (MQ=255)
 aagaaCGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATGTTc                                                                                                                                           <  1:1191554/149‑1 (MQ=255)
  agaaCGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATGTTcc                                                                                                                                          <  2:1322314/149‑1 (MQ=255)
   gaaCGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACCCGGCCCTGAGTCATGTTccg                                                                                                                                         >  1:1151535/1‑148 (MQ=255)
     acgacgCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATGTTCCAGa                                                                                                                                       <  1:31443/149‑1 (MQ=255)
      cgacgCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATGTTCCAGAt                                                                                                                                      >  2:1004495/1‑149 (MQ=255)
        acgCCCTTCTGGCCGCATCTCGCGCGGAATATCCAGCCCTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGGTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATGTTCCAGATaa                                                                                                                                    <  2:1258780/149‑1 (MQ=255)
           cACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATGTTCCAGATAAGcc                                                                                                                                 <  1:503437/149‑1 (MQ=255)
             cTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATGTTCCAGATAAGCCAc                                                                                                                               >  2:248587/1‑149 (MQ=255)
              ttCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTCAGCGGTCACATGGACCCGGCCCTGAGCCATGTTCCAGATAAGcccca                                                                                                                              >  1:26303/1‑146 (MQ=255)
              ttCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATGTTCCAGCTAAGCCCCg                                                                                                                              >  1:348743/1‑149 (MQ=255)
              ttCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATGTTCCAGATAAGCCACg                                                                                                                              >  2:820611/1‑149 (MQ=255)
              ttCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATGGTCAAGATAGGCCACg                                                                                                                              >  1:581357/1‑149 (MQ=255)
               tCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATGTTCCAGATAAGCCACGt                                                                                                                             >  1:268564/1‑149 (MQ=255)
                  ggCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGACGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATGTTCCAGATAAGCCACGTATc                                                                                                                          >  1:915994/1‑149 (MQ=255)
                    cagcaTCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATGTTCCAGATAAGCCACGTATCaa                                                                                                                        <  2:1147182/149‑1 (MQ=255)
                     agcaTCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATGTTCCAGATAAGCCACGTATCAAc                                                                                                                       <  1:1231999/149‑1 (MQ=255)
                       caTCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCAGGTTCCAGATAAGCCGCGTATCAACCg                                                                                                                     >  1:1042064/1‑149 (MQ=255)
                       caTCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTAGAACAACATAGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATGTTCCAGATAAGCCACGTATCAACCg                                                                                                                     >  1:1084419/1‑149 (MQ=255)
                         tctcGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATGTTCCAGATAAGCCACGTATCAACCGTa                                                                                                                   >  2:244608/1‑149 (MQ=255)
                           gcgcgcgGAATATCCAGCGCTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATGTTCCAGATAAGCCACGTATCAACCGTAcc                                                                                                                 <  1:1017892/148‑1 (MQ=255)
                              cgcgGAATATCCACCACTTCCAGCACTTTGTCGGCCCAGGCCAGGGTATGGATGTTGAACAAAATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATGTTCCAGATAAGCCACGTATCAACCGTACCaaa                                                                                                              <  1:916322/149‑1 (MQ=255)
                                     tatCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGTTGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATGTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAAt                                                                                                       <  2:435102/149‑1 (MQ=255)
                                        ccAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATGTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTca                                                                                                    >  2:114884/1‑149 (MQ=255)
                                         cAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATGTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTcac                                                                                                   >  1:1230672/1‑149 (MQ=255)
                                         cAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGAATGGATGTTGAACAACATGATACGAGAGGCGTGGGTGTAATCGGTAACATGGACACGGCCCTGAGTCATGTTCCAGAAAAGCCACGAAGCAACCGTACCAAACAGCAATTcac                                                                                                   >  2:1103983/1‑149 (MQ=255)
                                           gCACTTCCAGCATTTTGCCGTCCCAGTCCAGGCTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCAGGTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTcacca                                                                                                 <  1:340400/149‑1 (MQ=255)
                                               ttCCAGCATTTTGTCGTCCCAGTCCATGGTATGGTTGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTGACATGGACACGGCCCTGAGTCATGTTCCAGGTAAGCAACGTACACACCGTACCAAGCAGGAATTGACCccgag                                                                                             >  2:266039/1‑148 (MQ=255)
                                                tCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGACCAACATGGTACGAGAGGCGTTGGTGTAAACGGTCACATGGCCACGGCCCTGAGTCATGTTCCAGACAACCCACGTATCAACCGTACCAAACAGCAATTCACCacgacg                                                                                            >  2:515233/1‑149 (MQ=255)
                                                 ccAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATGTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGt                                                                                           <  2:132676/149‑1 (MQ=255)
                                                  cAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATGTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTg                                                                                          >  2:1008814/1‑149 (MQ=255)
                                                  cAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATGGTGCGAATAAGCCGCGTATCAACCGTACCACACAGCAATTCACCaagacgcg                                                                                          >  2:1140813/1‑147 (MQ=255)
                                                  cAGCATTGTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATGTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTg                                                                                          <  1:217281/149‑1 (MQ=255)
                                                   aGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATGTTCCAGATCAGCCCCGTATCAACCGTACCAAACAGCCATTCACCACGCCGTGc                                                                                         >  2:333514/1‑149 (MQ=255)
                                                        tttGTCGTCCCAGTCCAGGGTATGGATGTTCAACAACATGGCACGCGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATGTTCCAGACAAGCCACGTATCAACCGTACCAAACACAACTTCACCACGACGTGcctgcc                                                                                    >  2:946825/1‑144 (MQ=255)
                                                               tcccctgccAGTGTATGGACGTTGAACAACATGGTACGAGCGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATGTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAAGTCACCACGACGTGCACGCTCGCgaga                                                                             <  2:7878/142‑1 (MQ=255)
                                                               tcCCAGTCCAGGGTATGGATGTTGACCAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATGTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCgaga                                                                             >  1:675449/1‑149 (MQ=255)
                                                               tcCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGTACACGGCCCTGAGTCATGTTCCAGGTAGGCCCAGTTTTTGCTGGACCAAACAGCAATTCACACCGCCGTGCACGCCAGCgagc                                                                             >  2:478616/1‑148 (MQ=255)
                                                               tcCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATGTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTcgccgga                                                                             >  2:238298/1‑145 (MQ=255)
                                                               tcCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGGCATGGTCCAGATAAGCCACGGATCAACCGTACCAAACAGCACTTGACCACGACGTGCACGCTCGCgagc                                                                             >  1:95373/1‑148 (MQ=255)
                                                                 ccAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATGTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCCCGACGTGCACGCTCGCGAGAGc                                                                           >  1:76109/1‑149 (MQ=255)
                                                                 ccAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATGTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGc                                                                           >  1:1091982/1‑149 (MQ=255)
                                                                   aGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATGTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCt                                                                         >  1:389544/1‑149 (MQ=255)
                                                                     tCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATGTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCTTc                                                                       >  1:862487/1‑149 (MQ=255)
                                                                      ccAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATGTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCTTcc                                                                      <  1:670862/149‑1 (MQ=255)
                                                                       cAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGATCACATGGACACGGCCCTGAGTCATGTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTGCCCACGACGGGCACGCTCGCGAGAGCCTTCca                                                                     >  1:966611/1‑149 (MQ=255)
                                                                          ggTATGGATGTTGAACAACATGGTACAAGAGGCGTAGGAGTAATCGGTTACATGGACACGGCCCTGAGTCATGTTCCAGATAAGCCACGTATCAAACGTACCAAAAAGCACTTCCCCACGACGTGCACGCTCGCGGAAGCCTTCCCCAt                                                                  >  2:297323/1‑149 (MQ=255)
                                                                             aTGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATGTTCCAGATAAGCCACGTATCAACCGTACCAAACCGCAATTCACCACGACGTGCACGCTCGCGAGAGCCTTACACATGGt                                                               >  2:637610/1‑149 (MQ=255)
                                                                              tGGATTTTGAACAACGCGGTACGAGAGGCGTTGGTGTATTCGGTCACATGGCCACGGCCCTGATTCATGTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACTCTCGCGAGAGCCTTCCACATGGTc                                                              <  2:356857/149‑1 (MQ=255)
                                                                               ggATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATGTTCCAGATAAGCCACGTATCCACCGTACCAAACGGCAATTCACCACGACGTGCACGCTCGCGAGAGCCTTCCACATGGTCg                                                             >  1:646876/1‑149 (MQ=255)
                                                                                     tGATCAACGTGGTACGAGAGGCGTTTGTGGAATCGGTCACATCGTCACGGCCCTGAGTCATGTTCCAGATAAGCCACGTACCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGCTTTACACATTGTCGAGGATc                                                       <  2:358497/149‑1 (MQ=255)
                                                                                     tGAACCACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATGTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATc                                                       <  2:358025/149‑1 (MQ=255)
                                                                                       aacaacATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATGTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCca                                                     <  2:802361/149‑1 (MQ=255)
                                                                                       aacaacATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATGTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCca                                                     <  1:1289746/149‑1 (MQ=255)
                                                                                        acaacaGGGTCCGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATGTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCcac                                                    <  1:174206/149‑1 (MQ=255)
                                                                                         caacaTGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATGTTCCAGATAAGCCCCGTATCAACCGTACCAAACAGCAATTCACCACGCCGTGCACGCTCGCGAGAGACTTCCAAATAGTCGAGGATCcact                                                   >  1:589374/1‑149 (MQ=255)
                                                                                            caTGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACGCGGCCCTGAGTCATGTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCTTCCTCATGGTCGAGGATCcacttca                                                <  1:958028/149‑1 (MQ=255)
                                                                                              tGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATGTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGCGCCTTCCACATGGTCGAGGATCcccttcact                                              >  2:783618/1‑149 (MQ=255)
                                                                                               ggTACGAGAGGCGTTGGTGTAATCGGCCACATGCACACGGCCCTGAGTCATGTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCcacttcactt                                             <  2:1098008/149‑1 (MQ=255)
                                                                                                gTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATGTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTACCttt                                            >  2:1206864/1‑149 (MQ=255)
                                                                                                   cGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATGTTCCAGATAAGCCACGTATCGACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGt                                         >  1:1032802/1‑149 (MQ=255)
                                                                                                   cGAGAGCGGTGGGCGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCGCCACGGCGTGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGt                                         <  2:581357/149‑1 (MQ=255)
                                                                                                    gagaGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATGTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTg                                        >  1:297446/1‑149 (MQ=255)
                                                                                                     agagGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATGTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGc                                       <  2:1128090/149‑1 (MQ=255)
                                                                                                            ttGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATGTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGaaaa                                <  1:98249/149‑1 (MQ=255)
                                                                                                            ttGGTGTAAGCGGTCACATGCACACGGCCCTGAGTCATGTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGaaaa                                <  1:414423/149‑1 (MQ=255)
                                                                                                            gtGGTGTAATCGGCCACAGGGACACGGCCCTGAGTCATGTGCCAGATAAGCCGCGTATCAACCGTACCAAACAGCAATTGGCCACGACGTTCACGCTCGCGAGAGCCTTCTACATGGTCTAGGATCCACTTCACTTTGGTGCCAGaaaa                                <  2:51051/148‑1 (MQ=255)
                                                                                                                tgtAATCGGTCACATGGCCACGGCCCTGATTCATGTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTAc                            <  2:1151535/149‑1 (MQ=255)
                                                                                                                 gtAATCGGTCACATGGACACGGCCCTGAGTCATGTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACg                           <  1:248587/149‑1 (MQ=255)
                                                                                                                  tAATCGGTCACATGGACACGGCCCTGAGTCATGTTCCAGATAAGCCACGTATCAACCGTATTAATCAGCAAGTCGCCACGACGTGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACgg                          >  2:289404/1‑149 (MQ=255)
                                                                                                                  tAATCGGTCACATGGACACGGCCCTGAGTCATGTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACgg                          <  1:244608/149‑1 (MQ=255)
                                                                                                                  tAATCGGTCACATGGACACGGCCCTGAGTCATGTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACgg                          >  1:1214631/1‑149 (MQ=255)
                                                                                                                  tAATCGGTCACATGGACACGGCCCTGAGTCATGTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACACTCGCGAGAGCCTTCCACATGATCGAGGATCCACTTCACTTTGGTGCCCGAAAAGTACgg                          >  2:148242/1‑149 (MQ=255)
                                                                                                                  gaaTCGGTTACATGGACACGGCCCTGAGTCATGTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACgg                          <  1:478616/148‑1 (MQ=255)
                                                                                                                   aaTCGGTCACATGGACACGGCCCTGAGTCATGTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGCAAAGTACggg                         >  1:618505/1‑149 (MQ=255)
                                                                                                                   aaTCGGTCACATGGACACGGCCCTGAGTCATGTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACggg                         >  1:790212/1‑149 (MQ=255)
                                                                                                                        gTCACATGGACACGGCCCTGAGTCATGTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACGGGTCAAt                    >  1:786191/1‑149 (MQ=255)
                                                                                                                         tCACATGGACACGGCCCTGAGTCATGTTCCAGATAAGCCACGTATCAACCGTACCAAACGGCAATTCACCACGACGTGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACGGGTCAATc                   <  1:1317048/149‑1 (MQ=255)
                                                                                                                         acacaTGGACGCGGCCCGGAGTCATGTTCCAGATAAGCCACGTATCAACCGTCCCAAACAGCAATTCACCACGAAGTGCACGCCCGCGAGAGCCTTCAACATGGTCGAGGATCCACTCCTCTTTGGTGCCAGAAAAGTACGGGTCAATc                   <  2:26303/148‑1 (MQ=255)
                                                                                                                          cacaTGGACACGGCCCTGAGTCATGTTCCAGATAAGCCACGTATCAACCGCACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTATGGGTCAATCa                  <  2:1159699/149‑1 (MQ=255)
                                                                                                                           acttGGACACGGCCCTGAGTCATGTTCCAGATAGGCCACGTATCAACCGTACTAAACAGCAATTCACCACGACGTGCACGCTTGCGAGAGCCTTCCACGTGGTCGGGGATTCACTTCACTTTGGTGCCGGAAAAGTACGGGTCATTCAc                 >  2:481926/4‑149 (MQ=255)
                                                                                                                           acaTGGACACCGCCCTGAGTCATGTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACGGGTCAATCAc                 <  1:168157/149‑1 (MQ=255)
                                                                                                                              tGGACACGGCCCTGAGTCATGTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACGGGTCAATCACCAg              <  2:389544/149‑1 (MQ=255)
                                                                                                                                     ggCCCTGAGTCAGGTCCCAGATAAGCCACGTATCAACCGTCCAAAGCCGCAACTCACCACGACGTGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTa       <  2:1292364/149‑1 (MQ=255)
                                                                                                                                        ccTGAGTCATGTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGCTCCACGTCCCTTTGGTGCCAGAAAAGTAGGGGTCAATCACCCGACCGGTATtg    >  2:395752/1‑149 (MQ=255)
                                                                                                                                          tGAGTCATGTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATtgct  >  2:1190760/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                            
GAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATTGCT  >  NC_000913/4116557‑4116843

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