Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,971,446 G→T 5.1% E63* (GAA→TAA) ‡ 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,4460GT 5.1% 494.1 / 7.9 143E63* (GAA→TAA) ‡wecCUDP‑N‑acetyl‑D‑mannosamine dehydrogenase
Reads supporting (aligned to +/- strand):  ref base G (64/71);  new base T (2/6);  total (66/77)
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 144 total aligned reads displayed.

ATATCGGGCTGCCAACGGCAGCAGCGTTTGCCTCACGGCAAAAACAGGTAATTGGTGTCGATATCAACCAACATGCGGTTGATACCATCAATCGTGGCGAAATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTACGTTGAATCGGCTGCTCGCTCCATTGCGCCAGTGCTGAAAAAAGG  >  NC_000913/3971297‑3971594
                                                                                                                                                     |                                                                                                                                                    
atatCGGGCTGCCAACGGCAGCAGCGTTTGCCTCACGGCAAAAACAGGTAATTGGTGTCGATATCAACCAACATGCGGTTGATACCATCAATCGTGGCGAAATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAg                                                                                                                                                      >  2:898294/1‑150 (MQ=255)
 tatCGGGCTGCCAACGGCAGCAGCGTTTGCCTCACGTCAAAAACAGGTAATTGGTGTCGATATCAACCAACATGCGGTTGATACCATCAATCGTGGCGAAATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGa                                                                                                                                                     <  1:863788/150‑1 (MQ=255)
     gggCTGCTAACGGCAGCAGCGTTTGCCTCACGGCAAAAACAGGTAATTGGTGTCGATATCAACCAACATGCGGTTGATACCATCAATCGTGGCGAAATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGc                                                                                                                                                 >  2:737699/1‑150 (MQ=255)
     gggCTGCCAACGGCAGCAGCGTTTGCCTCACGGCAAAAACAGGTAATTGGTGTCGATATCAACCAACATGCGGTTGATACCATCAATCGTGGCGAAATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGc                                                                                                                                                 >  1:1449601/1‑150 (MQ=255)
      ggCTGCCAACGGCAGCAGCGTTTGCCTCACGGCAAAAACAGGTAATTGGTGTCGATATCAACCAACATGCGGTTGATACCATCAATCGTGGCGAAATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCg                                                                                                                                                >  1:457332/1‑150 (MQ=255)
        cTGCCAACGGCAGCAGCGTTTGCCTCACGGCAAAAACAGGTAATTGGTGTCGATATCAACCAACATGCGGTTGATACCATCAATCGTGGCGAAATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGt                                                                                                                                              >  1:1569419/1‑150 (MQ=255)
        cTGCCAACGGCAGCAGCGGTTGCCTCACGGCAAAAACAGGTAATTGGTGTCGATATCAACCAACATGCGGTTGATACCATCAATCGTGGCGAAATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGt                                                                                                                                              <  2:627504/150‑1 (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                                                                                                                                        <  1:326310/150‑1 (MQ=255)
              aCGGCAGCAGCGTTTGCCTCACGGCAAAAACAGGTAATTGGTGTCGATATCAACCAACATGCGGTTGATACCATCAATCGTGGCGAAATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTa                                                                                                                                        <  2:1267317/150‑1 (MQ=255)
              aCGGCAGCAGCGTTTGCCTCACGGCAAAAACAGGTAATTGGTGTCGATATCAACCAACATGCGGTTGATACCATCAATCGTGGCGAAATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACCGCCTTCTACGGCAGtttttac                                                                                                                                        <  1:353301/150‑3 (MQ=255)
                ggCAGCAGCGTTTGCCTCACGGCAAAAACAGGTAATTGGTGTCGATATCAACCAACATGCGGTTGATACCATCAATCGTGGCGAAATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACCGCCTTCTACGGCAGtttttacta                                                                                                                                      <  1:90699/150‑5 (MQ=255)
                 gcagcaGCGTTTGCCTCACGGCAAAAACAGGTAATTGGTGTCGATATCAACCAACATGCGGTTGATACCATCAATCGTGGCGAAATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTAcga                                                                                                                                     <  2:1181661/150‑1 (MQ=255)
                  cagcagCGTTTGCCTCACGGCAAAAACAGGTAATTGGTGTCGATATCAACCAACATGCGGTTGATACCATCAATCGTGGCGAAATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTAcgag                                                                                                                                    >  2:1306315/1‑150 (MQ=255)
                  cagcagCGTTTGCCTCACGGCAAAAACAGGTAATTGGTGTCGATATCAACCAACATGCGGTTGATACCATCAATCGTGGCGAAATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTAcgag                                                                                                                                    >  1:1577845/1‑150 (MQ=255)
                    gcagcGTTTGCCTCACGGCAAAAACAGGTAATTGGTGTCGATATCAACCAACATGCGGTTGATACCATCAATCGTGGCGAAATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTAcgagcg                                                                                                                                  <  1:1385772/150‑1 (MQ=255)
                       gcGTTTGCCTCACGGCAAAAACAGGTAATTGGTGTCGATATCAACCAACATGCGGTTGATACCATCAATCGTGGCGAAATCCATATCGTCGATCCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACgagcgagc                                                                                                                               >  2:388686/1‑150 (MQ=255)
                          tttGCCTCACGGCAAAAACAGGTAATTGGTGTCGATATCAACCAACATGCGGTTGATACCATCAATCGTGGCGAAATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCacg                                                                                                                            >  1:1295857/1‑150 (MQ=255)
                          tttGCCTCACGGCAAAAACAGGTAATTGGTGTCGATATCAACCAACATGCGGTTGATACCATCAATCGTGGCGAAATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCacg                                                                                                                            >  1:352394/1‑150 (MQ=255)
                             gCCTCACGGCAAAAACAGGTAATTGGTGTCGATATCAACCAACATGCGGTTGATACCATCAATCGTGGCGAAATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCacgacg                                                                                                                         >  2:956041/1‑150 (MQ=255)
                              ccTCACGGCAAAAACAGGTAATTGGTGTCGATATCAACCAACATGCGGTTGATACCATCAATCGTGGCGAAATCAATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGc                                                                                                                        >  2:242911/1‑150 (MQ=255)
                               cTCACGGCAAAAACAGGTAATTGGTGTCGATATCAACCAACATGCGGTTGATACCATCAATCGTGGCGAAATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGcc                                                                                                                       >  1:202259/1‑150 (MQ=255)
                               cTCACGGCAAAAACAGGTAATTGGTGTCGATATCAACCAACATGCGGTTGATACCATCAATCGTGGCGAAATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACCGCCTTCTACGGCAGtttttac                                                                                                                                        >  2:1329088/1‑131 (MQ=255)
                               cTCACGGCAAAAACAGGTAATTGGTGTCGATATCAACCAACATGCGGTTGATACCATCAATCGTGGCGAAATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACCGCCTTCTACGGCAGtttttac                                                                                                                                        <  1:1329088/133‑3 (MQ=255)
                                tCACGGCAAAAACAGGTAATTGGTGTCGATATCAACCAACATGCGGTTGATACCATCAATCGTGGCGAAATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCa                                                                                                                      >  1:1427337/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)
                                  aCGGCAAAAACAGGTAATTGGTGTCGATATCAACCAACATGCGGTTGATACCATCAATCGTGGCGAAATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGt                                                                                                                    <  2:1025866/150‑1 (MQ=255)
                                    ggCAAAAACAGGTAATTGGTGTCGATATCAACCAACATGCGGTTGATACCATCAATCGTGGCGAAATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTg                                                                                                                  <  2:274614/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)
                                                 ttttGGTGTCGATTTCAACCAACATGCGGTTTATACCATCATTCTTGGCGAAATCCATTTCTTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACTAGCTAGCACGACGCCAGTTGAAGCGGATGCCTg                                                                                                     <  2:156893/148‑1 (MQ=255)
                                                 aaTTGGTGTCGATATCAACCAACATGCGGTTGATACCATCAATCGTGGCGAAATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTg                                                                                                     <  2:687615/150‑1 (MQ=255)
                                                   ttGGTGTCGATATCAACCAACATGCGGTTGATACCATCAATCGTGGCGAAATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGc                                                                                                   <  2:561452/150‑1 (MQ=255)
                                                     ggTGTCGATATCAACCAACATGCGGTTGATACCATCAATCGTGGCGAAATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTg                                                                                                 <  2:345445/150‑1 (MQ=255)
                                                     ggTGTCGATATCAACCAACATGCGGTTGATACCATCAATCGTGGCGAAATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTg                                                                                                 <  2:339391/150‑1 (MQ=255)
                                                     ggTGTCAATATCAACCAACATGCGGTTGATACCATCAATCGTGGCGAAATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTg                                                                                                 >  1:681765/1‑150 (MQ=255)
                                                        gtCGATATCAACCAACATGCGGTTGATACCATCAATCGTGGCGAAATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGAtt                                                                                              >  1:1150905/1‑150 (MQ=255)
                                                           gATATCAACCAACATGCGGTTGATACCATCAATCGTGGCGAAATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCt                                                                                           >  1:421488/1‑150 (MQ=255)
                                                            atatCAACCAACATGCGGTTGATACCATCAATCGTGGCGAAATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTg                                                                                          <  1:322365/150‑1 (MQ=255)
                                                            atatCAACCAACATGCGGTTGATACCATCAATCGTGGCGAAATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTg                                                                                          >  1:842609/1‑150 (MQ=255)
                                                               tCATCCAAAATGCGGTTGATAGCATCAATCGTGGCGAAATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTTGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGCTGAAGCGGATGCCTGGCTGATTGCTGTAc                                                                                       <  2:845278/150‑1 (MQ=255)
                                                                 aaccaacATGCGGTTGATACCATCAATCGTGGCGAAATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTAccc                                                                                     >  1:1558654/1‑150 (MQ=255)
                                                                 aacaaaaATGCGCTTTATACCATCAAAAGTCGCGAAATCCATATCGTCTAACCTGATTTGGCGAGTGTAGTAAAAACTACCGTAGAAGGCGGTTTTTAAAGAGCGAGAACGACGACACATGAAGAGGATGCCTGGCTGATTGCTGTAccc                                                                                     >  2:329169/1‑150 (MQ=255)
                                                                    caacATGCGGTTGATACCATCAATCGTGGCGAAATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACg                                                                                  >  1:902425/1‑150 (MQ=255)
                                                                    caacATGCGGTTCATACCATCAATCGTGGCGAAATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCTGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACg                                                                                  <  1:1223658/150‑1 (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)
                                                                           cGGTTGATACCATCAATCGTGGCGAAATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCTTTTTTTTACGAGCTAGCACGACGCCATTTGAATCTGATGCCTTGCTGATTGCTGTACCCACGGCGTTTa                                                                           <  2:438759/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)
                                                                            ggTTGATACCATCAATCGTGGCGAAATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTaa                                                                          >  1:572566/1‑150 (MQ=255)
                                                                             gTTGATACCATCAATCGTGGCGAAATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAg                                                                         <  2:1160910/150‑1 (MQ=255)
                                                                              ttGATACGATCAATCGTGGTGTAATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGTCGTTTTTTTAGGAGCGATCACGTCGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAgg                                                                        <  1:1053793/150‑1 (MQ=255)
                                                                                 aTACCATCAATCGTGGCGAAATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCg                                                                     >  1:171774/1‑150 (MQ=255)
                                                                                 aTACCATCAATCGTGGCGAAATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCg                                                                     <  2:681765/150‑1 (MQ=255)
                                                                                  tACCATCAATCGTGGCGAAATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGa                                                                    >  1:485405/1‑150 (MQ=255)
                                                                                    ccATCAATCGTGGCGAAATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATc                                                                  <  1:513848/150‑1 (MQ=255)
                                                                                    ccATCAATCGTGGCGAAATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATc                                                                  >  2:297521/1‑150 (MQ=255)
                                                                                     cATCAATCGTGGCGAAATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCa                                                                 >  2:932891/1‑150 (MQ=255)
                                                                                        cAATCGTGGCGAAATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGa                                                              >  2:766550/1‑150 (MQ=255)
                                                                                            cGTGGCGAAATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCa                                                          >  1:1091265/1‑150 (MQ=255)
                                                                                               ggCGAAATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGat                                                       <  2:749544/150‑1 (MQ=255)
                                                                                                  gAAATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATg                                                    <  2:565398/150‑1 (MQ=255)
                                                                                                    aaTCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGAc                                                  <  2:1091265/150‑1 (MQ=255)
                                                                                                           atCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTACGtt                                           <  2:457332/150‑1 (MQ=255)
                                                                                                            tcgtcgAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTACGTTg                                          <  2:1577845/150‑1 (MQ=255)
                                                                                                            tcgtcgAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTACGTTg                                          <  1:766550/150‑1 (MQ=255)
                                                                                                            tcgtcgAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTACGTTg                                          <  1:152287/150‑1 (MQ=255)
                                                                                                                   aCCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTACGTTGAATCGgc                                   <  2:700845/150‑1 (MQ=255)
                                                                                                                   aCCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTACGTTGAATCGgc                                   >  2:88125/1‑150 (MQ=255)
                                                                                                                    ccTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTACGTTGAATCGgct                                  >  1:623938/1‑150 (MQ=255)
                                                                                                                         tttGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTACGTTGAATCGGCTgctcg                             >  1:340309/1‑150 (MQ=255)
                                                                                                                           tGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTACGTTGAATCGGCTgctcgct                           >  1:1548299/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)
                                                                                                                                    tagtaAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTACGTTGAATCGGCTGCTCGCTCCATTgcgc                  >  1:691202/1‑150 (MQ=255)
                                                                                                                                     agtaAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTACGTTGAATCGGCTGCTCGCTCCATTGCGcc                 <  1:1318131/150‑1 (MQ=255)
                                                                                                                                      gtaAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTACGTTGAATCGGCTGCTCGCTCCATTGCGccc                >  2:649327/1‑149 (MQ=255)
                                                                                                                                      gtaAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGAACGACGCCAGTTGAAGAGGATGCATGGATGATTGATGTAACCACGCCGTTTAATGGCGATCATAAGCAAGATATGACCTACGTTGAATCGGCTGCTCGCTCCATTGCGCCa                >  2:1434151/1‑150 (MQ=255)
                                                                                                                                        aaaaaCTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTACGTTGAATCGGCTGCTCGCTCCATTGCGCCAGt              >  2:1267433/1‑150 (MQ=255)
                                                                                                                                            aCTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTACGTTGAATCGGCTGCTCGCTCCATTGCGCCAGTGCTg          <  2:572566/150‑1 (MQ=255)
                                                                                                                                              tGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTACGTTGAATCGGCTGCTCGCTCCATTGCGCCAGTGCTGaa        >  2:1103956/1‑150 (MQ=255)
                                                                                                                                                ccGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGTGATCATGAGCCAGATATGACCTACGTTGAATCGGCTGCTCGCTCCATTGCGCCAGTGCTGaaaa      <  2:1199362/150‑1 (MQ=255)
                                                                                                                                                 cGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTACGTTGAATCGGCTGCTCGCTCCATTGCGCCAGTGCTGaaaaa     >  2:1149436/1‑150 (MQ=255)
                                                                                                                                                 cGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTACGTTGAATCGGCTGCTCGCTCCATTGCGCCAGTGCTGaaaaa     >  1:831974/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)
                                                                                                                                                     |                                                                                                                                                    
ATATCGGGCTGCCAACGGCAGCAGCGTTTGCCTCACGGCAAAAACAGGTAATTGGTGTCGATATCAACCAACATGCGGTTGATACCATCAATCGTGGCGAAATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTACGTTGAATCGGCTGCTCGCTCCATTGCGCCAGTGCTGAAAAAAGG  >  NC_000913/3971297‑3971594

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: