Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,504,205 T→A 5.4% Q94L (CAG→CTG)  fryA ← putative PTS multiphosphoryl transfer protein FryA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,504,2050TA 5.4% 301.0 / 5.4 93Q94L (CAG→CTG) fryAputative PTS multiphosphoryl transfer protein FryA
Reads supporting (aligned to +/- strand):  ref base T (45/43);  new base A (2/3);  total (47/46)
Fisher's exact test for biased strand distribution p-value = 6.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.57e-01

CCCGATAACTGTCGAGGCTGTCGCTTTGTAACAGGGTCAGCGTACCCACGCCGACGCCGCTTGCCAGCACATTGCCGTACAGTAAATCCGGGTTAAGGCGGCTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAACGCTGTCGCTGTCGATAAAGCGCACCTGGATGTACTCTTCCAGCACGCGCCGCGCCTGCTCTTCATCGCTGCCGCTAATGTTCAGGCTGCAACTGTCATTAAATAGGGTGCCGGTGCCAATCAGCGCCAGCGAGCTTTTCGCATC  >  NC_000913/2504056‑2504353
                                                                                                                                                     |                                                                                                                                                    
cccGATAACTGTCGAGGCTGTCGCTTTGTAACAGGGTCAGCGTACCCACGCCGACGCCGCTTGCCAGCACATTGCCGTACAGTAAATCCGGGTTAAGGCGGCTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCt                                                                                                                                                      <  2:140265/150‑1 (MQ=255)
 ccGATAACTGTCGAGGCTGTCGCTTTGTAACAGGGTCAGCGTACCCACGCCGACGCCGCTTGCCAGCACATTGCCGTACAGTAAATCCGGGTTAAGGCGGCTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTg                                                                                                                                                     >  2:434848/1‑150 (MQ=255)
    aTAACTGTCGAGGCTGTCGCTTTGTAACAGGGTCAGCGTACCCACGCCGACGCCGCTTGCCAGCACATTGCCGTACAGTAAATCCGGGTTAAGGCGGCTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCAGCCTACGCAGGCAGAAc                                                                                                                                                  <  1:81157/150‑1 (MQ=255)
     tAACTGTCGAGGCTGTCGCTTTGTAACAGGGTCAGCGTACCCACGCCGACGCCGCTTGCCAGCACATTGCCGTACAGTAAATCCGGGTTAAGGCGGCTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCAGCCTACGCAGGCAGAAct                                                                                                                                                 <  1:1257317/150‑2 (MQ=255)
      aaCTGTCGAGGCTGTCGCTTTGTAACAGGGTCAGCGTACCCACGCCGACGCCGCTTGCCAGCACATTGCCGTACAGTAAATCCGGGTTAAGGCGGCTGAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAAcgc                                                                                                                                                >  2:1522023/1‑150 (MQ=255)
       aCTGTCGAGGCTGTCGCTTTGTAACAGGGTCAGCGTACCCACGCCGACGCCGCTTGCCAGCACATTGCCGTACAGTAAATCCGGGTTAAGGCGGCTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAAcgct                                                                                                                                               >  1:1422795/1‑150 (MQ=255)
         tGTCGAGGCTGTCGCTTTGTAACAGGGTCAGCGTACCCACGCCGACGCCGCTTGCCAGCACATTGCCGTACAGTAAATCCGGGTTAAGGCGGCTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCAGCCTACGCAGGCAGAActgacg                                                                                                                                             >  2:1049438/1‑145 (MQ=255)
          gTCGAGGCTGTCGCTTTGTAACAGGGTCAGCGTACCCACGCCGACGCCGCTTGCCAGCACATTGCCGTACAGTAAATCCGGGTTAAGGCGGCTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAAcgctgtc                                                                                                                                            >  2:175381/1‑150 (MQ=255)
          gTCGAGGCTGTCGCTTTGTAACAGGGTCAGCGTACCCACGCCGACGCCGCTTGCCAGCACATTGCCGTACAGTAAATCCGGGTTAAGGCGGCTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCAGCCTACGCAGGCAGAActgacgg                                                                                                                                            <  1:1049438/150‑7 (MQ=255)
            cGAGGCTGTCGCTTTGTAACAGGGTCAGCGTACCCACGCCGACGCCGCTTGCCAGCACATTGCCGTACAGTAAATCCGGGTTAAGGCGGCTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAAcgctgtcgc                                                                                                                                          <  1:1515620/150‑1 (MQ=255)
            cGAGGCTGTCGCTTTGTAACAGGGTCAGCGTACCCACGCCGACGCCGCTTGCCAGCACATTGCCGTACAGTAAATCCGGGTTAAGGCGGCTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAAcgctgtcgc                                                                                                                                          >  1:630363/1‑150 (MQ=255)
             gAGGCTGTCGCTTTGTAACAGGGTCAGCGTACCCACGCCGACGCCGCTTGCCAGCACATTGCCGTACAGTAAATCCGGGTTAAGGCGGCTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAAcgctgtcgct                                                                                                                                         <  1:697599/150‑1 (MQ=255)
                gCTGTCGCTTTGTAACAGGGTCAGCGTACCCACGCCGACGCCGCTTGCCAGCACATTGCCGTACAGTAAATCCGGGTTAAGGCGGCTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAACgctgtcgctgtc                                                                                                                                      <  1:161015/150‑1 (MQ=255)
                gCTGTCGCTTTGTAACAGGGTCAGCGTACCCACGCCGACGCCGCTTGCCAGCACATTGCCGTACAGTAAATCCGGGTTAAGGCGGCTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAACgctgtcgctgtc                                                                                                                                      >  1:118139/1‑150 (MQ=255)
                   gTCGCTTTGTAACAGGGTCAGCGTACCCACGCCGACGCCGCTTGCCAGCACATTGCCGTACAGTAAATCCGGGTTAAGGCGGCTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAACGCTGTCGCTGTCGAt                                                                                                                                   >  1:607808/1‑150 (MQ=255)
                   gTCGCTTTGTAACAGGGTCAGCGTACCCACGCCGACGCCGCTTGCCAGCACATTGCCGTACAGTAAATCCGGGTTAAGGCGGCTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAACGCTGTCGCTGTCGAt                                                                                                                                   <  2:1062547/150‑1 (MQ=255)
                   gTCGCTTTGTAACAGGGTCAGCGTACCCACGCCGACGCCGCTTGCCAGCACATTGCCGTACAGTAAATCCGGGTTAAGGCGGCTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAACGCTGTCGCTGTCGAt                                                                                                                                   >  1:314847/1‑150 (MQ=255)
                        tttGTAACAGGGTCAGCGTACCCACGCCGACGCCGCTTGCCAGCACATTGCCGTACAGTAAATCCGGGTTAAGGCGGCTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAACGCTGTCGCTGTCGATAAAgc                                                                                                                              <  1:1239962/150‑1 (MQ=255)
                        tttGTAACAGGGTCAGCGTACCCACGCCGACGCCGCTTGCCAGCACATTGCCGTACAGTAAATCCGGGTTAAGGCGGCTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAACGCTGTCGCTGTCGATAAAgc                                                                                                                              <  1:670778/150‑1 (MQ=255)
                        tttGTAACAGGGTCAGCGTACCCACGCCGACGCCGCTTGCCAGCACATTGCCGTACAGTAAATCCGGGTTAAGGCGGCTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAACGCTGTCGCTGTCGATAAAgc                                                                                                                              <  1:675072/150‑1 (MQ=255)
                             aaCAGGGTCAGCGTACCCACGCCGACGCCGCTTGCCAGCACATTGCCGTACAGTAAATCCGGGTTAAGGCGGCTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAACGCTGTCGCTGTCGATAAAGCGCAcc                                                                                                                         >  1:1256162/1‑150 (MQ=255)
                             aaCAGGGTCAGCGTACCCACGCCGACGCCGCTTGCCAGCACATTGCCGTACAGTAAATCCGGGTTAAGGCGGCTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAACGCTGTCGCTGTCGATAAAGCGCAcc                                                                                                                         >  1:328068/1‑150 (MQ=255)
                               cAGGGTCAGCGTACCCACGCCGACGCCGCTTGCCAGCACATTGCCGTACAGTAAATCCGGGTTAAGGCGGCTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAACGCTGTCGCTGTCGATAAAGCGCACCTg                                                                                                                       >  2:238841/1‑150 (MQ=255)
                                  ggTCAGCGTACCCACGCCGACGCCGCTTGCCAGCACATTGCCGTACAGTAAATCCGGGTTAAGGCGGCTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAACGCTGTCGCTGTCGATAAAGCGCACCTGGAt                                                                                                                    >  2:1187351/1‑150 (MQ=255)
                                   gTCAGCGTACCCACGCCGACGCCGCTTGCCAGCACATTGCCGTACAGTAAATCCGGGTTAAGGCGGCTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAACgctgtcgctgtc                                                                                                                                      >  1:531595/1‑131 (MQ=255)
                                   gTCAGCGTACCCACGCCGACGCCGCTTGCCAGCACATTGCCGTACAGTAAATCCGGGTTAAGGCGGCTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAACgctgtcgctgtc                                                                                                                                      <  2:531595/131‑1 (MQ=255)
                                    tCAGCGTACCCACGCCGACGCCGCTTGCCAGCACATTGCCGTACAGTAAATCCGGGTTAAGGCGGCTGAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAACGCTGTCGCTGTCGATAAAGCGCACCTGGATGt                                                                                                                  <  1:1522023/150‑1 (MQ=255)
                                     cAGCGTACCCACGCCGACGCCGCTTGCTAGCACATTGCCGTACAGTAAATCCGGGTTAAGGCGGCTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAACGCTGTCGCTGTCGATAAAGCGCACCTGGATGTa                                                                                                                 <  2:65992/150‑1 (MQ=255)
                                      aGCGTACCCACGCCGACGCCGCTTGCCAGCACATTGCCGTACAGTAAATCCGGGTTAAGGCGGCTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAACGCTGTCGCTGTCGATAAAGCGCACCTGGATGTAc                                                                                                                <  1:109216/150‑1 (MQ=255)
                                      aGCGTACCCACGCCGACGCCGCTTGCCAGCACATTGCCGTACAGTAAATCCGGGTTAAGGCGGCTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAACGCTGTCGCTGTCGATAAAGCGCACCTGGATGTAc                                                                                                                >  2:640896/1‑150 (MQ=255)
                                      aGCGTACCCACGCCGACGCCGCTTGCCAGCACATTGCCGTACAGTAAATCCGGGTTAAGGCGGCTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAACGCTGTCGCTGTCGATAAAGCGCACCTGGATGTAc                                                                                                                >  2:61876/1‑150 (MQ=255)
                                          tACCCACGCCGACGCCGCTTGCCAGCACATTGCCGTACAGTAAATCCGGGTTAAGGCGGCTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAACGCTGTCGCTGTCGATAAAGCGCACCTGGATGTACTCtt                                                                                                            >  2:1603500/1‑150 (MQ=255)
                                            cccACGCCGACGCCGCTTGCCAGCACATTGCCGTACAGTAAATCCGGGTTAAGGCGGCTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAACGCTGTCGCTGTCGATAAAGCGCACCTGGATGTACTCTTcc                                                                                                          >  1:598789/1‑150 (MQ=255)
                                            cccACGCCGACGCCGCTTGCCAGCACATTGCCGTACAGTAAATCCGGGTTAAGGCGGCTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAACGCTGTCGCTGTCGATAAAGCGCACCTGGATGTACTCTTcc                                                                                                          >  2:800309/1‑150 (MQ=255)
                                            cccACGCCGACGCCGCTTGCCAGCACATTGCCGTACAGTAAATCCGGGTTAAGGCGGCTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAACGCTGTCGCTGTCGATAAAGCGCACCTGGATGTACTCTTcc                                                                                                          >  2:20076/1‑150 (MQ=255)
                                            cccACGCCGACGCCGCTTGCCAGCACATTGCCGTACAGTAAATCCGGGTTAAGGCGGCTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAACGCTGTCGCTGTCGATAAAGCGCACCTGGATGTACTCTTcc                                                                                                          <  2:1197633/150‑1 (MQ=255)
                                              cACGCCGACGCCGCTTGCCAGCACATTGCCGTACAGTAAATCCGGGTTAAGGCGGCTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAACGCTGTCGCTGTCGATAAAGCGCACCTGGATGTACTCTTCCAg                                                                                                        >  2:552657/1‑150 (MQ=255)
                                              cACGCCGACGCCGCTTGCCAGCACATTGCCGTACAGTAAATCCGGGTTAAGGCGGCTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAACGCTGTCGCTGTCGATAAAGCGCACCTGGATGTACTCTTCCAg                                                                                                        >  1:250598/1‑150 (MQ=255)
                                                cgccgacgccgCTTGCCAGCACATTGCCGTACAGTAAATCCGGGTTAAGGCGGCTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAACGCTGTCGCTGTCGATAAAGCGCACCTGGATGTACTCTTCCAGCa                                                                                                      >  1:320219/1‑150 (MQ=255)
                                                cgccgacgccgCTTGCCAGCACATTGCCGTACAGTAAATCCGGGTTAAGGCGGCTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAACGCTGTCGCTGTCGATAAAGCGCACCTGGATGTACTCTTCCAGCa                                                                                                      <  1:839255/150‑1 (MQ=255)
                                                         cgcTTGCCAGCACATTGCCGTACAGTAAATCCGGGTTAAGGCGGCTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAACGCTGTCGCTGTCGATAAAGCGCACCTGGATGTACTCTTCCAGCAcgcgccgcg                                                                                             <  1:434848/150‑1 (MQ=255)
                                                           cTTGCCAGCACATTGCCGTACAGTAAATCCGGGTTAAGGCGGCTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAACGCTGTCGCTGTCGATAAAGCGCACCTGGATGTACTCTTCCAGCACgcgccgcgcc                                                                                           <  1:600882/150‑1 (MQ=255)
                                                           cTTGCCAGCACATTGCCGTACAGTAAATCCGGGTTAAGGCGGCTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAACGCTGTCGCTGTCGATAAAGCGCACCTGGATGTACTCTTCCAGCACgcgccgcgcc                                                                                           <  1:44194/150‑1 (MQ=255)
                                                           cTTGCCAGCACATTGCCGTACAGTAAATCCGGGTTAAGGCGGCTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAACGCTGTCGCTGTCGATAAAGCGCACCTGGATGTACTCTTCCAGCACgcgccgcgcc                                                                                           <  1:20076/150‑1 (MQ=255)
                                                             tGCCAGCACATTGCCGTACAGTAAATCCGGGTTAAGGCGGCTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAACGCTGTCGCTGTCGATAAAGCGCACCTGGATGTACTCTTCCAGCACGCGCCGCGCCTg                                                                                         <  2:1524147/150‑1 (MQ=255)
                                                             tGCCAGCACATTGCCGTACAGTAAATCCGGGTTAAGGCGGCTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAACGCTGTCGCTGTCGATAAAGCGCACCTGGATGTACTCTTCCAGCACGCGCCGCGCCTg                                                                                         >  1:165431/1‑150 (MQ=255)
                                                             tGCCAGCACATTGCCGTACAGTAAATCCGGGTTAAGGCGGCTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAACGCTGTCGCTGTCGATAAAGCGCACCTGGATGTACTCTTCCAGCACGCGCCGCGCCTg                                                                                         <  2:1191277/150‑1 (MQ=255)
                                                              gCCAGCACATTGCCGTACAGTAAATCCGGGTTAAGGCGGCTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAACGCTGTCGCTGTCGATAAAGCGCACCTGGATGTACTCTTCCAGCACGCGCCGCGCCTGc                                                                                        >  2:4437/1‑150 (MQ=255)
                                                                        tGCCGTACAGTAAATCCGGGTTAAGGCGGCTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAACGCTGTCGCTGTCGATAAAGCGCACCTGGATGTACTCTTCCAGCACGCGCCGCGCCTGCTCTTCATCGc                                                                              <  2:475635/150‑1 (MQ=255)
                                                                          ccGTACAGTAAATCCGGGTTAAGGCGGCTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAACGCTGTCGCTGTCGATAAAGCGCACCTGGATGTACTCTTCCAGCACGCGCCGCGCCTGCTCTTCATCGCTg                                                                            <  2:784997/150‑1 (MQ=255)
                                                                          ccGTACAGTAAATCCGGGTTAAGGCGGCTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAACGCTGTCGCTGTCGATAAAGCGCACCTGGATGTACTCTTCCAGCACGCGCCGCGCCTGCTCTTCATCGCTg                                                                            <  2:1596234/150‑1 (MQ=255)
                                                                             tACAGTAAATCCGGGTTAAGGCGGCTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAACGCTGTCGCTGTCGATAAAGCGCACCTGGATGTACTCTTCCAGCACGCGCCGCGCCTGCTCTTCATCGCTGCCg                                                                         >  2:386239/1‑150 (MQ=255)
                                                                              aCAGTAAATCCGGGTTAAGGCGGCTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAACGCTGTCGCTGTCGATAAAGCGCACCTGGATGTACTCTTCCAGCACGCGCCGCGCCTGCTCTTCATCGCTGCCGc                                                                        <  1:238841/150‑1 (MQ=255)
                                                                              aCAGTAAAGCCGGGGTAAGGCGGCTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAACGCTGTCGCTGTCGATAAAGCGCACCTGGATGTACTCTTCCAGCACGCGCCGCGCCGGCTCTTCATCGCTGccta                                                                        <  2:157601/150‑3 (MQ=255)
                                                                               cAGTAAATCCGGGTTAAGGCGGCTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAACGCTGTCGCTGTCGATAAAGCGCACCTGGATGTACTCTTCCAGCACGCGCCGCGCCTGCTCTTCATCGCTGCCGCt                                                                       <  2:724679/150‑1 (MQ=255)
                                                                                aGTAAATCCGGGTTAAGGCGGCTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAACGCTGTCGCTGTCGATAAAGCGCACCTGGATGTACTCTTCCAGCACGCGCCGCGCCTGCTCTTCATCGCTGCCGCTa                                                                      <  1:61876/150‑1 (MQ=255)
                                                                                 gTAAATCCGGGTTAAGGCGGCTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAACGCTGTCGCTGTCGATAAAGCGCACCTGGATGTACTCTTCCAGCACGCGCCGCGCCTGCTCTTCATCGCTGCCGCaaa                                                                     >  2:1168755/1‑150 (MQ=255)
                                                                                   aaaTCCGGGTTAAGGCGGCTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAACGCTGTCGCTGTCGATAAAGCGCACCTGGATGTACTCTTCCAGCACGCGCCGCGCCTGCTCTTCATCGCTGCCGCTAATg                                                                   <  2:1256162/150‑1 (MQ=255)
                                                                                      tCCGGGTTAAGGCGGCTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAACGCTGTCGCTGTCGATAAAGCGCACCTGGATGTACTCTTCCAGCACGCGCCGCGCCTGCTCTTCATCGCTGCCGCTAATGTTc                                                                <  1:175381/150‑1 (MQ=255)
                                                                                          ggTTAAGGCGGCTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAACGCTGTCGCTGTCGATAAAGCGCACCTGGATGTACTCTTCCAGCACGCGCCGCGCCTGCTCTTCATCGCTGCCGCTAATGTTCAGGc                                                            <  2:978538/150‑1 (MQ=255)
                                                                                                     cTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAACGCTGTCGCTGTCGATAAAGCGCACCTGGATGTACTCTTCCAGCACGCGCCGCGCCTGCTCTTCATCGCTGCCGCTAATGTTCAGGCTGCAACTGTCa                                                 >  2:936650/1‑150 (MQ=255)
                                                                                                     cTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAACGCTGTCGCTGTCGATAAAGCGCACCTGGATGTACTCTTCCAGCACGCGCCGCGCCTGCTCTTCATCGCTGCCGCTAATGTTCAGGCTGCAACTGTCa                                                 <  2:607808/150‑1 (MQ=255)
                                                                                                       tAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAACGCTGTCGCTGTCGATAAAGCGCACCTGGATGTACTCTTCCAGCACGCGCCGCGCCTGCTCTTCATCGCTGCCGCTAATGTTCAGGCTGCAACTGTCAtt                                               <  1:918415/150‑1 (MQ=255)
                                                                                                        aaTGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAACGCTGTCGCTGTCGATAAAGCGCACCTGGATGTACTCTTCCAGCACGCGCCGCGCCTGCTCTTCATCGCTGCCGCTAATGTTCAGGCTGCAACTGTCATTa                                              <  1:277500/150‑1 (MQ=255)
                                                                                                         aTGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAACGCTGTCGCTGTCGATAAAGCGCACCTGGATGTACTCTTCCAGCACGCGCCGCGCCTGCTCTTCATCGCTGCCGCTAATGTTCAGGCTGCAACTGTCATTaa                                             <  1:1603500/150‑1 (MQ=255)
                                                                                                           gAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAACGCTGTCGCTGTCGATAAAGCGCACCTGGATGTACTCTTCCAGCACGCGCCGCGCCTGCTCTTCATCGCTGCCGCTAATGTTCAGGCTGCAACTGTCATTAAAt                                           >  2:56185/1‑150 (MQ=255)
                                                                                                           gAACGCGGCAGAGGATGCGCCGTCAGTTCTGCCTTCGTAGTCTGAACGCTGTCGCTGTCGATAAAGCGCACCTCGATGTACTCATCCAGCACGCGCCGCGCCTGCTCTTCATCGCTGCCGCTAATGTTCAGGATGCAACTGTCATTAAAt                                           >  2:833814/1‑150 (MQ=255)
                                                                                                               gcgGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAACGCTGTCGCTGTCGATAAAGCGCACCTGGATGTACTCTTCCAGCACGCGCCGCGCCTGCTCTTCATCGCTGCCGCTAATGTTCAGGCTGCAACTGTCATTAAATAggg                                       <  1:4437/150‑1 (MQ=255)
                                                                                                               gcgGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAACGCTGTCGCTGTCGATAAAGCGCACCTGGATGTACTCTTCCAGCACGCGCCGCGCCTGCTCTTCATCGCTGCCGCTAATGTTCAGGCTGCAACTGTCATTAAATAggg                                       <  2:118139/150‑1 (MQ=255)
                                                                                                                cgGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAACGCTGTCGCTGTCGATAAAGCGCACCTGGATGTACTCTTCCAGCACGCGCCGCGCCTGCTCTTCATCGCTGCCGCTAATGTTCAGGCTGCAACTGTCATTAAATAGggt                                      >  1:1663401/1‑150 (MQ=255)
                                                                                                                    aGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAACGCTGTCGCTGTCGATAAAGCGCACCTGGATGTACTCTTCCAGCACGCGCCGCGCCTGCTCTTCATCGCTGCCGCTAATGTTCAGGCTGCAACTGTCATTAAATAGggtgccg                                  >  2:845291/1‑150 (MQ=255)
                                                                                                                      cGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAACGCTGTCGCTGTCGATAAAGCGCACCTGGATGTACTCTTCCAGCACGCGCCGCGCCTGCTCTTCATCGCTGCCGCTAATGTTCAGGCTGCAACTGTCATTAAATAGggtgccggt                                <  1:974535/150‑1 (MQ=255)
                                                                                                                         aTGCGCCGTGAGTTCTGCCTGCGTAGGCTGAACGCTGTCGCTGTCGATAAAGCGCACCTGGATGTACTCTTCCAGCACGCGCCGCGCCTGCTCTTCATCGCTGCCGCTAATGTTCAGGCTGCAACTGTCATTAAATAGggtgccggtgcc                             <  1:1454215/150‑1 (MQ=255)
                                                                                                                              ccGTCAGTTCTGCCTGCGTAGGCTGAACGCTGTCGCTGTCGATAAAGCGCACCTGGATGTACTCTTCCAGCACGCGCCGCGCCTGCTCTTCATCGCTGCCGCTAATGTTCAGGCTGCAACTGTCATTAAATAGGGTGCCGGTGCCAATCa                        >  1:1428100/1‑150 (MQ=255)
                                                                                                                               acagcgTTCAGCCTACGCAGGCAGAACGCTGTCGCTGTCGATAAAGCGCACCTGGATGTACTCTTCCAGCACGCGCCGCGCCTGCTCTTCATCGCTGCCGCTAATGTTCAGGCTGCAACTGTCATTAAATAGGGTGCCGGTGCCAATCAg                       >  2:141925/6‑150 (MQ=255)
                                                                                                                                gTCAGTTCTGCCTGCGTAGGCTGAACGCTGTCGCTGTCGATAAAGCGCACCTGGATGTACTCTTCCAGCACGCGCCGCGCCTGCTCTTCATCGCTGCCGCTAATGTTCAGGCTGCAACTGTCATTAAATAGggtg                                     >  1:148277/1‑135 (MQ=255)
                                                                                                                                gTCAGTTCTGCCTGCGTAGGCTGAACGCTGTCGCTGTCGATAAAGCGCACCTGGATGTACTCTTCCAGCACGCGCCGCGCCTGCTCTTCATCGCTGCCGCTAATGTTCAGGCTGCAACTGTCATTAAATAGggtg                                     <  2:148277/135‑1 (MQ=255)
                                                                                                                                  cAGTTCTGCCTGCGTAGGCTGAACGCTGTCGCTGTCGATAAAGCGCACCTGGATGTACTCTTCCAGCACGCGCCGCGCCTGCTCTTCATCGCTGCCGCTAATGTTCAGGCTGCAACTGTCATTAAATAGGGTGCCGGTGCCAATCGgcgc                    >  2:971096/1‑150 (MQ=255)
                                                                                                                                   aGTTCTGCCTGCGTAGGCTGAACGCTGTCGCTGTCGATAAAGCGCACCTGGATGTACTCTTCCAGCACGCGCCGCGCCTGCTCTTCATCGCTGCCGCTAATGTTCAGGCTGCAACTGTCATTAAATAGGGTGCCGGTGCCAATCAGCGcc                   >  1:1128586/1‑150 (MQ=255)
                                                                                                                                   aGTTCTGCCTGCGTAGGCTGAACGCTGTCGCTGTCGATAAAGCGCACCTGGATGTACTCTTCCAGCACGCGCCGCGCCTGCTCTTCATCGCTGCCGCTAATGTTCAGGCTGCAACTGTCATTAAATAGGGTGCCGGTGCCAATCAGCGcc                   <  2:328068/150‑1 (MQ=255)
                                                                                                                                   aGTTCTGCCTGCGTAGGCTGAACGCTGTCGCTGTCGATAAAGCGCACCTGGATGTACTCTTCCAGCACGCGCCGCGCCTGCTCTTCATCGCTGCCGCTAATGTTCAGGCTGCAACTGTCATTAAATAGGGTGCCGGTGCCAATCAGCGcc                   >  1:1502105/1‑150 (MQ=255)
                                                                                                                                        tgcctgcGTAGGCTGAACGCTGTCGCTGTCGATAAAGCGCACCTGGATGTACTCTTCCAGCACGCGCCGCGCCTGCTCTTCATCGCTGCCGCTAATGTTCAGGCTGCAACTGTCATTAAATAGGGTGCCGGTGCCAATCAGCGCCAGCGa              >  2:479359/1‑150 (MQ=255)
                                                                                                                                         gcctgcGTAGGCTGAACGCTGTCGCTGTCGATAAAGCGCACCTGGATGTACTCTTCCAGCACGCGCCGCGCCTGCTCTTCATCGCTGCCGCTAATGTTCAGGCTGCAACTGTCATTAAATAGGGTGCCGGTGCCAATCAGCGCCAGCGAg             >  1:663734/1‑150 (MQ=255)
                                                                                                                                         gcctgcGTAGGCTGAACGCTGTCGCTGTCGATAAAGCGCACCTGGATGTACTCTTCCAGCACGCGCCGCGCCTGCTCTTCATCGCTGCCGCTAATGTTCAGGCTGCAACTGTCATTAAATAGGGTGCCGGTGCCAATCAGCGCCAGCGAg             >  1:230749/1‑150 (MQ=255)
                                                                                                                                           ctgcGTAGGCTGAACGCTGTCGCTGTCGATAAAGCGCACCTGGATGTACTCTTCCAGCACGCGCCGCGCCTGCTCTTCATCGCTGCCGCTAATGTTCAGGCTGCAACTGTCATTAAATAGGGTGTCGGTGCCAATCAGCGCCAGCGAGCt           >  2:1619452/1‑150 (MQ=255)
                                                                                                                                           ctgcGTAGGCTGAACGCTGTCGCTGTCGATAAAGCGCACCTGGATGTACTCTTCCAGCACGCGCCGCGCCTGCTCTTCATCGCTGCCGCTAATGTTCAGGCTGCAACTGTCATTAAATAGGGTGCCGGTGCCAATCAGCGCCAGCGAGCt           <  2:1125767/150‑1 (MQ=255)
                                                                                                                                                 aGGCTGAACGCTGTCGCTGTCGATAAAGCGCACCTGGATGTACTCTTCCAGCACGCGCCGCGCCTGCTCTTCATCGCTGCCGCTAATGTTCAGGCTGCAACTGTCATTAAATAGGGTGCCGGTGCCAATCAGCGCCAGCGAGCTTTTCGc     <  2:38711/150‑1 (MQ=255)
                                                                                                                                                  ggCTGAACGCTGTCGCTGTCGATAAAGCGCACCTGGATGTACTCTTCCAGCACGCGCCGCGCCTGCTCTTCATCGCTGCCGCTAATGTTCAGGCTGCAACTGTCATTAAATAGGGTGCCGGTGCCAATCAGCGCCAGCGAGCTTTTCGCa    >  2:397661/1‑150 (MQ=255)
                                                                                                                                                  ggCTGAACGCTGTCGCTGTCGATAAAGCGCACCTGGATGTACTCTTCCAGCACGCGCCGCGCCTGCTCTTCATCGCTGCCGCTAATGTTCAGGCTGCAACTGTCATTAAATAGGGTGCCGGTGCCAATCAGCGCCAGCGAGCTTTTCGCa    <  1:1247896/150‑1 (MQ=255)
                                                                                                                                                  ggCTGAACGCTGTCGCTGTCGATAAAGCGCACCTGGATGTACTCTTCCAGCACGCGCCGCGCCTGCTCTTCATCGCTGCCGCTAATGTTCAGGCTGCAACTGTCATTAAATAGGGTGCCGGTGCCAATCAGCGCCAGCGAGCTTTTCGCa    >  1:77817/1‑150 (MQ=255)
                                                                                                                                                  ggCTGAACGCTGTCGCTGTCGATAAAGCGCACCTGGATGTACTCTTCCAGCACGCGCCGCGCCTGCTCTTCATCGCTGCCGCTAATGTTCAGGCTGCAACTGTCATTAAATAGGGTGCCGGTGCCAATCAGCGCCAGCGAGCTTTTCGCa    <  2:1422795/150‑1 (MQ=255)
                                                                                                                                                  ggCTGAACGCTGTCGCTGTCGATAAAGCGCACCTGGATGTACTCTTCCAGCACGCGCCGCGCCTGCTCTTCATCGCTGCCGCTAATGTTCAGGCTGCAACTGTCATTAAATAGGGTGCCGGTGCCAATCAGCGCCAGCGAGCTTTTCGCa    >  2:1399561/1‑150 (MQ=255)
                                                                                                                                                   gCTGAACGCTGTCGCTGTCGATAAAGCGCACCTGGATGTACTCTTCCAGCACGCGCCGCGCCTGCTCTTCATCGCTGCCGCTAATGTTCAGGCTGCAACTGTCATTAAATAGGGTGCCGGTGCCAATCAGCGCCAGAGAGCTATTCGCAt   >  2:1514254/1‑150 (MQ=255)
                                                                                                                                                    cTGAAAGCTGTCGCTGTCGATAAAGCGCACCTGGATGTACTCTTCCAGCACGCGCCGCGCCTGCTCTTCATCGCTGCCGCTAATCTTCAGGCTGCAACTGTCATTAAATAGGGTGCCGGTGCCAATCAGCGCCAGCGAGCTTTTCGCATc  >  2:496588/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CCCGATAACTGTCGAGGCTGTCGCTTTGTAACAGGGTCAGCGTACCCACGCCGACGCCGCTTGCCAGCACATTGCCGTACAGTAAATCCGGGTTAAGGCGGCTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAACGCTGTCGCTGTCGATAAAGCGCACCTGGATGTACTCTTCCAGCACGCGCCGCGCCTGCTCTTCATCGCTGCCGCTAATGTTCAGGCTGCAACTGTCATTAAATAGGGTGCCGGTGCCAATCAGCGCCAGCGAGCTTTTCGCATC  >  NC_000913/2504056‑2504353

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