Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,971,439 T→C 5.6% T60T (ACT→ACC 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,4390TC 5.6% 437.0 / 8.9 140T60T (ACT→ACCwecCUDP‑N‑acetyl‑D‑mannosamine dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (60/72);  new base C (3/5);  total (63/77)
Fisher's exact test for biased strand distribution p-value = 7.30e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 142 total aligned reads displayed.

CTGGGTTATATCGGGCTGCCAACGGCAGCAGCGTTTGCCTCACGGCAAAAACAGGTAATTGGTGTCGATATCAACCAACATGCGGTTGATACCATCAATCGTGGCGAAATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTACGTTGAATCGGCTGCTCGCTCCATTGCGCCAGTGCTGAA  >  NC_000913/3971290‑3971588
                                                                                                                                                     |                                                                                                                                                     
cTGGGTTATATCGGGCTGCCAACGGCAGCAGCGTTTGCCTCACGGCAAAAACAGGTAATTGGTGTCGATATCAACCAACATGCGGTTGATACCATCAATCGTGGCGAAATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACt                                                                                                                                                       >  2:20423/1‑150 (MQ=255)
  gggTTATATCGGGCTGCCAACGGCAGCAGCGTTTGCCTCACGGCAAAAACAGGTAATTGGTGTCGATATCAACCAACATGCGGTTGATACCATCAATCGTGGCGAAATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACCGc                                                                                                                                                     >  1:1182937/1‑150 (MQ=255)
      tataTCGGGCTGCCAACGGCAGCAGCGTTTGCCTCACGGCAAAAACAGGTAATTGGTGTCGATATCAACCAACATGCGGTTGATACCATCAATCGTGGCGAAATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTa                                                                                                                                                 <  2:842504/150‑1 (MQ=255)
      tataTCGGGCTGCCAACGGCAGCAGCGTTTGCCTCACGGCAAAAACAGGTAATTGGTGTCGATATCAACCAACATGCGGTTGATACCATCAATCGTGGCGAAATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTa                                                                                                                                                 <  2:389592/150‑1 (MQ=255)
      tataTCGGGCTGCCAACGGCAGCAGCGTTTGCCTCACGGCAAAAACAGGTAATTGGTGTCGATATCAACCAACATGCGGTTGATACCATCAATCGTGGCGAAATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTa                                                                                                                                                 <  2:1549510/150‑1 (MQ=255)
       atatCGGGCTGCCAACGGCAGCAGCGTTTGCCTCACGGCAAAAACAGGTAATTGGTGTCGATATCAACCAACATGCGGTTGATACCATCAATCGTGGCGAAATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAg                                                                                                                                                >  2:898294/1‑150 (MQ=255)
        tatCGGGCTGCCAACGGCAGCAGCGTTTGCCTCACGTCAAAAACAGGTAATTGGTGTCGATATCAACCAACATGCGGTTGATACCATCAATCGTGGCGAAATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGa                                                                                                                                               <  1:863788/150‑1 (MQ=255)
        tatCGGGCTGCCAACGGCAGCAGCGTTTGCCTCACGGCAAAAACAGGTAATTGGTGTCGATATCAACCAACATGCGGTTGATACCATCAATCGTGGCGAAATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGa                                                                                                                                               >  1:869616/1‑150 (MQ=255)
            gggCTGCCAACGGCAGCAGCGTTTGCCTCACGGCAAAAACAGGTAATTGGTGTCGATATCAACCAACATGCGGTTGATACCATCAATCGTGGCGAAATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGc                                                                                                                                           >  1:1449601/1‑150 (MQ=255)
             ggCTGCCAACGGCAGCAGCGTTTGCCTCACGGCAAAAACAGGTAATTGGTGTCGATATCAACCAACATGCGGTTGATACCATCAATCGTGGCGAAATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCg                                                                                                                                          >  1:457332/1‑150 (MQ=255)
               cTGCCAACGGCAGCAGCGGTTGCCTCACGGCAAAAACAGGTAATTGGTGTCGATATCAACCAACATGCGGTTGATACCATCAATCGTGGCGAAATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGt                                                                                                                                        <  2:627504/150‑1 (MQ=255)
                tGCCAACGGCAGCAGCGTTTGCCTCACGGCAAAAACAGGTAATTGGTGTCGATATCAACCAACATGCGGTTGATACCATCAATCGTGGCGAAATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGtt                                                                                                                                       >  1:477293/1‑150 (MQ=255)
                 gCCAACGGCAGCAGCGTTTGCCTCACGGCAAAAACAGGTAATTGGTGTCGATATCAACCAACATGCGGTTGATACCATCAATCGTGGCGAAATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGttt                                                                                                                                      >  2:1385772/1‑150 (MQ=255)
                 gCCAACGGCAGCAGCGTTTGCCTCACGGCAAAAACAGGTAATTGGTGTCGATATCAACCAACATGCGGTTGATACCATCAATCGTGGCGAAATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGttt                                                                                                                                      >  2:1076179/1‑150 (MQ=255)
                 gCCAACGGCAGCAGCGTTTGCCTCACGGCAAAAACAGGTAATTGGTGTCGATATCAACCAACATGCGGTTGATACCATCAATCGTGGCGAAATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGttt                                                                                                                                      >  2:322365/1‑150 (MQ=255)
                   cAACGGCAGCAGCGTTTGCCTCACGGCAAAAACAGGTAATTGGTGTCGATATCAACCAACATGCGGTTGATACCATCAATCGTGGCGAAATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGttttt                                                                                                                                    >  1:700845/1‑150 (MQ=255)
                    aaCGGCAGCAGCGTTTGCCTCACGGCAAAAACAGGTAATTGGTGTCGATATCAACCAACATGCGGTTGATACCATCAATCGTGGCGAAATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGtttttt                                                                                                                                   >  1:1038904/1‑150 (MQ=255)
                     aCGGCAGCAGCGTTTGCCTCACGGCAAAAACAGGTAATTGGTGTCGATATCAACCAACATGCGGTTGATACCATCAATCGTGGCGAAATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTa                                                                                                                                  <  2:1267317/150‑1 (MQ=255)
                     aCGGCAGCAGCGTTTGCCTCACGGCAAAAACAGGTAATTGGTGTCGATATCAACCAACATGCGGTTGATACCATCAATCGTGGCGAAATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTa                                                                                                                                  <  1:326310/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)
                        gcagcaGCGTTTGCCTCACGGCAAAAACAGGTAATTGGTGTCGATATCAACCAACATGCGGTTGATACCATCAATCGTGGCGAAATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACCGCCTTCTACGGCAGtttttactac                                                                                                                               <  2:1182937/150‑6 (MQ=255)
                         cagcagCGTTTGCCTCACGGCAAAAACAGGTAATTGGTGTCGATATCAACCAACATGCGGTTGATACCATCAATCGTGGCGAAATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTAcgag                                                                                                                              >  1:1577845/1‑150 (MQ=255)
                         cagcagCGTTTGCCTCACGGCAAAAACAGGTAATTGGTGTCGATATCAACCAACATGCGGTTGATACCATCAATCGTGGCGAAATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTAcgag                                                                                                                              >  2:1306315/1‑150 (MQ=255)
                           gcagcGTTTGCCTCACGGCAAAAACAGGTAATTGGTGTCGATATCAACCAACATGCGGTTGATACCATCAATCGTGGCGAAATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTAcgagcg                                                                                                                            <  1:1385772/150‑1 (MQ=255)
                                    gCCTCACGGCAAAAACAGGTAATTGGTGTCGATATCAACCAACATGCGGTTGATACCATCAATCGTGGCGAAATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCacgacg                                                                                                                   >  2:956041/1‑150 (MQ=255)
                                      cTCACGGCAAAAACAGGTAATTGGTGTCGATATCAACCAACATGCGGTTGATACCATCAATCGTGGCGAAATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGcc                                                                                                                 >  1:202259/1‑150 (MQ=255)
                                        cACGGCAAAAACAGGTAATTGGTGTCGATATCAACCAACATGCGGTTGATACCATCAATCGTGGCGAAATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAg                                                                                                               >  1:749544/1‑150 (MQ=255)
                                         aCGGCAAAAACAGGTAATTGGTGTCGATATCAACCAACATGCGGTTGATACCATCAATCGTGGCGAAATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGcc                                                                                                                 >  1:1195920/1‑147 (MQ=255)
                                         aCGGCAAAAACAGGTAATTGGTGTCGATATCAACCAACATGCGGTTGATACCATCAATCGTGGCGAAATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGcc                                                                                                                 <  2:1195920/147‑1 (MQ=255)
                                         aCGGCAAAAACAGGTAATTGGTGTCGATATCAACCAACATGCGGTTGATACCATCAATCGTGGCGAAATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGt                                                                                                              <  2:1025866/150‑1 (MQ=255)
                                          cGGCAAAAACAGGTAATTGGTGTCGATATCAACCAACATGCGGTTGATACCATCAATCGTGGCGAAATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGtt                                                                                                             >  2:1457938/1‑150 (MQ=255)
                                           ggCAAAAACAGGTAATTGGTGTCGATATCAACCAACATGCGGTTGATACCATCAATCGTGGCGAAATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTg                                                                                                            <  2:274614/150‑1 (MQ=255)
                                                aaaCAGGTAATTGGTGTCGATATCAACCAACATGCGGTTGATACCATCAATCGTGGCGAAATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCg                                                                                                       <  1:1510869/150‑1 (MQ=255)
                                                aaaCAGGTAATTGGTGTCGATATCAACCAACATGCGGTTGATACCATCAATCGTGGCGAAATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCg                                                                                                       <  1:1362000/150‑1 (MQ=255)
                                                  aCAGTTATTTGGTGTGGATTTAAACCAAGATGCGGTTTTTACGATCTATCGTGGGGAAATCCATTTCGTGGAACCTTATTTGGCGAGTGTTGTAAAAACTGCCTTATACGGCTGTTTTTTACTAGCGAGCACGACGCCAGTTGAAGCGGa                                                                                                     <  2:202259/150‑1 (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)
                                                          ttGGTGTCGATATCAACCAACATGCGGTTGATACCATCAATCGTGGCGAAATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGc                                                                                             <  2:561452/150‑1 (MQ=255)
                                                          ttGGTGTCGATATCAACCAACATGCGGTTGATACCATCAATCGAGGCGAAATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGTCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGc                                                                                             <  1:331717/150‑1 (MQ=255)
                                                            ggTGTCGATATCAACCAACATGCGGTTGATACCATCAATCGTGGCGAAATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTg                                                                                           <  2:345445/150‑1 (MQ=255)
                                                            ggTGTCGATATCAACCAACATGCGGTTGATACCATCAATCGTGGCGAAATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTg                                                                                           <  2:339391/150‑1 (MQ=255)
                                                               gtCGATATCAACCAACATGCGGTTGATACCATCAATCGTGGCGAAATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGAtt                                                                                        >  1:1150905/1‑150 (MQ=255)
                                                                   atatCAACCAACATGCGGTTGATACCATCAATCGTGGCGAAATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTg                                                                                    <  2:1591192/150‑1 (MQ=255)
                                                                   atatCAACCAACATGCGGTTGATACCATCAATCGTGGCGAAATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTg                                                                                    <  1:1594626/150‑1 (MQ=255)
                                                                     atCAACCAACATGCGGTTGATACCATCAATCGTGGCGAAATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTa                                                                                  <  2:293187/150‑1 (MQ=255)
                                                                      tCATCCAAAATGCGGTTGATAGCATCAATCGTGGCGAAATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTTGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGCTGAAGCGGATGCCTGGCTGATTGCTGTAc                                                                                 <  2:845278/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)
                                                                           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)
                                                                                     ttGATACGATCAATCGTGGTGTAATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGTCGTTTTTTTAGGAGCGATCACGTCGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAgg                                                                  <  1:1053793/150‑1 (MQ=255)
                                                                                      tGATACCATCAATCGTGGCGAAATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAggg                                                                 <  1:682017/150‑1 (MQ=255)
                                                                                        aTACCATCAATCGTGGCGAAATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCg                                                               <  2:240218/150‑1 (MQ=255)
                                                                                        aTACCATCAATCGTGGCGAAATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCg                                                               <  2:681765/150‑1 (MQ=255)
                                                                                         tACCATCAATCGTGGCGAAATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGa                                                              >  1:485405/1‑150 (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                                                    <  2:1158963/150‑1 (MQ=255)
                                                                                                   cGTGGCGAAATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCa                                                    >  1:1091265/1‑150 (MQ=255)
                                                                                                   cGTGGCGAAATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCa                                                    <  1:1558769/150‑1 (MQ=255)
                                                                                                      ggCGAAATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGat                                                 <  2:749544/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:766550/150‑1 (MQ=255)
                                                                                                                          aCCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTACGTTGAATCGgc                             <  2:700845/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)
                                                                                                                               aTTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTACGTTGAATCGGCTgctc                        <  1:1306315/150‑1 (MQ=255)
                                                                                                                                tttGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTACGTTGAATCGGCTgctcg                       >  1:340309/1‑150 (MQ=255)
                                                                                                                                  tGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTACGTTGAATCGGCTgctcgct                     >  1:1548299/1‑150 (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)
                                                                                                                                                   aCTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTACGTTGAATCGGCTGCTCGCTCCATTGCGCCAGTGCTg    <  2:572566/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CTGGGTTATATCGGGCTGCCAACGGCAGCAGCGTTTGCCTCACGGCAAAAACAGGTAATTGGTGTCGATATCAACCAACATGCGGTTGATACCATCAATCGTGGCGAAATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTACGTTGAATCGGCTGCTCGCTCCATTGCGCCAGTGCTGAA  >  NC_000913/3971290‑3971588

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: