Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,055,407 C→T 10.6% S400L (TCG→TTG)  scpA → methylmalonyl‑CoA mutase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,055,4070CT10.6% 802.6 / 61.9 246S400L (TCG→TTG) scpAmethylmalonyl‑CoA mutase
Reads supporting (aligned to +/- strand):  ref base C (105/114);  new base T (11/15);  total (116/130)
Fisher's exact test for biased strand distribution p-value = 6.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 247 total aligned reads displayed.

CACTGCATACCAACGCCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCATCATCCAGGAAGAATCAGAACTCTGCCGCACCGTCGATCCACTGGCCGGATCCTATTACATTGAGTCGCTGACCGATCAAATCGTCAAACAAGCCAGAGCTATTATCCAACAGATCGACGAAGCCGGTGGCATGGCGAAAGCGATCGAAGCAGGTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCA  >  CP009273/3055260‑3055556
                                                                                                                                                   |                                                                                                                                                     
cACTGCATACCAACGCCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCATCATCCAGGAAGAATCAGAACTCTGCCGCACCGTCGATCCACTGGCCGGATCCTATTACATTGAGTCGc                                                                                                                                                     <  1:4119067/150‑1 (MQ=255)
 aCTGCATACCAACGCCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCATCATCCAGGAAGAATCAGAACTCTGCCGCACCGTCGATCCACTGGCCGGATCCTATTACATTGAGTCGCt                                                                                                                                                    <  2:3565615/150‑1 (MQ=255)
   tGCATACCAACGCCTTTGACGAAGCGCTTGGTTTGCCTACCGATTACTCAGCACGCATTGCCCGCAACACCCAGATCATCATCCAGGAAGAATCAGAACTCTGCCGCACCGTCGATCCACTGGCCGGATCCTATTACATTGAGTCGCTGa                                                                                                                                                  <  2:3884092/150‑1 (MQ=255)
    gCATACCAACGCCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCATCATCCAGGAAGAATCAGAACTCTGCCGCACCGTCGATCCACTGGCCGGATCCTATTACATTGAGTCGCTGAc                                                                                                                                                 >  2:2127583/1‑150 (MQ=255)
       tACCAACGCCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCATCATCCAGGAAGAATCAGAACTCTGCCGCACCGTCGATCCACTGGCCGGATCCTATTACATTGAGTCGCTGACCGa                                                                                                                                              <  2:306079/150‑1 (MQ=255)
       tACCAACGCCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCATCATCCAGGAAGAATCAGAACTCTGCCGCACCGTCGATCCACTGGCCGGATCCTATTACATTGAGTCGCTGACCGa                                                                                                                                              <  2:2878571/150‑1 (MQ=255)
        aCCAACGCCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCATCATCCAGGAAGAATCAGAACTCTGCCGCACCGTCGATCCACTGGCCGGATCCTATTACATTGAGTCGCTGACCGAt                                                                                                                                             <  1:2645684/150‑1 (MQ=255)
        aCCAACGCCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCATCATCCAGGAAGAATCAGAACTCTGCCGCACCGTCGATCCACTGGCCGGATCCTATTACATTGAGTCGCTGACCGAt                                                                                                                                             <  1:2576893/150‑1 (MQ=255)
        aCCAACGCCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCATCATCCAGGAAGAATCAGAACTCTGCCGCACCGTCGATCCACTGGCCGGATCCTATTACATTGAGTCGCTGACCGAt                                                                                                                                             >  1:1632859/1‑150 (MQ=255)
           aaCGCCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCATCATCCAGGAAGAATCAGAACTCTGCCGCACCGTCGATCCACTGGCCGGATCCTATTACATTGAGTCGCTGACCGATCaa                                                                                                                                          <  2:1492597/150‑1 (MQ=255)
             cGCCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCATCATCCAGGAAGAATCAGAACTCTGCCGCACCGTCGATCCACTGGCCGGATCCTATTACATTGAGTTGCTGACCGATCAAAt                                                                                                                                        <  2:3500705/150‑1 (MQ=255)
              gCCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCATCATCCAGGAAGAATCAGAACTCTGCCGCACCGTCGATCCACTGGCCGGATCCTATTACATTGAGTCGCTGACCGATCAAATc                                                                                                                                       >  2:2771213/1‑150 (MQ=255)
               ccTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCATCATCCAGGAAGAATCAGAACTCTGCCGCACCGTCGATCCACTGGCCGGATCCTATTACATTGAGTCGCTGACCGATCAAATCg                                                                                                                                      >  1:1773664/1‑150 (MQ=255)
                     aCGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCATCATCCAGGAAGAATCAGAACTCTGCCGCACCGTCGATCCACTGGCCGGATCCTATTACATTGAGTCGCTGACCGATCAAATCGTCAAAc                                                                                                                                >  1:1361013/1‑150 (MQ=255)
                     aCGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCATCATCCAGGAAGAATCAGAACTCTGCCGCACCGTCGATCCACTGGCCGGATCCTATTACATTGAGTCGCTGACCGATCAAATCGTCAAAc                                                                                                                                <  2:50212/150‑1 (MQ=255)
                     aCGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCATCATCCAGGAAGAATCAGAACTCTGCCGCACCGTCGATCCACTGGCCGGATCCTATTACATTGAGTCGCTGACCGATCAAATCGTCAAAc                                                                                                                                <  1:2032303/150‑1 (MQ=255)
                       gAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCATCATCCAGGAAGAATCAGAACTCTGCCGCACCGTCGATCCACTGGCCGGATCCTATTACATTGAGTCGCTGACCGATCAAATCGTCAAACaa                                                                                                                              >  1:638282/1‑150 (MQ=255)
                       gAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCATCATCCAGGAAGAATCAGAACTCTGCCGCACCGTCGATCCACTGGCCGGATCCTATTACATTGAGTCGCTGACCGATCAAATCGTCAAACaa                                                                                                                              >  2:2749878/1‑150 (MQ=255)
                         aGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCATCATCCAGGAAGAATCAGAACTCTGCCGCACCGTCGATCCACTGGCCGGATCCTATTACATTGAGTCGCTGACCGATCAAATCGTCAAACAAGc                                                                                                                            <  2:3130299/150‑1 (MQ=255)
                           cgcTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCATCATCCAGGAAGAATCAGAACTCTGCCGCACCGTCGATCCACTGGCCGGATCCTATTACATTGAGTCGCTGACCGATCAAATCGTCAAACAAGCCa                                                                                                                          <  1:1497356/150‑1 (MQ=255)
                            gcTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCATCATCCAGGAAGAATCAGAACTCTGCCGCACCGTCGATCCACTGGCCGGATCCTATTACATTGAGTCGCTGACCGATCAAATCGTCAAACAAGCCag                                                                                                                         >  1:3969902/1‑150 (MQ=255)
                            gcTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCATCATCCAGGAAGAATCAGAACTCTGCCGCACCGTCGATCCACTGGCCGGATCCTATTACATTGAGTCGCTGACCGATCAAATCGTCAAACAAGCCag                                                                                                                         >  1:594542/1‑150 (MQ=255)
                             cTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCATCATCCAGGAAGAATCAGAACTCTGCCGCACCGTCGATCCACTGGCCGGATCCTATTACATTGAGTCGCTGACCGATCAAATCGTCAAACAAGCCaga                                                                                                                        >  2:1624794/1‑150 (MQ=255)
                                ggTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCATCATCCAGGAAGAATCAGAACTCTGCCGCACCGTCGATCCACTGGCCGGATCCTATTACATTGAGTCGCTGACCGATCAAATCGTCAAACAAGCCAGAGCt                                                                                                                     >  2:3929662/1‑150 (MQ=255)
                                ggTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCATCATCCAGGAAGAATCAGAACTCTGCCGCACCGTCGATCCACTGGCCGGATCCTATTACATTGAGTCGCTGACCGATCAAATCGTCAAACAAGCCAGAGCt                                                                                                                     <  1:3483675/150‑1 (MQ=255)
                                ggTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCATCATCCAGGAAGAATCAGAACTCTGCCGCACCGTCGATCCACTGGCCGGATCCTATTACATTGAGTCGCTGACCGATCAAATCGTCAAACAAGCCAGAGCt                                                                                                                     <  1:3412481/150‑1 (MQ=255)
                                ggTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCATCATCCAGGAAGAATCAGAACTCTGCCGCACCGTCGATCCACTGGCCGGATCCTATTACATTGAGTCGCTGACCGATCAAATCGTCAAACAAGCCAGAGCt                                                                                                                     >  2:1831208/1‑150 (MQ=255)
                                ggTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCATCATCCAGGAAGAATCAGAACTCTGCCGCACCGTCGATCCACTGGCCGGATCCTATTACATTGAGTCGCTGACCGATCAAATCGTCAAACAAGCCAGAGCt                                                                                                                     >  2:1042397/1‑150 (MQ=255)
                                   ttGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCATCATCCAGGAAGAATCAGAACTCTGCCGCACCGTCGATCCACTGGCCGGATCCTATTACATTGAGTCGCTGACCGATCAAATCGTCAAACAAGCCAGAGCtatt                                                                                                                  <  2:3089547/150‑1 (MQ=255)
                                   ttGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCATCATCCAGGAAGAATCAGAACTCTGCCGCACCGTCGATCCACTGGCCGGATCCTATTACATTGAGTCGCTGACCGATCAAATCGTCAAACAAGCCAGAGCtatt                                                                                                                  <  1:3925009/150‑1 (MQ=255)
                                        tACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCATCATCCAGGAAGAATCAGAACTCTGCCGCACCGTCGATCCACTGGCCGGATCCTATTACATTGAGTCGCTGACCGATCAAATCGTCAAACAAGCCAGAGCTATTATCCa                                                                                                             >  2:424684/1‑150 (MQ=255)
                                           cGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCATCATCCAGGAAGAATCAGAACTCTGCCGCACCGTCGATCCACTGGCCGGATCCTATTACATTGAGTCGCTGACCGATCAAATCGTCAAACAAGCCAGAGCTATTATCCAACa                                                                                                          <  1:3504503/150‑1 (MQ=255)
                                               ttCTCAGCACGCATTGCCCGCAACACCCAGATCATCATCCAGGAAGAATCAGAACTCTGCCGCACCGTCGATCCACTGGCCGGATCCTATTACATTGAGTCGCTGACCGATCAAATCGTCAAACAAGCCAGAGCTATTATCCAACAGATc                                                                                                      <  2:21902/150‑1 (MQ=255)
                                                     gCACGCATTGCCCGCAACACCCAGATCATCATCCAGGAAGAATCAGAACTCTGCCGCACCGTCGATCCACTGGCCGGATCCTATTACATTGAGTCGCTGACCGATCAAATCGTCAAACAAGCCAGAGCTATTATCCAACAGATCGACGaa                                                                                                >  1:2975568/1‑150 (MQ=255)
                                                      cACGCATTGCCCGCAACACCCAGATCATCATCCATGTAGAATCATAACTCTGCCGCACCGTCGATCCACTGGCCGGATCCTATTACATTGAGTCGCTGACCGATCAAATCGTCAAACAAGCCAGAGCTATTATCCAACAGATCGACGAAg                                                                                               <  1:4006574/150‑1 (MQ=255)
                                                      cACGCATTGCCCGCAACACCCAGATCATCATCCAGGAAGAATCAGAACTCTGCCGCACCGTCGATCCACTGGCCGGATCCTATTACATTGAGTCGCTGACCGATCAAATCGTCAAACAAGCCAGAGCTATTATCCAACAGATCGACGAAg                                                                                               >  2:649027/1‑150 (MQ=255)
                                                          cATTGACCGCAACACCCTGATCATCATCCTGGAAGAATCAGAACTCTGCCGCACCGTCGATCCACTGGCCGGATCCTATTTAATTGTGTTGCTGACCGATCAAATCGTCATACTAGCAAGAGCTATTATCCTACATATCGACGAAGCCgg                                                                                           <  2:1013678/150‑1 (MQ=255)
                                                                ccGCAACACCCAGATCATCATCCAGGAAGAATCAGAACTCTGCCGCACCGTCGATCCACTGGCCGGATCCTATTACATTGAGTCGCTGACCGATCAAATCGTCAAACAAGCCAGAGCTATTATCCAACAGATCGACGAAGCCGGTGGCAt                                                                                     <  2:3722459/150‑1 (MQ=255)
                                                                ccGCAACACCCAGATCATCATCCAGGAAGAATCAGAACTCTGCCGCACCGTCGATCCACTGGCCGGATCCTATTACATTGAGTCGCTGACCGATCAAATCGTCAAACAAGCCAGAGCTATTATCCAACAGATCGACGAAGCCGGTGGCAt                                                                                     <  1:472540/150‑1 (MQ=255)
                                                                   cAACACCCAGATCATCCTCCAGGAAGAATCAGAACTCTGCCGCACCGTCGATCCACTGGCCGGATCCTATTACATTGAGTCGCTGACCGATCAAATCGTCAAACAAGCCAGAGCTATTATCCAACAGATCGACGAAGCCGGTGGCATGGc                                                                                  <  2:82216/150‑1 (MQ=255)
                                                                      cacCCAGATCATCATCCAGGAAGAATCAGAACTCTGCCGCACCGTCGATCCACTGGCCGGATCCTATTACATTGAGTCGCTGACCGATCAAATCGTCAAACAAGCCAGAGCTATTATCCAACAGATCGACGAAGCCGGTGGCATGGCGaa                                                                               >  1:1837519/1‑150 (MQ=255)
                                                                       acCCAGATCATCATCCAGGAAGAATCAGAACTCTGCCGCACCGTCGATCCACTGGCCGGATCCTATTACATTGAGTCGCTGACCGATCAAATCGTCAAACAAGCCAGAGCTATTATCCAACAGATCGACGAAGCCGGTGGCATGGCGaaa                                                                              >  1:824583/1‑150 (MQ=255)
                                                                         ccAGATCATCATCCAGGAAGAATCAGAACTCTGCCGCACCGTCGATCCACTGGCCGGATCCTATTACATTGAGTCGCTGACCGATCAAATCGTCAAACAAGCCAGAGCTATTATCCAACAGATCGACGAAGCCGGTGGCATGGCGAAAGc                                                                            <  2:4141323/150‑1 (MQ=255)
                                                                         ccAGATCATCATCCAGGAAGAATCAGAACTCTGCCGCACCGTCGATCCACTGGCCGGATCCTATTACATTGAGTCGCTGACCGATCAAATCGTCAAACAAGCCAGAGCTATTATCCAAAAGATCGACGAAGCCGGTGGCATGTCGAAAGc                                                                            >  2:3367506/1‑150 (MQ=255)
                                                                          cAGATCATCATCCAGGAAGAATCAGAACTCTGCCGCACCGTCGATCCACTGGCCGGATCCTATTACATTGAGTCGCTGACCGATCAAATCGTCAAACAAGCCAGAGCTATTATCCAACAGATCGACGAAGCCGGTGGCATGGCGAAAGCg                                                                           <  2:739608/150‑1 (MQ=255)
                                                                            gATCATCATCCAGGAAGAATCAGAACTCTGCCGCACCGTCGATCCACTGGCCGGATCCTATTACATTGAGTCGCTGACCGATCAAATCGTCAAACAAGCCAGAGCTATTATCCAACAGATCGACGAAGCCGGTGGCATGGCGAAAGCGAt                                                                         <  2:1713659/150‑1 (MQ=255)
                                                                              tcatcatcCAGGAAGAATCAGAACTCTGCCGCACCGTCGATCCACTGGCCGGATCCTATTACATTGAGTTGCTGACCGATCAAATCGTCAAACAAGCCAGAGCTATTATCCAACAGATCGACGAAGCCGGTGGCATGGCGAAAGCGATCg                                                                       <  1:664903/150‑1 (MQ=255)
                                                                              tcatcatcCAGGAAGAATCAGAACTCTGCCGCACCGTCGATCCACTGGCCGGATCCTATTACATTGAGTCGCTGACCGATCAAATCGTCAAACAAGCCAGAGCTATTATCCAACAGATCGACGAAGCCGGTGGCATGGCGAAAGCGATCg                                                                       <  1:2311767/150‑1 (MQ=255)
                                                                              tcatcatcCAGGAAGAATCAGAACTCTGCCGCACCGTCGATCCACTGGCCGGATCCTATTACATTGAGTCGCTGACCGATCAAATCGTCAAACAAGCCAGAGCTATTATCCAACAGATCGACGAAGCCGGTGGCATGGCGAAAGCGATCg                                                                       <  1:1146480/150‑1 (MQ=255)
                                                                                 tcatcCAGGAAGAATCAGAACTCTGCCGCACAGTCGATCCACTGGCCGGATCCTATTCCATTGAGTCGCTGACCGATCAAATCGTCAAACCAGCCAGAGCTATTATCCAACAGATCGACGAAGCCGGTGGCATGGCGAAAGCGATCGAAg                                                                    <  1:2495775/150‑1 (MQ=255)
                                                                                   atcCAGGAAGAATCAGAACTCTGCCGCACCGTCGATCCACTGGCCGGATCCTATTACATTGAGTCGCTGACCGATCAAATCGTCAAACAAGCCAGAGCTATTATCCAACAGATCGACGAAGCCGGTGGCATGGCGAAAGCGATCGAAGCa                                                                  >  1:291653/1‑150 (MQ=255)
                                                                                    tcCAGGAAGAATCAGAACTCTGCCGCACCGTCGATCCACTGGCCGGATCCTATTACATTGAGTCGCTGACCGATCAAATCGTCAAACAAGCCAGAGCTATTATCCAACAGATCGACGAAGCCGGTGGCATGGCGAAAGCGATCGAAGCAg                                                                 <  1:1579852/150‑1 (MQ=255)
                                                                                     ccAGGAAGAATCAGAACTCTGCCGCACCGTCGATCCACTGGCCGGATCCTATTACATTGAGTCGCTGACCGATCAAATCGTCAAACAAGCCAGAGCTATTATCCAACAGATCGACGAAGCCGGTGGCATGGCGAAAGCGATCGAAGCAgg                                                                >  1:824258/1‑150 (MQ=255)
                                                                                        ggAAGAATCAGAACTCTGCCGCACCGTCGATCCACTGGCCGGATCCTATTACATTGAGTCGCTGACCGATCAAATCGTCAAACAAGCCAGAGCTATTATCCAACAGATCGACGAAGCCGGTGGCATGGCGAAAGCGATCGAAGCAGGTCt                                                             <  2:1489281/150‑1 (MQ=255)
                                                                                         gaagaaTCAGAACTCTGCCGCACCGTCGATCCACTGGCCGGATCCTATTACATTGAGTTGCTGACCGATCAAATCGTCAAACAAGCCAGAGCTATTATCCAACAGATCGACGAAGCCGGTGGCATGGCGAAAGCGATCGAAGCAGGTCTg                                                            <  1:3849441/150‑1 (MQ=255)
                                                                                         gaagaaTCAGAACTCTGCCGCACCGTCGATCCACTGGCCGGATCCTATTACATTGAGTTGCTGACCGATCAAATCGTCAAACAAGCCAGAGCTATTATCCAACAGATCGACGAAGCCGGTGGCATGGCGAAAGCGATCGAAGCAGGTCTg                                                            <  1:1459580/150‑1 (MQ=255)
                                                                                           agaaTCAGAACTCTGCCGCACCGTCGATCCACTGGCCGGATCCTATTACATTGAGTCGCTGACCGATCAAATCGTCAAACAAGCCAGAGCTATTATCCAACAGATCGACGAAGCCGGTGGCATGGCGAAAGCGATCGAAGCAGGTCTGcc                                                          >  1:1300669/1‑150 (MQ=255)
                                                                                            gaaTCAGAACTCTGCCGCACCGTCGATCCACTGGCCGGATCCTATTACATTGAGTTGCTGACCGATCAAATCGTCAAACAAGCCAGAGCTATTATCCAACAGATCGACGAAGCCGGTGGCATGGCGAAAGCGATCGAAGCAGGTCTGCCa                                                         >  1:1702948/1‑150 (MQ=255)
                                                                                              aTCAGAACTCTGCCGCACCGTCGATCCACTGGCCGGATCCTATTACATTGAGTCGCTGACCGATCAAATCGTCAAACAAGCCAGAGCTATTATCCAACAGATCGACGAAGCCGGTGGCATGGCGAAAGCGATCGAAGCAGGTCTGCCaaa                                                       >  2:867341/1‑150 (MQ=255)
                                                                                                cAGAACTCTGCCGCACCGTCGATCCACTGGCCGGATCCTATTACATTGAGTCGCTGACCGATCAAATCGTCAAACAAGCCAGAGCTATTATCCAACAGATCGACGAAGCCGGTGGCATGGCGAAAGCGATCGAAGCAGGTCTGCCAAAAc                                                     >  2:3350834/1‑150 (MQ=255)
                                                                                                cAGAACTCTGCCGCACCGTCGATCCACTGGCCGGATCCTATTACATTGAGTCGCTGACCGATCAAATCGTCAAACAAGCCAGAGCTATTATCCAACAGATCGACGAAGCCGGTGGCATGGCGAAAGCGATCGAAGCAGGTCTGCCAAAAc                                                     >  2:2331040/1‑150 (MQ=255)
                                                                                                cAGAACTCAGCCGCACCGTCGATCCACTGGCCGGATCCTATTACATTGAGTCGCTGACCGATCAAATCGTCAAACAAGCCATAGCTATTATCCAACAGATCGACGAAGCCGGTGGCATGGCGAAAGCGATCGAAGCAGGTCTGCCAAAAc                                                     <  2:1965578/150‑1 (MQ=255)
                                                                                                 aGAACTCTGCCGCACCGTCGATCCACTGGCCGGATCCTATTACATTGAGTCGCTGACCGATCAAATCGTCAAACAAGCCAGAGCTATTATCCAACAGATCGACGAAGCCGGTGGCATGGCGAAAGCGATCGAAGCAGGTCTGCCAAAACg                                                    >  2:3687270/1‑150 (MQ=255)
                                                                                                 aGAACTCTGCCGCACCGTCGATCCACTGGCCGGATCCTATTACATTGAGTCGCTGACCGATCAAATCGTCAAACAAGCCAGAGCTATTATCCAACAGATCGACGAAGCCGGTGGCATGGCGAAAGCGATCGAAGCAGGTCTGCCAAAACg                                                    >  2:1230945/1‑150 (MQ=255)
                                                                                                 aGAACTCTGCCGCACCGTCGATCCACTGGCCGGATCCTATTACATTGAGTCGCTGACCGATCAAATCGTCAAACAAGCCAGAGCTATTATCCAACAGATCGACGAAGCCGGTGGCATGGCGAAAGCGATCGAAGCAGGTCTGCCAAAACg                                                    >  2:1230931/1‑150 (MQ=255)
                                                                                                 aGAACTCTGCCGCACCGTCGATCCACTGGCCGGATCCTATTACATTGAGTCGCTGACCGATCAAATCGTCAAACAAGCCAGAGCTATTATCCAACAGATCGACGAAGCCGGTGGCATGGCGAAAGCGATCGAAGCAGGTCTGCCAAAACg                                                    >  2:3687378/1‑150 (MQ=255)
                                                                                                    aCTCTGCCGCACCGTCGATCCACTGGCCGGATCCTATTACATTGAGTCGCTGACCGATCAAATCGTCAAACAAGCCAGAGCTATTATCCAACAGATCGACGAAGCCGGTGGCATGGCGAAAGCGATCGAAGCAGGTCTGCCAAAACGAAt                                                 <  2:264683/150‑1 (MQ=255)
                                                                                                        tGCCGCACCGTCGATCCACTGGCCGGATCCTATTACATTGAGTCGCTGACCGATCAAATCGTCAAACAAGCCAGAGCTATTATCCAACAGATCGACGAAGCCGGTGGCATGGCGAAAGCGATCGAAGCAGGTCTGCCAAAACGAATGATc                                             <  2:248447/150‑1 (MQ=255)
                                                                                                        tGCCGCACCGTCGATCCACTGGCCGGATCCTATTACATTGAGTCGCTGACCGATCAAATCGTCAAACAAGCCAGAGCTATTATCCAACAGATCGACGAAGCCGGTGGCATGGCGAAAGCGATCGAAGCAGGTCTGCCAAAACGAATGATc                                             >  1:2630840/1‑150 (MQ=255)
                                                                                                        tGCCGCACCGTCGATCCACTGGCCGGATCCTATTACATTGAGTCGCTGACCGATCAAATCGTCAAACAAGCCAGAGCTATTATCCAACAGATCGACGAAGCCGGTGGCATGGCGAAAGCGATCGAAGCAGGTCTGCCAAAACGAATGATc                                             <  1:58952/150‑1 (MQ=255)
                                                                                                        tGCCGCACCGTCGATCCACTGGCCGGATCCTATTACATTGAGTCGCTGACCGATCAAATCGTCAAACAAGCCAGAGCTATTATCCAACAGATCGACGAAGCCGGTGGCATGGCGAAAGCGATCGAAGCAGGTCTGCCAAAACGAATGATc                                             <  1:1319686/150‑1 (MQ=255)
                                                                                                          ccGCACCGTCGATCCACTGGCCGGATCCTATTACATTGAGTCGCTGACCGATCAAATCGTCAAACAAGCCAGAGCTATTATCCAACAGATCGACGAAGCCGGTGGCATGGCGAAAGCGATCGAAGCAGGTCTGCCAAAACGAATGATCGa                                           <  2:1340801/150‑1 (MQ=255)
                                                                                                           cGCACCGTCGATCCACTGGCCGGATCCTATTACATTGAGTCGCTGACCGATCAAATCGTCAAACAAGCCAGAGCTATTATCCAACAGATCGACGAAGCCGGTGGCATGGCGAAAGCGATCGAAGCAGGTCTGCCAAAACGAATGATCGaa                                          >  1:600833/1‑150 (MQ=255)
                                                                                                              aCCGTCGATCCACTGGCCGGATCCTATTACATTGAGTTGCTGACCGATCAAATCGTCAAACAAGCCAGAGCTATTATCCAACAGATCGACGAAGCCGGTGGCATGGCGAAAGCGATCGAAGCAGGTCTGCCAAAACGAATGATCGAagag                                       >  1:740008/1‑150 (MQ=255)
                                                                                                              aCCGTCGATCCACTGGCCGGATCCTATTACATTGAGTCGCTGACCGATCAAATCGTCAAACAAGCCAGAGCTATTATCCAACAGATCGACGAAGCCGGTGGCATGGCGAAAGCGATCGAAGCAGGTCTGCCAAAACGAATGATCGAagag                                       <  2:2351275/150‑1 (MQ=255)
                                                                                                              aCCGTCGATCCACTGGCCGGATCCTATTACATTGAGTCGCTGACCGATCAAATCGTCAAACAAGCCAGAGCTATTATCCAACAGATCGACGAAGCCGGTGGCATGGCGAAAGCGATCGAAGCAGGTCTGCCAAAACGAATGATCGAagag                                       >  2:2278992/1‑150 (MQ=255)
                                                                                                               ccGTCGATCCACTGGCCGGATCCTATTACATTGAGTCGCTGACCGATCAAATCGTCAAACAAGCCAGAGCTATTATCCAACAGATCGACGAAGCCGGTGGCATGGCGAAAGCGATCGAAGCAGGTCTGCCAAAACGAATGATCGAAGAgg                                      >  1:1014616/1‑150 (MQ=255)
                                                                                                                 gTCGATCCACTGGCCGGATCCTATTACATTGAGTCGCTGACCGATCAAATCGTCAAACAAGCCAGAGCTATTATCCAACAGATCGACGAAGCCGGTGGCATGGCGAAAGCGATCGAAGCAGGTCTGCCAAAACGAATGATCGAAGAGGcc                                    >  2:3751844/1‑150 (MQ=255)
                                                                                                                    gATCCACTGGCCGGATCCTATTACATTGAGTCGCTGACCGATCAAATCGTCAAACAAGCCAGAGCTATTATCCAACAGATCGACGAAGCCGGTGGCATGGCGAAAGCGATCGAAGCAGGTCTGCCAAAACGAATGATCGAAGAGGCCTCa                                 >  2:1378629/1‑150 (MQ=255)
                                                                                                                    gATCCACTGGCCGGATCCTATTACATTGAGTCGCTGACCGATCAAATCGTCAAACAAGCCAGAGCTATTATCCAACAGATCGACGAAGCCGGTGGCATGGCGAAAGCGATCGAAGCAGGTCTGCCAAAACGAATGATCGAAGAGGCCTCa                                 <  1:1097670/150‑1 (MQ=255)
                                                                                                                     aTCCACTGGCCGGATCCTATTACATTGAGTCGCTGACCGATCAAATCGTCAAACAAGCCAGAGCTATTATCCAACAGATCGACGAAGCCGGTGGCATGGCGAAAGCGATCGAAGCAGGTCTGCCAAAACGAATGATCGAAGAGGCCTCAg                                <  1:2264203/150‑1 (MQ=255)
                                                                                                                      tCCACTGGCCGGATCCTATTACATTGAGTCGCTGACCGATCAAATCGTCAAACAAGCCAGAGCTATTATCCAACAGATCGACGAAGCCGGTGGCATGGCGAAAGCGATCTAAGCAGGTCTGCCAAAACGAATGATCGAAGAGGCCTCAgc                               >  1:1975582/1‑150 (MQ=255)
                                                                                                                      tCCACTGGCCGGATCCTATTACATTGAGTCGCTGACCGATCAAATCGTCAAACAAGCCAGAGCTATTATCCAACAGATCGACGAAGCCGGTGGCATGGCGAAAGCGATCGAAGCAGGTCTGCCAAAACGAATGATCGAAGAGGCCTCAgc                               >  2:1087053/1‑150 (MQ=255)
                                                                                                                      tCCACTGGCCGGATCCTATTACATTGAGTCGCTGACCGATCAAATCGTCAAACAAGCCAGAGCTATTATCCAACAGATCGACGAAGCCGGTGGCATGGCGAAAGCGATCGAAGCAGGTCTGCCAAAACGAATGATCGAAGAGGCCTCAgc                               <  1:1944123/150‑1 (MQ=255)
                                                                                                                        cACTGGCCGGATCCTATTACATTGAGTCGCTGACCGATCAAATCGTCAAACAAGCCAGAGCTATTATCCAACAGATCGACGAAGCCGGTGGCATGGCGAAAGCGATCGAAGCAGGTCTGCCAAAACGAATGATCGAAGAGGCCTCAgcgc                             >  1:2822287/1‑150 (MQ=255)
                                                                                                                            ggCCGGATCCTATTACATTGAGTCGCTGACCGATCAAATCGTCAAACAATCCAGAGCTATTATCCAACATATCGACGAAGCCGGTGGCAAGGCGAAAGCGATCGAAGCAGGTCAGCCAAAACGAATGATAGAAGAGGCCACAGCGCGCGa                         >  1:3279592/1‑150 (MQ=255)
                                                                                                                              ccGGATCCTATTACATTGAGTCGCTGACCGATCAAATCGTCAAACAAGCCAGAGCTATTATCCAACAGATCGACGAAGCCGGTGGCATGGCGAAAGCGATCGAAGCAGGTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAAc                       >  2:2940751/1‑150 (MQ=255)
                                                                                                                               cGGATCCTATTACATTGAGTCGCTGACCGATCAAATCGTCAAACAAGCCAGAGCTATTATCCAACAGATCGACGAAGCCGGTGGCATGGCGAAAGCGATCGAAGCAGGTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACa                      <  1:2127583/150‑1 (MQ=255)
                                                                                                                                ggATCCTATTACATTGAGTCGCTGACCGATCAAATCGTCAAACAAGCCAGAGCTATTATCCAACAGATCGACGAAGCCGGTGGCATGGCGAAAGCGATCGAAGCAGGTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAg                     <  2:2926847/150‑1 (MQ=255)
                                                                                                                                ggAGCCTATTACATTGATTAGCTGACCGATCAAATCGTCAAACAAGCCAGAGCTATTTTCCAACAGATCGACGAAGCCGGTGGCATGGCGAAAGCGATCGAAGCAGGTCTGCCAAAACGAATGATAGAAGAGGCCTCAGAGCGCGAagag                     <  2:1647574/150‑1 (MQ=255)
                                                                                                                                 gATCCTATTACATTGAGTCGCTGACCGATCAAATCGTCAAACAAGCCAGAGCTATTATCCAACAGATCGACGAAGCCGGTGGCATGGCGAAAGCGATCGAAGCAGGTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGt                    >  1:1988725/1‑150 (MQ=255)
                                                                                                                                   tCCTATTACATTGAGTCGCTGACCGATCAAATCGTCAAACAAGCCAGAGCTATTATCCAACAGATCGACGAAGCCGGTGGCATGGCGAAAGCGATCGAAGCAGGTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCg                  >  2:20043/1‑150 (MQ=255)
                                                                                                                                     cTATTACATTGAGTCGCTGACCGATCAAATCGTCAAACAAGCCAGAGCTATTATCCAACAGATCGACGAAGCCGGTGGCATGGCGAAAGCGATCGAAGCAGGTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCt                <  1:3199277/150‑1 (MQ=255)
                                                                                                                                      tATTACATTGAGTCGCTGACCGATCAAATCGTCAAACAAGCCAGAGCTATTATCCAACAGATCGACGAAGCCGGTGGCATGGCGAAAGCGATCGAAGCAGGTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTg               <  2:3223699/150‑1 (MQ=255)
                                                                                                                                       aTTACATTGAGTCGCTGACCGATCAAATCGTCAAACAAGCCAGAGCTATTATCCAACAGATCGACGAAGCCGGTGGCATGGCGAAAGCGATCGAAGCAGGTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGa              >  1:2453903/1‑150 (MQ=255)
                                                                                                                                       aTTACATTGAGTCGCTGACCGATCAAATCGTCAAACAAGCCAGAGCTATTATCCAACAGATCGACGAAGCCGGTGGCATGGCGAAAGCGATCGAAGCAGGTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGa              <  2:3575760/150‑1 (MQ=255)
                                                                                                                                         tACATTGAGTCGCTGACCGATCAAATCGTCAAACAAGCCAGAGCTATTATCCAACAGATCGACGAAGCCGGTGGCATGGCGAAAGCGATCGAAGCAGGTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATc            <  1:3071082/150‑1 (MQ=255)
                                                                                                                                               gAGTCGCTGACCGATCAAATCGTCAAACAAGCCAGAGCTATTATCCAACAGATCGACGAAGCCGGTGGCATGGCGAAAGCGATCGAAGCAGGTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAg      >  2:4112024/1‑150 (MQ=255)
                                                                                                                                                aGTCGCTGACCGATCAAATCGTCAAACAAGCCAGAGCTATTATCCAACAGATCGACGAAGCCGGTGGCATGGCGAAAGCGATCGAAGCAGGTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAgg     >  2:3290005/1‑150 (MQ=255)
                                                                                                                                                 gTCGCTGACCGATCAAATAGTCAAACAAGCCAGAGCTATTATCCAACAGATCGACGAAGCCGGTGGCATGGCGAAAGCGATCGAAGCAGGTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAggg    <  1:1232795/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
CACTGCATACCAACGCCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCATCATCCAGGAAGAATCAGAACTCTGCCGCACCGTCGATCCACTGGCCGGATCCTATTACATTGAGTCGCTGACCGATCAAATCGTCAAACAAGCCAGAGCTATTATCCAACAGATCGACGAAGCCGGTGGCATGGCGAAAGCGATCGAAGCAGGTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCA  >  CP009273/3055260‑3055556

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: ‑