Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,055,546 C→A 25.6% I446I (ATC→ATA)  scpA → methylmalonyl‑CoA mutase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,055,5460CA25.6% 529.1 / 173.8 238I446I (ATC→ATA) scpAmethylmalonyl‑CoA mutase
Reads supporting (aligned to +/- strand):  ref base C (108/69);  new base A (25/36);  total (133/105)
Fisher's exact test for biased strand distribution p-value = 7.35e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

Only 100 of 241 total aligned reads displayed.

TACATTGAGTCGCTGACCGATCAAATCGTCAAACAAGCCAGAGCTATTATCCAACAGATCGACGAAGCCGGTGGCATGGCGAAAGCGATCGAAGCAGGTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGTACTTGAGATCGACAACGTGATGGTGCGTAACGAGCAAATTGCTTCGCTGGAACGCATTCGCGCCACCCGTGATGATGCCG  >  CP009273/3055397‑3055694
                                                                                                                                                     |                                                                                                                                                    
tACATTGAGTCGCTGACCGATCAAATCGTCAAACAAGCCAGAGCTATTATCCAACAGATCGACGAAGCCGGTGGCATGGCGAAAGCGATCGAAGCAGGTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATc                                                                                                                                                      <  2:2200151/150‑1 (MQ=255)
tACATTGAGTCGCTGACCGATCAAATCGTCAAACAAGCCAGAGCTATTATCCAACAGATCGACGAAGCCGGTGGCATGGCGAAAGCGATCGAAGCAGGTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATc                                                                                                                                                      <  2:1506472/150‑1 (MQ=255)
 aCATTGAGTCGCTGACCGATCAAATCGTCAAACAAGCCAGAGCTATTATCCAACAGATCGACGAAGCCGGTGGCATGGCGAAAGCGATCGAAGCAGGTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCg                                                                                                                                                     >  1:4632210/1‑150 (MQ=255)
  cATTGAGTCGCTGACCGATCAAATCGTCAAACAAGCCAGAGCTATTATCCAACAGATCGACGAAGCCGGTGGCATGGCGAAAGCGATCGAAGCAGGTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGa                                                                                                                                                    >  1:4772771/1‑150 (MQ=255)
       aGTCGCTGACCGATCAAATCGTCAAACAAGCCAGAGCTATTATCCAACAGATCGACGAAGCCGGTGGCATGGCGAAAGCGATCGAAGCAGGTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAgg                                                                                                                                               <  2:1992092/150‑1 (MQ=255)
       aGTCGCTGACCGATCAAATCGTCAAACAAGCCAGAGCTATTATCCAACAGATCGACGAAGCCGGTGGCATGGCGAAAGCGATCGAAGCAGGTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATAGACCAgg                                                                                                                                               <  1:4003078/150‑1 (MQ=255)
           gCTGACCGATCAAATCGTCAAACAAGCCAGAGCTATTATCCAACAGATCGACGAAGCCGGTGGCATGGCGAAAGCGATCGAAGCAGGTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATAGACCAGGGCaa                                                                                                                                           >  2:2510820/1‑150 (MQ=255)
             tGACCGATCAAATCGTCAAACAAGCCAGAGCTATTATCCAACAGATCGACGAAGCCGGTGGCATGGCGAAAGCGATCGAAGCAGGTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATAGACCAGGGCAAGc                                                                                                                                         >  1:1698995/1‑150 (MQ=255)
              gACCGATCAAATCGTCAAACAAGCCAGAGCTATTATCCAACAGATCGACGAAGCCGGTGGCATGGCGAAAGCGATCGAAGCAGGTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATAGACCAGGGCAAGCg                                                                                                                                        <  1:1208065/150‑1 (MQ=255)
                 cGATCAAATCGTCAAACAAGCCAGAGCTATTATCCAACAGATCGACGAAGCCGGTGGCATGGCGAAAGCGATCGAAGCAGGTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCgtgt                                                                                                                                     >  2:5046774/1‑150 (MQ=255)
                   aTCAAATCGTCAAACAAGCCAGAGCTATTATCCAACAGATCGACGACGCCGGTGGCATGGCGAAAGCGATCGAAGCAGGTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATAGACCAGGGCAAGCGTGTCa                                                                                                                                   <  2:379081/150‑1 (MQ=255)
                    tCAAATCGTCAAACAAGCCAGAGCTATTATCCAACAGATCGACGAAGCCGGTGGCATGGCGAAAGCGATCGAAGCAGGTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCAt                                                                                                                                  >  1:4094064/1‑150 (MQ=255)
                        aTCGTCAAACACGCCAGAGCTATTATCCAACAGATCGACGAAGCCGGTGGCATGGCGAAAGCGATCGAAGCAGGTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATAGACCAGGGCAAGCGTGTCATCGtt                                                                                                                              <  2:4732635/150‑1 (MQ=255)
                           gTCAAACAAGCCAGAGCTATTATCCAACAGATCGACGAAGCCGGTGGCATGGCGAAAGCGATCGAAGCAGGTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGgt                                                                                                                           <  1:3148622/150‑1 (MQ=255)
                           gTCAAACAAGCCAGAGCTATTATCCAACAGATCGACGAAGCCGGTGGCATGGCGAAAGCGATCGAAGCAGGTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATAGACCAGGGCAAGCGTGTCATCGTTGgt                                                                                                                           <  1:2243158/150‑1 (MQ=255)
                              aaaCAAGCCAGAGCTATTATCCAACAGATCGACGAAGCCGGTGGCATGGCGAAAGCGATCGAAGCAGGTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTc                                                                                                                        >  1:2173199/1‑150 (MQ=255)
                                  aaGCCAGAGCTATTATCCAACAGATCGACGAAGCCGGTGGCATGGCGAAAGCGATCGAAGCAGGTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTcaaca                                                                                                                    >  2:1138589/1‑150 (MQ=255)
                                   aGCCAGAGCTATTATCCAACAGATCGACGAAGCCGGTGGCATGGCGAAAGCGATCGAAGCAGGTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTcaacaa                                                                                                                   >  1:4419066/1‑150 (MQ=255)
                                    gCCAGAGCTATTATCCAACAGATCGACGAAGCCGGTGGCATGGCGAAAGCGATCGAAGCAGGTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAg                                                                                                                  <  2:1622555/150‑1 (MQ=255)
                                      cAGAGCTATTATCCAACAGATCGACGAAGCCGGTGGCATGGCGAAAGCGATCGAAGCAGGTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCAACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTa                                                                                                                <  1:4307201/150‑1 (MQ=255)
                                        gagCTATTATCCAACAGATCGACGAAGCCGGTGGCATGGCGAAAGCGATCGAAGCAGGTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACa                                                                                                              >  1:103813/1‑150 (MQ=255)
                                        gagCTATTATCCAACAGATCGACGAAGCCGGTGGCATGGCGAAAGCGATCGAAGCAGGTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATAGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACa                                                                                                              >  2:3511107/1‑150 (MQ=255)
                                        gagCTATTATCCAACAGATCGACGAAGCCGGTGGCATGGCGAAAGCGATCGAAGCAGGTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATAGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACa                                                                                                              <  2:1920016/150‑1 (MQ=255)
                                           cTATTATCCAACAGATCGACGAAGCCGGTGGCATGGCGAAAGCGATCGAAGCAGGTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAAc                                                                                                           >  1:1311902/1‑150 (MQ=255)
                                            tattatCCAACAGATCGACGAAGCCGGTGGCATGGCGAAAGCGATCGAAGCAGGTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACt                                                                                                          >  1:1890933/1‑150 (MQ=255)
                                              ttatCCAACAGATCGACGAAGCCGGTGGCATGGCGAAAGCGATCGAAGCAGGTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATAGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTgg                                                                                                        >  2:901602/1‑150 (MQ=255)
                                                 tCCAACAGATCGACGAAGCCGGTGGCATGGCGAAAGCGATCGAAGCAGGTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGGTGGTGTCAACAAGTACAAACTGGATc                                                                                                     >  2:276721/1‑150 (MQ=255)
                                                       aGATCGACGAAGCCGGTGGCATGGCGAAAGCGATCGAAGCAGGTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAg                                                                                               <  1:2381115/150‑1 (MQ=255)
                                                         aTCGACGAAGCCGGTGGCATGGCGAAAGCGATCGAAGCAGGTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGAc                                                                                             <  2:2983977/150‑1 (MQ=255)
                                                         aTCGACGAAGCCGGTGGCATGGCGAAAGCGATCGAAGCAGGTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATAGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGAc                                                                                             <  2:525074/150‑1 (MQ=255)
                                                         aTCGACGAAGCCGGTGGCATGGCGAAAGCGATCGAAGCAGGTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATAGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGAc                                                                                             <  2:3133278/150‑1 (MQ=255)
                                                             acgaAGCCGGTGGCATGGCGAAAGCGATCGAAGCAGGTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGaaa                                                                                         >  2:418891/1‑150 (MQ=255)
                                                               gaAGCCGGTGGCATGGCGAAAGCGATCGAAGCAGGTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAAcc                                                                                       <  1:1709312/150‑1 (MQ=255)
                                                               gaAGCCGGTGGCATGGCGAAAGCGATCGAAGCAGGTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATAGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAAcc                                                                                       <  1:2180372/150‑1 (MQ=255)
                                                               gaAGCCGGTGGCATGGCGAAAGCGATCGAAGCAGGTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATAGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAAcc                                                                                       <  1:3996018/150‑1 (MQ=255)
                                                                aaGCCGGTGGCATGGCGAAAGCGATCGAAGCAGGTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCg                                                                                      >  1:522578/1‑150 (MQ=255)
                                                                aaGCCGGTGGCATGGCGAAAGCGATCGAAGCAGGTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCg                                                                                      >  2:4892992/1‑150 (MQ=255)
                                                                aaGCCGGTGGCATGGCGAAAGCGATCGAAGCAGGTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATAGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCg                                                                                      >  1:609616/1‑150 (MQ=255)
                                                                    cGGTGGCATGGCGAAAGCGATCGAAGCAGGTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTATAAACTGGATCACGAAGACGAAACCGATGt                                                                                  <  1:3969210/150‑1 (MQ=255)
                                                                     ggTGGCATGGCGAAAGCGATCGAAGCAGGTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGTa                                                                                 <  2:4947583/150‑1 (MQ=255)
                                                                       tGGCATGGCGAAAGCGATCGAAGCAGGTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGTACt                                                                               <  2:1221730/150‑1 (MQ=255)
                                                                       tGGCATGGCGAAAGCGATCGAAGCAGGTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATAGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGTACt                                                                               <  2:2964449/150‑1 (MQ=255)
                                                                        ggCATGGCGAAAGCGATCGAAGCAGGTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATAGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGTACtt                                                                              <  1:484055/150‑1 (MQ=255)
                                                                           aTGGCGAAAGCGATCGAAGCAGGTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGTACTTgag                                                                           >  1:312711/1‑150 (MQ=255)
                                                                           aTGGCGAAAGCGATCGAAGCAGGTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGTACTTgag                                                                           >  1:5150715/1‑150 (MQ=255)
                                                                           aTGGCGAAAGCGATCGAAGCAGGTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGTACTTgag                                                                           >  1:2732468/1‑150 (MQ=255)
                                                                           aTGGCGAAAGCGATCGAAGCAGGTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGTACTGgag                                                                           >  1:116591/1‑150 (MQ=255)
                                                                           aTGGCGAAAGCGATCGAAGCAGGTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGTACTGgag                                                                           >  1:3856322/1‑150 (MQ=255)
                                                                           aTGGCGAAAGCGATCGAAGCAGGTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATAGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGTACTTgag                                                                           >  2:161280/1‑150 (MQ=255)
                                                                             ggCGAAAGCGATCGAAGCAGGTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGTACTTGAGAt                                                                         >  1:2410908/1‑150 (MQ=255)
                                                                               cGAAAGCGATCGAAGCAGGTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGTACTTGAGATCg                                                                       <  2:1001568/150‑1 (MQ=255)
                                                                                gAAAGCGATCGAAGCAGGTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGTACTTGAGATCGa                                                                      >  2:3594331/1‑150 (MQ=255)
                                                                                 aaaGCGATCGAAGCAGGTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGTACTGGAGATCGAc                                                                     >  1:868024/1‑150 (MQ=255)
                                                                                   aGCGATCGAAGCAGGTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGTACTTGAGATCGACaa                                                                   <  1:133752/150‑1 (MQ=255)
                                                                                    gCGATCGAAGCAGGTCTGCCACAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGTACTTGAGATCGACAAc                                                                  <  2:1431170/150‑1 (MQ=255)
                                                                                    gCGATCGAAGCAGGTCTGCCAAAACGAATGTTCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCTAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGTACTTGAGATCGACAAc                                                                  <  2:3848619/150‑1 (MQ=255)
                                                                                       aTCGAAGCAGGTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGTACTTGAGATCGACAACGTg                                                               >  1:1547182/1‑150 (MQ=255)
                                                                                        tCGAAGCAGGTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATAGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGTACTTGAGATCGACAACGTGa                                                              >  1:2673186/1‑150 (MQ=255)
                                                                                            aGCAGGTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATAGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGTACTTGAGATCGACAACGTGATGGt                                                          >  2:3085216/1‑150 (MQ=255)
                                                                                             gCAGGTCTGCCAACACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGTACTTGAGATCGACAACGTGATGGTg                                                         <  1:3434169/150‑1 (MQ=255)
                                                                                             gCAGGTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGTACTTGAGATCGACAACGTGATGGTg                                                         >  2:4158909/1‑150 (MQ=255)
                                                                                               aGGTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATAGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGTACTTGAGATCGACAACGTGATGGTGCg                                                       >  1:3089122/1‑150 (MQ=255)
                                                                                                 gTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGTACTTGAGATCGACAACGTGATGGTGCGTa                                                     <  1:1821093/150‑1 (MQ=255)
                                                                                                   cTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGTACTTGAGATCGACAACGTGATGGTGCGTAAc                                                   <  1:4302188/150‑1 (MQ=255)
                                                                                                   cTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGTACTGGAGATCGACAACGTGATGGTGCGTAAc                                                   >  1:2526146/1‑150 (MQ=255)
                                                                                                     gCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACACACTGGATCACGAAGACGAAACCGATGTACTTGAGATCGACAACGTGATGGTGCGTAACGa                                                 >  2:4511001/1‑150 (MQ=255)
                                                                                                       cAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGTACTTGAGATCGACAACGTGATGGTGCGTAACGAGc                                               >  2:1220011/1‑150 (MQ=255)
                                                                                                       cAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATAGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGTACTGGAGATCGACAACGTGATGGTGCGTAACGAGc                                               >  1:2462461/1‑150 (MQ=255)
                                                                                                                tGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGTACTTGAGATCGACAACGTGATGGTGCGTAACGAGCAAATTGCtt                                      <  1:264457/150‑1 (MQ=255)
                                                                                                                 gATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGTACTTGAGATCGACAACGTGATGGTGCGTAACGAGCAAATTGCTTc                                     >  1:2015748/1‑150 (MQ=255)
                                                                                                                 gATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATAGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGTACTTGAGATCGACAACGTGATGGTGCGTAACGAGCAAATTGCTTc                                     <  1:4803203/150‑1 (MQ=255)
                                                                                                                  aTCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGTACTTGAGATCGACAACGTGATGGTGCGTAACGAGCAAATTGCTTCg                                    >  2:2700902/1‑150 (MQ=255)
                                                                                                                  aTCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGTACTTGAGATCGACAACGTGATGGTGCGTAACGAGCAAATTGCTTCg                                    <  2:3331899/150‑1 (MQ=255)
                                                                                                                  aTCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGTACTTGAGATCGACAACGTGATGGTGCGTAACGAGCAAATTGCTTCg                                    >  2:3296221/1‑150 (MQ=255)
                                                                                                                    cGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATAGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGCTGTACTTGAGATCGACAACGTGATGGTGCGTAACGAGCAAATTGCTTCGCt                                  >  2:1842343/1‑150 (MQ=255)
                                                                                                                     gAAGAGGCCTCAGCGCGCGAACAGTCGCTGATAGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGTACTTGAGATCGACAACGTGATGGTGCGTAACGAGCAAATTGCTTCGCTg                                 <  1:226259/150‑1 (MQ=255)
                                                                                                                      aaGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGTACTTGAGATCGACAACGTGATGGTGCGTAACGAGCAAATTGCTTCGCTgg                                <  2:4579478/150‑1 (MQ=255)
                                                                                                                      aaGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGTACTTGAGATCGACAACGTGATGGTGCGTAACGAGCAAATTGCTTCGCTgg                                >  2:3923582/1‑150 (MQ=255)
                                                                                                                       agagGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGTACTTGAGATCGACAACGTGATGGTGCGTAACGAGCAAATTGCTTCGCTGGa                               >  1:431423/1‑150 (MQ=255)
                                                                                                                       agagGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGTACTGGAGATCGACAACGTGATGGTGCGTAACGAGCAAATTGCTTCGCTGGa                               >  1:847278/1‑150 (MQ=255)
                                                                                                                         agGCCTCAGCGCGCGAACAGTCGCTGATAGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGTACTTGAGATCGACAACGTGATGGTGCGTAACGAGCAAATTGCTTCGCTGGAAc                             <  2:41249/150‑1 (MQ=255)
                                                                                                                         agGCCTCAGCGCGCGAACAGTCGCTGATAGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGTACTTGAGATCGACAACGTGATGGTGCGTAACGAGCAAATTGCTTCGCTGGAAc                             <  2:41302/150‑1 (MQ=255)
                                                                                                                          ggCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGTACTTGAGATCGACAACGTGATGGTGCGTAACGAGCAAATTGCTTCGCTGGAACg                            >  1:2082416/1‑150 (MQ=255)
                                                                                                                             cTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGTACTTGAGATCGACAACGTGATGGTGCGTAACGAGCAAATTGCTTCGCTGGAACGCAt                         >  2:1850550/1‑150 (MQ=255)
                                                                                                                              tCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGTACTTGAGATCGACAACGTGATGGTGCGTAACGAGCAAATTGCTTCGCTGGAACGCagt                        >  2:686098/1‑148 (MQ=255)
                                                                                                                              tCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGTACTTGAGATCGACAACGTGATGGTGCGTAACGAGCAAATTGCTTCGCTGGAACGCAtt                        >  2:3291586/1‑150 (MQ=255)
                                                                                                                              tCAGCGCGCGAACAGTCGCTGATAGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGTACTTGAGATCGACAACGTGATGGTGCGTAACGAGCAAATTGCTTCGCTGGAACGCAtt                        >  1:4548469/1‑150 (MQ=255)
                                                                                                                                aGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGTACTTGAGATCGACAACGTGATGGTGCGTACCGAGCAAATTGCTTCGCTGGAACGCATTcg                      >  2:4323491/1‑150 (MQ=255)
                                                                                                                                  cgcgcgAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGTACTTGAGATCGACAACGTGATGGTGCGTAACGAGCAAATTGCTTCGCTGGAACGCATTcgcg                    >  2:2985066/1‑150 (MQ=255)
                                                                                                                                    cgcgAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGTACTTGAGATCGACAACGTGATGGTGCGTAACGAGCAAATTGCTTCGCTGGAACGCATTCGCGcc                  >  2:2379027/1‑150 (MQ=255)
                                                                                                                                     gcgAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGTACTGGAGATCGACAACGTGATGGTGCGTAACGAGCAAATTGCTTCGCTGGAACGCATTCGCGCCa                 >  1:4943452/1‑150 (MQ=255)
                                                                                                                                      cgAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAAATGGATCACGAAGACGAAACCGATGTACTTGAGATCGACAACGTGATGGTGCGTAACTAGCAAATTGCTTCGCTGGAACGCATTCGCGCCAc                >  1:3725313/1‑150 (MQ=255)
                                                                                                                                       gAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGTACTTGAGATCGACAACGTGATGGTGCGTAACGAGCAAATTGCTTCGCTGGAACGCATTCGCGCCAcc               >  2:780008/1‑150 (MQ=255)
                                                                                                                                        aaCAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGTACTTGAGATCGACAACGTGATGGTGCGTAACGAGCAAATTGCTTCGCTGGAACGCATTCGCGCCAccc              <  1:3056369/150‑1 (MQ=255)
                                                                                                                                         aCAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGTACTTGAGATCGACAACGTGATGGTGCGTAACGAGCAAATTGCTTCGCTGGAACGCATTCGCGCCACCCg             <  1:1951114/150‑1 (MQ=255)
                                                                                                                                           aGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGTACTTGAGATCGACAACGTGATGGTGCGTAACGAGCAAATTGCTTCGCTGGAACGCATTCGCGCCACCCGtg           >  1:2262843/1‑150 (MQ=255)
                                                                                                                                             tCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGTACTTGAGATCGACAACGTGATGGTGCGTAACGAGCAAATTGCTTCGCTGGAACGCATTCGCGCCACCCGtgat         <  2:5151157/150‑1 (MQ=255)
                                                                                                                                                cTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGTACTTGAGATCGACAACGTGATGGTGCGTAACGAGCAAATTGCTTCGCTGGAACGCATTCGCGCCACCCGTgatgat      >  2:1088330/1‑150 (MQ=255)
                                                                                                                                                cTGATAGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGTACTTGAGATCGACAACGTGATGGTGCGTAACGAGCAAATTGCTTCGCTGGAACGCATTCGCGCCACCCGTgatgat      <  1:4191948/150‑1 (MQ=255)
                                                                                                                                                   aTCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGTACTTGAGATCGACAACGTGATGGTGCGTAACGAGCAAATTGCTTCGCTGGAACGCATTCGCGCCACCCGTGATGATgcc   <  2:2570915/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
TACATTGAGTCGCTGACCGATCAAATCGTCAAACAAGCCAGAGCTATTATCCAACAGATCGACGAAGCCGGTGGCATGGCGAAAGCGATCGAAGCAGGTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGTACTTGAGATCGACAACGTGATGGTGCGTAACGAGCAAATTGCTTCGCTGGAACGCATTCGCGCCACCCGTGATGATGCCG  >  CP009273/3055397‑3055694

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