Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,503,930 A→C 7.6% S186A (TCG→GCG)  fryA ← putative PTS enzyme: Hpr, enzyme I and IIA components

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,503,9300AC 7.6% 201.2 / 3.9 74S186A (TCG→GCG) fryAputative PTS enzyme: Hpr, enzyme I and IIA components
Reads supporting (aligned to +/- strand):  ref base A (41/27);  new base C (3/3);  total (44/30)
Fisher's exact test for biased strand distribution p-value = 6.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGCGCAGATAATCGCTGGCAGAGGCAGAAAGTTTGGCGCAAACCTGCTCCATATTGCTGATGATCGCCGCCCCCAGCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCTGAATCAGCGACAAATGGGCGCTGAGGATAGTTTTGCTTTCGCCGTCACGCTCACGCAATTGCTGGTTCAGTTGCTCGGCAAGCGTTGCCAGGCTGTGCTCCAGCCGGGTGGAATCTTGCGCACTGGCGGGGATTGCCCGATAACTGTCGAGGCTGTCG  >  NC_000913/2503789‑2504078
                                                                                                                                             |                                                                                                                                                    
cgcgCAGATAATCGCTGGCAGAGGCAGAAAGTTTGGCGCAAACCTGCTCCATATTGCTGATGATCGCCGCCCCCAGCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCTGAATCAGCGACAAATggg                                                                                                                                              <  1:850939/150‑1 (MQ=255)
     aGATAATCGCTGGCAGAGGCAGAAAGTTTGGCGCAAACCTGCTCCATATTGCTGATGATCGCCGCCCCCAGCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCTGAATCAGCGACAAATGGGCGCTg                                                                                                                                         <  2:628258/150‑1 (MQ=255)
      gATAATCGCAGGCAGAGGCAGAAAGTATGGCGCAAACCTGCTCCATATTGCTGATGATCGCCGCCCCCAGCCCATGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCATGAATCAGCTACAAATGGGCGATGa                                                                                                                                        >  2:1348545/1‑150 (MQ=255)
       ttAATCGCTGGCAGAGGCAGAAAGTTTGGCGCAAACCTGCTCTATATTGCTGATGTTCGCCGCCCCCATCCCCTGATGCTGTTCTGTCATCAGGCTACGGATATTGCCTGCAAATTCATCATCCTGAATCAGCGACAAATGGGCGCTGAg                                                                                                                                       <  1:1092609/149‑1 (MQ=255)
         aaTCGCTGGCAGAGGCAGAAAGTTTGGCGCAAACCTGCTCCATATTGCTGATGATCGCCGCCCCCAGCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCTGAATCAGCGACAAATGGGCGCTGAGGa                                                                                                                                     <  1:1396225/150‑1 (MQ=255)
          aTCGCTGGCAGAGGCAGAAAGTTTGGCGCAAACCTGCTCCATATTGCTGATGATCGCCGCCCCCAGCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCTGAATCAGCGCCCATTTGTCGCTGattca                                                                                                                                    <  1:1510677/150‑5 (MQ=255)
          aTCGCTGGCAGAGGCAGAAAGTTTGGCGCAAACCTGCTCCATATTGCTGATGATCGCCGCCCCCAGCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCTGAATCAGCGCCCATTTGTCGCTGattaa                                                                                                                                    <  2:1504220/150‑5 (MQ=255)
           tCGCTGGCAGAGGCAGAAAGTTTGGCGCAAACCTGCTCCATATTGCTGATGATCGCCGCCCCCAGCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCTGAATCAGCGCCCATTTGTCGCTGattcag                                                                                                                                   <  1:360422/150‑6 (MQ=255)
           tCGCTGGCAGAGGCAGAAAGTTTGGCGCAAACCTGCTCCATATTGCTGATGATCGCCGCCCCCAGCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCTGAATCAGCGACAAATGGGCGCTGAGGATa                                                                                                                                   <  1:1690340/150‑1 (MQ=255)
             gCTGGCAGAGGCAGAAAGTTTGGCGCAAACCTGCTCCATATTGCTGATGATCGCCGCCCCCAGCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCTGAATCAGCGACAAATGGGCGCTGAGGATAGt                                                                                                                                 >  2:431698/1‑150 (MQ=255)
                ggcagaggcagaAAGTTTGGCGCAAACCTGCTCCATATTGCTGATGATCGCCGCCCCCAGCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCTGAATCAGCGACAAATGGGCGCTGAGGATAGtttt                                                                                                                              <  1:736897/150‑1 (MQ=255)
                   agaggcagaAAGTTTGGCGCAAACCTGCTCCATATTGCTGATGATCGCCGCCCCCAGCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCTGAATCAGCGACAAATGGGCGCTGAGGATAGTTTTGCt                                                                                                                           >  1:551837/1‑150 (MQ=255)
                   agaggcagaAAGTTTGGCGCAAACCTGCTCCATATTGCTGATGATCGCCGCCCCCAGCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCTGAATCAGCGACAAATGGGCGCTGAGGATAGTTTTGCt                                                                                                                           >  1:1667313/1‑150 (MQ=255)
                      ggcagaAAGTTTGGCGCAAACCTGCTCCATATTGCTGATGATCGCCGCCCCCAGCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCTGAATCAGCGACAAATGGGCGCTGAGGATAGTTTTGCTTTc                                                                                                                        >  1:29734/1‑150 (MQ=255)
                       gcagaAAGTTTGGCGCAAACCTGCTCCATATTGCTGATGATCGCCGCCCCCAGCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCTGAATCAGCGACAAATGGGCGCTGAGGATAGTTTTGCTTTCg                                                                                                                       <  1:996283/150‑1 (MQ=255)
                       gcagaAAGTTTGGCGCAAACCTGCTCCATATTGCTGATGATCGCCGCCCCCAGCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCTGAATCAGCGACAAATGGGCGCTGAGGATAGTTTTGCTTTCg                                                                                                                       <  1:187316/150‑1 (MQ=255)
                                 tGGCGCAAACCTGCTCCATATTGCTGATGATCGCCGCCCCCAGCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCTGAATCAGCGACAAATGGGCGCTGAGGATAGTTTTGCTTTCGCCGtcacgct                                                                                                             >  2:823834/1‑150 (MQ=255)
                                   gcgcAAACCTGCTCCATATTGCTGATGATCGCCGCCCCCAGCCCCTGCTGCTGGTCTGTCATCCGGCGACGGATATTGCCTGCAACTTCATCCTCCTGAATCAGCGCCAAATGGGCGCGGAGGCTAGTTTTGCTTTCGCCGtcacgctcc                                                                                                           >  1:174516/1‑149 (MQ=255)
                                      cAAACCTGCTCCATATTGCTGATGATCGCCGCCCCCAGCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCTGAATCAGCGACAAATGGGCGCTGAGGATAGTTTTGCTTTCGCCGtcacgctcacgc                                                                                                        >  2:507203/1‑150 (MQ=255)
                                          ccTGCTCCATATTGCTGATGATCGCCGCCCCCAGCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCTGAATCAGCGACAAATGGGCGCTGAGGATAGTTTTGCTTTCGCCGTCACGCTCACGCAAtt                                                                                                    >  2:892689/1‑150 (MQ=255)
                                          ccTGCTCCATATTGCTGATGATCGCCGCCACCAGCCACTGAAGCTGTTCTGTCATCAGGCGCAGGATATTGCCTGCAAATTCATCCTCCTGAATCAGCGACAAATGGGCGCTGAGGATAGATTTGCTTTCGCCGTCACGCTCAAGCAAtt                                                                                                    >  1:943403/1‑150 (MQ=255)
                                            tGCTCCATATTGCTGATGATCGCCGCCCCCAGCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCTGAATCAGCGACAAATGGGCGCTGAGGATAGTTTTGCTTTCGCCGTCACGCTCACGCAATTGc                                                                                                  <  2:377742/150‑1 (MQ=255)
                                              cTCCATATTGCTGATGATCGCCGCCCCCAGCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCTGAATCAGCGACAAATGGGCGCTGAGGATAGTTTTGCTTTCGCCGTCACGCTCACGCAATTGCTg                                                                                                >  1:435123/1‑150 (MQ=255)
                                              cTCCATATTGCTGATGAGAGCCGCCCCCAGCCCCTGCTGCTGTTCTGTCATCAGGCGCCGGAGATTGCCTGCAAATTCATCATCCTGAATCAGCGACCAATGGGCGCTGAGGATAGTTTTGCTTTCGCCGTCACGCTCACGCAATTGCTg                                                                                                <  1:320252/150‑1 (MQ=255)
                                                ccATATTGCTGATGATCGCCGCCCCCAGCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCTGAATCAGCGACAAATGGGCGCTGAGGATAGTTTTGCTTTCGCCGTCACGCTCACGCAATTGCTGGt                                                                                              >  2:1518664/1‑150 (MQ=255)
                                                        cTGATGATCGCCGCCCCCAGCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCTGAATCAGCGACAAATGGGCGCTGAGGATAGTTTTGCTTTCGCCGTCACGCTCACGCAATTGCTGGTTCAGttac                                                                                      >  2:1138766/1‑148 (MQ=255)
                                                        cTGATGATCGCCGCCCCAGCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCTGAATCAGCGACAAATGGGCGCTGAGGATAGTTTTGCTTTCGCCGTCACGCTCACGCAATTGCTGGTTCAGTTGCt                                                                                     >  1:481941/1‑150 (MQ=255)
                                                         tgatgaTCGCCGCCCCCAGCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCTGAATCAGCGACAAATGGGCGCTGAGGATAGTTTTGCTTTCGCCGTCACGCTCACGCAATTGCTGGTTCAGTTGCt                                                                                     >  1:1695111/1‑150 (MQ=255)
                                                           atgatCGCCGCCCCCAGCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCTGAATCAGCGACAAATGGGCGCTGAGGATAGTTTTGCTTTCGCCGTCACGCTCACGCAATTGCTGGTTCAGTTGCTCg                                                                                   <  1:1175210/150‑1 (MQ=255)
                                                            tgatCGCCGCCCCCAGCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCTGAATCAGCGACAAATGGGCGCTGAGGATAGTTTTGCTTTCGCCGTCACGCTCACGCAATTGCTGGTTCAGTTGCTCgg                                                                                  >  1:599037/1‑150 (MQ=255)
                                                              atCGCCGCCCCCAGCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGGCTGCAAATTCATCATCCTGAATCAGCGACAAATTGGCTCTGAGGATAGTTTTGCTTTCGCAGTCACGCTCACGCAATTGCTGGTTCAGTTGCTCGGCa                                                                                >  1:545304/1‑150 (MQ=255)
                                                               tCGCCGCCCCCAGCCCCTGATGCTGTTCTGTCATCAGGCGTCGGATTTTGCCTGCAATTTAATCAGCGTCAATCAGGGCCAATTGGGGGGTGGGGATAGTTTTGCGTTCACCTTCACCGGCACGCAATTGCTGCTACAGTTGCTCGGCag                                                                               >  2:350537/1‑149 (MQ=255)
                                                                 gccgccCCCAGCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCTGAATCAGCGCCAAATGGGCGCTGAGGATAGTTTTGCTTTCGCCGTCACTCTCACGCAATTGCTGGTTCAGTTGCTCGGCAAGc                                                                             >  2:1062267/1‑150 (MQ=255)
                                                                 gccgccCCCAGCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCTGAATCAGCGACAAATGGGCGCTGAGGAGAGTTTTGCTTTCGCCGTCACGCTCACGCAATTGCTGGTTCAGTTGCTCGGCAAGc                                                                             >  2:1290170/1‑150 (MQ=255)
                                                                             cccTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCTGAATCAGCGACAAATGGGCGCTGAGGATAGTTTTGCTTTCGCCGTCACGCTCACGCAATTGCTGGTTCAGTTGCTCGGCAAGCGTTGCCAGGCtg                                                                 >  2:1141306/1‑150 (MQ=255)
                                                                             cccTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGAAAATTCATCATCCGGAATCAGCGACAAATGGGCGCTGAGGATAGTTTTGCTTTCGCCGTCACGCTCACGCAATTGCTGGTTCAGTTGCTCGGCAAGCGTTGCCAGGCtg                                                                 <  1:1348545/150‑1 (MQ=255)
                                                                               cTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCTGAATCAGCGACAAATGGGCGCTGAGGATAGTTTTGCTTTCGCCGTCACGCTCACGCAATTGCTGGTTCAGTTGCTCGGCAAGCGTTGCCAGGCtgtg                                                               >  1:1629540/1‑150 (MQ=255)
                                                                                tGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCTGAATCAGCGACAAATGGGCGCTGAGGATAGTTTTGCTTTCGCCGTCACGCTCACGCAATTGCTGGTTCAGTTGCTCGGCAAGCGTTGCCAGGCTGTGc                                                              <  2:1265258/150‑1 (MQ=255)
                                                                                  aTGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCTGAATCAGCGACAAATGGGCGCTGAGGATAGTTTTGCTTTCGCCGTCACGCTCACGCAATTGCTGGTTCAGTTGCTCGGCAAGCGTTGCCAGGCTGTGCTc                                                            >  1:827867/1‑150 (MQ=255)
                                                                                   tGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCTGAATCAGCGACAAATGGGCGCTGAGGATAGTTTTGCTTTCGCCGTCACGCTCACGCAATTGCTGGTTCAGTTGCTCGGCAAGCGTTGCCAGGCTGTGCTcc                                                           >  1:1371266/1‑150 (MQ=255)
                                                                                    gCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCTGAATCAGCGACAAATGGGCGCTGAGGATAGTTTTGCTTTCGCCGTCACGCTCACGCAATTGCTGGTTCAGTTGCTCGGCAAGCGTTGCCAGGCTGTGCTCCa                                                          <  1:1518664/150‑1 (MQ=255)
                                                                                     cTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCTGAATCAGCGACAAATGGGCGCTGAGGATAGTTTTGCTTTCGCCGTCACGCTCACGCAATTGCTGGTTCAGTTGCTCGGCAAGCGTTGCCAGGCTGTGCTCCAg                                                         >  2:276467/1‑150 (MQ=255)
                                                                                         tCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCTGAATCAGCGACAAATGGGCGCTGAGGATAGTTTTGCTTTCGCCGTCACGCTCACGCAATTGCTGGTTCAGTTACTCGGCAAGCGTTGCCAGGCTGTGCTCCAGCCgg                                                     <  1:1138766/150‑1 (MQ=255)
                                                                                          cTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCTGAATCAGCGACAAATGGGCGCTGAGGATAGTTTTGCTTTCGCCGTCACGCTCACGCAATTGCTGGTTCAGTTGCTCGGCAAGCGTTGCCAGGCTGTGCTCCAGCCggg                                                    <  2:1695111/150‑1 (MQ=255)
                                                                                           tGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCTGAATCAGCGACAAATGGGCGCTGAGGATAGTTTTGCTTTCGCCGTCACGCTCACGCAATTGCTGGTTCAGTTGCTCGGCAAGCGTTGCCAGGCTGTGCTCCAGCCGGGt                                                   >  1:1149386/1‑150 (MQ=255)
                                                                                                  aGGCGACGGATATTGCCTGCAAATTCATCATCCTGAATCAGCGACAAATGGGCGCTGAGGATAGTTTTGCTTTCGCCGTCACGCTCACGCAATTGCTGGTTCAGTTGCTCGGCAAGCGTTGCCAGGCTGTGCTCCAGCCGGGTGGAATCt                                            <  1:823834/150‑1 (MQ=255)
                                                                                                    gCGACGGATATTGCCTGCAAATTCATCATCCTGAATCAGCGACAAATGGGCGCTGAGGATAGTTTTGCTTTCGCCGTCACGCTCACGCAATTGCTGGTTCAGTTGCTCGGCAAGCGTTGCCAGGCTGTGCTCCAGCCGGGTGGAATCTTg                                          >  1:939645/1‑150 (MQ=255)
                                                                                                    gCGACGGATATTGCCTGCAAATTCATCATCCTGAATCAGCGACAAATGGGCGCTGAGGATAGTTTTGCTTTCGCCGTCACGCTCACGCAATTGCTGGTTCAGTTGCTCGGCAAGCGTTGCCAGGCTGTGCTCCAGCCGGGTGGAATCTTg                                          <  1:1253200/150‑1 (MQ=255)
                                                                                                        cGGATATTGCCTGCAAATTCATCATCCTGAATCAGCGACAAATGGGCGCTGAGGATAGTTTTGCTTTCGCCGTCACGCTCACGCAATTGCTGGTTCAGTTGCTCGGCAAGCGTTGCCAGGCTGTGCTCCAGCCGGGTGGAATCTTGCGCa                                      >  1:1350071/1‑150 (MQ=255)
                                                                                                          gATATTGCCTGCAAATTCATCATCCTGAATCAGCGACAAATGGGCGCTGAGGATAGTTTTGCTTTCGCCGTCACGCTCACGCAATTGCTGGTTCAGTTGCTCGGCAAGCGTTGCCAGGCTGTGCTCCAGCCGGGTGGAATCTTGCGCACt                                    >  1:1287464/1‑150 (MQ=255)
                                                                                                             atTGCCTGCAAATTCATCATCCTGAATCAGCGACAAATGGGCGCTGAGGATAGTTTTGCTTTCGCCGTCACGCTCACGCAATTGCTGGTTCAGTTGCTCGGCAAGCGTTGCCAGGCTGTGCTCCAGCCGGGTGGAATCTTGCGCACTGGc                                 >  2:884410/1‑150 (MQ=255)
                                                                                                                gCCTGCAAATTCATCATCCTGAATCAGCGACAAATGGGCGCTGAGGATAGTTTTGCTTTCGCCGTCACGCTCACGCAATTGCTGGTTCAGTTGCTCGGCAAGCGTTGCCAGGCTGTGCTCCAGCCGGGTGGAATCTTGCGCACTGGCggg                              >  2:651513/1‑150 (MQ=255)
                                                                                                                 ccTGCAAATTCATCATCCTGAATCAGCGACAAATGGGCGCTGAGGATAGTTTTGCTTTCGCCGTCACGCTCACGCAATTGCTGGTTCAGTTGCTCGGCAAGCGTTGCCAGGCTGTGCTCCAGCCGGGTGGAATCTTGCGCACTGGCgggg                             >  1:1205862/1‑150 (MQ=255)
                                                                                                                 ccTGCAAATTCATCATCCTGAATCAGCGACAAATGGGCGCTGAGGATAGTTTTGCTTTCGCCGTCACGCTCACGCAATTGCTGGTTCAGTTGCTCGGCAAGCGTTGCCAGGCTGTGCTCCAGCCGGGTGGAATCTTGCGCACTGGCgggg                             >  2:1240388/1‑150 (MQ=255)
                                                                                                                 ccTGCAAATTCATCATCCTGAATCAGCGACAAATGGGCGCTGAGGATAGTTTTGCTTTCGCCGTCACGCTCACGCAATTGCTGGTTCAGTTGCGCGGCAAGCGTTGCCAGGCTGTGCTCCAGCCGGGTGGAATCTTGCGCACTGGCgggg                             <  1:77110/150‑1 (MQ=255)
                                                                                                                   tGCAAATTCATCATCCTGAATCAGCGACAAATGGGCGCTGAGGATAGTTTTGCTTTCGCCGTCACGCTCACGCAATTGCTGGTTCAGTTGCTCGGCAAGCGTTGCCAGGCTGTGCTCCAGCCGGGTGGAATCTTGCGCACTGGCGGGGAt                           >  2:743448/1‑150 (MQ=255)
                                                                                                                   tGCAAATTCATCATCCTGAATCAGCGACAAATGGGCGCTGAGGATAGTTTTGCTTTCGCCGTCACGCTCACGCAATTGCTGGTTCAGTTGCTCGGCAAGCGTTGCCAGGCTGTGCTCCAGCCGGGTGGAATCTTGCGCACTGGCGGGGAt                           <  2:435123/150‑1 (MQ=255)
                                                                                                                      aaaTTCATCATCCTGAATCAGCGACAAATGGGCGCTGAGGATAGTTTTGCTTTCGCCGTCACGCTCACGCAATTGCTGGTTCAGTTGCTCGGCAAGCGTTGCCAGGCTGTGCTCCAGCCGGGTGGAATCTTGCGCACTGGCGGGGATTGc                        >  2:1328930/1‑150 (MQ=255)
                                                                                                                         ttCATCATCCTGAATCAGCGACAAATGGGCGCTGAGGATAGTTTTGCTTTCGCCGTCACGCTCACGCAATTGCTGGTTCAGTTGCTCGGCAAGCGTTGCCAGGCTGTGCTCCAGCCGGGTGGAATCTTGCGCACTGGCGGGGATTGCCCg                     <  2:1149386/150‑1 (MQ=255)
                                                                                                                           catcatCCTGAATCAGCGACAAATGGGCGCTGAGGATAGTTTTGCTTTCGCCGTCACGCTCACGCAATTGCTGGTTCAGTTGCTCGGCAAGCGTTGCCAGGCTGTGCTCCAGCCGGGTGGAATCTTGCGCACTGGCg                                <  2:755579/137‑1 (MQ=255)
                                                                                                                           catcatCCTGAATCAGCGACAAATGGGCGCTGAGGATAGTTTTGCTTTCGCCGTCACGCTCACGCAATTGCTGGTTCAGTTGCTCGGCAAGCGTTGCCAGGCTGTGCTCCAGCCGGGTGGAATCTTGCGCACTGGCg                                >  1:755579/1‑137 (MQ=255)
                                                                                                                           catcatCCTGAATCAGCGACAAATGGGCGCTGAGGATAGTTTTGCTTTCGCCGTCACGCTCACGCAATTGCTGGTTCAGTTGCTCGGCAAGCGTTGCCAGGCTGTGCTCCAGCCGGGTGGAATCTTGCGCACTGGCGGGGATTGCCCGAt                   <  2:1026999/150‑1 (MQ=255)
                                                                                                                            atcatcCTGAATCAGCGACAAATGGGCGCTGAGGATAGTTTTGCTTTCGCCGTCACGCTCACGCAATTGCTGGTTCAGTTGCTCGGCAAGCGTTGCCAGGCTGTGCTCCAGCCGGGTGGAATCTTGCGCACTGGCGGGGATTGCCCGATa                  <  1:507203/150‑1 (MQ=255)
                                                                                                                              catcCTGACTCAGCGACAAATGGGCGCTGAGGATAGTTTTGCTTTCGCCGTCACGCTCACGCAATTGCTGGTTCAGTTGCTCGGCAAGCGTTGCCAGGCTGTGCTCCAGCCGGGTGGAATCTTGCGCACTGGCGGGGATTGCCCGATAAc                <  1:196857/150‑1 (MQ=255)
                                                                                                                              catcCTGAATCAGCGACAAATGGGCGCTGAGGATAGTTTTGCTTTCGCCATCACGCTCACGCAATTGCTGGTTCAGTTGCTCGGCAAGCGTTGCCAGGCTGTGCTCCAGCCGGGTGGAATCTTGCGCACTGGCGGGGATTGCCCGATAAc                >  1:549234/1‑150 (MQ=255)
                                                                                                                               atcCTGAATCAGCGACAAATGGGCGCTGAGGATAGTTTTGCTTTCGCCGTCACGCTCACGCAATTGCTGGTTCAGTTGCTCGGCAAGCGTTGCCAGGCTGTGCTCCAGCCGGGTGGAATCTTGCGCACTGGCGGGGATTGCCCGATAACt               >  2:354436/1‑150 (MQ=255)
                                                                                                                                      aTCAGCGACAAATGGGCGCTGAGGATAGTTTTGCTTTCGCCGTCACGCTCACGCAATTGCTGGTTCAGTTGCTCGGCAAGCGTTGCCAGGCTGTGCTCCAGCCGGGTGGAATCTTGCGCACTGGCGGGGATTGCCCGATAACTGTCGAgg        >  2:323479/1‑150 (MQ=255)
                                                                                                                                       tCAGCGACAAATGGGCGCTGAGGATAGTTTTGCTTTCGCCGTCACGCTCACGCAATTGCTGGTTCAGTTGCTCGGCTAGCGTTGCCAGGCTGTGCTCCAGCCGGGTGGAATCTTGCGCACTGGCGGGGATTGCCCGATAACTGTCGAGGc       <  1:276467/150‑1 (MQ=255)
                                                                                                                                       tCAGCGACAAATGGGCGCTGAGGATAGTTTTGCTTTCGCCGTCACGCTCACGCAATTGCTGGTTCAGTTGCTCGGCAAGCGTTGCCAGGCTGTGCTCCAGCCGGGTGGAATCTTGCGCACTGGCGGGGATTGCCCGATAACTGTCGAGGc       <  2:481941/150‑1 (MQ=255)
                                                                                                                                          gCGACAAATGGGCGCTGAGGATAGTTTTGCTTTCGCCGTCACGCTCACGCAATTGCTGGTTCAGTTGCTCGGCAAGCGTTGCCAGGCTGTGCTCCAGCCGGGTGGAATCTTGCGCACTGGCGGGGATTGCCCGATAACTGTCGAGGCTGt    >  1:1076881/1‑150 (MQ=255)
                                                                                                                                           cGACAAATGGGCGCTGAGGATAGTTTTGCTTTCGCCGTCACGCTCACGCAATTGCTGGTTCAGTTGCTCGGCAAGCGTTGCCAGGCTGTGCTCCAGCCGGGTGGAATCTTGCGCACTGGCGGGGATTGCCCGATAACTGTCGAGGCTGTc   >  1:464136/1‑150 (MQ=255)
                                                                                                                                           cGACAAATGGGCGCTGAGGATAGTTTTGCTTTCGCCGTCACGCTCACGCAATTGCTGGTTCAGTTGCTCGGCAAGCGTTGCCAGGCTGTGCTCCAGCCGGGTGGAATCTTGCGCACTGGCGGGGATTGCCCGATAACTGTCGAGGCTGTc   >  2:1136357/1‑150 (MQ=255)
                                                                                                                                            gACAAATGGGCGCTGAGGATAGTTTTGCTTTCGCCGTCACGCTCACGCAATTGCTGGTTCAGTTGCTCGGCAAGCGTTGCCAGGCTGTGCTCCAGCCGGGTGGAATCTTGCGCACTGGCGGGGATTGCCCGATAACTGTCGAGGCTGTCg  >  2:1575509/1‑150 (MQ=255)
                                                                                                                                            gACAAATGGGCGCTGAGGATAGTTTTGCTTTCGCCGTCACGCTCACGCAATTGCTGGTTCAGTTGCTCGGCAAGCGTTGCCAGGCTGTGCTCCAGCCGGGTGGAATCTTGCGCACTGGCGGGGATTGCCCGATAACTGTCGAGGCTGTCg  >  2:961936/1‑150 (MQ=255)
                                                                                                                                             |                                                                                                                                                    
CGCGCAGATAATCGCTGGCAGAGGCAGAAAGTTTGGCGCAAACCTGCTCCATATTGCTGATGATCGCCGCCCCCAGCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCTGAATCAGCGACAAATGGGCGCTGAGGATAGTTTTGCTTTCGCCGTCACGCTCACGCAATTGCTGGTTCAGTTGCTCGGCAAGCGTTGCCAGGCTGTGCTCCAGCCGGGTGGAATCTTGCGCACTGGCGGGGATTGCCCGATAACTGTCGAGGCTGTCG  >  NC_000913/2503789‑2504078

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