Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,055,069 C→A 25.7% V287V (GTC→GTA)  scpA → methylmalonyl‑CoA mutase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,055,0690CA25.7% 850.9 / 249.1 361V287V (GTC→GTA) scpAmethylmalonyl‑CoA mutase
Reads supporting (aligned to +/- strand):  ref base C (124/144);  new base A (55/38);  total (179/182)
Fisher's exact test for biased strand distribution p-value = 4.05e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 363 total aligned reads displayed.

GGTGCCAACTGCGTGCAGCAGGTAGCATTTACGCTCGCTGATGGGATTGAGTACATCAAAGCAGCAATCTCTGCCGGACTGAAAATTGATGACTTCGCTCCTCGCCTGTCGTTCTTCTTCGGCATCGGCATGGATCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCAC  >  CP009273/3054920‑3055218
                                                                                                                                                     |                                                                                                                                                     
   gCCAACTGCGTGCAGCAGGTAGCATTTACGCTCGCTGATGGGATTGAGTACATCAAAGCAGCAATCTCTGCCGGACTGAAAATTGATGACTTCGCTCCTCGCCTGTCGTTCTTCTTCGGCATCGGCATGGATCTGTTTATGAACGTCGcc                                                                                                                                                    <  2:6930/150‑1 (MQ=255)
    ccAACTGCGTGCAGCAGGTAGCATTTACGCTCGCTGATGGGATTGAGTACATCAAAGCAGCAATCTCTGCCGGACTGAAAATTGATGACTTCGCTCCTCGCCTGTCGTTCTTCTTCGGCATCGGCATGGATCTGTTTATGAACGTCGCCa                                                                                                                                                   <  2:1924105/150‑1 (MQ=255)
     cAACTGCGTGCAGCAGGTAGCATTTACGCTCGCTGATGGGATTGAGTACATCAAAGCAGCAATCTCTGCCGGACTGAAAATTGATGACTTCGCTCCTCGCCTGTCGTTCTTCTTCGGCATCGGCATGGATCTGTTTATGAACGTCGCCAt                                                                                                                                                  <  1:1213751/150‑1 (MQ=255)
       aCTGCGTGCAGCAGGTAGCATTTACGCTCGCTGATGGGATTGAGTACATCAAAGCAGCAATCTCTGCCGGACTGAAAATTGATGACTTCGCTCCTCGCCTGTCGTTCTTCTTCGGCATCGGCATGGATCTGTTTATGAACGTCGCCATGt                                                                                                                                                <  1:2852502/150‑1 (MQ=255)
         tGCGTGCAGCAGGTAGCATTTACGCTCGCTGATGGGATTGAGTACATCAAAGCAGCAATCTCTGCCGGACTGAAAATTGATGACTTCGCTCCTCGCCTGTCGTTCTTCTTCGGCATCGGCATGGATCTGTTTATGAACGTCGCCATGTtg                                                                                                                                              >  1:2398073/1‑150 (MQ=255)
         tGCGTGCAGCAGGTAGCATTTACGCTCGCTGATGGGATTGAGTACATCAAAGCAGCAATCTCTGCCGGACTGAAAATTGATGACTTCGCTCCTCGCCTGTCGTTCTTCTTCGGCATCGGCATGGATCTGTTTATGAACGTCGCCATGTtg                                                                                                                                              >  2:4664440/1‑150 (MQ=255)
            gTGCAGCAGGTAGCATTTACGCTCGCTGATGGGATTGAGTACATCAAAGCAGCAATCTCTGCCGGACTGAAAATTGATGACTTCGCTCCTCGCCTGTCGTTCTTCTTCGGCATCGGCATGGATCTGTTTATGAACGTAGCCATGTtgcgt                                                                                                                                           <  1:3476331/150‑1 (MQ=255)
             tGCAGCAGGTAGCATTTACGCTCGCTGATGGGATTGAGTACATCAAAGCAGCAATCTCTGCCGGACTGAAAATTGATGACTTCGCTCCTCGCCTGTCGTTCTTCTTCGGCATCGGCATGGATCTGTTTATGAACGTCGCCATGTtgcgtg                                                                                                                                          <  1:3282223/150‑1 (MQ=255)
                    ggTAGCATTTACGCTCGCTGATGGGATTGAGTACATCAAAGCAGCAATCTCTGCCGGACTGAAAATTGATGACTTCGCTCCTCGCCTGTCGTTCTTCTTCGGCATCGGCATGGATCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACg                                                                                                                                   <  1:2133844/150‑1 (MQ=255)
                     gTAGCATTTACGCTCGCTGATGGGATTGAGTACATCAAAGCAGCAATCTCTGCCGGACTGAAAATTGATGACTTCGCTCCTCGCCTGTCGTTCTTCTTCGGCATCGGCATGGATCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGt                                                                                                                                  >  2:2537340/1‑150 (MQ=255)
                          aTTTACGCTCGCTGATGGGATTGAGTACATCAAAGCAGCAATCTCTGCCGGACTGAAAATTGATGACTTCGCTCCTCGCCTGTCGTTCTTCTTCGGCATCGGCATGGATCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGttattt                                                                                                                             >  2:603940/1‑150 (MQ=255)
                           tttACGCTCGCTGATGGGATTGAGTACATCAAAGCAGCAATCTCTGCCGGACTGAAAATTGATGACTTCGCTCCTCGCCTGTCGTTCTTCTTCGGCATCGGCATGGATCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGttattta                                                                                                                            >  2:3105069/1‑150 (MQ=255)
                               cgctcgctGATGGGATTGAGTACATCAAAGCAGCAATCTCTGCCGGACTGAAAATTGATGACTTCGCTCCTCGCCTGTCGTTCTTCTTCGGCATCGGCATGGATCTGTTTATGAACGTAGCCATGTTGCGTGCGGCACGTTATTTATGGa                                                                                                                        <  1:2156826/150‑1 (MQ=255)
                                  tcgctGATGGGATTGAGTACATCAAAGCAGCAATCTCTGCCGGACTGAAAATTGATGACTTCGCTCCTCGCCTGTCGTTCTTCTTCGGCATCGGCATGGATCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCg                                                                                                                     >  1:4690320/1‑150 (MQ=255)
                                  tcgctGATGGGATTGAGTACATCAAAGCAGCAATCTCTGCCGGACTGAAAATTGATGACTTCGCTCCTCGCCTGTCGTTCTTCTTCGGCATCGGCATGGATCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCg                                                                                                                     <  1:1777873/150‑1 (MQ=255)
                                  tcgctGATGGGATTGAGTACATCAAAGCAGCAATCTCTGCCGGACTGAAAATTGATGACTTCGCTCCTCGCCTGTCGTTCTTCTTCGGCATCGGAATGGATCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCg                                                                                                                     >  1:5140705/1‑150 (MQ=255)
                                     ctGATGGGATTGAGTACATCAAAGCAGCAATCTCTGCCGGACTGAAAATTGATGACTTCGCTCCTCGCCTGTCGTTCTTCTTCGGCATCGGCATGGATCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAg                                                                                                                  <  1:4648102/150‑1 (MQ=255)
                                      tGATGGGATTGAGTACATCAAAGCAGCAATCTCTGCCGGACTGAAAATTGATGACTTCGCTCCTCGCCTGTCGTTCTTCTTCGGCATCGGCATGGATCTGTTTATGAACGTAGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGc                                                                                                                 <  1:3201610/150‑1 (MQ=255)
                                      tGATGGGATTGAGTACATCAAAGCAGCAATCTCTGCCGGACTGAAAATTGATGACTTCGCTCCTCGCCTGTCGTTCTTCTTCGGCATCGGCATGGATCTGTTTATGAACGTAGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGc                                                                                                                 >  2:2140347/1‑150 (MQ=255)
                                          gggATTGAGTACATCAAAGCAGCAATCTCTGCCGGACTGAAAATTGATGACTTCGCTCCTCGCCTGTCGTTCTTCTTCGGCATCGGCATGGATCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTc                                                                                                             <  2:3699423/150‑1 (MQ=255)
                                             aTTGAGTACATCAAAGCAGCAATCTCTGCCGGACTGAAAATTGATGACTTCGCTCCTCGCCTGTCGTTCTTCTTCGGCATCGGCATGGATCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGt                                                                                                          <  1:1455299/150‑1 (MQ=255)
                                             aTTGAGTACATCAAAGCAGCAATCTCTGCCGGACTGAAAATTGATGACTTCGCTCCTCGCCTGTCGTTCTTCTTCGGCATCGGCATGGATCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGt                                                                                                          <  1:2965199/150‑1 (MQ=255)
                                                gAGTACATCAAAGCAGCAATCTCTGCCGGACTGAAAATTGATGACTTCGCTCCTCGCCTGTCGTTCTTCTTCGGCATCGGCATGGATCTGTTTATGAACGTAGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGa                                                                                                       >  1:1163632/1‑150 (MQ=255)
                                                 aGTACATCAAAGCAGCAATCTCTGCCGGACTGAAAATTGATGACTTCGCTCCTCGCCTGTCGTTCTTCTTCGGCATCGGCATGGATCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGAt                                                                                                      <  1:4341090/150‑1 (MQ=255)
                                                 aGTACATCAAAGCAGCAATCTCTGCCGGACTGAAAATTGATGACTTCGCTCCTCGCCTGTCGTTCTTCTTCGGCATCGGCATGGATCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGAt                                                                                                      >  2:1411340/1‑150 (MQ=255)
                                                   tACATCAAAGCAGCAATCTCTGCCGGACTGAAAATTGATGACTTCGCTCCTCGCCTGTCGTTCTTCTTCGGCATCGGCATGGATCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGAttt                                                                                                    <  2:4529319/150‑1 (MQ=255)
                                                   tACATCAAAGCAGCAATCTCTGCCGGACTGAAAATTGATGACTTCGCTCCTCGCCTGTCGTTCTTCTTCGGCATCGGCATGGATCTGTTTATGAACGTAGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGAttt                                                                                                    <  1:756785/150‑1 (MQ=255)
                                                   tACATCAAAGCAGCAATCTCTGCCGGACTGAAAATTGATGACTTCGCTCCTCGCCTGTCGTTCTTCTTCGGCATCGGCATGGATCTGTTTATGAACGTAGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGAttt                                                                                                    <  1:3630233/150‑1 (MQ=255)
                                                     cATCAAAGCAGCAATCTCTGCCGGACTGAAAATTGATGACTTCGCTCCTCGCCTGTCGTTCTTCTTCGGCATCGGCATGGATCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTgg                                                                                                  <  1:4773718/150‑1 (MQ=255)
                                                     cATCAAAGCAGCAATCTCTGCCGGACTGAAAATTGATGACTTCGCTCCTCGCCTGTCGTTCTTCTTCGGCATCGGCATGGATCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTgg                                                                                                  <  2:4086415/150‑1 (MQ=255)
                                                     cATCAAAGCAGCAATCTCTGCCGGACTGAAAATTGATGACTTCGCTCCTCGCCTGTCGTTCTTCTTCGGCATCGGCATGGATCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTgg                                                                                                  >  1:438044/1‑150 (MQ=255)
                                                       tCAAAGCAGCAATCTCTGCCGGACTGAAAATTGATGACTTCGCTCCTCGCCTGTCGTTCTTCTTCGGCATCGGCATGGATCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGgcg                                                                                                >  1:3762853/1‑150 (MQ=255)
                                                          aaGCAGCAATCTCTGCCGGACTGAAAATTGATGACTTCGCTCCTCGCCTGTCGTTCTTCTTCGGCATCGGCATGGATCTGTTTATGAACGTAGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGcac                                                                                             >  1:926041/1‑150 (MQ=255)
                                                            gcagcaATCTCTGCCGGACTGAAAATTGATGACTTCGCTCCTCGCCTGTCGTTCTTCTTCGGCATCGGCATGGATCTGTTTATGAACGTAGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAg                                                                                           >  2:920067/1‑150 (MQ=255)
                                                               gcaATCTCTGCCGGACTGAAAATTGATGACTTCGCTCCTCGCCTGTCGTTCTTCTTCGGCATCGGCATGGATCTGTTTATGAACGTAGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGAc                                                                                        >  1:897940/1‑150 (MQ=255)
                                                                 aaTCTCTGCCGGACTGAAAATTGATGACTTCGCTCCTCGCCTGTCGTTCTTCTTCGGCATCGGCATGGATCTGTTTATGAACGTCGCCCTGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGAccc                                                                                      >  1:1002648/1‑150 (MQ=255)
                                                                 aaTCTCTGCCGGACTGAAAATTGATGACTTCGCTCCTCGCCTGTCGTTCTTCTTCGGCATCGGCATGGATCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGAccc                                                                                      >  2:1795816/1‑150 (MQ=255)
                                                                 aaTCTCTGCCGGACTGAAAATTGATGACTTCGCTCCTCGCCTGTCGTTCTTCTTCGGCATCGGCATGGATCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGAccc                                                                                      >  1:2656620/1‑150 (MQ=255)
                                                                  aTCTCTGCCGGACTGAAAATTGATGACTTCGCTCCTCGCCTGTCGTTCTTCTTCGGCATCGGCATGGATCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCg                                                                                     >  2:1123764/1‑150 (MQ=255)
                                                                  aTCTCTGCCGGACTGAAAATTGATGACTTCGCTCCTCGCCTGTCGTTCTTCTTCGGCATCGGCATGGATCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCg                                                                                     >  1:1352467/1‑150 (MQ=255)
                                                                    ctctGCCGGACTGAAAATTGATGACTTCGCTCCTCGCCTGTCGTTCTTCTTCGGCATCGGCATGGATCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGaa                                                                                   <  2:3158618/150‑1 (MQ=255)
                                                                    ctctGCCGGACTGAAAATTGATGACTTCGCTCCTCGCCTGTCGTTCTTCTTCGGCATCGGCATGGATCTGTTTATGAACGTAGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGaa                                                                                   <  2:1272581/150‑1 (MQ=255)
                                                                    ctctGCCGGACTGAAAATTGATGACTTCGCTCCTCGCCTGTCGTTCTTCTTCGGCATCGGCATGGATCTGTTTATGAACGTAGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGaa                                                                                   <  2:4688240/150‑1 (MQ=255)
                                                                                gAAAATTGATGACTTCGCTCCTCGCCTGTCGTTCTTCTTCGGCATCGGCATGGATCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCt                                                                       <  1:4992245/150‑1 (MQ=255)
                                                                                 aaaaTTGATGACTTCGCTCCTCGCCTGTCGTTCTTCTTCGGCATCGGCATGGATCTGTTTATGAACGTAGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTg                                                                      <  1:1699680/150‑1 (MQ=255)
                                                                                 aaaaTTGATGACTTCGCTCCTCGCCTGTCGTTCTTCTTCGGCATCGGCATGGATCTGTTTATGAACGTAGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTg                                                                      >  2:4984180/1‑150 (MQ=255)
                                                                                  aaaTTGATGACTTCGCTCCTCGCCTGTCGTTCTTCTTCGGCATCGGCATGGATCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGc                                                                     <  2:4264430/150‑1 (MQ=255)
                                                                                    aTTGATGACTTCGCTCCTCGCCTGTCGTTCTTCTTCGGCATCGGCATGGATCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGt                                                                   >  1:1263572/1‑150 (MQ=255)
                                                                                       gatgaCTTCGCTCCTCGCCTGTCGTTCTTCTTCGGCATCGGCATGGATCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTAcc                                                                <  1:3341440/150‑1 (MQ=255)
                                                                                        atgaCTTCGCTCCTCGCCTGTCGTTCTTCTTCGGCATCGGCATGGATCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTAccc                                                               >  1:260198/1‑150 (MQ=255)
                                                                                        atgaCTTCGCTCCTCGCCTGTCGTTCTTCTTCGGCATCGGCATGGATCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTAccc                                                               <  1:4312474/150‑1 (MQ=255)
                                                                                        atgaCTTCGCTCCTCGCCTGTCGTTCTTCTTCGGCATCGGCATGGATCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTAccc                                                               >  2:543501/1‑150 (MQ=255)
                                                                                         tgaCTTCGCTCCTCGCCTGTCGTTCTTCTTCGGCATCGGCATGGATCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCa                                                              <  2:416683/150‑1 (MQ=255)
                                                                                         tgaCTTCGCTCCTCGCCTGTCGTTCTTCTTCGGCATCGGCATGGATCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCa                                                              <  1:742892/150‑1 (MQ=255)
                                                                                            cTTCGCTCCTCGCCTGTCGTTCTTCTTCGGCATCGGCATGGATCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTg                                                           <  1:1540827/150‑1 (MQ=255)
                                                                                             ttCGCTCCTCGCCTGTCGTTCTTCTTCGGCATCGGCATGGATCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGc                                                          >  2:844740/1‑150 (MQ=255)
                                                                                             ttCGCTCCTCGCCTGTCGTTCTTCTTCGGCATCGGCATGGATCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGc                                                          <  1:400095/150‑1 (MQ=255)
                                                                                                   cctcGCCTGTCGTTCTTCTTCGGCATCGGCATGGATCTGTTTATGAACGTAGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGAcc                                                    >  2:307924/1‑150 (MQ=255)
                                                                                                     tcGCCTGTCGTTCTTCTTCGGCATCGGCATGGATCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTc                                                  <  1:796307/150‑1 (MQ=255)
                                                                                                     tcGCCTGTCGTTCTTCTTCGGCATCGGCATGGATCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTc                                                  >  1:1857892/1‑150 (MQ=255)
                                                                                                      cGCCTGTCGTTCTTCTTCGGCATCGGCATGGATCTGTTTATGAACGTAGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCa                                                 <  1:4600456/150‑1 (MQ=255)
                                                                                                       gCCTGTCGTTCTTCTTCGGCATCGGCATGGATCTGTTTATGAACGTAGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAg                                                >  1:4926080/1‑150 (MQ=255)
                                                                                                       gCCTGTCGTTCTTCTTAGGCATCGGCATGGATCTTTTTATGAACGTAGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACAGCCAGACCTCAg                                                >  1:2212044/1‑150 (MQ=255)
                                                                                                        ccTGTCGTTCTTCTTCGGCATCGGCATGGATCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAgg                                               <  2:1849257/150‑1 (MQ=255)
                                                                                                        ccTGTCGTTCTTCTTCGGCATAGGCATGGATCTGTTTATGAACGTAGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAgg                                               >  2:2326035/1‑150 (MQ=255)
                                                                                                          tGTCGTTCTTCTTCGGCATCGGCATGGATCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCt                                             >  2:4868944/1‑150 (MQ=255)
                                                                                                           gTCGTTCTTCTTCGGCATCGGCATGGATCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTg                                            >  1:175585/1‑150 (MQ=255)
                                                                                                            tCGTTCTTCTTCGGCATCGGCATGGATCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTgg                                           >  1:2877971/1‑150 (MQ=255)
                                                                                                            tCGTTCTTCTTCGGCATCGGCATGGATCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTgg                                           <  2:3389955/150‑1 (MQ=255)
                                                                                                              gTTCTTCTTCGGCATCGGCATGGATCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAg                                         >  2:2259533/1‑150 (MQ=255)
                                                                                                               ttcttcttcGGCATCGGCATGGATCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGc                                        >  1:1755099/1‑150 (MQ=255)
                                                                                                                tcttcttcGGCATCGGCATGGATCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCTAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGcc                                       >  1:4213228/1‑150 (MQ=255)
                                                                                                                  ttcttcGGCATCGGCATGGATCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCctg                                     <  2:1477918/150‑1 (MQ=255)
                                                                                                                     ttcGGCATCGGCATGGATCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCctgact                                  >  2:739323/1‑150 (MQ=255)
                                                                                                                        ggcatcggcatGGATCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGaa                               <  1:4664444/150‑1 (MQ=255)
                                                                                                                        ggcatcggcatGGATCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGaa                               >  2:752741/1‑150 (MQ=255)
                                                                                                                        ggcatcggcatGGATCTGTTTATGAACGTAGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGaa                               <  1:4931190/150‑1 (MQ=255)
                                                                                                                          catcggcatGGATCTGTTTATGAACGTAGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACa                             >  1:4847484/1‑150 (MQ=255)
                                                                                                                           atcggcatGGATCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAg                            <  1:394665/150‑1 (MQ=255)
                                                                                                                           atcggcatGGATCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAg                            <  1:2990606/150‑1 (MQ=255)
                                                                                                                           atcggcatGGATCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAg                            <  1:2447780/150‑1 (MQ=255)
                                                                                                                             cggcatGGATCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGa                          <  1:4583233/150‑1 (MQ=255)
                                                                                                                             cggcatGGATCTGTTTATGAACGTAGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGa                          >  2:4386771/1‑150 (MQ=255)
                                                                                                                                catGGATCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATcc                       >  2:3554781/1‑150 (MQ=255)
                                                                                                                                 atGGATCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCg                      <  1:2178447/150‑1 (MQ=255)
                                                                                                                                  tGGATCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGt                     <  1:1076578/150‑1 (MQ=255)
                                                                                                                                   ggATCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGta                    >  2:3882308/1‑150 (MQ=255)
                                                                                                                                   ggATCTGTTTATGAACGTAGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGta                    <  1:1967534/150‑1 (MQ=255)
                                                                                                                                    gATCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGtat                   <  1:4167150/150‑1 (MQ=255)
                                                                                                                                     aTCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGtata                  >  2:2632358/1‑150 (MQ=255)
                                                                                                                                       cTGTTTATGAACGTAGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATaac                >  2:5090617/1‑150 (MQ=255)
                                                                                                                                        tGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATaaca               <  2:398601/150‑1 (MQ=255)
                                                                                                                                         gTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATaacaa              >  2:1740162/1‑150 (MQ=255)
                                                                                                                                           ttATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACg            <  2:2910264/150‑1 (MQ=255)
                                                                                                                                             aTGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGtt          >  1:672394/1‑150 (MQ=255)
                                                                                                                                              tGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTa         <  2:2799060/150‑1 (MQ=255)
                                                                                                                                              tGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTa         <  2:1638948/150‑1 (MQ=255)
                                                                                                                                                aaCGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATc       <  1:2058549/150‑1 (MQ=255)
                                                                                                                                                 aCGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATcc      <  2:834230/150‑1 (MQ=255)
                                                                                                                                                   gTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGc    >  1:4884144/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GGTGCCAACTGCGTGCAGCAGGTAGCATTTACGCTCGCTGATGGGATTGAGTACATCAAAGCAGCAATCTCTGCCGGACTGAAAATTGATGACTTCGCTCCTCGCCTGTCGTTCTTCTTCGGCATCGGCATGGATCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCAC  >  CP009273/3054920‑3055218

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