Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,074,024 G→A 7.4% A96T (GCC→ACC)  yihV → 6‑deoxy‑6‑sulfofructose kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,074,0240GA 7.4% 203.9 / 4.3 67A96T (GCC→ACC) yihV6‑deoxy‑6‑sulfofructose kinase
Reads supporting (aligned to +/- strand):  ref base G (30/32);  new base A (2/3);  total (32/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGCGGCAAGGCTGGGGGCGCAGGTCGATTTTATTGGTCGCGTAGGTGATGACGACACGGGCAACAGCCTGCTGGCAGAGCTGGAATCCTGGGGCGTTAACACCCGTTACACCAAACGGTATAACCAGGCGAAATCTTCGCAATCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCAGCCCGGATCTGCTGCCTGACGCAGAGTGGTTGGAGGAAATTGATTTCTCTCAGTGGGATGTTGTGCTGGCAGATGTACGCTGGCACGACGGCGCTA  >  NC_000913/4073879‑4074168
                                                                                                                                                 |                                                                                                                                                
ggcggcAAGGCTGGGGGCGCAGGTCGATTTTATTGGTCGCGTAGGTGATGACGACACGGGCAACAGCCTGCTGGCAGAGCTGGAATCCTGGGGCGTTAACACCCGTTACACCAAACGGTATAACCAGGCGAAATCTTCGCAATCCGCcat                                                                                                                                              <  2:98703/150‑1 (MQ=255)
ggcggcAAGGCTGGGGGCGCAGGTCGATTTTATTGGTCGCGTAGGTGATGACGACACGGGCAACAGCCTGCTGGCAGAGCTGGAATCCTGGGGCGTTAACACCCGTTACACCAAACGGTATAACCAGGCGAAATCTTCGCAATCCGCcat                                                                                                                                              >  2:424010/1‑150 (MQ=255)
  cggcAAGGCTGGGGGCGCAGGTCGATTTTATTGGTCGCGTAGGTGATGACGACACGGGCAACAGCCTGCTGGCAGAGCTGGAATCCTGGGGCGTTAACACCCGTTACACCAAACGGTATAACCAGGCGAAATCTTCGCAATCCGCcatca                                                                                                                                            >  2:1266654/1‑150 (MQ=255)
   ggcAAGGCTGGGGGCGCAGGTCGATTTTATTGGTCGCGTAGGTGATGACGCCACGGGCAACAGCCTGCTGGCAGAGCTGGAATCCTGGGGCGTTAACACCCGTTACACCAAACGGTATAACCAGGCGAAAACTTCGCAATCCGCcatcat                                                                                                                                           >  2:495815/1‑150 (MQ=255)
      aaGGCTGGGGGCGCAGGTCGATTTTATTGGTCGCGTAGGTGATGACGACACGGGCAACAGCCTGCTGGCAGAGCTGGAATCCTGGGGCGTTAACACCCGTTACACCAAACGGTATAACCAGGCGAAATCTTCGCAATCCGCCATCAtggt                                                                                                                                        <  2:1085614/150‑1 (MQ=255)
       aGGCTGGGGGCGCAGGTCGATTTTATTGGTCGCGTAGGTGATGACGACACGGGCAACAGCCTGCTGGCAGAGCTGGAATCCTGGGGCGTTAACACCCGTTACACCAAACGGTATAACCAGGCGAAATCTTCGCAATCCACCATGAtggcg                                                                                                                                       <  2:1518792/150‑3 (MQ=255)
        ggCTGGGGGCGCAGGTCGATTTTATTGGTCGCGTAGGTGATGACGACACGGGCAACAGCCTGCTGGCAGAGCTGGAATCCTGGGGCGTTAACACCCGTTACACCAAACGGTATAACCAGGCGAAATCTTCGCAATCCGCCATCAtggtgg                                                                                                                                      >  1:114361/1‑150 (MQ=255)
        ggCGGGGGGCGCAGGTCGATTTTATTGGTCGCGTAGGGGATGACGACACGGGCAACAGCCTGCTGGCAGAGCTGGAATCCTGGGGCGTTAACACCCGTTACACCAAACGGTATAACCAGGCGAAATCTTCGCAATCCGCCATCAtggtgg                                                                                                                                      <  1:1519899/150‑1 (MQ=255)
          cTGGGGGCGCAGGTCGATTTTATTGGTCGCGTAGGTGATGACGACACGGGCAACAGCCTGCTGGCAGAGCTGGAATCCTGGGGCGTTAACACCCGTTACACCAAACGGTATAACCAGGCGAAATCTTCGCAATCCGCCATCATGGTGGAt                                                                                                                                    <  2:1510597/150‑1 (MQ=255)
            gggggCGCAGGTCGATTTTATTGGTCGCGTAGGTGATGACGACACGGGCAACAGCCTGCTGGCAGAGCTGGAATCCTGGGGCGTTAACACCCGTTACACCAAACGGTATAACCAGGCGAAATCTTCGCAATCCACCATGATGGCGGAttg                                                                                                                                  <  2:6666/150‑3 (MQ=255)
             ggggCGCAGGTCGATTTTATTGGTCGCGTAGGTGATGACGACACGGGCAACAGCCTGCTGGCAGAGCTGGAATCCTGGGGCGTTAACACCCGTTACACCAAACGGTATAACCAGGCGAAATCTTCGCAATCCACCATGATGGCGGAttgc                                                                                                                                 <  1:356958/150‑4 (MQ=255)
               ggCGCAGGTCGATTTTATTGGTCGCATAGGTGATGACGACACGGGCAACAGCCTGCTGGCAGAGCTGGAATCCTGGGGCGTTAACACCCGTTACACCAAACGGTATAACCAGGCGAAATCTTCGCAATCCGCCATCATGGTGGATAc                                                                                                                                  <  2:974310/147‑1 (MQ=255)
               ggCGCAGGTCGATTTTATTGGTCGCATAGGTGATGACGACACGGGCAACAGCCTGCTGGCAGAGCTGGAATCCTGGGGCGTTAACACCCGTTACACCAAACGGTATAACCAGGCGAAATCTTCGCAATCCGCCATCATGGTGGATAc                                                                                                                                  >  1:974310/1‑147 (MQ=255)
                 cgcAGGTCGATTTTATTGGTCGCGTAGGTGATGACGACACGGGCAACAGCCTGCTGGCAGAGCTGGAATCCTGGGGCGTTAACACCCGTTACACCAAACGGTATAACCAGGCGAAATCTTCGCAATCCGCCATCATGGTGGATACCAAAg                                                                                                                             >  1:976644/1‑150 (MQ=255)
                    aGGTCGATTTTATTGGTCGCGTAGGTGATGACGACACGGGCAACAGCCTGCTGGCAGAGCTGGAATCCTGGGGCGTTAACACCCGTTACACCAAACGGTATAACGAGGCGAAATCTTCGCAATCCGCCATCATGGTGGATACCAAAGGCg                                                                                                                          <  2:1393917/150‑1 (MQ=255)
                        cGATTTTATTGGTCGCGTAGGTGATGACGAGACGGGCAACAGCCTGCTGGCAGAGCTGGAATCCTGGGGCGTTAACACCCGTTACATCAAACGGTATAACCAGGCGAAATCTTCGCAATCCGCCATCATGGTGGATACCAAAGGCGAGCg                                                                                                                      <  2:1383777/150‑1 (MQ=255)
                        cGATTTTATTGGTCGCGTAGGTGATGACGACACGGGCAACAGCCTGCTGGCAGAGCTGGAATCCTGGGGCGTTAACACCCGTTACACCAAACGGTATAACCAGGCGAAATCTTCGCAATCCGCCATCATGGTGGATACCAAAGGCGAGCg                                                                                                                      <  2:550120/150‑1 (MQ=255)
                        cGATTTTATTGGTCGCGTAGGTGATGACGACACGGGCAACAGCCTGCTGGCAGAGCTGGAATCCTGGGGCGTTAACACCCGTTACACCAAACGGTATAACCAGGCGAAATCTTCGCAATCCGCCATCATGGTGGATACCAAAGGCGAGCg                                                                                                                      <  1:495070/150‑1 (MQ=255)
                         gATTTTATTGGTCGCGTAGGTGATGACGACACGGGCAACAGCCTGCTGGCAGAGCTGGAATCCTGGGGCGTTAACACCCGTTACACCAAACGGTATAAACAGGCGAAATCTTCGCAATCCGCCATCATGGTGGATACCAAAGGCGAGCgg                                                                                                                     >  2:163290/1‑150 (MQ=255)
                           ttttATTGGTCGCGTAGGTGATGACGACACGGGCAACAGCCTGCTGGCAGAGCTGGAATCCTGGGGCGTTAACACCCGTTACACCAAACGGTATAACCAGGCGAAATCTTCGCAATCCGCCATCATGGTGGATACCAAAGGCGAGCGGAt                                                                                                                   <  2:494454/150‑1 (MQ=255)
                              tATTGGTCGCGTAGGTGATGACGACACGGGCAACAGCCTGCTGGCAGAGCTGGAATCCTGGGGCGTTAACACCCGTTACACCAAACGGTATAACCAGGCGAAATCTTCGCAATCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAAt                                                                                                                <  2:1590188/150‑1 (MQ=255)
                                  ggTCGCGTAGGTGATGACGACACGGGCAACAGCCTGCTGGCAGAGCTGGAATCCTGGGGCGTTAACACCCGTTACACCAAACGGTATAACCAGGCGAAATCTTCGCAATCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCAtt                                                                                                            <  2:1503396/150‑1 (MQ=255)
                                   gTCGCGTAGGTGATGACGACACGGGCAACAGCCTGCCGGCAGAGCTAGAATCCTGGGGCGTTAACACCCGTTACACCAAACGGTATAACCAGGCGAAATCTTCGCAATCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTa                                                                                                           <  2:1587089/150‑1 (MQ=255)
                                           ggTGATGACGACACGGGCAACAGCCTGCTGGCAGAGCTGGAATCCTGGGGCGTTAACACCCGTTACACCAAACGGTATAACCAGGCGAAATCTTCGCAATCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTAACTAcccc                                                                                                   >  2:627092/1‑150 (MQ=255)
                                               atgaCGACACGGGCAACAGCCTGCTGGCAGAGCTGGAATCCTGGGGCGTTAACACCCGTTACACCAAACGGTATAACCAGGCGAAATCTTCGCAATCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCAGcc                                                                                               <  1:644849/150‑1 (MQ=255)
                                                 gacgacACGGGCAACAGCCTGCTGGCAGAGCTGGAATCCTGGGGCGTTAACACCCGTTACACCAAACGGTATAACCAGGCGAAATCTTCGCAATCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCAGCCCg                                                                                             >  1:1573813/1‑150 (MQ=255)
                                                      cacGGGCAACAGCCTGCTGGCAGAGCTGGAATCCTGGGGCGTTAACACCCGTTACACCAAACGGTATAACCAGGCGAAATCTTCGCAATCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCAGCCCGGATct                                                                                        >  1:1138932/1‑150 (MQ=255)
                                                      cacGGGCAACAGCCTGCTGGCAGAGCTGGAATCCTGGGGCGTTAACACCCGTTACACCAAACGGTATAACCAGGCGAAATCTTCGCAATCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCAGCCCGGATct                                                                                        <  2:1623952/150‑1 (MQ=255)
                                                           gCAACAGCCTGCTGGCAGAGCTGGAATCCTGGGGCGTTAACACCCGTTACACCAAACGGTATAACCAGGCGAAATCTTCGCAATCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCAGCCCGGATCtgctgc                                                                                   >  2:1254970/1‑150 (MQ=255)
                                                                   ctgctgGCAGAGCTGGAATCCTGGGGCGTTAACACCCGTTACACCAAACGGTATAACCAGGCGAAATCTTCGCAATCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCAGCCCGGATCTGCTGCCTGACGCa                                                                           <  2:1414915/150‑1 (MQ=255)
                                                                        ggCAGAGCTGGAATCCTGGGGCGTTAACACCCGTTACACCAAACGGTATAACCAGGCGAAATCTTCGCAATCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCAGCCCGGATCTGCTGCCTGACGCAGAGTg                                                                      >  2:1133256/1‑150 (MQ=255)
                                                                         gCAGAGCTGGAATCCTGGGGCGTTAACACCCGTTACACCAAACGGTATAACCAGGCGAAATCTTCGCAATCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCAGCCCGGATCTGCTGCCTGACGCAGAGTgg                                                                     >  2:78682/1‑150 (MQ=255)
                                                                             agCTGGAATCCTGGGGCGTTAACACCCGTTACACCAAACGGTATAACCAGGCGAAATCTTCGCAATCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCAGCCCGGATCTGCTGCCTGCCGCAGAGTGGTTgg                                                                 >  2:1433631/1‑150 (MQ=255)
                                                                               cTGGAATCCTGGGGCGTTAACACCCGTTACACCAAACGGTATAACCAGGCGAAATCTTCGCAATCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCAGCCCGGATCTGCTGCCTGACGCAGAGTGGTTggag                                                               <  1:929244/150‑1 (MQ=255)
                                                                                       cTGGGGCGTTAACACCCGTTACACCAAACGGTATAACCAGGCGAAATCTTCGCAATCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCAGCCCGGATCTGCTGCCTGACGCAGAGTGGTTGGAGGAAATTGa                                                       <  2:347239/150‑1 (MQ=255)
                                                                                        tGGGGCGTTAACACCCGTTACACCAAACGGTATAACCAGGCGAAATCTTCGCAATCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCAGCCCGGATCTGCTGCCTGACGCAGAGTGGTTGGAGGAAATTGAt                                                      <  1:424010/150‑1 (MQ=255)
                                                                                         ggggCGTTAACACCCGTTACACCAAACGGTATAACCAGGCGAAATCTTCGCAATCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCAGCCCGGATCTGCTGCCTGACGCAGAGTGGTTGGAGGAAATTGAtt                                                     >  1:101473/1‑150 (MQ=255)
                                                                                             cGTTAACACCCGTTACACCAAACGGTATAACCAGGCGAAATCTTCGCAATCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCAGCCCGGATCTGCTGCCTGACGCAGAGTGGTTGGAGGAAATTGATTtctc                                                 >  1:1034949/1‑150 (MQ=255)
                                                                                                 aaCACCCGTTACACCAAACGGTATAACCAGGCGAAATCTTCGCAATCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTa                                                                                                           <  2:1444264/88‑1 (MQ=255)
                                                                                                 aaCACCCGTTACACCAAACGGTATAACCAGGCGAAATCTTCGCAATCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTa                                                                                                           >  1:1444264/1‑88 (MQ=255)
                                                                                                 aaCACCCGTTACACCAAACGGTATAACCAGGCGAAATCTTCGCAATCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCCGCCCGGATCTGCTGCCTGACGCAGAGTGGTTGGAGGAAATTGATTTCTCTCAg                                             >  2:968227/1‑150 (MQ=255)
                                                                                                 aaCACCCGTTACACCAAACGGTATAACCAGGCGAAATCTTCGCAATCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCAGCCCGGATCTGCTGCCTGACGCAGAGTGGTTGGAGGAAATTGATTTCTCTCAg                                             >  2:412076/1‑150 (MQ=255)
                                                                                                  acacCCGTTACACCAAACGGTATAACCAGGCGAAATCTTCGCAATCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCAGCCCGGATCTGCTGCCTGACGCAGAGTGGTTGGAGGAAATTGATTTCTCTCAGt                                            <  2:1371986/150‑1 (MQ=255)
                                                                                                   cacCCGTTACACCAAACGGTATAACCAGGCGAAATCTTCGCAATCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCAGCCCGGATCTGCTGCCTGACGCAGAGTGGTTGGAGGAAATTGATTTCTCTCAGTg                                           <  2:475903/150‑1 (MQ=255)
                                                                                                   cacCCGTTACACCAAACGGTATAACCAGGCGAAATCTTCGCAATCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCAGCCCGGATCTGCTGCCTGACGCAGAGTGGTTGGAGGAAATTGATTTCTCTCAGTg                                           >  2:1409089/1‑150 (MQ=255)
                                                                                                          tACACCAAACGGTATAACCAGGCGAAATCTTCGCAATCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCAGCCCGGATCTGCTGCCTGACGCAGAGTGGTTGGAGGAAATTGATTTCTCTCAGTGGGAtgtt                                    >  1:947857/1‑150 (MQ=255)
                                                                                                           acacCAAACGGTATAACCAGGCGAAATCTTCGCAATCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCAGCCCGGATCTGCTGCCTGACGCAGAGTGGTTGGAGGAAATTGATTTCTCTCAGTGGGAtgttg                                   >  2:433318/1‑150 (MQ=255)
                                                                                                               cAAACGGTATAACCAGGCGAAATCTTCGCAATCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCAGCCCGGATCTGCTGCCTGACGCAGAGTGGTTGGAGGAAATTGATTTCTCTCAGTGGGATGTTGTGCt                               <  2:964659/150‑1 (MQ=255)
                                                                                                               cAAACGGTATAACCAGGCGAAATCTTCGCAATCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCAGCCCGGATCTGCTGCCTGACGCAGAGTGGTTGGAGGAAATTGATTTCTCTCAGTGGGATGTTGTGCt                               <  2:1189445/150‑1 (MQ=255)
                                                                                                                 aaCGGTATAACCAGGCGAAATCTTCGCAATCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCAGCCCGGATCTGCTGCCTGACGCAGAGTGGTTGGAGGAAATTGATTTCTCTCAGTGGGATGTTGTGCTgg                             >  2:649486/1‑150 (MQ=255)
                                                                                                                 aaCGGTATAACCAGGCGAAATCTTCGCAATCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCAGCCCGGATCTGCTGCCTGACGCAGAGTGGTTGGAGGAAATTGATTTCTCTCAGTGGGATGTTGTGCTgg                             <  2:520006/150‑1 (MQ=255)
                                                                                                                  aCGGTATAACCAGGCGAAATCTTCGCAATCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCAGCCCGGATCTGCTGCCTGACGCAGAGTGGTTGGAGGAAATTGATTTCTCTCAGTGGGATGTTGTGCTGGc                            <  2:1316119/150‑1 (MQ=255)
                                                                                                                  aCGGTATAACCAGGCGAAATCTTCGCAATCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCAGCCCGGATCTGCTGCCTGACGCAGAGTGGTTGGAGGAAATTGATTTCTCTCAGTGGGATGTTGTGCTGGc                            >  2:1436268/1‑150 (MQ=255)
                                                                                                                  aCGGTATAACCAGGCGAAATCTTCGCAATCCGCCATCATGGTGGAGACCAAAGGCGAGCGGATAATCATTAACTACCCCAGCCCGGATCTGCTGCCTGACGCAGAGTGGTTGGAGGAAATTGATTTCTCTCAGTGGGATGTTGTGCTGGc                            <  1:751956/150‑1 (MQ=255)
                                                                                                                   cGGAATAACCAGGCGAAATCTTCGCAATCCGCCATCATGCTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCAGCCCGGATCTGCTGCCTGACGCAGAGTGGTTGGAGGAAATTGATTTCTCTCAGTGGGATGTTGTGCTGGCa                           >  2:1193195/1‑150 (MQ=255)
                                                                                                                            cAGGCGAAATCTTCGCAACCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCAGCCCGGATCTGCTGCCTGACGCAGAGTGGTTGGAGGAAATTGATTTCTCTCAGTGGGATGTTGTGCTGGCAGATGTACGc                  >  1:561343/1‑150 (MQ=255)
                                                                                                                              ggCGAAATCTTCGCAATCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCAGCCCGGATCTGCTGCCTGACGCAGAGTGGTTGGAGGAAATTGATTTCTCTCAGTGGGATGTTGTGCTGGCAGATGTACGCTg                <  1:1160661/150‑1 (MQ=255)
                                                                                                                              ggCGAAATCTTCGCAATCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCAGCCCGGATCTGCTGCCTGACGCAGAGTGGTTGGAGGAAATTGATTTCTCTCAGTGGGATGTTGTGCTGGCAGATGTACGCTg                >  1:811918/1‑150 (MQ=255)
                                                                                                                               gCGAAATCTTCGCAATCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCAGCCCGGATCTGCTGCCTGACGCAGAGTGGTTGGAGGAAATTGATTTCTCTCAGTGGGATGTTGTGCTGGCAGATGTACGCTgg               <  1:1266654/150‑1 (MQ=255)
                                                                                                                               gCGAAATCTTCGCAATCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCAGCCCGGATCTGCTGCCTGACGCAGAGTGGTTGGAGGAAATTGATTTCTCTCAGTGGGATGTTGTGCTGGCAGATGTACGCTgg               >  2:613634/1‑150 (MQ=255)
                                                                                                                                cGAAATCTTCGCAATCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCAGCCCGGATCTGCTGCCTGACGCAGAGTGGTTGGAGGAAATTGATTTCTCTCAGTGGGATGTTGTGCTGGCAGATGTACGCTgcc              >  2:946365/1‑148 (MQ=255)
                                                                                                                                      cTTCGCAATCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCAGCCCGGATCTGCTGCCTGACGCAGAGTGGTTGGAGGAAATTGATTTCTCTCAGTGGGATGTTGTGCTGGCAGATGTACGCTGGCacgacg        <  1:1361931/150‑1 (MQ=255)
                                                                                                                                          gCAATCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCAGCCCGGATCTGCTGCCTGACGCAGAGTGGTTGGAGGAAATTGATTTCTCTCAGTGGGATGTTGTGCTGGCa                           <  1:1487069/127‑1 (MQ=255)
                                                                                                                                          gCAATCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCAGCCCGGATCTGCTGCCTGACGCAGAGTGGTTGGAGGAAATTGATTTCTCTCAGTGGGATGTTGTGCTGGCa                           >  2:1487069/1‑127 (MQ=255)
                                                                                                                                           gtaTCCACCATGATGGCGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCAGCCCGGATCTGCTGCCTGACGCAGAGTGGTTGGAGGAAATTGATTTCTCTCAGTGGGATGTTGTGCTGGCAGATGTACGCTGGCACGACGGCGCt   >  2:1004132/3‑150 (MQ=255)
                                                                                                                                           gtaTCCACCATGATGGCGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCAGCCCGGATCTGCTGCCTGACGCAGAGTGGTTGGAGGAAATTGATTTCTCTCAGTGGGATGTTGTGCTGGCAGATGTACGCTGGCACGACGGCGCt   >  1:1496030/3‑150 (MQ=255)
                                                                                                                                            aaTCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCAGCCCGGATCTGCTGCCTGACGCAGAGTGGTTGGAGGAAATTGATTTCTCTCAGTGGGATGTTGTGCTGGCAGATGTACGCTGGCACGACGGCGCTa  <  2:1138932/150‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
GGCGGCAAGGCTGGGGGCGCAGGTCGATTTTATTGGTCGCGTAGGTGATGACGACACGGGCAACAGCCTGCTGGCAGAGCTGGAATCCTGGGGCGTTAACACCCGTTACACCAAACGGTATAACCAGGCGAAATCTTCGCAATCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCAGCCCGGATCTGCTGCCTGACGCAGAGTGGTTGGAGGAAATTGATTTCTCTCAGTGGGATGTTGTGCTGGCAGATGTACGCTGGCACGACGGCGCTA  >  NC_000913/4073879‑4074168

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