Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,971,448 A→C 5.3% E63D (GAA→GAC) ‡ wecC → UDP‑N‑acetyl‑D‑mannosamine dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,971,4480AC 5.3% 434.2 / 10.2 141E63D (GAA→GAC) ‡wecCUDP‑N‑acetyl‑D‑mannosamine dehydrogenase
Reads supporting (aligned to +/- strand):  ref base A (63/70);  new base C (2/6);  total (65/76)
Fisher's exact test for biased strand distribution p-value = 2.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 141 total aligned reads displayed.

GGGCTGCCAACGGCAGCAGCGTTTGCCTCACGGCAAAAACAGGTAATTGGTGTCGATATCAACCAACATGCGGTTGATACCATCAATCGTGGCGAAATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTACGTTGAATCGGCTGCTCGCTCCATTGCGCCAGTGCTGAAAAAAGG  >  NC_000913/3971302‑3971594
                                                                                                                                                  |                                                                                                                                                  
gggCTGCTAACGGCAGCAGCGTTTGCCTCACGGCAAAAACAGGTAATTGGTGTCGATATCAACCAACATGCGGTTGATACCATCAATCGTGGCGAAATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGc                                                                                                                                                 >  2:737699/1‑150 (MQ=255)
gggCTGCCAACGGCAGCAGCGTTTGCCTCACGGCAAAAACAGGTAATTGGTGTCGATATCAACCAACATGCGGTTGATACCATCAATCGTGGCGAAATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGc                                                                                                                                                 >  1:1449601/1‑150 (MQ=255)
   cTGCCAACGGCAGCAGCGTTTGCCTCACGGCAAAAACAGGTAATTGGTGTCGATATCAACCAACATGCGGTTGATACCATCAATCGTGGCGAAATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGt                                                                                                                                              >  2:513848/1‑150 (MQ=255)
    tGCCAACGGCAGCAGCGTTTGCCTCACGGCAAAAACAGGTAATTGGTGTCGATATCAACCAACATGCGGTTGATACCATCAATCGTGGCGAAATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGtt                                                                                                                                             >  1:477293/1‑150 (MQ=255)
     gCCAACGGCAGCAGCGTTTGCCTCACGGCAAAAACAGGTAATTGGTGTCGATATCAACCAACATGCGGTTGATACCATCAATCGTGGCGAAATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGttt                                                                                                                                            >  2:322365/1‑150 (MQ=255)
     gCCAACGGCAGCAGCGTTTGCCTCACGGCAAAAACAGGTAATTGGTGTCGATATCAACCAACATGCGGTTGATACCATCAATCGTGGCGAAATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGttt                                                                                                                                            >  2:1076179/1‑150 (MQ=255)
      ccAACGGCAGCAGCGTTTGCCTCACGGCAAAAACAGGTAATTGGTGTCGATATCAACCAACATGCGGTTGATACCATCAATCGTGGCGAAATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGtttt                                                                                                                                           <  1:480147/150‑1 (MQ=255)
       cAACGGCAGCAGCGTTTGCCTCACGGCAAAAACAGGTAATTGGTGTCGATATCAACCAACATGCGGTTGATACCATCAATCGTGGCGAAATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGttttt                                                                                                                                          >  1:700845/1‑150 (MQ=255)
        aaCGGCAGCAGCGTTTGCCTCACGGCAAAAACAGGTAATTGGTGTCGATATCAACCAACATGCGGTTGATACCATCAATCGTGGCGAAATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGtttttt                                                                                                                                         >  1:1038904/1‑150 (MQ=255)
         aCGGCAGCAGCGTTTGCCTCACGGCAAAAACAGGTAATTGGTGTCGATATCAACCAACATGCGGTTGATACCATCAATCGTGGCGAAATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTa                                                                                                                                        <  2:1267317/150‑1 (MQ=255)
         aCGGCAGCAGCGTTTGCCTCACGGCAAAAACAGGTAATTGGTGTCGATATCAACCAACATGCGGTTGATACCATCAATCGTGGCGAAATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACCGCCTTCTACGGCAGtttttac                                                                                                                                        <  2:71669/150‑3 (MQ=255)
         aCGGCAGCAGCGTTTGCCTCACGGCAAAAACAGGTAATTGGTGTCGATATCAACCAACATGCGGTTGATACCATCAATCGTGGCGAAATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACCGCCTTCTACGGCAGtttttac                                                                                                                                        <  1:353301/150‑3 (MQ=255)
            gcagcaGCGTTTGCCTCACGGCAAAAACAGGTAATTGGTGTCGATATCAACCAACATGCGGTTGATACCATCAATCGTGGCGAAATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTAcga                                                                                                                                     <  2:1181661/150‑1 (MQ=255)
             cagcagCGTTTGCCTCACGGCAAAAACAGGTAATTGGTGTCGATATCAACCAACATGCGGTTGATACCATCAATCGTGGCGAAATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTAcgag                                                                                                                                    >  1:1577845/1‑150 (MQ=255)
             cagcagCGTTTGCCTCACGGCAAAAACAGGTAATTGGTGTCGATATCAACCAACATGCGGTTGATACCATCAATCGTGGCGAAATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTAcgag                                                                                                                                    >  2:1306315/1‑150 (MQ=255)
                  gcGTTTGCCTCACGGCAAAAACAGGTAATTGGTGTCGATATCAACCAACATGCGGTTGATACCATCAATCGTGGCGAAATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACgagcgagc                                                                                                                               <  2:66133/150‑1 (MQ=255)
                     tttGCCTCACGGCAAAAACAGGTAATTGGTGTCGATATCAACCAACATGCGGTTGATACCATCAATCGTGGCGAAATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCacg                                                                                                                            >  1:1295857/1‑150 (MQ=255)
                        gCCTCACGGCAAAAACAGGTAATTGGTGTCGATATCAACCAACATGCGGTTGATACCATCAATCGTGGCGAAATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCacgacg                                                                                                                         >  2:956041/1‑150 (MQ=255)
                         ccTCACGGCAAAAACAGGTAATTGGTGTCGATATCAACCAACATGCGGTTGATACCATCAATCGTGGCGAAATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGc                                                                                                                        <  2:871493/150‑1 (MQ=255)
                         ccTCACGGCAAAAACAGGTAATTGGTGTCGATATCAACCAACATGCGGTTGATACCATCAATCGTGGCGAAATCAATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGc                                                                                                                        >  2:242911/1‑150 (MQ=255)
                          cTCACGGCAAAAACAGGTAATTGGTGTCGATATCAACCAACATGCGGTTGATACCATCAATCGTGGCGAAATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGcc                                                                                                                       >  1:202259/1‑150 (MQ=255)
                          cTCACGGCAAAAACAGGTAATTGGTGTCGATATCAACCAACATGCGGTTGATACCATCAATCGTGGCGAAATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACCGCCTTCTACGGCAGtttttac                                                                                                                                        <  1:1329088/133‑3 (MQ=255)
                          cTCACGGCAAAAACAGGTAATTGGTGTCGATATCAACCAACATGCGGTTGATACCATCAATCGTGGCGAAATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACCGCCTTCTACGGCAGtttttac                                                                                                                                        >  2:1329088/1‑131 (MQ=255)
                           tCACGGCAAAAACAGGTAATTGGTGTCGATATCCACCAACATGCGGTTGATACCATCAATCGTGGCGAAATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCa                                                                                                                      >  1:1199362/1‑150 (MQ=255)
                            cACGGCAAAAACAGGTAATTGGTGTCGATATCAACCAACATGCGGTTGATACCATCAATCGTGGCGAAATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAg                                                                                                                     >  1:749544/1‑150 (MQ=255)
                             aCGGCAAAAACAGGTAATTGGTGTCGATATCAACCAACATGCGGTTGATACCATCAATCGTGGCGAAATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGcc                                                                                                                       <  2:1195920/147‑1 (MQ=255)
                             aCGGCAAAAACAGGTAATTGGTGTCGATATCAACCAACATGCGGTTGATACCATCAATCGTGGCGAAATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGcc                                                                                                                       >  1:1195920/1‑147 (MQ=255)
                               ggCAAAAACAGGTAATTGGTGTCGATATCAACCAACATGCGGTTGATACCATCAATCGTGGCGAAATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTg                                                                                                                  <  2:274614/150‑1 (MQ=255)
                                    aaaCAGGTAATTGGTGTCGATATCAACCAACATGCGGTTGATACCATCAATCGTGGCGAAATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCg                                                                                                             <  1:1362000/150‑1 (MQ=255)
                                    aaaCAGGTAATTGGTGTCGATATCAACCAACATGCGGTTGATACCATCAATCGTGGCGAAATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCg                                                                                                             <  1:1510869/150‑1 (MQ=255)
                                      aCAGTTATTTGGTGTGGATTTAAACCAAGATGCGGTTTTTACGATCTATCGTGGGGAAATCCATTTCGTGGAACCTTATTTGGCGAGTGTTGTAAAAACTGCCTTATACGGCTGTTTTTTACTAGCGAGCACGACGCCAGTTGAAGCGGa                                                                                                           <  2:202259/150‑1 (MQ=255)
                                       cAGGTAATTGGTGTCGATATCAACCAACATGCGGTTGATACCATCAATCGTGGCGAAATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCacgac                                                                                                                          >  2:1193795/1‑134 (MQ=255)
                                       cAGGTAATTGGTGTCGATATCAACCAACATGCGGTTGATACCATCAATCGTGGCGAAATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCacgac                                                                                                                          <  1:1193795/134‑1 (MQ=255)
                                         ggTAATTGGTGTCGATATCAACCAACATGCGGTTGATACCATCAATCGTGGCGAAATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACTAGCGAGCACGACGCCAGTTGAAGCGGATGc                                                                                                        <  1:94590/150‑1 (MQ=255)
                                          gTAATTGGTGTCGATATCAACCAACATGCGGTTGATACCATCAATCGTGGCGAAATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGcc                                                                                                       <  1:313861/150‑1 (MQ=255)
                                            ttttGGTGTCGATTTCAACCAACATGCGGTTTATACCATCATTCTTGGCGAAATCCATTTCTTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACTAGCTAGCACGACGCCAGTTGAAGCGGATGCCTg                                                                                                     <  2:156893/148‑1 (MQ=255)
                                            aaTTGGTGTCGATATCAACCAACATGCGGTTGATACCATCAATCGTGGCGAAATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTg                                                                                                     <  2:687615/150‑1 (MQ=255)
                                              ttGGTGTCGATATCAACCAACATGCGGTTGATACCATCAATCGAGGCGAAATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGTCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGc                                                                                                   <  1:331717/150‑1 (MQ=255)
                                                ggTGTCGATATCAACCAACATGCGGTTGATACCATCAATCGTGGCGAAATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTg                                                                                                 <  2:339391/150‑1 (MQ=255)
                                                ggTGTCGATATCAACCAACATGCGGTTGATACCATCAATCGTGGCGAAATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTg                                                                                                 <  2:345445/150‑1 (MQ=255)
                                                 gtgtCGATATCAACCAACATGCGGTTGATACCATCAATCGTGGCGAAATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGa                                                                                                <  2:821606/150‑1 (MQ=255)
                                                   gtCGATATCAACCAACATGCGGTTGATACCATCAATCGTGGCGAAATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGAtt                                                                                              >  1:1150905/1‑150 (MQ=255)
                                                      gATATCAACCAACATGCGGTTGATACCATCAATCGTGGCGAAATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCt                                                                                           >  1:421488/1‑150 (MQ=255)
                                                       atatCAACCAACATGCGGTTGATACCATCAATCGTGGCGAAATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTg                                                                                          <  1:1594626/150‑1 (MQ=255)
                                                       atatCAACCAACATGCGGTTGATACCATCAATCGTGGCGAAATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTg                                                                                          >  1:842609/1‑150 (MQ=255)
                                                       atatCAACCAACATGCGGTTGATACCATCAATCGTGGCGAAATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTg                                                                                          <  1:322365/150‑1 (MQ=255)
                                                         atCAACCAACATGCGGTTGATACCATCAATCGTGGCGAAATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTa                                                                                        <  2:293187/150‑1 (MQ=255)
                                                           caaccaacATGCGGTTGATACCATCAATCGTGGCGAAATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTAcc                                                                                      >  2:1270399/1‑150 (MQ=255)
                                                            aaccaacATGCGGTTGATACCATCAATCGTGGCGAAATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTAccc                                                                                     >  1:1558654/1‑150 (MQ=255)
                                                              ccaacATGCGGTTGATACCATCAATCGTGGCGAAATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCAc                                                                                   >  2:1129330/1‑150 (MQ=255)
                                                               caacATGCGGTTGATACCATCAATCGTGGCGAAATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACg                                                                                  >  1:902425/1‑150 (MQ=255)
                                                                 acATGCGGTTGATACCATCAATCGTGGCGAAATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGcc                                                                                <  1:1407530/150‑1 (MQ=255)
                                                                  cATGCGGTTGATACCATCAATCGTGGCGAAATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCg                                                                               <  2:842609/150‑1 (MQ=255)
                                                                    tGCGGTTGATACCATAAATCGTGGCGAAATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGTCTGATTGCTGTACCAACGCCGtt                                                                             >  1:800104/1‑150 (MQ=255)
                                                                     gCGGTTGATACCATCAATCGTGGCGAAATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGttt                                                                            <  1:1238055/150‑1 (MQ=255)
                                                                       ggTTGATACCATCAATCGTGGCGAAATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTa                                                                                        <  2:113113/136‑1 (MQ=255)
                                                                       ggTTGATACCATCAATCGTGGCGAAATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTa                                                                                        >  1:113113/1‑136 (MQ=255)
                                                                       ggTTGATACCATCAATCGTGGCGAAATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTaa                                                                          >  2:682017/1‑150 (MQ=255)
                                                                        gTTGATACGATCAATCGTGGCGAAATCCCTATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGCAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGTAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAg                                                                         <  2:625219/150‑1 (MQ=255)
                                                                        gTTGATACCATCAATCGTGGCGAAATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAg                                                                         <  2:1160910/150‑1 (MQ=255)
                                                                         ttGATACGATCAATCGTGGTGTAATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGTCGTTTTTTTAGGAGCGATCACGTCGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAgg                                                                        <  1:1053793/150‑1 (MQ=255)
                                                                          tGATACCATCAATCGTGGCGAAATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAggg                                                                       <  1:682017/150‑1 (MQ=255)
                                                                            ttACAATAAATCTTGGCTAAATCAATATCGTGTAACCTTATTTGGGTAGTGTAGTAAAAATTGCCGTAGAAGGGGGTTTTTTACGAGGGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCg                                                                     <  2:421488/149‑1 (MQ=255)
                                                                               ccATCAATCGTGGCGAAATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATc                                                                  <  1:513848/150‑1 (MQ=255)
                                                                                cATCAATCGTGGCGAAATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCa                                                                 >  2:932891/1‑150 (MQ=255)
                                                                                cATCAATCGTGGCGAAATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCa                                                                 >  2:687276/1‑150 (MQ=255)
                                                                                   cAATCGTGGCGAAATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGa                                                              >  2:766550/1‑150 (MQ=255)
                                                                                     aTCGTGGCGAAATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGc                                                            >  2:388478/1‑150 (MQ=255)
                                                                                       cGTGGCGAAATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCa                                                          >  1:1091265/1‑150 (MQ=255)
                                                                                       cGTGGCGAAATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCa                                                          <  1:1558769/150‑1 (MQ=255)
                                                                                       cGTGGCGAAATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCa                                                          <  2:1158963/150‑1 (MQ=255)
                                                                                             gAAATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATg                                                    <  2:565398/150‑1 (MQ=255)
                                                                                               aaTCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGAc                                                  <  1:594851/150‑1 (MQ=255)
                                                                                               aaTCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGAc                                                  <  2:1091265/150‑1 (MQ=255)
                                                                                                aTCCATATCGTCGATCCTGATTTGGCGAGTGTAGTAAAAACTGCCTTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGAcc                                                 <  1:388686/150‑1 (MQ=255)
                                                                                                 tCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCt                                                >  1:836648/1‑150 (MQ=255)
                                                                                                      atCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTACGtt                                           <  2:457332/150‑1 (MQ=255)
                                                                                                       tcgtcgAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTACGTTg                                          <  2:1577845/150‑1 (MQ=255)
                                                                                                       tcgtcgAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTACGTTg                                          <  1:152287/150‑1 (MQ=255)
                                                                                                       tcgtcgAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTACGTTg                                          <  1:766550/150‑1 (MQ=255)
                                                                                                            gAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTACGTTGAATCg                                     <  2:1038904/150‑1 (MQ=255)
                                                                                                              aCCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTACGTTGAATCGgc                                   >  2:88125/1‑150 (MQ=255)
                                                                                                               ccTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTACGTTGAATCGgct                                  <  2:1265508/150‑1 (MQ=255)
                                                                                                               ccTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTACGTTGAATCGgct                                  >  1:623938/1‑150 (MQ=255)
                                                                                                                    tttGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTACGTTGAATCGGCTgctcg                             >  1:340309/1‑150 (MQ=255)
                                                                                                                      tGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTACGTTGAATCGGCTgctcgct                           >  2:201458/1‑150 (MQ=255)
                                                                                                                             ggtagtaAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTACGTTGAATCGGCTGCTCGCTCCATTgc                    <  2:1569419/149‑1 (MQ=255)
                                                                                                                             ctcgtaaaaaaCCGCCTTCTACGGCAGGTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTACGTTGAATCGGCTGCTCGATCCATTgc                    >  2:389767/7‑150 (MQ=255)
                                                                                                                              gtagtaAAACGTGGCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTACGTTGAATCGGCTGCTCGCTCCATTgcg                   <  2:1416265/150‑1 (MQ=255)
                                                                                                                              gtagtaAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTACGTTGAATCGGCTGCTCGCTCCATTgcg                   <  1:488472/150‑1 (MQ=255)
                                                                                                                               tagtaAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTACGTTGAATCGGCTGCTCGCTCCATTgcgc                  >  1:691202/1‑150 (MQ=255)
                                                                                                                                agtaAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTACGTTGAATCGGCTGCTCGCTCCATTGCGcc                 <  1:1318131/150‑1 (MQ=255)
                                                                                                                                 gtaAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTACGTTGAATCGGCTGCTCGCTCCATTGCGccc                >  2:649327/1‑149 (MQ=255)
                                                                                                                                 gtaAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTACGTTGAATCGGCTGCTCGCTCCATTGCGCCa                >  2:32587/1‑150 (MQ=255)
                                                                                                                                 gtaAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGAACGACGCCAGTTGAAGAGGATGCATGGATGATTGATGTAACCACGCCGTTTAATGGCGATCATAAGCAAGATATGACCTACGTTGAATCGGCTGCTCGCTCCATTGCGCCa                >  2:1434151/1‑150 (MQ=255)
                                                                                                                                       aCTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTACGTTGAATCGGCTGCTCGCTCCATTGCGCCAGTGCTg          <  2:572566/150‑1 (MQ=255)
                                                                                                                                             gTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTACGTTGAATCGGCTGCTCGCTCCATTGCGCCAGTGCTGaaaaaa    >  1:1413226/1‑150 (MQ=255)
                                                                                                                                             gTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTACGTCGAATCGGCTGCTCGCTCCATTGCGCCAGTGCTGaaaaaa    >  1:793885/1‑150 (MQ=255)
                                                                                                                                              tAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTACGTTGAATCGGCTGCTCGCTCCATTGCGCCAGTGCTGAAAAAAg   <  1:948688/150‑1 (MQ=255)
                                                                                                                                               aGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCTTTTAAGGGCGATCATGAGCCAGATATGACCTACGTTGAATCGGCTGCTCGCTCCATTGCGCCAGTGCTGAAAAAAgg  <  1:54897/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
GGGCTGCCAACGGCAGCAGCGTTTGCCTCACGGCAAAAACAGGTAATTGGTGTCGATATCAACCAACATGCGGTTGATACCATCAATCGTGGCGAAATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTACGTTGAATCGGCTGCTCGCTCCATTGCGCCAGTGCTGAAAAAAGG  >  NC_000913/3971302‑3971594

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: