Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,921,853 C→T 5.4% G83D (GGT→GAT) ‡ gudP ← galactarate/glucarate/glycerate transporter GudP

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,921,8530CT 5.4% 473.7 / 9.7 147G83D (GGT→GAT) ‡gudPgalactarate/glucarate/glycerate transporter GudP
Reads supporting (aligned to +/- strand):  ref base C (69/70);  new base T (4/4);  total (73/74)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.23e-01

Only 100 of 151 total aligned reads displayed.

GGAAAGATGGCGCTTCAGCAAGCCCGACCAGGAAGCGCAGCGTAAACAGGGCAACGATAATGCCGAATCCACTAAAGATATCGACGAAGCCTTGCAGCAAGGTAAACATCGACCAGATAAAGATCGACCAGAAGTAGACGCGTTTTGAACCAAAACGGTCCAGCAACCAGCCACCAGGGATCTGCCCGATAACATAAGCCCATGAGAAAGCAGAGAACACATAGCCCATTCCCACGGGATCAAGGCCGATATCTTTGGCCATTTCCGAACCGGCGATAGAGAGCGTAGCGCGGTCGCC  >  NC_000913/2921704‑2922001
                                                                                                                                                     |                                                                                                                                                    
ggAAAGATGGCGCTTCAGCAAGCCCGACCAGGAAGCGCAGCGTAAACAGGGCAACGATAATGCCGAATCCACTAAAGATATCGACGAAGCCTTGCAGCAAGGTAAACATCGACCAGATAAAGATCGACCAGAAGTAGACGCGTTTTGAAc                                                                                                                                                      <  1:1256685/150‑1 (MQ=255)
ggAAAGATGGCGCTTCAGCAAGCCCGACCAGGAAGCGCAGCGTAAACAGGGCAACGATAATGCCGAATCCACTAAAGATATCGACGAAGCCTTGCAGCAAGGTAAACATCGACCAGATAAAGATCGACCAGAAGTAGACGCGTTTTGAAc                                                                                                                                                      <  2:2058144/150‑1 (MQ=255)
ggAAAGATGGCGCTTCAGCAAGCCCGACCAGGAAGCGCAGCGTAAACAGGGCAACGATAATGCCGAATCCACTAAAGATATCGACGAAGCCTTGCAGCAAGGTAAACATCGACCAGATAAAGATCGACCAGAAGTAGACGCGTTTTGAAc                                                                                                                                                      <  2:17100/150‑1 (MQ=255)
 gAAAGATGGCGCTTCAGCAAGCCCGACCAGGAAGCGCAGCGTAAACAGGGCAACGATAATGCCGAATCCACTAAAGATATCGACGAAGCCTTGCAGCAAGGTAAACATCGACCAGATAAAGATCGACCAGAAGTAGACGCGTTTTGAAcc                                                                                                                                                     <  2:1867015/150‑1 (MQ=255)
   aaGATGGCGCTTCAGCAAGCCCGACCAGGAAGCGCAGCGTAAACAGGGCAACGATAATGCCGAATCCACTAAAGATATCGACGAAGCCTTGCAGCAAGGTAAACATCGACCAGATAAAGATCGACCAGAAGTAGACGCGTTTTGAACCaa                                                                                                                                                   <  1:2318596/150‑1 (MQ=255)
     gATGGCGCTTCAGCAAGCCCGACCAGGAAGCGCAGCGTAAACAGGGCAACGATAATGCCGAATCCACTAAAGATATCGACGAAGCCTTGCAGCAAGGTAAACATCGACCAGATAAAGATCGACCAGAAGTAGACGCGTTTTGAACCaaaa                                                                                                                                                 >  1:1129240/1‑150 (MQ=255)
     gATGGCGCTTCAGCAAGCCCGACCAGGAAGCGCAGCGTAAACAGGGCAACGATAATGCCGAATCCACTAAAGATATCGACGAAGCCTTGCAGCAAGGTAAACATCGACCAGATAAAGATCGACCAGAAGTAGACGCGTTTTGAACCaaaa                                                                                                                                                 >  1:1470012/1‑150 (MQ=255)
      aTGGCGCTTCAGCAAGCCCGACCAGGAAGCGCAGCGTAAACAGGGCAACGATAATGCCGAATCCACTAAAGATATCGACGAAGCCTTGCAGCAAGGTAAACATCGACCAGATAAAGATCGACCAGAAGTAGACGCGTTTTGAACCAAAAc                                                                                                                                                >  1:1212184/1‑150 (MQ=255)
        ggCGCTTCAGCAAGCCCGACCAGGAAGCGCAGCGTAAACAGGGCAACGATAATGCCGAATCCACTTAAGATATCGACGAAGCCTTGCAGCAAGGTAAACATCGACCAGATAAAGATCGACCAGAAGTAGACGCGTTTTGAACCAAAACgg                                                                                                                                              >  2:2388743/1‑150 (MQ=255)
        ggCGCTTCAGCAAGCCCGACCAGGAAGCGCAGCGTAAACAGGGCAACGATAATGCCGAATCCACTAAAGATATCGACGAAGCCTTGCAGCAAGGTAAACATCGACCAGATAAAGATCGACCAGAAGTAGACGCGTTTTGAACCAAAACgg                                                                                                                                              <  1:2407919/150‑1 (MQ=255)
         gcgcTTCAGCAAGCCCGACCAGGAAGCGCAGCGTAAACAGGGCAACGATAATGCCGAATCCACTAAAGATATCGACGAAGCCTTGCAGCAAGGTAAACATCGACCAGATAAAGATCGACCAGAAGTAGACGCGTTTTGAACCAAAACGGt                                                                                                                                             <  2:1782314/150‑1 (MQ=255)
         gcgcTTCAGCAAGCCCGACCAGGAAGCGCAGCGTAAACAGGGCAACGATAATGCCGAATCCACTAAAGATATCGACGAAGCCTTGCAGCAAGGTAAACATCGACCAGATAAAGATCGACCAGAAGTAGACGCGTTTTGAACCAAAACGGt                                                                                                                                             <  2:1189111/150‑1 (MQ=255)
         gcgcTTCAGCAAGCCCGACCAGGAAGCGCAGCGTAAACAGGGCAACGATAATGCCGAATCCACTAAAGATATCGACGAAGCCTTGCAGCAAGGTAAACATCGACCAGATAAAGATCGACCAGAAGTAGACGCGTTTTGAACCAAAACGGt                                                                                                                                             <  2:669705/150‑1 (MQ=255)
          cgcTTCAGCAAGCCCGACCAGGAAGCGCAGCGTAAACAGGGCAACGATAATGCCGAATCCACTAAAGATATCGACGAAGCCTTGCAGCAAGGTAAACATCGACCAGATAAAGATCGACCAGAAGTAGACGCGTTTTGGTTCAAAACgcgt                                                                                                                                            <  2:823145/150‑4 (MQ=255)
          cgcTTCAGCAAGCCCGACCAGGAAGCGCAGCGTAAACAGGGCAACGATAATGCCGAATCCACTAAAGATATCGACGAAGCCTTGCAGCAAGGTAAACATCGACCAGATAAAGATCGACCAGAAGTAGACGCGTTTTGAACCAAAACGGTc                                                                                                                                            <  2:677625/150‑1 (MQ=255)
           gcTTCAGCAAGCCCGACCAGGAAGCGCAGCGTAAACAGGGCAACGATAATGCCGAATCCACTAAAGATATCGACGAAGCCTTGCAGCAAGGTAAACATCGACCAGATAAAGATCGACCAGAAGTAGACGCGTTTTGGTTCAAAACgcgtc                                                                                                                                           <  1:2043299/150‑5 (MQ=255)
           gcTTCAGCAAGCCCGACCAGGAAGCGCAGCGTAAACAGGGCAACGATAATGCCGAATCCACTAAAGATATCGACGAAGCCTTGCAGCAAGGTAAACATCGACCAGATAAAGATCGACCAGAAGTAGACGCGTTTTGAACCAAAACGGTcc                                                                                                                                           <  1:1390911/150‑1 (MQ=255)
           gcTTCAGCAAGCCCGACCAGGAAGCGCAGCGTAAACAGGGCAACGATAATGCCGAATCCACTAAAGATATCGACGAAGCCTTGCAGCAAGGTAAACATCGACCAGATAAAGATCGACCAGAAGTAGACGCGTTTTGAACCAAAACGGTcc                                                                                                                                           >  2:575910/1‑150 (MQ=255)
            cTTCAGCAAGCCCGACCAGGAAGCGCAGCGTAAACAGGGCAACGATAATGCCGAATCCACTAAAGATATCGACGAAGCCTTGCAGCAAGGTAAACATCGACCAGATAAAGATCGACCAGAAGTAGACGCGTTTTGAACCAAAACGGTCCa                                                                                                                                          >  2:2103408/1‑150 (MQ=255)
            cTTCAGCAAGCCCGACCAGGAAGCGCAGCGTAAACAGGGCAACGATAATGCCGAATCCACTAAAGATATCGACGAAGCCTTGCAGCAAGGTAAACATCGACCAGATAAAGATCGACCAGAAGTAGACGCGTTTTGAACCAAAACGGTCCa                                                                                                                                          <  2:307214/150‑1 (MQ=255)
             ttCAGCAAGCCCGACCAGGAAGCGCAGCGTAAACAGGGCAACGATAATGCCGAATCCACTAAAGATATCGACGAAGCCTTGCAGCAAGGTAAACATCGACCAGATAAAGATCGACCAGAAGTAGACGCGTTTTGAACCAAAACGGTCCAg                                                                                                                                         <  1:1548252/150‑1 (MQ=255)
                 gCAAGCCCGACCAGGAAGCGCAGCGTAAACAGGGCAACGATAATGCCGAATCCACTAAAGATATCGACGAAGCCTTGCAGCAAGGTAAACATCGACCAGATAAAGATCGACCAGAAGTAGACGCGTTTTGAACCAAAACGGTCCAGCAAc                                                                                                                                     >  2:297637/1‑150 (MQ=255)
                  cAAGCCCGACCAGGAAGCGCAGCGTAAACAGGGCAACGATAATGCCGAATCCACTAAAGATATCGACGAAGCCTTGCAGCAAGGTAAACATCGACCAGATAAAGATCGACCAGAAGTAGACGCGTTTTGAACCAAAACGGTCCAGCAAcc                                                                                                                                    <  2:1958100/150‑1 (MQ=255)
                      cccGACCAGGAAGCGCAGCGTAAACAGGGCAACGATAATGCCGAATCCACTAAAGATATCGACGAAGCCTTGCAGCAAGGTAAACATCGACCAGATAAAGATCGACCAGAAGTAGACGCGTTTTGAACCAAAACGGTCCAGCAACCAGcc                                                                                                                                >  2:2052487/1‑150 (MQ=255)
                       ccGACCAGGAAGCGCAGCGTAAACAGGGCAACGATAATGCCGAATCCACTAAAGATATCGACGAAGCCTTGCAGCAAGGTAAACATCGACCAGATAAAGATCGACCAGAAGTAGACGCGTTTTGAACCAAAACGGTCCAGCAACCAGcca                                                                                                                               >  1:476972/1‑150 (MQ=255)
                           ccAGGAAGCGCAGCGTAAACAGGGCAACGATAATGCCGAATCCACTAAAGATATCGACGAAGCCTTGCAGCAAGGTAAACATCGACCAGATAAAGATCGACCAGAAGTAGACGCGTTTTGAACCAAAACGGTCCAGCAACCAGCCACCAg                                                                                                                           <  2:963783/150‑1 (MQ=255)
                            cAGGAAGCGCAGCGTAAACAGGGCAACGATAATGCCGAATCCACTAAAGATATCGACGAAGCCTTGCAGCAAGGTAAACATCGACCAGATAAAGATCGACCAGAAGTAGACGCGTTTTGAACCAAAACGGTCCAGCAACCAGCCACCAgg                                                                                                                          >  1:905698/1‑150 (MQ=255)
                              tgAATCGCAGAGTAAACAGTGCAAAGATAATGCCTAATCCACTAAAGATATCGACTAAGCCTTGCAGCAAGGTAAACATCGTACAGATAAAGATAGTCCAGAAGTAGACGCGTTTTGATCCAAAACGGTCCAGCAACCAGCCACCAGGGa                                                                                                                        <  1:2225022/149‑1 (MQ=255)
                              ggAAGCGCAGCGTAAACAGGGCAACGATAATGCCGAATCCACTAAAGATATCGACGAAGCCTTGCAGCAAGGTAAACATCGACCAGATAAAGATCGACCAGAAGTAGACGCGTTTTGAACCAAAACGGTCCAGCAACCAGCCACCAGGGa                                                                                                                        >  1:27685/1‑150 (MQ=255)
                              ggAAGCGCAGCGTAAACAGGGCAACGATAATGCCGAATCCACTAAAGATATCGACGAAGCCTTGCAGCAAGGTAAACATCGACCAGATAAAGATCGACCAGAAGTAGACGCGTTTTGAACCAAAACGGTCCAGCAACCAGCCACCAGGGa                                                                                                                        >  2:2207449/1‑150 (MQ=255)
                                       gCGTAAACAGGGCAACGATAATGCCGAATCCACTAAAGATATCGACGAAGCCTTGCAGCAAGGTAAACATCGACCAGATAAAGATCGACCAGAAGTAGACGCGTTTTGAACCAAAACGGTCCAGCAACCAGCCACCAGGGATCTGCCCGa                                                                                                               >  1:2243466/1‑150 (MQ=255)
                                       gCGTAAACAGGGCAACGATAATGCCGAATCCACTAAAGATATCGACGAAGCCTTGCAGCAAGGTAAACATCGACCAGATAAAGATCGACCAGAAGTAGACGCGTTTTGAACCAAAACGGTCCAGCAACCAGCCACCAGGGATCTGCCCGa                                                                                                               >  1:45899/1‑150 (MQ=255)
                                        cGTAAACAGGGCAACGATAATGCCGAATCCACTAAAGATATCGACGAAGCCTTGCAGCAAGGTAAACATCGACCAGATAAAGATCGACCAGAAGTAGACGCGTTTTGAACCAAAACGGTCCAGCAACCAGCCACCAGGGATCTGCCCGAt                                                                                                              <  2:2128696/150‑1 (MQ=255)
                                          tAAACAGGGCAACGATAATGCCGAATCCACTAAAGATATCGACGAAGCCTTGCAGCAAGGTAAACATCGCCCAGATAAAGATCGACCAGAAGTAGACGCGTTTTGAACCAAAACGGTCCAGCAACCAGCCCCCAGGGAACTGCCCGATaa                                                                                                            >  2:2463707/1‑150 (MQ=255)
                                          tAAACAGGGCAACGATAATGCCGAATCCACTAAAGATATCGACGAAGCCTTGCAGCAAGGTAAACATCGACCAGATAAAGATCGACCAGAAGTAGACGCGTTTTGAACCAAAACGGTCCAGCAACCAGCCACCAGGGATCTGCCCGATaa                                                                                                            <  1:1818907/150‑1 (MQ=255)
                                            aaCAGGGCAACTATAATGCCGAATCCACTAAAGATATCGACGAAGCCTTGCAGCAAGGTAAACATCGACCAGATAAAGATCGACCAGAAGTAGACGCGTTTTGAACCAAAACGGTCCAGCAACCAGCCACCAGGGATCTGCCCGATAACa                                                                                                          <  2:1974154/150‑1 (MQ=255)
                                             aCAGGGCAACGATAATGCCGAATCCACTAAAGATATCGACGAAGCCTTGCAGCAAGGTAAACATCGACCAGATAAAGATCGACCAGAAGTAGACGCGTTTTGAACCAAAACGGTCCAGCAACCAGCCACCAGGGATCTGCCCGATAACAt                                                                                                         <  1:2103408/150‑1 (MQ=255)
                                                 ggCAACGATAATGCCGAATCCACTAAAGATATCGACGAAGCCTTGCAGCAAGGTAAACATCGACCAGATAAAGCTCGACCAGAAGTAGACGCGTTTTGAACCAAAACGGTCCAGCAACCAGCCACCAGGGATCTGCCCGATAACATAAGc                                                                                                     >  2:2176210/1‑150 (MQ=255)
                                                  gCAACGATAATGCCGAATCCACTAAAGATATCGACGAAGCCTTGCAGCAAGGTAAACATCGACCAGATAAAGATCGACCAGAAGTAGACGCGTTTTGAACCAAAACGGTCCAGCAACCAGCCACCAGGGATCTGCCCGATAACATAAGcc                                                                                                    <  1:2003805/150‑1 (MQ=255)
                                                  gCAACGATAATGCCGAATCCACTAAAGATATCGACGAAGCCTTGCAGCAAGGTAAACATCGACCAGATAAAGATCGACCAGAAGTAGACGCGTTTTGAACCAAAACGGTCCAGCAACCAGCCACCAGGGATCTGCCCGATAACATAAGcc                                                                                                    >  2:116644/1‑150 (MQ=255)
                                                  gCAACGATAATGCCGAATCCACTAAAGATATCGACGAAGCCTTGCAGCAAGGTAAACATCGACCAGATAAAGATCGACCAGAAGTAGACGCGTTTTGAACCAAAACGGTCCAGCAACCAGCCACCAGGGATCTGCCCGATAACATAAGcc                                                                                                    >  1:1815864/1‑150 (MQ=255)
                                                       gATAATGCCGAATCCACTAAAGATATCGACGAAGCCTTGCAGCAAGGTAAACATCGACCAGATAAAGATCGACCAGAAGTAGACGCGTTTTGAACCAAAACGGTCCAGCAACCAGCCACCAGGGATCTGCCCGATAACATAAGCCCATga                                                                                               <  2:156959/150‑1 (MQ=255)
                                                             gCCGAATCCACTAAAGATATCGACGAAGCATTGCAGCAAGGTAAACATCCACCAGATAAAGATCGACCAGAAGTAGACGCGTTTTGAACCAAAACGGTCCAGCAACCAGCCACCAGGGATCTGCCCGAAAACAAAAGCCCATGAGAAAGc                                                                                         >  2:2221283/1‑150 (MQ=255)
                                                               cGAATCCACTAAAGATATCGACGAAGCCTTGCAGCAAGGTAAACATCGACCAGATAAAGATCGACCAGAAGTAGACGCGTTTTGAACCAAAACGGTCCAGCAACCAGCCACCAGGGATCTGCCCGATAACATAAGCCCATGAGAAAGCag                                                                                       >  1:1871160/1‑150 (MQ=255)
                                                               cGAATCCACTAAAGATATCGACGAAGCCTTGCAGCAAGGTAAACATCGACCAGATAAAGATCGACCAGAAGTAGACGCGTTTTGAACCAAAACGGTCCAGCAACCAGCCACCAGGGATCTGCCCGATAACATAAGCCCATGAGAAAGCag                                                                                       >  1:1049768/1‑150 (MQ=255)
                                                                    ccACTAAAGATATCGACGAAGCCTTGCAGCAAGGTAAACATCGACCAGATAAAGATCGACCAGAAGTAGACGCGTTTTGAACCAAAACGGTCCAGCAACCAGCCACCAGGGATCTGCCCGATAACATAAGCCCATGAGAAAGCAGAGAac                                                                                  <  2:827066/150‑1 (MQ=255)
                                                                       cTAAAGATATCGACGAAGCCTTGCAGCAAGGTAAACATCGACCAGATAAAGATCGACCAGAAGTAGACGCGTTTTGAACCAAAACGGTCCAGCAACCAGCCACCAGGGATCTGCCCGATAACATAAGCCCATGAGAAAGCAGAGAAcaca                                                                               <  2:730131/150‑1 (MQ=255)
                                                                        tAAAGATATCGACGAAGCCTTGCAGCAAGGTAAACATCGACCAGATAAAGATCGACCAGAAGTAGACGCGTTTTGAACCAAAACGGTCCAGCAACCAGCCACCAGGGATCTGCCCGATAACATAAGCCCATGAGAAAGCAGAGAACACAt                                                                              >  2:1078633/1‑150 (MQ=255)
                                                                                  gacgaAGCCTTGCAGCAAGGTAAACATCGACCAGATAAAGATCGACCAGAAGTAGACGCGTTTTGAACCAAAACGGTCCAGCAACCTGCCACCAGGGATCTGCCCGATAACATAAGCCCATGAGAAAGCAGAGAACACATAGCCCATTcc                                                                    >  2:2362014/1‑150 (MQ=255)
                                                                                    cgaAGCCTTGCAGCAAGGTAAACATCGACCAGATAAAGATCGACCAGAAGTAGACGCGTTTTGAACCAAAACGGTCCAGCAACCAGCCACCAGGGATCTGCCCGATAACATAAGCCCATGAGAAAGCAGAGAACACATAGCCCATTCCCa                                                                  >  2:211888/1‑150 (MQ=255)
                                                                                     gaAGCCTTGCAGCAAGGTAAACATCGACCAGATAAAGATCGACCAGAAGTAGACGCGTTTTGAACCAAAACGGTCCAGCAACCAGCCACCAGGGATCTGCCCGATAACATAAGCCCATGAGAAAGCAGAGAACACATAGCCCATTCCCAc                                                                 >  1:225484/1‑150 (MQ=255)
                                                                                      aaGCCTTGCAGCAAGGTAAACATCGACCAGATAAAGATCGACCAGAAGTAGACGCGTTTTGAACCAAAACGGTCCAGCAACCAGCCACCAGGGATCTGCCCGATAACATAAGCCCATGAGAAAGCAGAGAACACATAGCCCATTCCCACg                                                                <  2:368264/150‑1 (MQ=255)
                                                                                      aaGCCTTGCAGCAAGGTAAACATCGACCAGATAAAGATCGACCAGAAGTAGACGCGTTTTGAACCAAAACGGTCCAGCAACCAGCCACCAGGGATCTGCCCGATAACATAAGCCCATGAGAAAGCAGAGAACACATAGCCCATTCCCACg                                                                >  1:1330721/1‑150 (MQ=255)
                                                                                       aGCCTTGCAGCAATGTAAACATCTACCAGATAAAGATCGACCAGAAGTAGACGCGTTTTGAACCAAAACGGTCCAGCAACCAGCCACCAGGGATCTGCCCGATAACATAAGCCCATGAGAAAGCAGAGAACACATAGCCCATTCCCACgg                                                               <  1:310145/150‑1 (MQ=255)
                                                                                        gCCTTGCAGCAAGGTAAACATCGACCAGATAAAGATCGACCAGAAGTAGACGCGTTTTGAACCAAAACGGTCCAGCAACCAGCCACCAGGGATCTGCCCGATAACATAAGCCCATGAGAAAGCAGAGAACACATAGCCCATTCCCACggg                                                              >  2:1466670/1‑150 (MQ=255)
                                                                                         ccTTGCAGCAAGGTAAACATCGACCAGATAAAGATCGACCAGAAGTAGACGCGTTTTGAACCAAAACGGTCCAGCAACCAGCCACCAGGGATCTGCCCGATAACATAAGCCCATGAGAAAGCAGAGAACACATAGCCCATTCCCACGGGa                                                             >  2:279420/1‑150 (MQ=255)
                                                                                          cTTGCAGCAAGGTAAACATCGACCAGATAAAGATCGACCAGAAGTAGACGCGTTTTGAACCAAAACGGTCCAGCAACCAGCCACCAGGGATCTGCCCGATAACATAAGCCCATGAGAAAGCAGAGAACACATAGCCCATTCCCACGGGAt                                                            >  2:163826/1‑150 (MQ=255)
                                                                                          cTTGCAGCAAGGTAAACATCGACCAGATAAAGATCGACCAGAAGTAGACGCGTTTTGAACCAAAACGGTCCAGCAACCAGCCACCAGGGATCTGCCCGATAACATAAGCCCATGAGAAAGCAGAGAACACATAGCCCAATCCCACGGGAt                                                            >  2:2085598/1‑150 (MQ=255)
                                                                                           ttGCAGCAAGGTAAACATCGACCAGATAAAGATCGACCAGAAGTAGACGCGTTTTGAACCAAAACGGTCCAGCAACCAGCCACCAGGGATCTGCCCGATAACATAAGCCCATGAGAAAGCAGAGAACACATAGCCCATTCCCACGGGATc                                                           <  2:312694/150‑1 (MQ=255)
                                                                                            tGCAGCAAGGTAAACATCGACCAGATAAAGATCGACCAGAAGTAGACGCGTTTTGAACCAAAACGGTCCAGCAACCAGCCACCAGGGATCTGCCCGATAACATAAGCCCATGAGAAAGCAGAGAACACATAGCCCATTCCCACGGGATCa                                                          <  1:853804/150‑1 (MQ=255)
                                                                                            tGCAGCAAGGTAAACATCAACCAGATAAAGATCGACCAGAAGTAGACGCGTTTTGAACCAAAACGGTCCAGCAACCAGCCACCAGGGATCTGCCCGATAACATAAGCCCATGAGAAAGCAGAGAACACATAGCCCATTCCCACGGGATCa                                                          >  1:461907/1‑150 (MQ=255)
                                                                                              cagcaAGGTAAACATCGACCAGATAAAGATCGACCAGAAGTAGACGCGTTTTGAACCAAAACGGTCCAGCAACCAGCCACCAGGGATCTGCCCGATAACATAAGCCCATGAGAAAGCAGAGAACACATAGCCCATTCCCACGGGATCAAg                                                        <  2:7125/150‑1 (MQ=255)
                                                                                              cagcaAGGTAAACATCGACCAGATAAAGATCGACCAGAAGTAGACGCGTTTTGAACCAAAACGGTCCAGCAACCAGCCACCAGGGATCTGCCCGATAACATAAGCCCATGAGAAAGCAGAGAACACATAGCCCATTCCCACGGGATCAAg                                                        <  1:1433382/150‑1 (MQ=255)
                                                                                                gcaAGGTAAACATCGACCAGATAAAGATCGACCAGAAGTAGACGCGTTTTGAACCAAAACGGTCCAGCAACCAGCCACCAGGGATCTGCCCGATAACATAAGCCCATGAGAAAGCAGAGAACACATAGCCCATTCCCACGGGATCAAGGc                                                      >  1:1605981/1‑150 (MQ=255)
                                                                                                 caAGGTAAACATCGACCAGATAAAGATCGACCAGAAGTAGACGCGTTTTGAACCAAAACGGTCCAGCAACCAGCCACCAGGGATCTGCCCGATAACATAAGCCCATGAGAAAGCAGAGAACACATAGCCCATTCCCACGGGATCAAGGcc                                                     <  1:1492407/150‑1 (MQ=255)
                                                                                                 caAGGTAAACATCGACCAGATAAAGATCGACCAGAAGTAGACGCGTTTTGAACCAAAACGGTCCAGCAACCAGCCACCAGGGATCTGCCCGATAACATAAGCCCATGAGAAAGCAGAGAACACATAGCCCATTCCCACGGGATCAAGGcc                                                     <  2:29256/150‑1 (MQ=255)
                                                                                                  aaGGTAAACATCGACCAGATAAAGATCGAACAGAAGTAGACGCGTTTTGAACCAAAACGGTCCAGCAACCAGCCACCAGGGATCTGCCCGATAACATAAGCCCATGAGAAAGCAGAGAACACATAGCCCATTCCCACGGGATCAAGGCCg                                                    <  1:2476546/150‑1 (MQ=255)
                                                                                                   aGGTAAACATCGACCAGATAAAGATCGACCAGAAGTAGACGCGTTTTGAACCAAAACGGTCCAGCAACCAGCCACCAGGGATCTGCCCGATAACATAAGCCCATGAGAAAGCAGAGAACACATAGCCCATTCCCACGGGATCAAGGCCGa                                                   <  2:355713/150‑1 (MQ=255)
                                                                                                   aGGTAAACATCGACCAGATAAAGATCGACCAGAAGTAGACGCGTTTTGAACCAAAACGGTCCAGCAACCAGCCACCAGGGATCTGCCCGATAACATAAGCCCATGAGAAAGCAGAGAACACATAGCCCATTCCCACGGGATCAAGGCCGa                                                   <  2:1129240/150‑1 (MQ=255)
                                                                                                   aGGTAAACATCGACCAGATAAAGATCGACCAGAAGTAGACGCGTTTTGAACCAAAACGGTCCAGCAACCAGCCACCAGGGATCTGCCCGATAACATAAGCCCATGAGAAAGCAGAGAACACATAGCCCATTCCCACGGGATCAAGGCCGa                                                   <  2:1540886/150‑1 (MQ=255)
                                                                                                          cATCGACCAGATAAAGATCGACCAGAAGTAGACGCGTTTTGAACCAAAACGGTCCAGCAACCAGCCACCAGGGATCTGCCCGATAACATAAGCCCATGAGAAAGCAGAGAACACATAGCCCATTCCCACGGGATCAAGGCCGATATCttt                                            >  2:1486285/1‑150 (MQ=255)
                                                                                                            tCGACCAGATAAAGATCGACCAGAAGTAGACGCGTTTTGAACCAAAACGGTCCAGCAACCAGCCACCAGGGATCTGCCCGATAACATAAGCCCATGAGAAAGCAGAGAACACATAGCCCATTCCCACGGGATCAAGGCCGATATCTTTgg                                          <  1:108161/150‑1 (MQ=255)
                                                                                                            tCGACCAGATAAAGATCGACCAGAAGTAGACGCGTTTTGAACCAAAACGGTCCAGCAACCAGCCACCAGGGATCTGCCCGATAACATAAGCCCATGAGAAAGCAGAGAACACATAGCCCATTCCCACGGGATCAAGGCCGATATCTTTgg                                          <  2:1473842/150‑1 (MQ=255)
                                                                                                            tCGACCAGATAAAGATCGACCAGAAGTAGACGCGTTTTGAACCAAAACGGTCCAGCAACCAGCCACCAGGGATCTGCCCGATAACATAAGCCCATGAGAAAGCAGAGAACACATAGCCCATTCCCACGGGATCAAGGCCGATATCTTTgg                                          <  1:2207449/150‑1 (MQ=255)
                                                                                                             cGACCAGATAAAGATCGACCAGAAGTAGACGCGTTTTGAACCAAAACGGTCCAGCAACCAGCCACCAGGGATCTGCCCGATAACATAAGCCCATGAGAAAGCAGAGAACACATAGCCCATTCCCACGGGATCAAGGCCGATATCTTTGGc                                         <  1:860976/150‑1 (MQ=255)
                                                                                                             cGACCAGATAAAGATCGACCAGAAGTAGACGCGTTTTGAACCAAAACGGTCCAGCAACCAGCCACCAGGGATCTGCCCGATAACATAAGCCCATGAGAAAGCAGAGAACACATAGCCCATTCCCACGGGATCAAGGCCGATATCTTTGGc                                         <  1:2299572/150‑1 (MQ=255)
                                                                                                              gACCAGATAAAGATCGACCAGAAGTAGACGCGTTTTGAACCAAAACGGTCCAGCAACCAGCCACCAGGGATCTGCCCGATAACATAAGCCCATGAGAAAGCAGAGAACACATAGCCCATTCCCACGGGATCAAGGCCGATATCTTTGGcc                                        >  2:1635700/1‑150 (MQ=255)
                                                                                                                    aTAAAGATCGACCAGAAGTAGACGCGTTTTGAACCAAAACGGTCCAGCAACCAGCCACCAGGGATCTGCCCGATAACATAAGCCCATGAGAAAGCAGAGAACACATAGCCCATTCCCACGGGATCAAGGCCGATATCTTTGGCCATTTcc                                  <  1:2088885/150‑1 (MQ=255)
                                                                                                                     tAAAGATCGACCAGAAGTAGACGCGTTTTGAACCAAAACGGTCCAGCAACCAGCCACCAGGGATCTGCCCGAGAACATAAGCCCATGAGAACGCAGAGAACACCTCGCCCATTCCCACGGGATCAAGGCCGATAGCGTTGGCCATTTCCg                                 >  1:1688552/1‑150 (MQ=255)
                                                                                                                        aGATCGACCAGAAGTAGACGCGTTTTGAACCAAAACGGTCCAGCAACCAGCCACCAGGGATCTGCCCGATAACATAAGCCCATGAGAAAGCAGAGAACACATAGCCCATTCCCACGGGATCAAGGCCGATATCTTTGGCCATTTCCGAAc                              >  1:1508901/1‑150 (MQ=255)
                                                                                                                        aGATCGACCAGAAGTAGACGCGTTTTGAACCAAAACGGTCCAGCAACCAGCCACCAGGGATCTGCCCGATAACATAAGCCCATGAGAAAGCAGAGAACACATAGCCCATTCCCACGGGATCAAGGCCGATATCTTTGGCCATTTCCGAAc                              <  2:2452757/150‑1 (MQ=255)
                                                                                                                         gATCGACCAGAAGTAGACGCGTTTTGAACCAAAACGGTCCAGCAACCAGCCACCAGGGATCTGCCCGATAACATAAGCCCATGAGAAAGCAGAGAACACATAGCCCATTCCCACGGGATCAAGGCCGATATCTTTGGCCATTTCCGAAcc                             >  2:1077138/1‑150 (MQ=255)
                                                                                                                         gATCGACCAGAAGTAGACGCGTTTTGAACCAAAACGGTCCAGCAACCAGCCACCAGGGATCTGCCCGATAACATAAGCCCATGAGAAAGCAGAGAACACATAGCCCATTCCCACGGGATCAAGGCCGATATCTTTGGCCATTTCCGAAcc                             >  2:1069634/1‑150 (MQ=255)
                                                                                                                          aTCGACCAGAAGTAGACGCGTTTTGAACCAAAACGGTCCAGCAACCAGCCACCAGGGATCTGCCCGATAACATAAGCCCATGAGAAAGCAGAGAACACATAGCCCATTCCCACGGGATCAAGGCCGATATCTTTGGCCATTTCCGAACCg                            <  1:52442/150‑1 (MQ=255)
                                                                                                                            cGACCAGAAGTAGACGCGTTTTGAACCAAAACGGTCCAGCAACCAGCCACCAGGGATCTGCCCGATAACATAAGCCCATGAGAAAGCAGAGAACACATAGCCCATTCCCACGGGATCAAGGCCGATATCTTTGGCCATTTCCGAACCGGc                          <  1:256753/150‑1 (MQ=255)
                                                                                                                            cGACCAGAAGTAGACGCGTTTTGAACCAAAACGGTCCAGCAACCAGCCACCAGGGATCTGCCCGATAACATAAGCCCATGAGAAAGCAGAGAAACATAGCCCATTCCCACGGGATCAAGGCCGATATCTTTGGCCATTTCCGAACCGGCg                         <  1:297637/150‑1 (MQ=255)
                                                                                                                              aCCAGAAGTAGACGCGTTTTGAACCAAAACGGTCCAGCAACCAGCCACCAGGGATCTGCCCGATAACATAAGCCCATGAGAAAGCAGAGAACACATAGCCCATTCCCACGGGATCAAGGCCGATATCTTTGGCCATTTCCGAACCGGCGa                        >  2:26940/1‑150 (MQ=255)
                                                                                                                               ccAGAAGTAGACGCGTTTTGAACCAAAACGGTCCAGCAACCAGCCACCAGGGATCTGCCCGATAACATAAGCCCATGAGAAAGCAGAGAACACATAGCCCATTCCCACGGGATCAAGGCCGATATCTTTGGCCATTTCCGAACCGGCGAt                       <  1:1467773/150‑1 (MQ=255)
                                                                                                                                cAGAAGTAGACGCGTTTTGAACCAAAACGGTCCAGCAACCAGCCACCAGGGATCTGCCCGATAACATAAGCCCATGAGAAAGCAGAGAACACATAGCCCATTCCCACGGGATCAAGGCCGATATCTTTGGCCATTTCCGAACCGGCGATa                      >  1:479240/1‑150 (MQ=255)
                                                                                                                                        ggaccgTTTTGGTTCAAAACGGTCCAGCAACCAGCCACCAGGGATCTGCCCGATAACATAAGCCCATGAGAAAGCAGAGAACACATAGCCCATTCCCACGGGATCAAGGCCGATATCTTTGGCCATTTCCGAACCGGCGATAGAGAGCGt              >  1:2442405/5‑150 (MQ=255)
                                                                                                                                        gACGCGTTTTGAACCAAAACGGTCCAGCAACCAGCCACCAGGGATCTGCCCGATAACATAAGCCCATGAGAAAGCAGAGAACACATAGCCCATTCCCACGGGATCAAGGCCGATATCTTTGGCCATTTCCGAACCGGCGATAGAGAGCGt              >  1:1308835/1‑150 (MQ=255)
                                                                                                                                         gaccgTTTTGGTTCAAAACGGTCCAGCAACCAGCCACCAGGGATCTGCCCGATAACATAAGCCCATGAGAAAGCAGAGAACACATAGCCCATTCCCACGGGATCAAGGCCGATATCTTTGGCCATTTCCGAACCGGCGATAGAGAGCGTa             >  2:1872851/4‑150 (MQ=255)
                                                                                                                                         gaccgTTTTGGTTCAAAACGGTCCAGCAACCAGCCACCAGGGATCTGCCCGATAACATAAGCCCATGAGAAAGCAGAGAACACATAGCCCATTCCCACGGGATCAAGGCCGATATCTTTGGCCATTTCCGAACCGGCGATAGAGAGCGTa             >  2:583577/4‑150 (MQ=255)
                                                                                                                                               tttGAACCAAAACGGTCCAGCAACCAGCCACCAGGGATCTGCCCGATAACATAAGCCCATGAGAAAGCAGAGAACACATGGCCCATTCCCACGGGATCAAGGCCGATATCTTTGGCCATTTCCGAACCGGCGATAGAGAGCGTAGCGCgg       >  1:577967/1‑150 (MQ=255)
                                                                                                                                                ttGAACCAAAACGGTCCAGCAACCAGCCACCAGGGATCTGCCCGATTACATAAGCCCATGAGAAAGCAGAGAACACATAGCCCATTCCCACGGGATCAAGGCCGATATCTTTGGCCATTTCCGAACCGGCGATAGAGAGCGTAGCGCggg      >  2:1677943/1‑149 (MQ=255)
                                                                                                                                                ttGAACCAAAACGGTCCAGCAACCAGCCACCAGGGATCTGCCCGATAACATAAGCCCATGAGAAAGCAGAGAACACATAGCCCATTCCCACGGGATCAAGGCCGATATCTTTGGCCATTTCCGAACCGGCGATAGAGAGCGTAGCGCGGt      <  1:116644/150‑1 (MQ=255)
                                                                                                                                                   aaCCAAAACGGTCCAGCAACCAGCCACCAGGGATCTGCCCGATAACATAAGCCCATGAGAAAGCAGAGAACACATAGCCCATTCCCACGGGATCAAGGCCGATATCTTTGGCCATTTCCGAACCGGCGATAGAGAGCGTAGCGCGGTCGc   <  2:2402766/150‑1 (MQ=255)
                                                                                                                                                   aaCCAAAACGGTCCAGCAACCAGCCACCAGGGATCTGCCCGATAACATAAGCCCATGAGAAAGCAGAGAACACATAGCCCATTCCCACGGGATCAAGGCCGATATCTTTGGCCATTTCCGAACCGGCGATAGAGAGCGTAGCGCGGTCGc   <  2:1790744/150‑1 (MQ=255)
                                                                                                                                                    aCCAAAACGGTCCAGCAACCAGCCACCAGGGATCTGCCCGATAACATAAGCCCATGAGAAAGCAGAGAACACATAGCCCATTCCCACGGGATCAAGGCCGATATCTTTGGCCATTTCCGAACCGGCGATAGAGAGCGTAGCGCGGTCGcc  >  1:2482084/1‑150 (MQ=255)
                                                                                                                                                    aCCAAAACGGTCCAGCAACCAGCCACCAGGGATCTGCCCGATAACATAAGCCCATGAGAAAGCAGAGAACACATAGCCCATTCCCACGGGATCAAGGCCGATATCTTTGGCCATTTCCGAACCGGCGATAGAGAGCGTAGCGCGGTCGcc  <  2:144938/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GGAAAGATGGCGCTTCAGCAAGCCCGACCAGGAAGCGCAGCGTAAACAGGGCAACGATAATGCCGAATCCACTAAAGATATCGACGAAGCCTTGCAGCAAGGTAAACATCGACCAGATAAAGATCGACCAGAAGTAGACGCGTTTTGAACCAAAACGGTCCAGCAACCAGCCACCAGGGATCTGCCCGATAACATAAGCCCATGAGAAAGCAGAGAACACATAGCCCATTCCCACGGGATCAAGGCCGATATCTTTGGCCATTTCCGAACCGGCGATAGAGAGCGTAGCGCGGTCGCC  >  NC_000913/2921704‑2922001

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: