Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,921,852 A→T 6.0% G83G (GGT→GGA) ‡ 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,8520AT 6.0% 445.5 / 11.8 147G83G (GGT→GGA) ‡gudPgalactarate/glucarate/glycerate transporter GudP
Reads supporting (aligned to +/- strand):  ref base A (68/70);  new base T (4/5);  total (72/75)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.40e-01

Only 100 of 153 total aligned reads displayed.

GGGAAAGATGGCGCTTCAGCAAGCCCGACCAGGAAGCGCAGCGTAAACAGGGCAACGATAATGCCGAATCCACTAAAGATATCGACGAAGCCTTGCAGCAAGGTAAACATCGACCAGATAAAGATCGACCAGAAGTAGACGCGTTTTGAACCAAAACGGTCCAGCAACCAGCCACCAGGGATCTGCCCGATAACATAAGCCCATGAGAAAGCAGAGAACACATAGCCCATTCCCACGGGATCAAGGCCGATATCTTTGGCCATTTCCGAACCGGCGATAGAGAGCGTAGCGCGGTCGCC  >  NC_000913/2921703‑2922001
                                                                                                                                                     |                                                                                                                                                     
gggAAAGATGGCGCTTCAGCAAGCCCGACCAGGAAGCGCAGCGTAAACAGGGCAACGATAATGCCGAATCCACTAAAGATATCGACGAAGCCTTGCAGCAAGGTAAACATCGACCAGATAAAGATCGACCAGAAGTAGACGCGTTTTGaa                                                                                                                                                       <  1:1095452/150‑1 (MQ=255)
 ggAAAGATGGCGCTTCAGCAAGCCCGACCAGGAAGCGCAGCGTAAACAGGGCAACGATAATGCCGAATCCACTAAAGATATCGACGAAGCCTTGCAGCAAGGTAAACATCGACCAGATAAAGATCGACCAGAAGTAGACGCGTTTTGAAc                                                                                                                                                      <  2:2058144/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)
       aTGGCGCTTCAGCAAGCCCGCCCAGGAAGCGCAGCGTAAACAGGGCAACGATAATGCCGAATCCACTAAAGATATCGACGAAGCCTTGCAGCAAGGTAAACATCGACCAGATAAAGATCGACCAGAAGTAGACGCGTTTTGAACCAAAAc                                                                                                                                                >  1:144938/1‑150 (MQ=255)
       aTGGCGCTTCAGCAAGCCCGACCAGGAAGCGCAGCGTAAACAGGGCAACGATAATGCCGAATCCACTAAAGATATCGACGAAGCCTTGCAGCAAGGTAAACATCGACCAGATAAAGATCGACCAGAAGTAGACGCGTTTTGAACCAAAAc                                                                                                                                                >  1:1212184/1‑150 (MQ=255)
          gcgcTTCAGCAAGCCCGACCAGGAAGCGCAGCGTAAACAGGGCAACGATAATGCCGAATCCACTAAAGATATCGACGAAGCCTTGCAGCAAGGTAAACATCGACCAGATAAAGATCGACCAGAAGTAGACGCGTTTTGAACCAAAACGGt                                                                                                                                             <  2:669705/150‑1 (MQ=255)
          gcgcTTCAGCAAGCCCGACCAGGAAGCGCAGCGTAAACAGGGCAACGATAATGCCGAATCCACTAAAGATATCGACGAAGCCTTGCAGCAAGGTAAACATCGACCAGATAAAGATCGACCAGAAGTAGACGCGTTTTGAACCAAAACGGt                                                                                                                                             <  2:1782314/150‑1 (MQ=255)
           cgcTTCAGCAAGCCCGACCAGGAAGCGCAGCGTAAACAGGGCAACGATAATGCCGAATCCACTAAAGATATCGACGAAGCCTTGCAGCAAGGTAAACATCGACCAGATAAAGATCGACCAGAAGTAGACGCGTTTTGGTTCAAAACgcgt                                                                                                                                            <  2:823145/150‑4 (MQ=255)
           cgcTTCAGCAAGCCCGACCAGGAAGCGCAGCGTAAACAGGGCAACGATAATGCCGAATCCACTAAAGATATCGACGAAGCCTTGCAGCAAGGTAAACATCGACCAGATAAAGATCGACCAGAAGTAGACGCGTTTTGGTTCAAAACgcgt                                                                                                                                            <  2:2482831/150‑4 (MQ=255)
           cgcTTCAGCAAGCCCGACCAGGAAGCGCAGCGTAAACAGGGCAACGATAATGCCGAATCCACTAAAGATATCGACGAAGCCTTGCAGCAAGGTAAACATCGACCAGATAAAGATCGACCAGAAGTAGACGCGTTTTGAACCAAAACGGTc                                                                                                                                            <  2:677625/150‑1 (MQ=255)
            gcTTCAGCAAGCCCGACCAGGAAGCGCAGCGTAAACAGGGCAACGATAATGCCGAATCCACTAAAGATATCGACGAAGCCTTTCAGCAAGGTAAACATCGACCAGATAAAGATCGACCAGAAGTAGACGCGTTTTGAACCAAAACGGTcc                                                                                                                                           <  2:1750349/150‑1 (MQ=255)
            gcTTCAGCAAGCCCGACCAGGAAGCGCAGCGTAAACAGGGCAACGATAATGCCGAATCCACTAAAGATATCGACGAAGCCTTGCAGCAAGGTAAACATCGACCAGATAAAGATCGACCAGAAGTAGACGCGTTTTGGTTCAAAACgcgtc                                                                                                                                           <  2:1626957/150‑5 (MQ=255)
            gcTTCAGCAAGCCCGACCAGGAAGCGCAGCGTAAACAGGGCAACGATAATGCCGAATCCACTAAAGATATCGACGAAGCCTTGCAGCAAGGTAAACATCGACCAGATAAAGATCGACCAGAAGTAGACGCGTTTTGGTTCAAAACgcgtc                                                                                                                                           <  1:2043299/150‑5 (MQ=255)
            gcTTCAGCAAGCCCGACCAGGAAGCGCAGCGTAAACAGGGCAACGATAATGCCGAATCCACTAAAGATATCGACGAAGCCTTGCAGCAAGGTAAACATCGACCAGATAAAGATCGACCAGAAGTAGACGCGTTTTGAACCAAAACGGTcc                                                                                                                                           <  1:872197/150‑1 (MQ=255)
            gcTTCAGCAAGCCCGACCAGGAAGCGCAGCGTAAACAGGGCAACGATAATGCCGAATCCACTAAAGATATCGACGAAGCCTTGCAGCAAGGTAAACATCGACCAGATAAAGATCGACCAGAAGTAGACGCGTTTTGAACCAAAACGGTcc                                                                                                                                           <  1:1484114/150‑1 (MQ=255)
            gcTTCAGCAAGCCCGACCAGGAAGCGCAGCGTAAACAGGGCAACGATAATGCCGAATCCACTAAAGATATCGACGAAGCCTTGCAGCAAGGTAAACATCGACCAGATAAAGATCGACCAGAAGTAGACGCGTTTTGAACCAAAACGGTcc                                                                                                                                           <  1:1390911/150‑1 (MQ=255)
            gcTTCAGCAAGCCCGACCAGGAAGCGCAGCGTAAACAGGGCAACGATAATGCCGAATCCACTAAAGATATCGACGAAGCCTTGCAGCAAGGTAAACATCGACCAGATAAAGATCGACCAGAAGTAGACGCGTTTTGAACCAAAACGGTcc                                                                                                                                           >  2:575910/1‑150 (MQ=255)
             cTTCAGCAAGCCCGACCAGGAAGCGCAGCGTAAACAGGGCAACGATAATGCCGAATCCACTAAAGATATCGACGAAGCCTTGCAGCAAGGTAAACATCGACCAGATAAAGATCGACCAGAAGTAGACGCGTTTTGAACCAAAACGGTccc                                                                                                                                          >  2:2421862/1‑149 (MQ=255)
             cTTCAGCAAGCCCGACCAGGAAGCGCAGCGTAAACAGGGCAACGATAATGCCGAATCCACTAAAGATATCGACGAAGCCTTGCAGCAAGGTAAACATCGACCAGATAAAGATCGACCAGAAGTAGACGCGTTTTGAACCAAAACGGTCCa                                                                                                                                          >  2:2364883/1‑150 (MQ=255)
             cTTCAGCAAGCCCGACCAGGAAGCGCAGCGTAAACAGGGCAACGATAATGCCGAATCCACTAAAGATATCGACGAAGCCTTGCAGCAAGGTAAACATCGACCAGATAAAGATCGACCAGAAGTAGACGCGTTTTGAACCAAAACGGTCCa                                                                                                                                          >  2:2103408/1‑150 (MQ=255)
             cTTCAGCAAGCCCGACCAGGAAGCGCAGCGTAAACAGGGCAACGATAATGCCGAATCCACTAAAGATATCGACGAAGCCTTGCAGCAAGGTAAACATCGACCAGATAAAGATCGACCAGAAGTAGACGCGTTTTGAACCAAAACGGTCCa                                                                                                                                          <  2:307214/150‑1 (MQ=255)
             cTTCAGCAAGCCCGACCAGGAAGCGCAGCGTAAACAGGGCAACGATAATGCCGAATCCACTAAAGATATCGACGAAGCCTTGCAGCAAGGTAAACATCGACCAGATAAAGATCGACCAGAAGTAGACGCGTTTTGAACCAAAACGGTCCa                                                                                                                                          >  2:1084917/1‑150 (MQ=255)
              ttCAGCAAGCCCGACCAGGAAGCGCAGCGTAAACAGGGCAACGATAATGCCGAATCCACTAAAGATATCGACGAAGCCTTGCAGCAAGGTAAACATCGACCAGATAAAGATCGACCAGAAGTAGACGCGTTTTGAACCAAAACGGTCCAg                                                                                                                                         <  1:1548252/150‑1 (MQ=255)
              ttCAGCAAGCCCGACCAGGAAGCGCAGCGTAAACAGGGCAACGATAATGCCGAATCCACTAAAGATATCGACGAAGCCTTGCAGCAAGGTAAACATCGACCAGATAAAGATCGACCAGAAGTAGACGCGTTTTGAACCAAAACGGTCCAg                                                                                                                                         <  2:1672129/150‑1 (MQ=255)
                 aGCAAGCCCGACCAGGAAGCGCAGCGTAAACAGGGCAACGATAATGCCGAATCCACTAAAGATATCGACGAAGCCTTGCAGCAAGGTAAACATCGACCAGATAAAGATCGACCAGAAGTAGACGCGTTTTGAACCAAAACGGTCCAGCaa                                                                                                                                      >  1:228463/1‑150 (MQ=255)
                  gCAAGCCCGACCAGGAAGCGCAGCGTAAACAGGGCAACGATAATGCCGAATCCACTAAAGATATCGACGAAGCCTTGCAGCAAGGTAAACATCGACCAGATAAAGATCGACCAGAAGTAGACGCGTTTTGAACCAAAACGGTCCAGCAAc                                                                                                                                     >  2:297637/1‑150 (MQ=255)
                        ccGACCAGGAAGCGCAGCGTAAACAGGGCAACGATAATGCCGAATCCACTAAAGATATCGACGAAGCCTTGCAGCAAGGTAAACATCGACCAGATAAAGATCGACCAGAAGTAGACGCGTTTTGAACCAAAACGGTCCAGCAACCAGcca                                                                                                                               >  1:476972/1‑150 (MQ=255)
                            ccAGGAAGCGCAGCGTAAACAGGGCAACGATAATGCCGAATCCACTAAAGATATCGACGAAGCCTTGCAGCAAGGTAAACATCGACCAGATAAAGATCGACCAGAAGTAGACGCGTTTTGAACCAAAACGGTCCAGCAACCAGCCACCAg                                                                                                                           >  2:1430045/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)
                                  aGCGCAGCGTAAACAGGGCAACGATAATGCCGAATCCACTAAAGATATCGACGAAGCCTTGCAGCAAGGTAAACATCGACCAGATAAAGATCGACCAGAAGTAGACGCGTTTTGAACCAAAACGGTCCAGCAACCAGCCACCAGGGATCt                                                                                                                     >  1:368264/1‑150 (MQ=255)
                                     gcAGCGTAAACAGGGCAACGATAATGCCGAATCCACTAAAGATATCGACGAAGCCTTGCAGCAAGGTAAACATCGACCAGATAAAGATCGACCAGAAGTAGACGCGTTTTGAACCAAAACGGTCCAGCAACCAGCCACCAGGGATCTGcc                                                                                                                  >  1:2230729/1‑150 (MQ=255)
                                     gcAGCGTAAACAGGGCAACGATAATCCCGAATCCACTAAAGATATCGACGAAGCCTTGCAGCAAGGTAAACATCGACCAGATAAAGATCGACCAGAAGTAGACGCGTTTTGAACCAAAACGGTCCAGCAACCAGCCACCAGGGATCTGcc                                                                                                                  <  1:163990/150‑1 (MQ=255)
                                      cAGCGTAAACAGGGCAACGATAATGCCGAATCCACTAAAGATATCGACGAAGCCTTGCAGCAAGGTAAACATCGACCAGATAAAGATCGACCAGAAGTAGACGCGTTTTGAACCAAAACGGTCCAGCAACCAGCCACCAGGGATCTGccc                                                                                                                 <  2:1470012/150‑1 (MQ=255)
                                        gCGTAAACAGGGCAACGATAATGCCGAATCCACTAAAGATATCGACGAAGCCTTGCAGCAAGGTAAACATCGACCAGATAAAGATCGACCAGAAGTAGACGCGTTTTGAACCAAAACGGTCCAGCAACCAGCCACCAGGGATCTGCCCGa                                                                                                               >  1:45899/1‑150 (MQ=255)
                                         cGTAAACAGGGCAACGATAATGCCGAATCCACTAAAGATATCGACGAAGCCTTGCAGCAAGGTAAACATCGACCAGATAAAGATCGACCAGAAGTAGACGCGTTTTGAACCAAAACGGTCCAGCAACCAGCCACCAGGGATCTGCCCGAt                                                                                                              <  2:2128696/150‑1 (MQ=255)
                                          gTAAACAGGGCAACGATAATGCCGAATCCACTAAAGATATCGACGAAGCCTTGCAGCAAGGTAAACATCGACCAGATAAAGATCGACCAGAAGTAGACGCGTTTTGAACCAAAACGGTCCAGCAACCAGCCACCAGGGATCTGCCCGATa                                                                                                             <  1:1816665/150‑1 (MQ=255)
                                           tAAACAGGGCAACGATAATGCCGAATCCACTAAAGATATCGACGAAGCCTTGCAGCAAGGTAAACATCGCCCAGATAAAGATCGACCAGAAGTAGACGCGTTTTGAACCAAAACGGTCCAGCAACCAGCCCCCAGGGAACTGCCCGATaa                                                                                                            >  2:2463707/1‑150 (MQ=255)
                                            aaaCAGGGCAACGATAATGCCGAATCCACTAAAGATATCGACGAAGCCTTGCAGCAAGGTAAACATCGACCAGATAAAGATCGACCAGAAGTAGACGCGTTTTGAACCAAAACGGTCCAGCAACCAGCCACCAGGGATCTGCCCGATAAc                                                                                                           >  2:1492407/1‑150 (MQ=255)
                                            aaaCAGGGCAACGATAATGCCGAATCCACTAAAGATATCGACGAAGCCTTGCAGCAAGGTAAACATCGACCAGATAAAGATCGACCAGAAGTAGACGCGTTTTGAACCAAAACGGTCCAGCAACCAGCCACCAGGGATCTGCCCGATAAc                                                                                                           <  1:1527377/150‑1 (MQ=255)
                                            aaaCAGGGCAACGATAATGCCGAATCCACTAAAGATATCGACGAAGCCTTGCAGCAAGGTAAACATCGACCAGATAAAGATCGACCAGAAGTAGACGCGATTAGAACCAAACCGGTCCCGCAACCAGCCCCCAGGGATCTGCCCGATAAc                                                                                                           >  2:1072774/1‑150 (MQ=255)
                                             aaCAGGGCAACTATAATGCCGAATCCACTAAAGATATCGACGAAGCCTTGCAGCAAGGTAAACATCGACCAGATAAAGATCGACCAGAAGTAGACGCGTTTTGAACCAAAACGGTCCAGCAACCAGCCACCAGGGATCTGCCCGATAACa                                                                                                          <  2:1974154/150‑1 (MQ=255)
                                              aCAGGGCAACGATAATGCCGAATCCACTAAAGATATCGACGAAGCCTTGCAGCAAGGTAAACATCGACCAGATAAAGATCGACCAGAAGTAGACGCGTTTTGAACCAAAACGGTCCAGCAACCAGCCACCAGGGATCTGCCCGATAACAt                                                                                                         <  1:2103408/150‑1 (MQ=255)
                                                   gCAACGATAATGCCGAATCCACTAAAGATATCGACGAAGCCTTGCAGCAAGGTAAACATCGACCAGATAAAGATCGACCAGAAGTAGACGCGTTTTGAACCAAAACGGTCCAGCAACCAGCCACCAGGGATCTGCCCGATAACATAAGcc                                                                                                    >  1:1815864/1‑150 (MQ=255)
                                                   gCAACGATAATGCCGAATCCACTAAAGATATCGACGAAGCCTTGCAGCAAGGTAAACATCGACCAGATAAAGATCGACCAGAAGTAGACGCGTTTTGAACCAAAACGGTCCAGCAACCAGCCACCAGGGATCTGCCCGATAACATAAGcc                                                                                                    >  2:116644/1‑150 (MQ=255)
                                                        gATAATGCCGAATCCACTAAAGATATCGACGAAGCCTTGCAGCAAGGTAAACATCGACCAGATAAAGATCGACCAGAAGTAGCCGCGTTTTGAACCAAAACGGTCCAGCAACCAGCCACAAGGGATCTTCCCGATAACATAAGACCATga                                                                                               >  1:1389261/1‑150 (MQ=255)
                                                            aTGCCGAATCCACTAAAGATATCGACGAAGCCTTGCAGCAAGGTAAACATCGACCAGATAAAGATCGACCAGAAGTAGACGCGTTTTGAACCAAAACGGTCCAGCAACCAGCCACCAGGGATCTGCCCGATAACATAAGCCCATGAGaaa                                                                                           <  1:129747/150‑1 (MQ=255)
                                                                cGAATCCACTAAAGATATCGACGAAGCCTTGCAGCAAGGTAAACATCGACCAGATAAAGATCGACCAGAAGTAGACGCGTTTTGAACCAAAACGGTCCAGCAACCAGCCACCAGGGATCTGCCCGATAACATAAGCCCATGAGAAAGCag                                                                                       >  1:1871160/1‑150 (MQ=255)
                                                                cGAATCCACTAAAGATATCGACGAAGCCTTGCAGCAAGGTAAACATCGACCAGATAAAGATCGACCAGAAGTAGACGCGTTTTGAACCAAAACGGTCCAGCAACCAGCCACCAGGGATCTGCCCGATAACATAAGCCCATGAGAAAGCag                                                                                       <  2:1159167/150‑1 (MQ=255)
                                                                cGAATCCACTAAAGATATCGACGAAGCCTTGCAGCAAGGTAAACATCGACCAGATAAAGATCGACCAGAAGTAGACGCGTTTTGAACCAAAACGGTCCAGCAACCAGCCACCAGGGATCTGCCCGATAACATAAGCCCATGAGAAAGCag                                                                                       <  1:1191399/150‑1 (MQ=255)
                                                                cGAATCCACTAAAGATATCGACGAAGCCTTGCAGCAAGGTAAACATCGACCAGATAAAGATCGACCAGAAGTAGACGCGTTTTGAACCAAAACGGTCCAGCAACCAGCCACCAGGGATCTGCCCGATAACATAAGCCCATGAGAAAGCag                                                                                       >  1:1049768/1‑150 (MQ=255)
                                                                   aTCCACTAAAGATATCGACGAAGCCTTGCAGCAAGGTAAACATCGACCAGATAAAGATCGACCAGAAGTAGACGCGTTTTGAACCAAAACGGTCCAGCAACCAGCCACCAGGGATCTGCCCGATAACATAAGCCCATGAGAAAGCAgaga                                                                                    <  2:1953693/150‑1 (MQ=255)
                                                                     ccACTAAAGATATCGACGAAGCCTTGCAGCAAGGTAAACATCGACCAGATAAAGATCGACCAGAAGTAGACGCGTTTTGAACCAAAACGGTCCAGCAACCAGCCACCAGGGATCTGCCCGATAACATAAGCCCATGAGAAAGCAGAGAac                                                                                  <  2:827066/150‑1 (MQ=255)
                                                                      cACTAAAGATATCGACGAAGCCTTGCAGCAAGGTAAACATCGACCAGATAAAGATCGACCAGAAGTAGACGCGTTTTGAACCAAAACGGTCCAGCAACCAGCCACCAGGGATCTGCCCGATAACATAAGCCCATGAGAAAGCAGAGAaca                                                                                 >  2:2108748/1‑150 (MQ=255)
                                                                        cTAAAGATATCGACGAAGCCTTGCAGCAAGGTAAACATCGACCAGATAAAGATCGACCAGAAGTAGACGCGTTTTGAACCAAAACGGTCCAGCAACCAGCCACCAGGGATCTGCCCGATAACATAAGCCCATGAGAAAGCAGAGAAcaca                                                                               <  2:730131/150‑1 (MQ=255)
                                                                         tAAAGATATCGACGAAGCCTTGCAGCAAGGTAAACATCGACCAGATAAAGATCGACCAGAAGTAGACGCGTTTTGAACCAAAACGGTCCAGCAACCAGCCACCAGGGATCTGCCCGATAACATAAGCCCATGAGAAAGCAGAGAACACAt                                                                              <  2:2220300/150‑1 (MQ=255)
                                                                               tatCGACGAAGCCTTGCAGCAAGGTAAACATCGACCAGATAAAGATCGACCAGAAGTAGACGCGTTTTGAACCAAAACGGTCCAGCAACCAGCCACCAGGGATCTGCCCGATAACATAAGCCCATGAGAAAGCAGAGAACACATAGCCCa                                                                        <  1:47184/150‑1 (MQ=255)
                                                                               tatCGACGAAGCCTTGCAGCAAGGTAAACATCGACCAGATAAAGATCGACCAGAAGTAGACGCGTTTTGAACCAAAACGGTCCAGCAACCAGCCACCAGGGATCTGCCCGATAACATAAGCCCATGAGAAAGCAGAGAACACATAGCCCa                                                                        >  1:2170135/1‑150 (MQ=255)
                                                                                atCGACGAAGCCTTGCAGCAAGGTAAACATCGACCAGATAAAGATCGACCAGAAGTAGACGCGTTTTGAACCAAAACGGTCCAGCAACCAGCCACCAGGGATCTGCCCGATAACATAAGCCCATGAGAAAGCAGAGAACACATAGCCCAt                                                                       >  2:2003813/1‑150 (MQ=255)
                                                                                   gacgaAGCCTTGCAGCAAGGTAAACATCGACCAGATAAAGATCGACCAGAAGTAGACGCGTTTTGAACCAAAACGGTCCAGCAACCTGCCACCAGGGATCTGCCCGATAACATAAGCCCATGAGAAAGCAGAGAACACATAGCCCATTcc                                                                    >  2:2362014/1‑150 (MQ=255)
                                                                                   gacgaAGCCTTGCAGCAAGGTAAACATCGACCAGATAAAGATCGACCAGAAGTAGACGCGTTTTGAACCAAAACGGTCCAGCAACCAGCCACCAGGGATCTGCCCGATAACATAAGCCCATGAGAAAGCAGAGAACACATAGCCCATTcc                                                                    <  2:1858411/150‑1 (MQ=255)
                                                                                      gaAGCCTTGCAGCAAGGTAAACATCGACCAGATAAAGATCGACCAGAAGTAGACGCGTTTTGAACCAAAACGGTCCAGCAACCAGCCACCAGGGATCTGCCCGATAACATAAGCCCATGAGAAAGCAGAGAACACATAGCCCATTCCCAc                                                                 >  1:225484/1‑150 (MQ=255)
                                                                                       aaGCCTTGCAGCAAGGTAAACATCGACCAGATAAAGATCGACCAGAAGTAGACGCGTTTTGAACCAAAACGGTCCAGCAACCAGCCACCAGGGATCTGCCCGATAACATAAGCCCATGAGAAAGCAGAGAACACATAGCCCATTCCCACg                                                                <  2:368264/150‑1 (MQ=255)
                                                                                        aGCCTTGCAGCAATGTAAACATCTACCAGATAAAGATCGACCAGAAGTAGACGCGTTTTGAACCAAAACGGTCCAGCAACCAGCCACCAGGGATCTGCCCGATAACATAAGCCCATGAGAAAGCAGAGAACACATAGCCCATTCCCACgg                                                               <  1:310145/150‑1 (MQ=255)
                                                                                          ccTTGCAGCAAGGTAAACATCGACCAGATAAAGATCGACCAGAAGTAGACGCGTTTTGAACCAAAACGGTCCAGCAACCAGCCACCAGGGATCTGCCCGATAACATAAGCCCATGAGAAAGCAGAGAACACATAGCCCATTCCCACGGGa                                                             >  2:279420/1‑150 (MQ=255)
                                                                                            ttGCAGCAAGGTAAACATCGACCAGATAAAGATCGACCAGAAGTAGACGCGTTTTGAACCAAAACGGTCCAGCAACCAGCCACCAGGGATCTGCCCGATAACATAAGCCCATGAGAAAGCAGAGAACACATAGCCCATTCCCACGGGATc                                                           <  2:312694/150‑1 (MQ=255)
                                                                                             tGCAGCAAGGTAAACATCGACCAGATAAAGATCGACCAGAAGTAGACGCGTTTTGAACCAAAACGGTCCAGCAACCAGCCACCAGGGATCTGCCCGATAACATAAGCCCATGAGAAAGCAGAGAACACATAGCCCATTCCCACGGGATCa                                                          <  1:853804/150‑1 (MQ=255)
                                                                                               cagcaAGGTAAACATCGACCAGATAAAGATCGACCAGAAGTAGACGCGTTTTGAACCAAAACGGTCCAGCAACCAGCCACCAGGGATCTGCCCGATAACATAAGCCCATGAGAAAGCAGAGAACACATAGCCCATTCCCACGGGATCAAg                                                        <  2:7125/150‑1 (MQ=255)
                                                                                                  caAGGTAAACATCGACCAGATAAAGATCGACCAGAAGTAGACGCGTTTTGAACCAAAACGGTCCAGCAACCAGCCACCAGGGATCTGCCCGATAACATAAGCCCATGAGAAAGCAGAGAACACATAGCCCATTCCCACGGGATCAAGGcc                                                     <  1:1492407/150‑1 (MQ=255)
                                                                                                  caAGGTAAACATCGACCAGATAAAGATCGACCAGAAGTAGACGCGTTTTGAACCAAAACGGTCCAGCAACCAGCCACCAGGGATCTGCCCGATAACATAAGCCCATGAGAAAGCAGAGAACACATAGCCCATTCCCACGGGATCAAGGcc                                                     <  2:29256/150‑1 (MQ=255)
                                                                                                    aGGTAAACATCGACCAGATAAAGATCGACCAGAAGTAGACGCGTTTTGAACCAAAACGGTCCAGCAACCAGCCACCAGGGATCTGCCCGATAACATAAGCCCATGAGAAAGCAGAGAACACATAGCCCATTCCCACGGGATCAAGGCCGa                                                   <  2:355713/150‑1 (MQ=255)
                                                                                                    aGGTAAACATCGACCAGATAAAGATCGACCAGAAGTAGACGCGTTTTGAACCAAAACGGTCCAGCAACCAGCCACCAGGGATCTGCCCGATAACATAAGCCCATGAGAAAGCAGAGAACACATAGCCCATTCCCACGGGATCAAGGCCGa                                                   <  2:1540886/150‑1 (MQ=255)
                                                                                                             tCGACCAGATAAAGATCGACCAGAAGTAGACGCGTTTTGAACCAAAACGGTCCAGCAACCAGCCACCAGGGATCTGCCCGATAACATAAGCCCATGAGAAAGCAGAGAACACATAGCCCATTCCCACGGGATCAAGGCCGATATCTTTgg                                          <  2:1473842/150‑1 (MQ=255)
                                                                                                             tCGACCAGATAAAGATCGACCAGAAGTAGACGCGTTTTGAACCAAAACGGTCCAGCAACCAGCCACCAGGGATCTGCCCGATAACATAAGCCCATGAGAAAGCAGAGAACACATAGCCCATTCCCACGGGATCAAGGCCGATATCTTTgg                                          <  1:108161/150‑1 (MQ=255)
                                                                                                               gACCAGATAAAGATCGACCAGAAGTAGACGCGTTTTGAACCAAAACGGTCCAGCAACCAGCCACCAGGGATCTGCCCGATAACATAAGCCCATGAGAAAGCAGAGAACACATAGCCCATTCCCACGGGATCAAGGCCGATATCTTTGGcc                                        >  1:2402766/1‑150 (MQ=255)
                                                                                                               gACCAGATAAAGATCGACCAGAAGTAGACGCGTTTTGAACCAAAACGGTCCAGCAACCAGCCACCAGGGATCTGCCCGATAACATAAGCCCATGAGAAAGCAGAGAACACATAGCCCATTCCCACGGGATCAAGGCCGATATCTTTGGcc                                        >  2:1635700/1‑150 (MQ=255)
                                                                                                                          gATCGACCAGAAGTAGACGCGTTTTGAACCAAAACGGTCCAGCAACCAGCCACCAGGGATCTGCCCGATAACATAAGCCCATGAGAAAGCAGAGAACACATAGCCCATTCCCACGGGATCAAGGCCGATATCTTTGGCCATTTCCGAAcc                             >  2:1077138/1‑150 (MQ=255)
                                                                                                                           aTCGACCAGAAGTAGACGCGTTTTGAACCAAAACGGTCCAGCAACCAGCCACCAGGGATCTGCCCGATAACATAAGCCCATGAGAAAGCAGAGAACACATAGCCCATTCCCACGGGATCAAGGCCGATATCTTTGGCCATTTCCGAACCg                            <  1:52442/150‑1 (MQ=255)
                                                                                                                             cGACCAGAAGTAGACGCGTTTTGAACCAAAACGGTCCAGCAACCAGCCACCAGGGATCTGCCCGATAACATAAGCCCATGAGAAAGCAGAGAACACATAGCCCATTCCCACGGGATCAAGGCCGATATCTTTGGCCATTTCCGAACCGGc                          <  1:2471381/150‑1 (MQ=255)
                                                                                                                             cGACCAGAAGTAGACGCGTTTTGAACCAAAACGGTCCAGCAACCAGCCACCAGGGATCTGCCCGATAACATAAGCCCATGAGAAAGCAGAGAACACATAGCCCATTCCCACGGGATCAAGGCCGATATCTTTGGCCATTTCCGAACCGGc                          <  1:256753/150‑1 (MQ=255)
                                                                                                                             cGACCAGAAGTAGACGCGTTTTGAACCAAAACGGTCCAGCAACCAGCCACCAGGGATCTGCCCGATAACATAAGCCCATGAGAAAGCAGAGAAACATAGCCCATTCCCACGGGATCAAGGCCGATATCTTTGGCCATTTCCGAACCGGCg                         <  1:297637/150‑1 (MQ=255)
                                                                                                                                ccAGAAGTAGACGCGTTTTGAACCAAAACGGTCCAGCAACCAGCCACCAGGGATCTGCCCGATAACATAAGCCCATGAGAAAGCAGAGAACACATAGCCCATTCCCACGGGATCAAGGCCGATATCTTTGGCCATTTCCGAACCGGCGAt                       <  1:1467773/150‑1 (MQ=255)
                                                                                                                                  aGAAGTAGACGCGTTTTGAACCAAAACGGTCCAGCAACCAGCCACCAGGGATCTGCCCGATAACATAAGCCCATGAGAAAGCAGAGAACACATAGCCCATTCCCACGGGATCAAGGCCGATATCTTTGGCCATTTCCGaa                               <  2:2085272/140‑1 (MQ=255)
                                                                                                                                  aGAAGTAGACGCGTTTTGAACCAAAACGGTCCAGCAACCAGCCACCAGGGATCTGCCCGATAACATAAGCCCATGAGAAAGCAGAGAACACATAGCCCATTCCCACGGGATCAAGGCCGATATCTTTGGCCATTTCCGaa                               >  1:2085272/1‑140 (MQ=255)
                                                                                                                                  aGAAGTAGACGCGTTTTGAACCAAAACGGTCCAGCAACCAGCCACCAGGGATCTGCCCGATAACATAAGCCCATGAGAAAGCAGAGAACACATAGCCCATTCCCACGGGATCAAGGCCGATATCTTTGGCCATTTCCGAACCGGCGATag                     <  1:483733/150‑1 (MQ=255)
                                                                                                                                         ggaccgTTTTGGTTCAAAACGGTCCAGCAACCAGCCACCAGGGATCTGCCCGATAACATAAGCCCATGAGAAAGCAGAGAACACATAGCCCATTCCCACGGGATCAAGGCCGATATCTTTGGCCATTTCCGAACCGGCGATAGAGAGCGt              >  1:762737/5‑150 (MQ=255)
                                                                                                                                         ggaccgTTTTGGTTCAAAACGGTCCAGCAACCAGCCACCAGGGATCTGCCCGATAACATAAGCCCATGAGAAAGCAGAGAACACATAGCCCATTCCCACGGGATCAAGGCCGATATCTTTGGCCATTTCCGAACCGGCGATAGAGAGCGt              >  1:2442405/5‑150 (MQ=255)
                                                                                                                                         gACGCGTTTTGAACCAAAACGGTCCAGCAACCAGCCACCAGGGATCTGCCCGATAACATAAGCCCATGAGAAAGCAGAGAACACATAGCCCATTCCCACGGGATCAAGGCCGATATCTTTGGCCATTTCCGAACCGGCGATAGAGAGCGt              >  1:1308835/1‑150 (MQ=255)
                                                                                                                                          gaccgTTTTGGTTCAAAACGGTCCAGCAACCAGCCACCAGGGATCTGCCCGATAACATAAGCCCATGAGAAAGCAGAGAACACATAGCCCATTCCCACGGGATCAAGGCCGATATCTTTGGCCATTTCCGAACCGGCGATAGAGAGCGTa             >  2:583577/4‑150 (MQ=255)
                                                                                                                                              gTTTTGAACCAAAACGGTCCAGCAACCAGCCACCAGGGATCTGCCCGATAACATAAGCCCATGAGAAAGCAGAGAACACATAGCCCATTCCCACGGGATCAAGGCCGATATCTTTGGCCATTTCCGAACCGGCGATAGAGAGCGTAgcgc         >  1:1951313/1‑150 (MQ=255)
                                                                                                                                                 ttGAACCAAAACGGTCCAGCAACCAGCCACCAGGGATCTGCCCGATTACATAAGCCCATGAGAAAGCAGAGAACACATAGCCCATTCCCACGGGATCAAGGCCGATATCTTTGGCCATTTCCGAACCGGCGATAGAGAGCGTAGCGCggg      >  2:1677943/1‑149 (MQ=255)
                                                                                                                                                  tGAACCAAAACGGTCCAGCAACCAGCCACCAGGGATCTGCCCGATAACATAAGCCCATGAGAAAGCAGAGAACACATAGCCCATTCCCACGGGATCAAGGCCGATATCTTTGGCCATTTCCGAACCGGCGATAGAGAGCGTAGCGCGGTc     >  2:2151771/1‑150 (MQ=255)
                                                                                                                                                    aaCCAAAACGGTCCAGCAACCAGCCACCAGGGATCTGCCCGATAACATAAGCCCATGAGAAAGCAGAGAACACATAGCCCATTCCCACGGGATCAAGGCCGATATCTTTGGCCATTTCCGAACCGGCGATAGAGAGCGTAGCGCGGTCGc   <  2:2402766/150‑1 (MQ=255)
                                                                                                                                                     aCCAAAACGGTCCAGCAACCAGCCACCAGGGATCTGCCCGATAACATAAGCCCATGAGAAAGCAGAGAACACATAGCCCATTCCCACGGGATCAAGGCCGATATCTTTGGCCATTTCCGAACCGGCGATAGAGAGCGTAGCGCGGTCGcc  <  2:144938/150‑1 (MQ=255)
                                                                                                                                                     aCCAAAACGGTCCAGCAACCAGCCACCAGGGATCTGCCCGATAACATAAGCCCATGAGAAAGCAGAGAACACATAGCCCATTCCCACGGGATCAAGGCCGATATCTTTGGCCATTTCCGAACCGGCGATAGAGAGCGTAGCGCGGTCGcc  >  1:2482084/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GGGAAAGATGGCGCTTCAGCAAGCCCGACCAGGAAGCGCAGCGTAAACAGGGCAACGATAATGCCGAATCCACTAAAGATATCGACGAAGCCTTGCAGCAAGGTAAACATCGACCAGATAAAGATCGACCAGAAGTAGACGCGTTTTGAACCAAAACGGTCCAGCAACCAGCCACCAGGGATCTGCCCGATAACATAAGCCCATGAGAAAGCAGAGAACACATAGCCCATTCCCACGGGATCAAGGCCGATATCTTTGGCCATTTCCGAACCGGCGATAGAGAGCGTAGCGCGGTCGCC  >  NC_000913/2921703‑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: