Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,043,595 C→A 11.9% M144I (ATG→ATT gcvT ← aminomethyltransferase, tetrahydrofolate‑dependent, subunit (T protein) of glycine cleavage complex

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,043,5950CA11.9% 1059.2 / 78.3 311M144I (ATG→ATTgcvTaminomethyltransferase, tetrahydrofolate‑dependent, subunit (T protein) of glycine cleavage complex
Reads supporting (aligned to +/- strand):  ref base C (124/150);  new base A (19/18);  total (143/168)
Fisher's exact test for biased strand distribution p-value = 4.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 319 total aligned reads displayed.

CCCGCTTCACCGGTATAACCAGTGGTGGCAATAAACAGATCGCCCGCCTGCACGCCAAAGAACGGTTTCATCCCTTCCACCGCCTGACGCTGGGCGTCATTAAACAGTGTGGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTCAGCGTGTTGGGTAATCCAGGAGAGGTCTTTTTCGCGGGTGGCGGAGTTAACAACGAGGCGGAAGAAATCTTCAGTAAAGTAGTAGACGATGAGGTCATCTATCAC  >  CP009273/3043446‑3043744
                                                                                                                                                     |                                                                                                                                                     
cccGCTTCACCGGTATAACCAGTGGTGGCAATAAACAGATCGCCCGCCTGCACGCCAAAGAACGGTTTCATCCCTTCCACCGCCTGACGCTGGGCGTCATTAAACAGTGTGGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATc                                                                                                                                                       <  2:3180736/150‑1 (MQ=255)
cccGCTTCACCGGTATAACCAGTGGTGGCAATAAACAGATCGCCCGCCTGCACGCCAAAGAACGGTTTCATCCCTTCCACCGCCTGACGCTGGGCGTCATTAAACAGTGTGGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATc                                                                                                                                                       <  2:1292852/150‑1 (MQ=255)
 ccGCTTCACCGGTATAACCAGTGGTGGCAATAAACAGATCGCCCGCCTGCACGCCAAAGAACGGTTTCATCCCTTCCACCGCCTGACGCTGGGCGTCATTAAACAGTGTGGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCa                                                                                                                                                      <  2:1589132/150‑1 (MQ=255)
  cGCTTCACCGGTATAACCAGTGGTGGCAATAAACAGATCGCCCGCCTGCACGCCAAAGAACGGTTTCATCCCTTCCACCGCCTGACGCTGGGCGTCATTAAACAGTGTGGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCAt                                                                                                                                                     <  1:2310657/150‑1 (MQ=255)
  cGCTTCACCGGTATAACCAGTGGTGGCAATAAACAGATCGCCCGCCTGCACGCCAAAGAACGGTTTCATCCCTTCCACCGCCTGACGCTGGGCGTCATTAAACAGTGTGGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCAt                                                                                                                                                     <  2:1403751/150‑1 (MQ=255)
     ttCACCGGTATAACCAGTGGTGGCAATAAACAGATCGCCCGCCTGCACGCCAAAGAACGGTTTCATCCCTTCCACCGCCTGACGCTGGGCGTCATTAAACAGTGTGGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGa                                                                                                                                                  >  1:3880798/1‑150 (MQ=255)
      tCACCGGTATAACCAGTGGTGGCAATAAACAGATCGCCCGCCTGCACGCCAAAGAACGGTTTCATCCCTTCCACCGCCTGACGCTGGGCGTCATTAAACAGTGTGGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGaa                                                                                                                                                 >  1:1005418/1‑150 (MQ=255)
       cACCGGTATAACCAGTGGTGGCAATAAACAGATCGCCCGCCTGCACGCCAAAGAACGGTTTCATCCCTTCCACCGCCTGACGCTGGGCGTCATTAAACAGTGTGGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGaaa                                                                                                                                                >  1:2412529/1‑150 (MQ=255)
         ccGGTATAACCAGTGGTGGCAATAAACAGATCGCCCGCCTGCACGCCAAAGAACGGTTTCATCCCTTCCACCGCCTGACGCTGGGCGTCATTAAACAGTGTGGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGAAAgg                                                                                                                                              >  2:309558/1‑150 (MQ=255)
          cGGTATAACCAGTGGTGGCAATAAACAGATCGCCCGCCTGCACGCCAAAGAACGGTTTCATCCCTTCCACCGCCTGACGCTGGGCGTCATTAAACAGTGTGGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATAATGGAAAGGt                                                                                                                                             >  2:852969/1‑150 (MQ=255)
             tataACCAGTGGTGGCAATAAACAGATCGCCCGCCTGCACGCCAAAGAACGGTTTCATCCCTTCCACCGCCTGACGCTGGGCGTCATTAAACAGTGTGGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGtcat                                                                                                                                          <  2:3153799/150‑1 (MQ=255)
              ataACCAGTGGTGGCAATAAACAGATCGCCCGCCTGCACGCCAAAGAACGGTTTCATCCCTTCCACCGCCTGACGCTGGGCGTCATTAAACAGTGTGGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGtcatc                                                                                                                                         >  1:4221033/1‑150 (MQ=255)
                aaCCAGTGGTGGCAATAAACAGATCGCCCGCCTGCACGCCAAAGAACGGTTTCATCCCTTCCACCGCCTGACGCTGGGCGTCATTAAACAGTGTGGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGTCATCAc                                                                                                                                       >  2:3546012/1‑150 (MQ=255)
                   cAGTGGTGGCAATAAACAGATCGCCCGCCTGCACGCCAAAGAACGGTTTCATCCCTTCCACCGCCTGACGCTGGGCGTCATTAAACAGTGTGGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGTCATCACGaa                                                                                                                                    >  1:4531102/1‑150 (MQ=255)
                    aGTGGTGGCAATAAACAGATCGCCCGCCTGCACGCCAAAGAACGGTTTCATCCCTTCCACCGCCTGACGCTGGGCGTCATTAAACAGTGTGGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGTCATCACGAAc                                                                                                                                   >  1:1709937/1‑150 (MQ=255)
                    aGTGGTGGCAATAAACAGATCGCCCGCCTGCACGCCAAAGAACGGTTTCATCCCTTCCACCGCCTGACGCTGGGCGTCATTAAACAGTGTGGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGTCATCACGAAc                                                                                                                                   <  1:449299/150‑1 (MQ=255)
                       ggtggCAATAAACAGATCGCCCGCCTGCACGCCAAAGAACGGTTTCATCCCTTCCACCGCCTGACGCTGGGCGTCATTAAACAGTGTGGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGTCATCACGAACGGt                                                                                                                                >  2:2321259/1‑150 (MQ=255)
                          ggCAATAAACAGATCGCCCGCCTGCACGCCAAAGAACGGTTTCATCCCTTCCACCGCCTGACGCTGGGCGTCATTAAACAGTGTGGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGTCATCACGAACGGTAAt                                                                                                                             <  2:1439493/150‑1 (MQ=255)
                          ggCAATAAACAGATCGCCCGCCTGCACGCCAAAGAACGGTTTCATCCCTTCCACCGCCTGACGCCGGGCGTCATTAAACAGTGTGGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGTCATCACGAACGGTAAt                                                                                                                             <  2:1674762/150‑1 (MQ=255)
                           gCAATAAACAGATCGCCCGCCTGCACGCCAAAGAACGGTTTCATCCCTTCCACCGCCTGACGCTGGGCGTCATTAAACAGTGTGGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGTCATCACGAACGGTAAtt                                                                                                                            >  1:61255/1‑150 (MQ=255)
                           gCAATAAACAGATCGCCCGCCTGCACGCCAAAGAACGGTTTCATCCCTTCCACCGCCTGACGCTGGGCGTCATTAAACAGTGTGGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGTCATCACGAACGGTAAtt                                                                                                                            <  2:2212012/150‑1 (MQ=255)
                            cAATAAACAGATCGCCCGCCTGCACGCCAAAGAACGGTTTCATCCCTTCCACCGCCTGACGCTGGGCGTCATTAAACAGTGTGGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGTCATCACGAACGGTAAttt                                                                                                                           <  2:2826182/150‑1 (MQ=255)
                            cAATAAACAGATCGCCCGCCTGCACGCCAAAGAACGGTTTCATCCCTTCCACCGCCTGACGCTGGGCGTCATTAAACAGTGTGGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGTCATCACGAACGGTAAttt                                                                                                                           <  2:4073857/150‑1 (MQ=255)
                              aTAAACAGATCGCCCGCCTGCACGCCAAAGAACGGTTTCATCCCTTCCACCGCCTGACGCTGGGCGTCATTAAACAGTGTGGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGTCATCACGAACGGTAATTTCg                                                                                                                         <  1:1300038/150‑1 (MQ=255)
                                aaaCAGATCGCCCGCCTGCACGCCAAAGAACGGTTTCATCCCTTCCACCGCCTGACGCTGGGCGTCATTAAACAGTGTGGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGTCATCACGAACGGTAATTTCGAt                                                                                                                       >  1:830905/1‑150 (MQ=255)
                                aaaCAGATCGCCCGCCTGCACGCCAAAGAACGGTTTCATCCCTTCCACCGCCTGACGCTGGGCGTCATTAAACAGTGTGGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGTCATCACGAACGGTAATTTCGAt                                                                                                                       >  1:4342565/1‑150 (MQ=255)
                                   cAGATCGCCCGCCTGCACGCCAAAGAACGGTCTCATCCCTTCCACCGCCTGACGCTGGGCGTCATTAAACAGTGTGGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATAATGGAAAGGTCATCACGAACGGTAATTTCGATGcc                                                                                                                    >  1:3104321/1‑150 (MQ=255)
                                    aGATCGCCCGCCTGCACGCCAAAGAACGGTTTCATCCCTTCCACCGCCTGACGCTGGGCGTCATTAAACAGTGTGGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGTCATCACGAACGGTAATTTCGATGCCg                                                                                                                   <  2:1840152/150‑1 (MQ=255)
                                    aGATCGCCCGCCTGCACGCCAAAGAACGGTTTCATCCCTTCCACCGCCTGACGCTGGGCGTCATTAAACAGTGTGGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGTCATCACGAACGGTAATTTCGATGCCg                                                                                                                   <  2:806743/150‑1 (MQ=255)
                                      aTCGCCCGCCTGCACGCCAAAGAACGGTTTCATCCCTTCCACCGCCTGACGCTGGGCGTCATTAAACAGTGTGGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGaa                                                                                                                 >  2:2275274/1‑150 (MQ=255)
                                        cgcccgccTGCACGCCAAAGAACGGTTTCATCCCTTCCACCGCCTGACGCTGGGCGTCATTAAACAGTGTGGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAg                                                                                                               >  2:2929478/1‑150 (MQ=255)
                                         gcccgccTGCACGCCAAAGAACGGTTTCATCCCTTCCACCGCCTGACGCTGGGCGTCATTAAACAGTGTGGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAgg                                                                                                              <  1:1235287/150‑1 (MQ=255)
                                           ccgccTGCACGCCATAGAACGGTTTCATCCCTTCCAACGCCTGACGCTGGGCGTCATTAAACAGTGTGGCAGCTTTTGCCTGCGCATTCGGCCCTTTCACGGCAATCATGGAAAGGTCATCACGAACGGTATTTTCGATGCCGAAAGGtt                                                                                                            <  2:4305604/150‑1 (MQ=255)
                                            cgccTGCACGCCAAAGAACGGTTTCATCCCTTCCACCGCCTGACGCTGGGCGTCATTAAACAGTGTGGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTc                                                                                                           <  2:3104703/150‑1 (MQ=255)
                                             gccTGCACGCCAAAGAACGGTTTCATCCCTTCCACCGCCTGACGCTGGGCGTCATTAAACAGTGTGGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTCa                                                                                                          >  2:4135847/1‑150 (MQ=255)
                                             gccTGCACGCCAAAGAACGGTTTCATCCCTTCCACCGCCTGACGCTGGGCGTCATTAAACAGTGTGGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTCa                                                                                                          >  1:2460092/1‑150 (MQ=255)
                                              ccTGCACGCCAAATAACGGTTTCATCCCTTCCACCGCCTGACGCTGGGCGTCATTAAACAGTGTGGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTCAg                                                                                                         <  1:3013408/150‑1 (MQ=255)
                                                 gCACGCCAAAGAACGGTTTCATCCCTTCCACCGCCTGACGCTGGGCGTCATTAAACAGTGTGGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATAATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTCAGCgt                                                                                                      >  1:2277544/1‑150 (MQ=255)
                                                 gCACGCCAAAGAACGGTTTCATCCCTTCCACCGCCTGACGCTGGGCGTCATTAAACAGTGTGGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATAATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTCAGCgt                                                                                                      <  2:773649/150‑1 (MQ=255)
                                                     gCCAAAGAACGGTTTCATCCCTTCCACCGCCTGACGCTGGGCGTCATTAAACAGTGTGGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTCAGCGGGTTg                                                                                                  >  2:1431914/1‑150 (MQ=255)
                                                      ccAAAGAACGGTTTCATCCCTTCCACCGCCTGACGCTGGGCGTCAGTAAACAGTGTGGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGTCATCACGAACGGTCATTTCGCTGCCGACAGGTTCAGCGTGTTgg                                                                                                 >  2:4353960/1‑150 (MQ=255)
                                                         aaGAACGGTTTCATCCCTTCCACCGCCTGACGCTGGGCGTCATTAAACAGTGTGGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTCAGCGTGTTGGGTa                                                                                              >  1:1941663/1‑150 (MQ=255)
                                                         aaGAACGGTTTCATCCCTTCCACCGCCTGACGCTGGGCGTCATTAAACAGTGTGGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTCAGCGTGTTGGGTa                                                                                              <  2:3697063/150‑1 (MQ=255)
                                                         aaGAACGGTTTCATCCCTTCCACCGCCTGACGCTGGGCGTCATTAAACAGTGTGGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTCAGCGTGTTGGGTa                                                                                              >  2:702669/1‑150 (MQ=255)
                                                               ggTTTCATCCCTTCCACCGCCTGACGCTGGGCGTCATTAAACAGTGTGGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTCAGCGTGTTGGGTAATCCAg                                                                                        >  2:3624095/1‑150 (MQ=255)
                                                               ggTTTCATCCCTTCCACCGCCTGACGCTGGGCGTCATTAAACAGTGTGGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTCAGCGTGTTGGGTAATCCAg                                                                                        >  2:2113077/1‑150 (MQ=255)
                                                                  ttCATCCCTTCCACCGCCTGACGCTGGGCGTCATTAAACAGTGTGGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTCAGCGTGTTGGGTAATCCAGgag                                                                                     >  1:2618055/1‑150 (MQ=255)
                                                                  ttCATCCCTTCCACCGCCTGACGCTGGGCGTCATTAAACAGTGTGGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATAATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTCAGCGTGTTGGGTAATCAAGgag                                                                                     <  1:2726013/150‑1 (MQ=255)
                                                                   tCATCCCTTCCACCGCCTGACGCTGGGCGTCATTAAACAGTGTGGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTCAGCGTGTTGGGTAATCCAGgaga                                                                                    <  1:3079629/150‑1 (MQ=255)
                                                                     aTCCCTTCCACCGCCTGACGCTGGGCGTCATTAAACAGTGTGGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTCAGCGTGTTGGGTAATCCAGGAGAgg                                                                                  <  2:3715231/150‑1 (MQ=255)
                                                                      tCCCTTCCACCGCCTGACGCTGGGCGTCATTAAACAGTGTGGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTCAGCGTGTTGGGTAATCCAGGAGAGGt                                                                                 <  2:2594459/150‑1 (MQ=255)
                                                                      tCCCTTCCACCGCCTGACGCTGGGCGTCATTAAACAGTGTGGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTCAGCGTGTTGGGTAATCCAGGAGAGGt                                                                                 >  2:398412/1‑150 (MQ=255)
                                                                       cccTTCCACCGCCTGACGCTGGGCGTCATTAAACAGTGTGGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTCAGCGTGTTGGGTAATCCAGGAGAGGTc                                                                                <  1:1076932/150‑1 (MQ=255)
                                                                       cccTTCCACCGCCTGACGCTGGGCGTCATTAAACAGTGTGGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTCAGCGTGTTGGGTAATCCAGGAGAGGTc                                                                                <  1:1606600/150‑1 (MQ=255)
                                                                         cTTCCACCGCCTGACGCTGGGCGTCATTAAACAGTGTGGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTCAGCGTGTTGGGTAATCCAGGAGAGGTCtt                                                                              >  1:906389/1‑150 (MQ=255)
                                                                          ttCCACCGCCTGACGCTGGGCGTCATTAAACAGTGTGGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTCAGCGTGTTGGGTAATCCAGGAGAGGTCttt                                                                             <  1:762990/150‑1 (MQ=255)
                                                                           tCCACCGCCTGACGCTGGGCGTCATTAAACAGTGTGGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTCAGCGTGTTGGGTAATCCAGGAGAGGTCtttt                                                                            <  2:3867109/150‑1 (MQ=255)
                                                                             cACCGCCTGACGCTGGGCGTCATTAAACAGTGTGGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATAATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTCAGCGTGTTGGGTAATCCAGGAGAGGTCTTTTTc                                                                          <  1:4392476/150‑1 (MQ=255)
                                                                                    tGACGCTGGGCGTCATTAAACAGTGTGGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGTTTCAGCGTGTTGGGTAATCCAGGAGAGGTCTTTTTCGCGGGTg                                                                   >  2:1794108/1‑150 (MQ=255)
                                                                                    tGACGCTGGGCGTCATTAAACAGTGTGGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTCAGCGTGTTGGGTAATCCAGGAGAGGTCTTTTTCGCGGGTg                                                                   >  2:524812/1‑150 (MQ=255)
                                                                                    tGACGCTGGGCGTCATTAAACAGTGTGGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTCAGCGTGTTGGGTAATCCAGGAGAGGTCTTGTTCGCGGGTg                                                                   >  2:1794425/1‑150 (MQ=255)
                                                                                     gACGCTGGGCGTCATTAAACAGTGTGGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTCAGCGTGTTGGGTAATCCAGGAGAGGTCTTTTTCGCGGGTgg                                                                  >  1:2517947/1‑150 (MQ=255)
                                                                                         cTGGGCGTCATTAAACAGTGTGGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTCAGCGTGTTGGGTAATCCAGGAGAGGTCTTTTTCGCGGGTGGCGGa                                                              <  2:1485181/150‑1 (MQ=255)
                                                                                          tGGGCGTCATTAAACAGTGTGGCAGCTTTTGCCTGCGCATTCGGCCCGTGCACGGCAATCATGGAAAGGTCATCACGAACGGTAATGTCGATGCCGAAAGGTTCAGCGTGTTGGGTAATCCAGGAGAGGTCTTTTTCGCGGGTGGCGGAg                                                             <  2:2481625/150‑1 (MQ=255)
                                                                                              cGTCATTAAACAGTGTGGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTCAGCGTGTTGGGTAATCCAGGAGAGGTCTTTTTCGCGGGTGGCGGAGTTaa                                                         >  2:3007808/1‑150 (MQ=255)
                                                                                              cGTCATTAAACAGTGTGGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTCAGCGTGTTGGGTAATCCAGGAGAGGTCTTTTTCGCGGGTGGCGGAGTTaa                                                         >  1:1764913/1‑150 (MQ=255)
                                                                                               gTCATTAAACAGTGTGGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGTCATCACGAACGGTGATTTCGATGCCGAAAGGTTCAGCGTGTTGGGTAATCCAGGAGAGGTCTTTTTCGCGGGTGGCGGAGTTaac                                                        <  1:416160/150‑1 (MQ=255)
                                                                                                 cATTAAACAGTGTGGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTCAGCGTGTTGGGTAATCCAGGAGAGGTCTTTTTCGCGGGTGGCGGAGTTaacaa                                                      >  2:3539233/1‑150 (MQ=255)
                                                                                                  aTTAAACCGTGTGGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTCAGCGTGTTGGGTAATCCAGGAGAGGTCTTTTTCGCGGGTGGCGGAGTTaacaac                                                     <  2:3654262/150‑1 (MQ=255)
                                                                                                  aTTAAACAGTGTGGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTCAGCGTGTTGGGTAATCCAGGAGAGGTCTTTTTCGCGGGTGGCGGAGTTaacaac                                                     <  2:2713624/150‑1 (MQ=255)
                                                                                                  aTTAAACAGTGTGGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTCAGCGTGTTGGGTAATCCAGGAGAGGTCTTTTTCGCGGGTGGCGGAGTTaacaac                                                     <  2:1112604/150‑1 (MQ=255)
                                                                                                     aaaCAGTGTGGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTCAGCGTGTTGGGTAATCCAGGAGAGGTCTTTTTCGCGGGTGGCGGAGTTAACAACGAg                                                  <  2:1852865/150‑1 (MQ=255)
                                                                                                     aaaCAGTGTGGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTCAGCGTGTTGGGTAATCCAGGAGAGGTCTTTTTCGCGGGTGGCGGAGTTAACAACGAg                                                  <  2:4530295/150‑1 (MQ=255)
                                                                                                       aCAGTGTGGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATAATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTCAGCGTGTTGGGTAATCCAGGAGAGGTCTTTTTCGCGGGTGGCGGAGTTAACAACGAGGc                                                <  2:704793/150‑1 (MQ=255)
                                                                                                         aGTGTGGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTCAGCGTGTTGGGTAATCCAGGAGAGGTCTTTTTCGCGGGTGGCGGAGTTAACAACGAGGCgg                                              >  2:1718902/1‑150 (MQ=255)
                                                                                                           tgtgGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTCAGCGTGTTGGGTAATCCAGGAGAGGTCTTTTTCGCGGGTGGCGGAGTTAACAACGAGGCGgaa                                            >  1:3050191/1‑150 (MQ=255)
                                                                                                             tgGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTCAGCGTGTTGGGTAATCCAGGAGAGGTCTTTTTCGCGGGTGGCGGAGTTAACAACGAGGCGgaaga                                          <  1:1397051/150‑1 (MQ=255)
                                                                                                             tgGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTCAGCGTGTTGGGTAATCCAGGAGAGGTCTTTTTCGCGGGTGGCGGAGTTAACAACGAGGCGgaaga                                          <  1:1742164/150‑1 (MQ=255)
                                                                                                             tgGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATAATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTCAGCGTGTTGGGTAATCCAGGAGAGGTCTTTTTCGCGGGTGGCGGAGTTAACAACGAGGCGgaaga                                          >  1:1206359/1‑150 (MQ=255)
                                                                                                              ggCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATAATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTCAGCGTGTTGGGTAATCCAGGAGAGGTCTTTTTCGCGGGTGGCGGAGTTAACAACGAGGCGgaagaa                                         >  1:1918018/1‑150 (MQ=255)
                                                                                                                cAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTCAGCGTGTTGGGTAATCCAGGAGAGGTCTTTTTCGCGGGTGGCGGAGTTAACAACGAGGCGGAAGAAAt                                       <  1:3200552/150‑1 (MQ=255)
                                                                                                                cAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATAATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTCAGCGTGTTGGGTAATCCAGGAGAGGTCTTTTTCGCGGGTGGCGGAGTTAACAACGAGGCGGAAGAAAt                                       <  1:724365/150‑1 (MQ=255)
                                                                                                                    ttttGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTCAGCGTGTTGGGTAATCCAGGAGAGGTCTTTTTCGCGGGTGGCGGAGTTAACAACGAGGCGGAAGAAATCTTc                                   <  2:2412529/150‑1 (MQ=255)
                                                                                                                     tttGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTCAGCGTGTTGGGTAATCCAGGAGAGGTCTTTTTCGCGGGTGGCGGAGTTAACAACGAGGCGGAAGAAATCTTCa                                  <  1:2901401/150‑1 (MQ=255)
                                                                                                                       tGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTCAGCGTGTTGGGTAATCCAGGAGAGGTCTTTTTCGCGGGTGGCGGAGTTAACAACGAGGCGGAAGAAATCTTCAGt                                >  2:3910946/1‑150 (MQ=255)
                                                                                                                         ccTGCGCATTCGGCCCTTGCACGGCAATAATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTCAGCGTGTTGGGTAATCCAGGAGAGGTCTTTTTCGCGGGTGGCGGAGTTAACAACGAGGCGGAAGAAATCTTCAGTaa                              >  1:113949/1‑150 (MQ=255)
                                                                                                                          cTGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTCAGCGTGTTGGGTAATCCAGGAGAGGTCTTTTTCGCGGGTGGCGGAGTTAACAACGAGGCGGAAGAAATCTTCAGTaaa                             <  2:2670148/150‑1 (MQ=255)
                                                                                                                              gcATTCGGCCCTTGCACGGCAATCATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTCAGCGTGTTGGGTAATCCAGGAGAGGTCTTTTTCGCGGGTGGCGGAGTTAACAACGAGGCGGAAGAAATCTTCAGTAAagtag                         >  2:2788558/1‑150 (MQ=255)
                                                                                                                                aTTCGGCCCTTGCACGGCAATCATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTCAGCGTGTTGGGTAATCCAGGAGAGGTCTTTTTCGCGGGTGGCGGAGTTAACAACGAGGCGGAAGAAATCTTCAGTAAAgtagta                       <  1:4018310/150‑1 (MQ=255)
                                                                                                                                 ttCGGCCCTTGCACGGCAATAATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTCAGCGTGTTGGGTAATCCAGGAGAGGTCTTTTTCGCGGGTGGCGGAGTTAACAACGAGGCGGAAGAAATCTTCAGTAAAGtagtag                      >  1:524687/1‑150 (MQ=255)
                                                                                                                                     gCCCTTGCACGGCAATCATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTCAGCGTGTTGGGTAATCCAGGAGAGGTCTTTTTCGCGGGTGGCGGAGTTAACAACGAGGCGGAAGAAATCTTCAGTAAAGTAGTAGACGa                  >  1:700523/1‑150 (MQ=255)
                                                                                                                                        cTTGCACGGCAATCATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTCAGCGTGTTGGGTAATCCAGGAGAGGTCTTTTTCGCGGGTGGCGGAGTTAACAACGAGGCGGAAGAAATCTTCAGTAAAGTAGTAGACGATGa               <  2:3197264/150‑1 (MQ=255)
                                                                                                                                           gCACGGCAATAATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTCAGCGTGTTGGGTAATCCAGGAGAGGTCTTTTTCGCGGGTGGCGGAGTTAACAACGAGGCGGAAGAAATCTTCAGTAAAGTAGTAGACGATGAGGt            <  2:1206359/150‑1 (MQ=255)
                                                                                                                                             aCGGCAATCATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTCAGCGTGTTGGGTAATCCAGGAGAGGTCTTTTTCGCGGGTGGCGGAGTTAACAACGAGGCGGAAGAAATCTTCAGTAAAGTAGTAGACGATGAGGTCa          <  1:180846/150‑1 (MQ=255)
                                                                                                                                              cGGCAATCATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTCAGCGTGTTGGGTAATCCAGGAGAGGTCTTTTTCGCGGGTGGCGGAGTTAACAACGAGGCGGAAGAAATCTTCAGTAAAGTAGTAGACGATGAGGTCAt         >  2:3791619/1‑150 (MQ=255)
                                                                                                                                              cGGCAATCATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTCAGCGTGTTGGGTAATCCAGGAGAGGTCTTTTTCGCGGGTGGCGGAGTTAACAACGAGGCGGAAGAAATCTTCAGTAAAGTAGTAGACGATGAGGTCAt         <  2:33531/150‑1 (MQ=255)
                                                                                                                                              cGGCAATCATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTCAGCGTGTTGGGTAATCCAGGAGAGGTCTTTTTCGCGGGTGGCGGAGTTAACAACGAGGCGGAAGAAATCTTCAGTAAAGTAGTAGACGATGAGGTCAt         <  2:1671690/150‑1 (MQ=255)
                                                                                                                                                gCAATCATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTCAGCGTGTTGGGTAATCCAGGAGAGGTCTTTTTCGCGGGTGGCGGAGTTAACAACGAGGCGGAAGAAATCTTCAGTAAAGTAGTAGACGATGAGGTCATCt       >  1:3537288/1‑150 (MQ=255)
                                                                                                                                                 cAATCATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTCAGCGTGTTGGGTAATCCAGGAGAGGTCTTTTTCGCGGGTGGCGGAGTTAACAACGAGGCGGAAGAAATCTTCAGTAAAGTAGTAGACGATGAGGTCATCTa      >  1:4377195/1‑150 (MQ=255)
                                                                                                                                                     cATGGAAAGGTCATCACGACCGGTAATTTCGATGCCGAAAGGTTCAGCGTGTTGGGTAATCCAGGAGAGGTCTTTTTCGCGGGTGGCGGAGTTAACAACGAGGCGGAAGAAATCTTCAGTAAAGTAGTAGACGATGAGGTCATCTATcac  <  2:1005418/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CCCGCTTCACCGGTATAACCAGTGGTGGCAATAAACAGATCGCCCGCCTGCACGCCAAAGAACGGTTTCATCCCTTCCACCGCCTGACGCTGGGCGTCATTAAACAGTGTGGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTCAGCGTGTTGGGTAATCCAGGAGAGGTCTTTTTCGCGGGTGGCGGAGTTAACAACGAGGCGGAAGAAATCTTCAGTAAAGTAGTAGACGATGAGGTCATCTATCAC  >  CP009273/3043446‑3043744

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: