Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,116,274 G→T 8.2% L317M (CTG→ATG)  glpK ← glycerol kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,116,2740GT 8.2% 323.7 / 13.3 98L317M (CTG→ATG) glpKglycerol kinase
Reads supporting (aligned to +/- strand):  ref base G (36/54);  new base T (4/4);  total (40/58)
Fisher's exact test for biased strand distribution p-value = 7.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.01e-01

Only 100 of 101 total aligned reads displayed.

GAAAATCGCCCCGCGCGCATACGGGTCCCAGTACGGCGCACCCAGCCCGGTAAATGCCGGAACCACATACACACCATTGGTGTTTTGCACTTTGGTGGCGAAATATTCGGAATCGTAGGCGTCGTTAATCAACTTCATTTCATCGCGCAGCCACTGAATGGATGCGCCTGCCATAAACACCGCACCTTCCAACGCATAGTTCACTTCGCCAGTCGGGCCGCAGGCGATGGTGGTCAGCAGGCCGTTTTCTGATTTCACCGCTTTCTCGCCAGTGTTCATCAGCATAAAGCAGCCAGTGC  >  NC_000913/4116125‑4116423
                                                                                                                                                     |                                                                                                                                                     
gAAAATCGCCCCGCGCGCATACGGGTCCCAGTACGGCGCACCCAGCCCGGTAAATGCCGGAACCACATACACACCATTGGTGTTTTGCACTTTGGTGGCGAAATATTCGGAATCGTAGGCGTCGTTAATCAACTTCATTTCATCGCGCAg                                                                                                                                                       >  2:1310438/1‑150 (MQ=255)
gAAAATCGCCCCGCGCGCATACGGGTCCCAGTACGGCGCACCCAGCCCGGTAAATGCCGGAACCACATACACACCATTGGTGTTTTGCACTTTGGTGGCGAAATATTCGGAATCGTAGGCGTCGTTAATCAACTTCATTTCATCGCGCAg                                                                                                                                                       >  2:529571/1‑150 (MQ=255)
   aaTCGCCCCGCGCGCATACGGGTCCCAGTACGGCGCACCCAGCCCGGTAAATGCCGGAACCACATACACACCATTGGTGTTTTGCACTTTGGTGGCGAAATATTCGGAATCGTAGGCGTCGTTAATCAACTTCATTTCATCGCGCAGCCa                                                                                                                                                    >  2:1492880/1‑150 (MQ=255)
      cGCCCCGCTCGCATACGGGTCCCAGTACGGCGCACCCAGCCCGGTAAATGCCGGAACCACATACACACCATTGGTGTTTTGCACTTTTGTGGCGAAATATTCGGAATCGTAGGCGTCGTTAATCAACTTCATTTCATCGCGCAGCCACTg                                                                                                                                                 <  2:1404697/150‑1 (MQ=255)
      cGCCCCGCGCGCATACGGGTCCCAGTACGGCGCACCCAGCCCGGTAAATGCCGGAACCACATACACACCATTGGTGTTTTGCACTTTGGTGGCGAAATATTCGGAATCGTAGGCGTCGTTAATCAACTTCATTTCATCGCGCAGCCACTg                                                                                                                                                 <  1:1310438/150‑1 (MQ=255)
       gCCCCGCGCGCATACGGGTCCCAGTACGGCGCACCCAGCCCGGTAAATGCCGGAACCACATACACACCATTGGTGTTTTGCACTTTGGTGGCGAAATATTCGGAATCGTAGGCGTCGTTAATCAACTTCATTTCATCGCGCAGCCACTGa                                                                                                                                                >  1:1059592/1‑150 (MQ=255)
        ccccGCGCGCATACGGGTCCCAGTACGGCGCACCCAGCCCGGTAAATGCCGGAACCACATACACACCATTGGTGTTTTGCACTTTGGTGGCGAAATATTCGGAATCGTAGGCGTCGTTAATCAACTTCATTTCATCGCGCAGCCACTGaa                                                                                                                                               <  1:672907/150‑1 (MQ=255)
          ccGCGCGCATACGGGTCCCAGTACGGCGCACCCAGCCCGGTAAATGCCGGAACCACATACACACCATTGGTGTTTTGCACTTTGGTGGCGAAATATTCGGAATCGTAGGCGTCGTTAATCAACTTCATTTCATCGCGCAGCCACTGAATg                                                                                                                                             <  2:952906/150‑1 (MQ=255)
           cgcgcgCATACGGGTCCCAGTACGGCGCACCCAGCCCGGTAAATGCCGGAACCACATACACACCATTGGTGTTTTGCACTTTGGTGGCGAAATATTCGGAATCGTAGGCGTCGTTAATCAACTTGATTTCATCGCGCAGCCACTGAATgg                                                                                                                                            <  2:1605428/150‑1 (MQ=255)
                  aTACGGGTCCCAGTACGGCGCACCCAGCCCGGTAAATGCCGGAACCACATACACACCATTGGTGTTTTGCACTTTGGTGGCGAAATATTCGGAATCGTAGGCGTCGTTAATCAACTTCATTTCATCGCGCAGCCACTGAATGGATGCGcc                                                                                                                                     <  2:1091449/150‑1 (MQ=255)
                   tACGGGTCCCAGTACGGCGCACCCAGCCCGGTAAATGCCGGAACCACATACACACCATTGGTGTTTTGCACTTTGGTGGCGAAATATTCGGAATCGTAGGCGTCGTTAATCAACTTCATTTCATCGCGCAGCCACTGAATGGATGCGCCt                                                                                                                                    >  2:704548/1‑150 (MQ=255)
                   tACGGGTCCCAGTACGGCGCACCCAGCCCGGTAAATGCCGGAACCACATACACACCATTGGTGTTTTGCACTTTGGTGGCGAAATATTCGGAATCGTAGGCGTCGTTAATCAACTTCATTTCATCGCGCAGCCACTGAATGGATGCGCCt                                                                                                                                    >  2:734788/1‑150 (MQ=255)
                     cGGGTCCCAGTACGGCGCACCCAGCCCGGTAAATGCCGGAACCACATACACACCATTGGTGTTTTGCACTTTGGTGGCGAAATATTCGGAATCGTAGGCGTCGTTAATCAACTTCATTTCATCGCGCAGCCACTGAATGGATGCGCCTGc                                                                                                                                  >  2:534452/1‑150 (MQ=255)
                     cGGGTCCCAGTACGGCGCACCCAGCCCGGTAAATGCCGGAACCACATACACACCATTGGTGTTTTGCACTTTGGTGGCGAAATATTCGGAATCGTAGGCGTCGTTAATCAACTTCATTTCATCGCACATCCATTCAGTGGCTgcgcgatg                                                                                                                                  <  1:272919/150‑5 (MQ=255)
                       ggTCCCAGTACGGCGCACCCAGCCCGGTAAATGCCGGAACCACATACACACCATTGGTGTTTTGCACTTTGGTGGCGAAATATTCGGAATCGTAGGCGTCGTTAATCAACTTCATTTCATCGCGCAGCCACTGAATGGATGCGCCTGCCa                                                                                                                                <  1:197495/150‑1 (MQ=255)
                       ggTCACAGTACGGCGCACCCAGCCCGGTAAATGCCGGAACCACATACACACCATTGGTGTTTTGCACTTTGGTGGCGAAATATTCGGAATCGTAGGCGTCGTTAATCAACTTCATTTCATCGCGCAGCCACTGAATGGATGCGCCTGCCa                                                                                                                                >  2:1179528/1‑150 (MQ=255)
                        gTCCCAGTACGGCGCACCCAGCCCGGTAAATGCCGGAACCACATACACACCATTGGTGTTTTGCACTTTGGTGGCGAAATATTCGGAATCGTAGGCGTCGTTAATCAACTTCATTTCATCGCGCATCCATTCAGTGGCTgcgcgatgaaa                                                                                                                               <  1:1340086/150‑8 (MQ=255)
                        gTCCCAGTACGGCGCACCCAGCCCGGTAAATGCCGGAACCACATACACACCATTGGTGTTTTGCACTTTGGTGGCGAAATATTCGGAATCGTAGGCGTCGTTAATCAACTTCATTTCATCGCGCATCCATTCAGTGGCTgcgcgatgaaa                                                                                                                               <  1:1491983/150‑8 (MQ=255)
                        gTCCCAGTACGGCGCACCCAGCCCGGTAAATGCCGGAACCACATACACACCATTGGTGTTTTGCACTTTGGTGGCGAAATATTCGGAATCGTAGGCGTCGTTAATCAACTTCATTTCATCGCGCATCCATTCAGTGGCTgcgcgatgaaa                                                                                                                               <  2:1225808/150‑8 (MQ=255)
                          cccAGTACGGCGCACCCAGCCCGGTAAATGCCGGAACCACATACACACCATTGGTGTTTTGCACTTTGGTGGCGAAATATTCGGAATCGTAGGCGTCGTTAATCAACTTCATTTCATCGCGCAGCCACTGAATGGATGCGCCTGCCATaa                                                                                                                             >  1:77229/1‑150 (MQ=255)
                          cccAGTACGGCGCACCCAGCCCGGTAAATGCCGGAACCACATACACACCATTGGTGTTTTGCACTTTGGTGGCGAAATATTCGGAATCGTAGGCGTCGTTAATCAACTTCATTTCATCGCGCAGCCACTGAATGGATGCGCCTGCCATaa                                                                                                                             >  1:112050/1‑150 (MQ=255)
                                aCGGCGCACCCAGCCCGGTAAATGCCGGAACCACATACACACCATTGGTGTTTTGCACTTTGGTGGCGAAATATTCGGAATCGTAGGCGTCGTTAATCAACTTCATTTCATCGCGCAGCCACTGAATGGATGCGCCTGCCATAAACACCg                                                                                                                       <  1:802616/150‑1 (MQ=255)
                                 cGGCGCACCCAGCCCGGTAAATGCCGGAACCACATACACACCATTGGTGTTTTGCACTTTGGTGGCGAAATATTCGGAATCGTAGGCGTCGTTAATCAACTTCATTTCATCGCGCAGCCACTGAATGGATGCGCCTGCCATAAACACCGc                                                                                                                      <  2:669908/150‑1 (MQ=255)
                                 cGGCGCACCCAGCCCGGTAAATGCCGGAACCACATACACACCATTGGTGTTTTGCACTTTGGTGGCGAAATATTCGGAATCGTAGGCGTCGTTAATCAACTTCATTTCATCGCGCAGCCACTGAATGGATGCGCCTGCCATAAACACCGc                                                                                                                      <  2:283444/150‑1 (MQ=255)
                                  ggCGCACCCAGCCCGGTAAATGCCGGAACCACATACACACCATTGGTGTTTTGCACTTTGGTGGCGAAATATTCGGAATCGTAGGCGTCGTTAATCAACTTCATTTCATCGCGCAGCCACTGAATGGATGCGCCTGCCATAAACACCGcc                                                                                                                     >  1:419449/1‑149 (MQ=255)
                                   gcgcACCCAGCCCGGTAAATGCCGGAACCACATACACACCATTGGTGTTTTGCACTTTGGTGGCGAAATATTCGGAATCGTAGGCGTCGTTAATCAACTTCATTTCATCGCGCAGCCACTGAATGGATGCGCCTGCCATAAACACCGCAc                                                                                                                    <  1:1197865/150‑1 (MQ=255)
                                       aCCCAGCCCGGTAAATGCCGGAACCACATACACACCATTGGTGTTTTGCACTTTGGTGGCGAAATATTCGGAATCGTAGGCGTCGTTAATCAACTTCATTTCATCGCGCAGCCACTGAATGGATGCGCCTGCCATAAACACCGCACCTTc                                                                                                                >  2:228229/1‑150 (MQ=255)
                                        cccATCCCGGTAAATGCCGGAACCACATACACACCCTTGGTGTTTTGCACTTTGGTGGCGAAATATTCGGACTCGTAGGCGTCGTTAATCAACTTCATTTCCTCGCGCAGCCACTGAATGGATGCGCCTGCCATAAACACCGCACCTTcc                                                                                                               <  2:1085171/150‑1 (MQ=255)
                                        cccAGCCCGGTAAATGCCGGAACCACATACACACCATTGGTGTTTTGCACTTTGGTGGCGAAATATTCGGAATCGTAGGCGTCGTTAATCAACTTCATTTCATCGCGCAGCCACTGAATGGATGCGCCTGCCATAAACACCGCACCTTcc                                                                                                               <  1:734788/150‑1 (MQ=255)
                                         ccAGCCCGGTAAATGCCGGAACCACATACACACCATTGGTGTTTTGCACTTTGGTGGCGAAATATTCGGAATCGTAGGCGTCGTTAATCAACTTCATTTCATCGCGCAGCCACTGAATGGATGCGCCTGCCATAAACACCGCACCTTc                                                                                                                >  2:1140840/1‑148 (MQ=255)
                                         ccAGCCCGGTAAATGCCGGAACCACATACACACCATTGGTGTTTTGCACTTTGGTGGCGAAATATTCGGAATCGTAGGCGTCGTTAATCAACTTCATTTCATCGCGCAGCCACTGAATGGATGCGCCTGCCATAAACACCGCACCTTc                                                                                                                >  2:1135218/1‑148 (MQ=255)
                                         ccAGCCCGGTAAATGCCGGAACCACATACACACCATTGGTGTTTTGCACTTTGGTGGCGAAATATTCGGAATCGTAGGCGTCGTTAATCAACTTCATTTCATCGCGCAGCCACTGAATGGATGCGCCTGCCATAAACACCGCACCTTc                                                                                                                <  1:1140840/148‑1 (MQ=255)
                                         ccAGCCCGGTAAATGCCGGAACCACATACACACCATTGGTGTTTTGCACTTTGGTGGCGAAATATTCGGAATCGTAGGCGTCGTTAATCAACTTCATTTCATCGCGCAGCCACTGAATGGATGCGCCTGCCATAAACACCGCACCTTc                                                                                                                <  1:1135218/148‑1 (MQ=255)
                                            gCCCGGTAAATGCCGGAACCACATACACACCATTGGTGTTTTGCACTTTTGTGGCGAAATATTCGGAATCGTAGGCGTCGTTAATCAACTTCATTTCATCGCGCAGCCACTGAATGGATGCGCCTGCCATAAACACCGCACCTTCCAACg                                                                                                           <  1:1025111/150‑1 (MQ=255)
                                              ccGGTAAATGCCGGAACCACATACACACCATTGGTGTTTTGCACTTTGGTGGCGAAATATTCGGAATCGTAGGCGTCGTTAATCAACTTCATTTCATCGCGCAGCCACTGAATGGATGCGCCTGCCATAAACACCGCACCTTCCAACGCa                                                                                                         >  1:838651/1‑150 (MQ=255)
                                                    aaTGCCGGAACCACATACACACCATTGGTGTTTTGCACTTTGGTGGCGAAATATTCGGAATCGTAGGCGTCGTTAATCAACTTCATTTCATCGCGCAGCCACTGAATGGATGCGCCTGCCATAAACACCGCACCTTCCAACGCATAGTTc                                                                                                   <  2:486573/150‑1 (MQ=255)
                                                    aaTGCCGGAACCACATACACACCATTGGTGTTTTGCACTTTGGTGGCGAAATATTCGGAATCGTAGGCGTCGTTAATCAACTTCATTTCATCGCGCAGCCACTGAATGGATGCGCCTGCCATAAACACCGCACCTTCCAACGCATAGTTc                                                                                                   >  1:1014287/1‑150 (MQ=255)
                                                     aTGCCGGAACCACATACACACCATTGGTGTTTTGCACTTTGGTGGCGAAATATTCGGAATCGTAGGCGTCGTTAATCAACTTCATTTCATCGCGCAGCCACTGAATGGATGCGCCTGCCATAAACACCGCACCTTCCAACGCATAGTTCa                                                                                                  <  1:1505622/150‑1 (MQ=255)
                                                     aTGCCGGAACCACATACACACCATTGGTGTTTTGCACTTTGGTGGCGAAATATTCGGAATCGTAGGCGTCGTTAATCAACTTCATTTCATCGCGCAGCCACTGAATGGATGCGCCTGCCATAAACACCGCACCTTCCAACGCATAGTTCa                                                                                                  >  1:245966/1‑150 (MQ=255)
                                                       gCCGGAACCACATACACACCATTGGTGTTTTGCACTTTGGTGGCGAAATATTCGGAATCGTAGGCGTCGTTAATCAACTTCATTTCATCGCGCAGCCACTGAATGGATGCGCCTGCCATAAACACCGCACCTTCCAACGCATAGTTCACt                                                                                                >  1:605284/1‑150 (MQ=255)
                                                                  aTACACACCATTGGTGTTTTGCACTTTGGTGGCGAAATATGCGTAATCGTAGGCGTCGTTAATCAACTTCATTTCATCGCGCAGCCACTGAATGGATGCGCCTGCCATAAACACCGCACCTTCCAACGCATAGTTCACTTCGCCAGTCgg                                                                                     <  2:434709/150‑1 (MQ=255)
                                                                    acacacCATTGGTGTTTTGCACTTTGGTGGCGAAATATTCGGAATCGTAGGCGTCGTTAATCAACTTCATTTCATCGCGCAGCCACTGAATGGATGCGCCTGCCATAAACACCGCACCTTCCTACGCATAGTTCACTTCGCCAGTCGGGc                                                                                   <  1:624459/150‑1 (MQ=255)
                                                                       cacCATTGGTGTTTTGCACTTTGGTGGCGAAATATTCGGAATCGTAGGCGTCGTTAATCAACTTCATTTCATCGCGCAGCCACTGAATGGATGCGCCTGCCATAAACACCGCACCTTCCAACGCATAGTTCACTTCGCCAGTCGGGCCGc                                                                                <  2:567433/150‑1 (MQ=255)
                                                                       cacCATTGGTGTTTTGCACTTTGGTGGCGAAATATTCGGAATCGTAGGCGTCGTTAATCAACTTCATTTCATCGCGCAGCCACTGAATGGATGCGCCTGCCATAAACACCGCACCTTCCAACGCATAGTTCACTTCGCCAGTCGGGCCGc                                                                                >  2:298041/1‑150 (MQ=255)
                                                                        acCATTGGTGTTTTGCACTTTGGTGGCGAAATATTCGGAATCGTAGGCGTCGTTAATCAACTTCATTTCATCGCGCAGCCACTGAATGGATGCGCCTGCCATAAACACCGCACCTTCCAACGCATAGTTCACTTCGCCAGTCGGGCCGCa                                                                               <  1:1576945/150‑1 (MQ=255)
                                                                               gtgtTTTGCACTTTGGTGGCGAAATATTCGGAATCGTAGGCGTCGTTAATCAACTTCATTTCATCGCGCAGCCACTGAATGGATGCGCCTGCCATAAACACCGCACCTTCCAACGCATAGTTCACTTCGCCAGTCGGGCCGCAGGCGAtg                                                                        <  1:114424/150‑1 (MQ=255)
                                                                                ggtTTTGCACTTTGGTGGCGAAATATTCGGAATCGTAGGCGGCGTTAATCAACTTCATTTCATCGCGCAGCCACTGAATGGATGCGCCTGCCATAAACACCGCACCTTCCAACGCATAGTTCACTTCGCCAGTCGGGCCGCAGGCGAtgg                                                                       <  1:1594057/149‑1 (MQ=255)
                                                                                 gtTTTGCACTTTGGTGGCGAAATATTCGGAATCGTAGGCGTCGTTAATCAACTTCATTTCATCGCGCAGCCACTGAATGGATGCGCCTGCCATAAACACCGCACCTTCCAACGCATAGTTCACTTCGCCAGTCGGGCCGCAGGCGAtggt                                                                      <  1:1426505/150‑1 (MQ=255)
                                                                                  ttttGCACTTTGGTGGCGAAATATTCGGAATCGTAGGCGTCGTTAATCAACTTCATTTCATCGCGCAGCCACTGAATGGATGCGCCTGCCATAAACACCGCACCTTCCAACGCATAGTTCACTTCGCCAGTCGGGCCGCAGGCGAtggtg                                                                     >  1:443724/1‑150 (MQ=255)
                                                                                  ttttGCACTTTGGTGGCGAAATATTCGGAATCGTAGGCGTCGTTAATCAACTTCATTTCATCGCGCAGCCACTGAATGGATGCGCCTGCCATAAACACCGCACCTTCCAACGCATAGTTCACTTCGCCAGTCGGGCCGCAGGCGAtggtg                                                                     <  1:1492880/150‑1 (MQ=255)
                                                                                   tttGCACTTTGGTGGCGAAATATTCGGAATCGTAGGCGTCGTTAATCAACTTCATTTCATCGCGCAGCCACTGAATGGATGCGCCTGCCATAAACACCGCACCTTCCAACGCATAGTTCACTTCGCCAGTCGGGCCGCAGGCGAtggtgg                                                                    <  2:1014287/150‑1 (MQ=255)
                                                                                   tttGCACTTTGGTGGCGAAATATTCGGAATCGTAGGCGTCGTTAATCAACTTCATTTCATCGCGCAGCCACTGAATGGATGCGCCTGCCATAAACACCGCACCTTCCAACGCATAGTTCACTTCGCCAGTCGGGCCGCAGGCGAtggtgg                                                                    <  1:941880/150‑1 (MQ=255)
                                                                                    ttGCACTTTGGTGGCGAAATATTCGGAATCGTAGGCGTCGTTAATCAACTTCATTTCATCGCGCAGCCACTGAATGGATGCGCCTGCCATAAACACCGCACCTTCCAACGCATAGTTCACTTCGCCAGTCGGGCCGCAGGCGATggtggt                                                                   <  1:562057/150‑1 (MQ=255)
                                                                                      gCACTTTGGTGGCGAAATATTCGGAATCGTAGGCGTCGTTAATCAACTTCATTTCATCGCGCAGCCACTGAATGGATGCGCCTGCCATAAACACCGCACCTTCCAACGCATAGTTCACTTCGCCAGTCGGGCCGCAGGCGATGGTGGTca                                                                 <  2:1112289/150‑1 (MQ=255)
                                                                                         cTTTGGTGGCGAAATATTCGGAATCGTAGGCGTCGTTAATCAACTTCATTTCATCGCGCAGCCACTGAATGGATGCGCGTGCCATAAACACCGCACCTTCCAACGCATAGTTCACTTCGCCAGTCGGGCCCCAGGCGATGGTGGTcagcc                                                              >  2:846078/1‑149 (MQ=255)
                                                                                         cTTTGGTGGCGAAATATTCGGAATCGTAGGCGTCGTTAATCAACTTCATTTCATCGCGCAGCCACTGAATGGATGCGCCTGCCATAAACACCGCACCTTCCAACGCATAGTTCACTTCGCCAGTCGGGCCGCAGGCGATGGTGGTcagca                                                              >  2:1553419/1‑150 (MQ=255)
                                                                                         cTTTGGTGGCGAAATATTCGGAATCGTAGGCGTCGTTAATCAACTTCATTTCATCGCGCAGCCACTGAATGGATGCGCCTGCCATAAACACCGCACCTTCCAACGCATAGTTCACTTCGCCAGTCGGGCCGCAGGCGATGGTGGTcagca                                                              >  1:642974/1‑150 (MQ=255)
                                                                                         cTTTGGTGGCGAAATATTCGGAATCGTAGGCGTCGTTAATCAACTTCATTTCATCGCGCAGCCACTGAATGGATGCGCCTGCCATAAACACCGCACCTTCCAACGCATAGTTCACTTCGCCAGTCGGGCCGCAGGCGATGGTGGTcagca                                                              >  2:1239243/1‑150 (MQ=255)
                                                                                                  cGAAATATTCGGAATCGTAGGCGTCGTTAATCAACTTCATTTCATCGCGCAGCCACTGAATGGATGCGCCTGCCATAAACACCGCACCTTCCAACGCATAGTTCACTTCGCCAGTCGGGCCGCAGGCGATGGTGGTCAGCAGGCCGtttt                                                     >  1:882367/1‑150 (MQ=255)
                                                                                                  cGAAATATTCGGAATCGTAGGCGTCGTTAATCAACTTCATTTCATCGCGCAGCCACTGAATGGATGCGCCTGCCATAAACACCGCACCTTCCAACGCATAGTTCACTTCGCCAGTCGGGCCGCAGGCGATGGTGGTCAGCAGGCCGtttt                                                     <  1:84054/150‑1 (MQ=255)
                                                                                                    aaaTATTCGGAATCGTAGGCGTCGTTAATCAACTTCATTTCATCGCGCAGCCACTGAATGGATGCGCCTGCCATAAACACCGCACCTTCCAACGCATAGTTCACTTCGCCAGTCGGGCCGCAGGCGATGGTGGTCAGCAGGCCGTTTTCt                                                   <  2:975565/150‑1 (MQ=255)
                                                                                                    aaaTATTCGGAATCGTAGGCGTCGTTAATCAACTTCATTTCATCGCGCAGCCACTGAATGGATGCGCCTGCCATAAACACCGCACCTTCCAACGCATAGTTCACTTCGCCAGTCGGGCCGCAGGCGATGGTGGTCAGCAGGCCGTTTTCt                                                   <  1:1179528/150‑1 (MQ=255)
                                                                                                    aaaTATCCGGAATCGTAGGCGTCGTTAATCAACTTCATTTCATCGCGCAGCCACTGAATGGATGCGCCTGCCATAAACACCGCACCTTCCAACGCATAGTTCACTTCGCCAGTCGGGCCGCAGGCGATGGTGGTCAGCAGGCCGTTTTCt                                                   <  2:608005/150‑1 (MQ=255)
                                                                                                     aaTATTCGGAATCGTAGGCGTCGTTAATCAACTTCATTTCATCGCGCAGCCACTGAATGGATGCGCCTGCCATAAACACCGCACCTTCCAACGCATAGTTCACTTCGCCAGTCGGGCCGCAGGCGATGGTGGTCAGCAGGCCGTTTTCTg                                                  <  2:1059592/150‑1 (MQ=255)
                                                                                                       tatTCGGATTCGTAGGCGTGGTTAATCATCTTCATTTCATCGCGCAGCCACTGAATGGATGCGCCTGCCATAAACACCGCACCTTCCTACGCATAGTTCACTTCGCCAGTCGGGCCGCAGGCGATGGTGGTCAGCAGGCCGTTTTCTGAt                                                <  2:1454098/150‑1 (MQ=255)
                                                                                                       tatTCGGAATCGTAGGCGTCGTTAATCAACTTCATTTCATCGCGCAGCCACTGAATGGATGCGCCTGCCATAAACACCGCACCTTCCAACGCATAGTTCACTTCGCCAGTCGGGCCGCAGGCGATGGTGGTCAGCAGGCCGTTTTCTGAt                                                >  2:1445729/1‑150 (MQ=255)
                                                                                                         ttCGGAATCGTAGGCGTCGTTAATCAACTTCATTTCATCGCGCAGCCACTGAATGGATGCGCCTGCCATAAACACCGCACCTTCCAACGCATAGTTCACTTCGCCAGTCGGGCCGCAGGCGATGGTGGTCAGCAGGCCGTTTTCTGAttt                                              >  2:1436391/1‑150 (MQ=255)
                                                                                                          tCGGAATCGTAGGCGTCGTTAATCAACTTCATTTCATCGCGCAGCCACTGAATGGATGCGCCTGCCATAAACACCGCACCTTCCAACGCATAGTTCACTTCGCCAGTCGGGCCGCAGGCGATGGTGGTCAGCAGGCCGTTTTCTGATTTc                                             >  2:132258/1‑150 (MQ=255)
                                                                                                           cGGAATCGTAGGCGTCGTTAATCAACTTCATTTCCTCGCGCAGCCACTGAAAGGATGCGCCTGCCATAACCACCGCACCTTCCAACGCATAGTTCACTTCGCCAGTCGGGCCGCAGGCGATGGTGGTCAGCAGGGCCATTTCTGATTTCa                                            >  2:378933/1‑150 (MQ=255)
                                                                                                             ccaTCGTAGGCGTCGTTAATCAACTTCATTTCATCGCGCAGCCACTGAATGGATGCGCCTGCCATAAACACCGCACCTTCCAACGCATAGTTCACTTCGCCAGTCGGGCCGCAGGCGATGGTGGTCAGCAGGCCGTTTTCTGATTTCAcc                                          >  1:1123833/3‑150 (MQ=255)
                                                                                                              aaTCGTAGGCGTCGTTAATCAACTTCATTTCATCGCGCAGCCACTGAATGGATGCGCCTGACATAAACACCGCACCTTCCAACGCATAGTTCACTTCGCCAGTCGGGCCGCAGGCGATGGTGGTCAGCAGGCCGTTTTCTGATTTCACCg                                         >  1:357547/1‑150 (MQ=255)
                                                                                                               aTCGTAGGCGTCGTTAATCAACTTCATTTCATCGCGCAGCCACTGAATGGATGCGCCTGCCATAAACACCGCACCTTCCAACGCATAGTTCACTTCGCCAGTCGGGCCGCAGGCGATGGTGGTCAGCAGGCCGTTTTCTGATTTCACCGc                                        <  1:1514871/150‑1 (MQ=255)
                                                                                                                 cGTAGGCGTCGTTAATCAACTTCATTTCATCGCGCAGCCACTGAATGGATGCGCCTGTCATAAACACCGCACCTTCCAACGCATAGTTCACTTCGCCAGTCGGGCCGCAGGCGATGGTGGTCAGCAGGCCGTTTTCTGATTTCACCGCtt                                      <  2:357547/150‑1 (MQ=255)
                                                                                                                   tAGGCGTCGTTAATCAACTTCATTTCATCGCGCAGCCACTGAATGGATGCGCCTGCCAAAAACACCGCACCTTCCAACGCATAGTTCACTTCGCCAGTCGGGCCGCAGGCGATGGTGGTCAGCAGGCCGTTTTCTGATTTCACCGCTTtc                                    <  1:1137928/150‑1 (MQ=255)
                                                                                                                    aGGCGTCGTTAATCAACTTCATTTCATCGCGCAGCCACTGAATGGATGCGCCTGCCATAAACACCGCACCTTCCAACGCATAGTTCACTTCGCCAGTCGGGCCGCAGGCGATGGTGGTCAGCAGGCCGTTTTCTGATTTCACCGCTTtct                                   <  2:77229/150‑1 (MQ=255)
                                                                                                                     ggCGTCGTTAATCAACTTCATTTCATCGCGCAGCCACTGAATGGATGCGCCTGCCATAAACACCGCACCTTCCAACGCATAGTTCACTTCGCCAGTCGGGCCGCAGGCGATGGTGGTCAGCAGGCCGTTTTTTGATTTCACCGCTTtctc                                  <  2:245966/150‑1 (MQ=255)
                                                                                                                        gtcgtTAATCAACTTCATTTCATCGCGCAGCCACTGAATGGATGCGCCTGCCATAAACACCGCACCTTCCAACGCATAGTTCACTTCGCCAGTCGGGCCGCAGGCGATGGTGGTCAGCAGGCCGTTTTCTGATTTCACCGCTTTCTCGcc                               <  1:529571/150‑1 (MQ=255)
                                                                                                                        gtcgtTAATCAACTTCATTTCATCGCGCAGCCACTGAATGGATGCGCCTGCCATAAACACCGCACCTTCCAACGCATAGTTCACTTCGCCAGTCGGGCCGCAGGCGATGGTGGTCAGCAGGCCGTTTTCTGATTTCACCGCTTTCTCGcc                               <  2:1326655/150‑1 (MQ=255)
                                                                                                                          cgtTAATCAACTTCATTTCATCGCGCAGCCACTGAATGGATGCGCCTGCCATAAACACCGCACCTTCCAACGCATAGTTCACTTCGCCAGTCGGGCCGCAGGCGATGGTGGTCAGCAGGCCGTTTTCTGATTTCACCGCTTTCTCGCCAg                             <  1:132258/150‑1 (MQ=255)
                                                                                                                          cgtTAATCAACTTCATTTCATCGCGCAGCCACTGAATGGATGCGCCTGCCATAAACACCGCACCTTCCAACGCATAGTTCACTTCGCCAGTCGGGCCGCAGGCGATGGTGGTCAGCAGGCCGTTTTCTGATTTCACCGCTTTCTCGCCAg                             >  1:1123905/1‑150 (MQ=255)
                                                                                                                          cgtTAATCAACTTCATTTCATCGCGCAGCCACTGAATGGATGCGCCTGCCATAAACACCGCACCTTCCAACGCATAGTTCACTTCGCCAGTCGGGCCGCAGGCGATGGTGGTCAGCAGGCCGTTTTCTGATTTCACCGCTTTCTCGCCAg                             <  1:42571/150‑1 (MQ=255)
                                                                                                                            ttAATCAACTTCATTTCATCGCGCAGCCACTGAATGGATGCGCCTGCCATAAACACCGCACCTTCCAACTCATAGTTCACTTCGCCAGTCGGGCCGCAGGCGATGGTGGTCAGCAGGCCGTTTTCTGATTTCACCGCTTTCTCGCCAgtg                           <  2:882367/150‑1 (MQ=255)
                                                                                                                             tAATCAACTTCATTTCATCGCGCAGCCACTGAATGGATGCGCCTGCCATAAACACCGCACCTTCCAACGCATAGTTCACTTCGCCAGTCGGGCCGCAGGCGATGGTGGTCAGCAGGCCGTTTTCTGATTTCACCGCTTTCTCGCCAgtgt                          >  1:731806/1‑150 (MQ=255)
                                                                                                                             tAATCAACTTCATTTCATCGCGCAGCCACTGAATGGATGCGCCTGCCATAAACACCGCACCTTCCAACGCATAGTTCACTTCGCCAGTCGGGCCGCAGGCGATGGTGGTCAGCAGGCCGTTTTCTGATTTCACCGCTTTCTCGCCAgtgt                          >  1:890393/1‑150 (MQ=255)
                                                                                                                             tAATCAACTTCATTTCATCGCGCAGCCACTGAATGGATGCGCCTGCCATAAACACCGCACCTTCCAACGCATAGTTCACTTCGCCAGTCGGGCCGCAGGCGATGGTGGTCAGCAGGCCGTTTTCTGATTTCACCGCTTTCTCGCCAgtgt                          <  1:298041/150‑1 (MQ=255)
                                                                                                                              aaTCAACTTCATTTCATCGCGCAGCCACTGAATGGATGCGCCTGCCATAAACACCGCACCTTCCAACGCATAGTTCACTTCGCCAGTCGGGCCGCAGGCGATGGTGGTCAGCAGGCCGTTTTCTGATTTCACCGCTTTCTCGCCAGTGtt                         <  2:1123833/150‑1 (MQ=255)
                                                                                                                                 cAACTTCATTTCATCGCGCAGCCACTGAATGGATGCGCCTGCCATAAACACCGCACCTTCCAACGCATAGTTCACTTCGCCAGTCGGGCCGCAGGCGATGGTGGTCAGCAGGCCGTTTTCTGATTTCACCGCTTTCTCGCCAGTGTtcat                      <  1:1239243/150‑1 (MQ=255)
                                                                                                                                    cTTCATTTCATCGCGCAGCCACTGAATGGATGCGCCTGCCATAAACACCGCACCTTCCAACGCATAGTTCACTTCGCCAGTCGGGCCGCAGGCGATGGTGGTCAGCAGGCCGTTTTCTGATTTCACCGCTTTCTCGCCAGTGTTCATCAg                   >  1:1075522/1‑150 (MQ=255)
                                                                                                                                      tcatttcatCGCGCAGCCACTGAATGGATGCGCCTGCCATAAACACCGCACCTTCCAACGCATAGTTCACTTCGCCAGTCGGGCCGCAGGCGATGGTGGTCAGCAGGCCGTTTTCTGATTTCACCGCTTTCTCGCCAGTGTTCATCAGCa                 <  1:1203211/150‑1 (MQ=255)
                                                                                                                                      tcatttcatCGCGCAGCCACTGAATGGATGCGCCTGCCATAAACACCGCACCTTCCAACGCATAGTTCACTTCGCCAGTCGGGCCGCAGGCGATGGTGGTCAGCAGGCCGTTTTCTGATTTCACCGCTTTCTCGCCAGTGTTCATCAGCa                 <  1:704548/150‑1 (MQ=255)
                                                                                                                                       catttcatCGCGCAGCCACTGAATGGATGCGCCTGCCATAAACACCGCACCTTCCAACGCATAGTTCACTTCGCCAGTCGGGCCGCAGGCGATGGTGGTCAGCAGGCCGTTTTCTGATTTCACCGCTTTCTCGCCAGTGTTCATCAGCAt                <  2:644474/150‑1 (MQ=255)
                                                                                                                                        atttcatCGCGCAGCCACTGAATGGATGCGCCTGCCATAAACACCGCACCTTCCAACGCATAGTTCACTTCGCCAGTCGGGCCGCAGGCGATGGTGGTCAGCAGGCCGTTTTCTGATTTCACCGCTTTCTCGCCAGTGTTCATCAGCATa               <  1:527057/150‑1 (MQ=255)
                                                                                                                                          ttcatCGCGCAGCCACTGAATGGATGCGCCTGCCATAAACACCGCACCTTCCAACGCATAGTTCACTTCGCCAGTCGGGCCGCAGGCGATGGTGGTCAGCAGGCCGTTTTCTGATTTCACCGCTTTCTCGCCAGTGTTCATCAGCATaaa             >  2:1313133/1‑150 (MQ=255)
                                                                                                                                           ggcaggcgcATCCATTCAGTGGCTGCGCCTGCCATAAACACCGCACCTTCCAACGCATAGTTCACTTCGCCAGTCGGGCCACAGGCGATGGTGGTCAGCAGGCCGTTTTCTGATTTCACCGCTTTCTCGCCAGTGTTCATCAGCATAAag            >  1:570266/6‑150 (MQ=255)
                                                                                                                                            catCGCGCAGCCACTGAATGGATGCGCCTGCCATAAACACCGCACCTTCCAACGCATAGTTCACTTCGCCAGTCGGGCCGCAGGCGATGGTGGTCAGCAGGCCGTTTTCTGATTTCACCGCTTTCTCGCCAGTGTTCATCAGCATAAagc           >  2:560860/1‑150 (MQ=255)
                                                                                                                                             caggcgcATCCATTCAGTGGCTGCGCCTGCCATAAACACCGCACCTTCCAACGCATAGTTCACTTCGCCAGTCGGGCCGCAGGCGATGGTGGTCAGCAGGCCGTTTTCTGATTTCACCGCTTTCTCGCCAGTGTTCATCAGCATAAagca          >  1:1463378/4‑150 (MQ=255)
                                                                                                                                             caggcgcATCCATTCAGTGGCTGCGCCTGCCATAAACACCGCACCTTCCAACGCATAGTTCACTTCGCCAGTCGGGCCGCAGGCGATGGTGGTCAGCAGGCCGTTTTCTGATTTCACCGCTTTCTCGCCAGTGTTCATCAGCATAAagca          >  1:342441/4‑150 (MQ=255)
                                                                                                                                               cgcgCAGCCACTGAATGGATGCGCCTGCCATAAACACCGCACCTTCCAACGCATAGTTCACTTCGCCAGTCGGGCCGCAGGCGATGGTGGTCAGCAGGCCGTTTTCTGATTTCACCGCTTTCTCGCCAGTGTTCATCAGCATAAagcagc        >  1:1175658/1‑150 (MQ=255)
                                                                                                                                                gcgcATCCATTCAGTGGCTGCGCCTGCCATAAACACCGCACCTTCCAACGCATAGTTCACTTCGCCAGTCGGGCCGCAGGCGATGGTGGTCAGCAGGCCGTTTTCTGATTTCACCGCTTTCTCGCCAGTGTTCATCAGCATAAAGCAGcc       >  2:426986/1‑150 (MQ=255)
                                                                                                                                                     gCCACTGAATGGATGCGCCTGCCATAAACACCGCACCTTCCAACGCATAGTTCACTTCGCCAGTCGGGCCGCAGGCGATGGTGGTCAGCAGGCCGTTTTCTGATTTCACCGCTTTCTCGCCAGTGTTCATCAGCATAAAGCAGCCAGTGc  <  1:560860/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GAAAATCGCCCCGCGCGCATACGGGTCCCAGTACGGCGCACCCAGCCCGGTAAATGCCGGAACCACATACACACCATTGGTGTTTTGCACTTTGGTGGCGAAATATTCGGAATCGTAGGCGTCGTTAATCAACTTCATTTCATCGCGCAGCCACTGAATGGATGCGCCTGCCATAAACACCGCACCTTCCAACGCATAGTTCACTTCGCCAGTCGGGCCGCAGGCGATGGTGGTCAGCAGGCCGTTTTCTGATTTCACCGCTTTCTCGCCAGTGTTCATCAGCATAAAGCAGCCAGTGC  >  NC_000913/4116125‑4116423

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