Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,055,576 C→A 9.3% V456V (GTC→GTA scpA → methylmalonyl‑CoA mutase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,055,5760CA 9.3% 565.4 / 25.3 161V456V (GTC→GTAscpAmethylmalonyl‑CoA mutase
Reads supporting (aligned to +/- strand):  ref base C (69/77);  new base A (9/6);  total (78/83)
Fisher's exact test for biased strand distribution p-value = 4.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-01

Only 100 of 174 total aligned reads displayed.

CAAGCCAGAGCTATTATCCAACAGATCGACGAAGCCGGTGGCATGGCGAAAGCGATCGAAGCAGGTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGTACTTGAGATCGACAACGTGATGGTGCGTAACGAGCAAATTGCTTCGCTGGAACGCATTCGCGCCACCCGTGATGATGCCGCCGTAACCGCCGCGTTGAACGCCCTGACT  >  CP009273/3055430‑3055723
                                                                                                                                                  |                                                                                                                                                   
cAAGCCAGAGCTATTATCCAACAGATCGACGAAGCCGGTGGCATGGCGAAAGCGATCGAAGCAGGTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTca                                                                                                                                                    >  1:2861431/1‑148 (MQ=255)
cAAGCCAGAGCTATTATCCAACAGATCGACGAAGCCGGTGGCATGGCGAAAGCGATCGAAGCAGGTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTca                                                                                                                                                    >  2:3025692/1‑148 (MQ=255)
 aaGCCAGAGCTATTATCCAACAGATCGACGAAGCCGGTGGCATGGCGAAAGCGATCGAAGCAGGTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTcaa                                                                                                                                                   <  1:2517071/148‑1 (MQ=255)
   gCCAGAGCTAGTATCCAACAGATCGACGAAGCCGGTGGCATGGCGAAAGCGATCGAAGCAGGTCTGCCATAACGAAGGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCACGCGGGTCCTCGTTGGTGGcacca                                                                                                                                                 >  1:297040/1‑148 (MQ=255)
    ccAGAGCTATTATCCAACAGATCGACGAAGCCGGTGGCATGGCGAAAGCGATCGAAGCAGGTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTcaacaa                                                                                                                                                >  2:327031/1‑148 (MQ=255)
          cTATTATCCAACAGATCGACGAAGCCGGTGGCATGGCGAAAGCGATCGAAGCAGGTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACaa                                                                                                                                          >  1:277141/1‑148 (MQ=255)
          cTATTATCCAACAGATCGACGAAGCCGGTGGCATGGCGAAAGCGATCGAAGCAGGTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACaa                                                                                                                                          >  2:1222956/1‑148 (MQ=255)
            attatCCAACAGATCGACGAAGCCGGTGGCATGGCGAAAGCGATCGAAGCAGGTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAAc                                                                                                                                        >  1:2456799/1‑148 (MQ=255)
            attatCCAACAGATCGACGAAGCCGGTGGCATGGCGAAAGCGATCGAAGCAGGTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTAAACAAGTACAAAc                                                                                                                                        >  2:2822903/1‑148 (MQ=255)
            attatCCAACAGATCGACGAAGCCGGTGGCATGGCGAAAGCGATCGAAGCAGGTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTAAACAAGTACAAAc                                                                                                                                        >  1:2895620/1‑148 (MQ=255)
             ttatCCAACAGATCGACGAAGCCGGTGGCATGGCGAAAGCGATCGAAGCAGGTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTAAACAAGTACAAACt                                                                                                                                       >  1:1903840/1‑148 (MQ=255)
              tatCCAACAGATCGACGAAGCCGGTGGCATGGCGAAAGCGATCGAAGCAGGTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTg                                                                                                                                      >  1:2030857/1‑148 (MQ=255)
                   aaCAGATCGACGAAGCCGGTGGCATGGCGAAAGCGATCGAAGCAGGTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCa                                                                                                                                 <  2:2331302/148‑1 (MQ=255)
                       gATCGACGAAGCCGGTGGCATGGCGAAAGCGATCGAAGCAGGTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGaa                                                                                                                             >  2:168481/1‑148 (MQ=255)
                       gATCGACGAAGCCGGTGGCATGGCGAAAGCGATCGAAGCAGGTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTAAACAAGTACAAACTGGATCACGaa                                                                                                                             >  1:1372775/1‑148 (MQ=255)
                        aTCGACGAAGCCGGTGGCATGGCGAAAGCGATCGAAGCAGGTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAg                                                                                                                            <  2:268882/148‑1 (MQ=255)
                        aTCGACGAAGCCGGTGGCATGGCGAAAGCGATCGAAGCAGGTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAg                                                                                                                            >  1:1882358/1‑148 (MQ=255)
                        aTCGACGAAGCCGGTGGCATGGCGAAAGCGATCGAAGCAGGTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAg                                                                                                                            <  1:291461/148‑1 (MQ=255)
                        aTCGACGAAGCCGGTGGCATGGCGAAAGCGATCGAAGCAGGTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAg                                                                                                                            <  2:1918204/148‑1 (MQ=255)
                         tCGACGAAGCCGGTGGCATGGCGAAAGCGATCGAAGCAGGTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGa                                                                                                                           <  1:1873240/148‑1 (MQ=255)
                            acgaAGCCGGTGGCATGGCGAAAGCGATCGAAGCAGGTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGa                                                                                                                        >  2:2531562/1‑148 (MQ=255)
                            acgaAGCCGGTGGCATGGCGAAAGCGATCGAAGCAGGTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGa                                                                                                                        >  2:824147/1‑148 (MQ=255)
                              gaAGCCGGTGGCATGGCGAAAGCGATCGAAGCAGGTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGaaa                                                                                                                      <  2:497472/148‑1 (MQ=255)
                              gaAGCCGGTGGCATGGCGAAAGCGATCGAAGCAGGTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGaaa                                                                                                                      >  1:1235700/1‑148 (MQ=255)
                               aaGCCGGTGGCATGGCGAAAGCGATCGAAGCAGGTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAAc                                                                                                                     >  1:2944503/1‑148 (MQ=255)
                               aaGCCGGTGGCATGGCGAAAGCGATCGAAGCAGGTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTAAAAACTGGATCACGAAGACGAAAc                                                                                                                     >  1:2197918/1‑148 (MQ=255)
                                aGCCGGTGGCATGGCGAAAGCGATCGAAGCAGGTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAAcc                                                                                                                    >  2:1388260/1‑148 (MQ=255)
                                  ccGGTGGCATGGCGAAAGCGATCGAAGCAGGTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGa                                                                                                                  <  2:2279814/148‑1 (MQ=255)
                                   cGGTGGCATGGCGAAAGCGATCGAAGCAGGTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGAt                                                                                                                 <  1:1007828/148‑1 (MQ=255)
                                   cGGTGGCATGGCGAAAGCGATCGAAGCAGGTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGAt                                                                                                                 <  1:1007924/148‑1 (MQ=255)
                                        gCATGGCGAAAGCGATCGAAGCAGGTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGTACt                                                                                                            <  1:2636421/148‑1 (MQ=255)
                                        gCATGGCGAAAGCGATCGAAGCAGGTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGTACt                                                                                                            <  1:981135/148‑1 (MQ=255)
                                             gCGAAAGCGATCGAAGCAGGTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGTACTTgaga                                                                                                       <  2:3084031/148‑1 (MQ=255)
                                                 aaGCGATCGAAGCAGGTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGTACTTGAGATCGa                                                                                                   >  1:2435003/1‑148 (MQ=255)
                                                  aGCGATCGAAGCAGGTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGTACTTGAGATCGAc                                                                                                  >  2:880205/1‑148 (MQ=255)
                                                       tCGAAGCAGGTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGTACTTGAGATCGACAACGt                                                                                             >  1:2559805/1‑148 (MQ=255)
                                                       tCGAAGCAGGTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGTACTTGAGATCGACAACGt                                                                                             >  1:2559510/1‑148 (MQ=255)
                                                        cGAAGCAGGTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGTACTTGAGATCGACAACGTg                                                                                            <  2:2748314/148‑1 (MQ=255)
                                                        cGAAGCAGGTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGTACTTGAGATCGACAACGTg                                                                                            <  2:2745528/148‑1 (MQ=255)
                                                         gAAGCAGGTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGTACTTGAGATCGACAACGTGa                                                                                           <  1:2002895/148‑1 (MQ=255)
                                                          aaGCAGGTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGTACTTGAGATCGACAACGTGAt                                                                                          <  2:2913808/148‑1 (MQ=255)
                                                           aGCAGGTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGTACTTGAGATCGACAACGTGATg                                                                                         >  2:2205911/1‑148 (MQ=255)
                                                           aGCAGGTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGTACTTGAGATCGACAACGTGATg                                                                                         >  2:1410138/1‑148 (MQ=255)
                                                           aGCAGGTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGTACTTGAGATCGACAACGTGATg                                                                                         <  1:1316906/148‑1 (MQ=255)
                                                            gCAGGTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGTACTTGAGATCGACAACGTGATgg                                                                                        >  2:3180902/1‑148 (MQ=255)
                                                              aGGTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGTACTTGAGATCGACAACGTGATGGTg                                                                                      >  2:1471703/1‑148 (MQ=255)
                                                                gTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGTACTTGAGATCGACAACGTGATGGTGCg                                                                                    <  2:277141/148‑1 (MQ=255)
                                                                gTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGTACTTGAGATCGACAACGTGATGGTGCg                                                                                    >  1:1988006/1‑148 (MQ=255)
                                                                gTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGTACTTGAGATCGACAACGTGATGGTGCg                                                                                    <  1:3180902/148‑1 (MQ=255)
                                                                   tGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGTACTTGAGATCGACAACGTGATGGTGCGTaa                                                                                 <  1:574337/148‑1 (MQ=255)
                                                                     ccAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTAAACAAGTACAAACTGGATCACGAAGACGAAACAGATGTACTTGAGATCGACAACGTTATGGTGCGTAACg                                                                               >  1:2661636/1‑148 (MQ=255)
                                                                       aaaaCGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGTACTTGAGATCGACAACGTGATGGTGCGTAACGAg                                                                             <  2:1485183/148‑1 (MQ=255)
                                                                       aaaaCGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGTACTTGAGATCGACAACGTGATGGTGCGTAACGAg                                                                             <  2:836250/148‑1 (MQ=255)
                                                                        aaaCGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGTACTTGAGATCGACAACGTGATGGTGCGTAACGAGc                                                                            >  1:1097855/1‑148 (MQ=255)
                                                                         aaCGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGTACTTGAGATCGACAACGTGATGGTGCGTAACGAGCa                                                                           >  1:2694064/1‑148 (MQ=255)
                                                                          aCGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGTACTTGAGATCGACAACGTGATGGTGCGTAACGAGCaa                                                                          <  2:582758/148‑1 (MQ=255)
                                                                          aCGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGTACTTGAGATCGACAACGTGATGGTGCGTAACGAGCaa                                                                          <  1:1290247/148‑1 (MQ=255)
                                                                          aCGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGTACTTGAGATCGACAACGTGATGGTGCGTAACGAGCaa                                                                          >  2:2297864/1‑148 (MQ=255)
                                                                          aCGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGTACTTGAGATCGACAACGTGATGGTGCGTAACGAGCaa                                                                          >  2:637970/1‑148 (MQ=255)
                                                                           cGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGTACTTGAGATCGACAACGTGATGGTGCGTAACGAGCaaa                                                                         <  2:1163074/148‑1 (MQ=255)
                                                                           cGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGTACTTGAGATCGACAACGTGATGGTGCGTAACGAGCaaa                                                                         <  1:2214887/148‑1 (MQ=255)
                                                                              aTGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGTACTTGAGATCGACAACGTGATGGTGCGTAACGAGCAAATTg                                                                      <  1:3125220/148‑1 (MQ=255)
                                                                               tGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGTACTTGAGATCGACAACGTGATGGTGCGTAACGAGCAAATTGc                                                                     <  1:2218889/148‑1 (MQ=255)
                                                                                 aTCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTAAACAAGTACAAACTGGATCACGAAGACGAAACCGATGTACTTGAGATCGACAACGTGATGGTGCGTAACGAGCAAATTGCtt                                                                   <  2:1903840/148‑1 (MQ=255)
                                                                                 aTCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTAAACAAGTACAAACTGGATCACGAAGACGAAACCGATGTACTTGAGATCGACAACGTGATGGTGCGTAACGAGCAAATTGCtt                                                                   <  2:1905439/148‑1 (MQ=255)
                                                                                  tCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGTACTTGAGATCGACAACGTGATGGTGCGTAACGAGCAAATTGCTTc                                                                  >  2:256824/1‑148 (MQ=255)
                                                                                     aaGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGTACTTGAGATCGACAACGTGATGGTGCGTAACGAGCAAATTGCTTCGCt                                                               <  1:1469843/148‑1 (MQ=255)
                                                                                      agagGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGTACTTGAGATCGACAACGTGATGGTGCGTAACGAGCAAATTGCTTCGCTg                                                              >  2:2630744/1‑148 (MQ=255)
                                                                                       gagGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGTACTTGAGATCGACAACGTGATGGTGCGTAACGAGCAAATTGCTTCGCTgg                                                             >  1:3218197/1‑148 (MQ=255)
                                                                                         ggCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGTACTTGAGATCGACAACGTGATGGTGCGTAACGAGCAAATTGCTTCGCTGGaa                                                           <  2:3143830/148‑1 (MQ=255)
                                                                                         ggCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTAAACAAGTACAAACTGGATCACGAAGACGAAACCGATGTACTTGAGATCGACAACGTGATGGTGCGTAACGAGCAAATTGCTTCGCTGGaa                                                           <  2:1337498/148‑1 (MQ=255)
                                                                                         ggCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTAAACAAGTACAAACTGGATCACGAAGACGAAACCGATGTACTTGAGATCGACAACGTGATGGTGCGTAACGAGCAAATTGCTTCGCTGGaa                                                           <  1:2510877/148‑1 (MQ=255)
                                                                                           ccTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTAAACAAGTACAAACTGGATCACGAAGACGAAACCGATGTACTTGAGATCGACAACGTGATGGTGCGTAACGAGCAAATTGCTTCGCTGGAACg                                                         <  2:1372775/148‑1 (MQ=255)
                                                                                             tCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGTACTTGAGATCGACAACGTGATGGTGCGTAACGAGCAAATTGCTTCGCTGGAACGCa                                                       >  2:1603418/1‑148 (MQ=255)
                                                                                             tCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGTACTTGAGATCGACAACGTGATGGTGCGTAACGAGCAAATTGCTTCGCTGGAACGCa                                                       <  2:2796917/148‑1 (MQ=255)
                                                                                              cAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGTACTTGAGATCGACAACGTGATGGTGCGTAACGAGCAAATTGCTTCGCTGGAACGCAt                                                      <  2:2091520/148‑1 (MQ=255)
                                                                                              cAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGTACTTGAGATCGACAACGTGATGGTGCGTAACGAGCAAATTGCTTCGCTGGAACGCAt                                                      >  1:1790821/1‑148 (MQ=255)
                                                                                               aGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGTACTTGAGATCGACAACGTGATGGTGCGTAACGAGCAAATTGCTTCGCTGGAACGCAtt                                                     <  1:796628/148‑1 (MQ=255)
                                                                                                 cgcgcgAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGTACTTGAGATCGACAACGTGATGGTGCGTAACGAGCAAATTGCTTCGCTGGAACGCATTcg                                                   >  2:3245221/1‑148 (MQ=255)
                                                                                                       aaCAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGTACTTGAGATCGACAACGTGATGGTGCGTAACGAGCAAATTGCTTCGCTGGAACGCATTCGCGCCAc                                             >  2:1304906/1‑148 (MQ=255)
                                                                                                       aaCAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGTACTTGAGATCGACAACGTGATGGTGCGTAACGAGCAAATTGCTTCGCTGGAACGCATTCGCGCCAc                                             >  2:734602/1‑148 (MQ=255)
                                                                                                       aaCAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGTACTTGAGATCGACAACGTGATGGTGCGTAACGAGCAAATTGCTTCGCTGGAACGCATTCGCGCCAc                                             >  2:29002/1‑148 (MQ=255)
                                                                                                        aCAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGTACTTGAGATCGACAACGTGATGGTGCGTAACGAGCAAATTGCTTCGCTGGAACGCATTCGCGCCAcc                                            >  1:2063547/1‑148 (MQ=255)
                                                                                                        aCAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGTACTTGAGATCGACAACGTGATGGTGCGTAACGAGCAAATTGCTTCGCTGGAACGCATTCGCGCCAcc                                            >  2:2887905/1‑148 (MQ=255)
                                                                                                         cAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGTACTTGAGATCGACAACGTGATGGTGCGTAACGAGCAAATTGCTTCGCTGGAACGCATTCGCGCCAccc                                           >  2:1936207/1‑148 (MQ=255)
                                                                                                           gTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGTACTTGAGATCGACAACGTGATGGTGCGTAACGAGCAAATTGCTTCGCTGGAACGCATTCGCGCCACCCGt                                         >  1:687393/1‑148 (MQ=255)
                                                                                                                           ggCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGTACTTGAGATCGACAACGTGATGGTGCGTAACGAGCAAATTGCTTCGCTGGAACGCATTCGCGCCACCCGTGATGATGCCGCCGTaa                         <  1:1077607/148‑1 (MQ=255)
                                                                                                                             cAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGTACTTGAGATCGACAACGTGATGGTGCGTAACGAGCAAATTGCTTCGCTGGAACGCATTCGCGCCACCCGTGATGATGCCGCCGTAAcc                       >  1:120989/1‑148 (MQ=255)
                                                                                                                               aGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGTACTTGAGATCGACAACGTGATGGTGCGTAACGAGCAAATTGCTTCGCTGGAACGCATTCGCGCCACCCGTGATGATGCCGCCGTAAccgc                     >  2:2250412/1‑148 (MQ=255)
                                                                                                                                   tgtCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGTACTTGAGATCGACAACGTGATGGTGCGTAACGAGCAAATTGCTTCGCTGGAACGCATTCGCGCCACCCGTGATGATGCCGCCGTAACCGCcgcg                 <  1:327031/148‑1 (MQ=255)
                                                                                                                                        tCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGTACTTGAGATCGACAACGTGATGGTGCGTAACGAGCAAATTGCTTCGCTGGAACGCATTCGCGCCACCCGTGATGATGCCGCCGTAACCGCCGCGTTGaa            <  1:2029164/148‑1 (MQ=255)
                                                                                                                                            tGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGTACTTGAGATCGACAACGTGATGGTGCGTAACGAGCAAATTGCTTCGCTGGAACGCATTCGCGCCACCCGTGATGATGCCGCCGTAACCGCCGCGTTGAACGcc        >  2:2891346/1‑148 (MQ=255)
                                                                                                                                            tGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGTACTTGAGATCGACAACGTGATGGTGCGTAACGAGCAAATTGCTTCGCTGGAACGCATTCGCGCCACCCGTGATGATGCCGCCGTAACCGCCGCGTTGAACGcc        <  1:390698/148‑1 (MQ=255)
                                                                                                                                             ggTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGTACTTGAGATCGACAACGTGATGGTGCGTAACGAGCAAATTGCTTCGCTGGAACGCATTCGCGCCACCCGTGATGATGCCGCCGTAACCGCCGCGTTGAACGccc       <  2:3218197/148‑1 (MQ=255)
                                                                                                                                               tgtCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGTACTTGAGATCGACAACGTGATGGTGCGTAACGAGCAAATTGCTTCGCTGGAACGCATTCGCGCCACCCGTGATGATGCCGCCGTAACCGCCGCGTTGAACGCCCTg     <  1:1866719/148‑1 (MQ=255)
                                                                                                                                               tgtCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGTACTTGAGATCGACAACGTGATGGTGCGTAACGAGCAAATTGCTTCGCTGGAACGCATTCGCGCCACCCGTGATGATGCCGCCGTAACCGCCGCGTTGAACGCCCTg     >  2:95558/1‑148 (MQ=255)
                                                                                                                                                 tCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGTACTTGAGATCGACAACGTGATGGTGCGTAACGAGCAAATTGCTTCGCTGGACCGCATTCGCGCCACCCGTGATGATGCCGCCGTAACCGCCGCGTTGAACGCCCTGAc   >  1:1724230/1‑148 (MQ=255)
                                                                                                                                                 tCAACAAGTACAAACTGGATCACGAAGACCAAACCGATGTAAATTAGATCTAAAAAGTGATGTTGCGTAACAAGCAAATTGCTACGATTGAAAGCATTAGCGCCACCCGTGATAATGACTCCGTAACAGACGATATGAACGACCTGAc   >  1:1725358/1‑148 (MQ=255)
                                                                                                                                                  caacaaGTACAAACTGGATCACGAAGACGAAACCGATGTACTTGAGATCGACAACGTGATGGTGCGTAACGAGCAAATTGCTTCGCTGGAACGCATTCGCGCCACCCGTGATGATGCCGCCGTAACCGCCGCGTTGAACGCCCTGACt  <  2:392983/148‑1 (MQ=255)
                                                                                                                                                  caacaaGTACAAACTGGATCACGAAGACGAAACCGATGTACTTGAGATCGACAACGTGATGGTGCGTAACGAGCAAATTGCTTCGCTGGAACGCATTCGCGCCACCCGTGATGATGCCGCCGTAACCGCCGCGTTGAACGCCCTGACt  <  1:168481/148‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
CAAGCCAGAGCTATTATCCAACAGATCGACGAAGCCGGTGGCATGGCGAAAGCGATCGAAGCAGGTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGTACTTGAGATCGACAACGTGATGGTGCGTAACGAGCAAATTGCTTCGCTGGAACGCATTCGCGCCACCCGTGATGATGCCGCCGTAACCGCCGCGTTGAACGCCCTGACT  >  CP009273/3055430‑3055723

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: