Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,214,112 C→A 7.0% T119T (ACC→ACA metH → homocysteine‑N5‑methyltetrahydrofolate transmethylase, B12‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,214,1120CA 7.0% 1225.7 / 40.4 324T119T (ACC→ACAmetHhomocysteine‑N5‑methyltetrahydrofolate transmethylase, B12‑dependent
Reads supporting (aligned to +/- strand):  ref base C (135/164);  new base A (15/10);  total (150/174)
Fisher's exact test for biased strand distribution p-value = 2.10e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 349 total aligned reads displayed.

GCCTACTTTGAAGCGGGCGCGGATATCATCGAAACCAACACCTTCAACTCCACGACCATTGCGATGGCGGATTACCAGATGGAATCCCTGTCGGCGGAAATCAACTTTGCGGCGGCGAAACTGGCGCGAGCTTGTGCTGACGAGTGGACCGCGCGCACGCCAGAGAAACCGCGCTACGTTGCCGGTGTTCTCGGCCCGACCAACCGCACGGCGTCTATTTCTCCGGACGTCAACGATCCGGCATTTCGTAATATCACTTTTGACGGGCTGGTGGCGGCTTATCGAGAGTCCACCAAAGC  >  CP009273/4213963‑4214261
                                                                                                                                                     |                                                                                                                                                     
gCCTACTTTGAAGCGGGCGCGGATATCATCGAAACCAACACCTTCAACTCCACGACCATTGCGATGGCGGATTACCAGATGGAATCCCTGTCGGCGGAAATCAACTTTGCGGCGGCGAAACTGGCGCGAGCTTGTGCTGACGAGTGGAcc                                                                                                                                                       <  1:480852/150‑1 (MQ=255)
     cTTTGAAGCGGGCGCGGATATCATCGAAACCAACACCTTCAACTCCACGACCATTGCGATGGCGGATTACCAGATGGAATCCCTGTCGGCGGAAATCAACTTTGCGGCGGCGAAACTGGCGCGAGCTTGTGCTGACGAGTGGACcgcgcg                                                                                                                                                  <  1:2663769/150‑1 (MQ=255)
        tGAAGCGGGCGCGGATATCATCGAAACCAACACCTTCAACTCCACGACCATTGCGATGGCGGATTACCAGATGGAATCCCTGTCGGCGGAAATCAACTTTGCGGCGGCGAAACTGGCGCGAGCTTGTGCTGACGAGTGGACCGCGCGCAc                                                                                                                                               >  2:5032518/1‑150 (MQ=255)
         gAAGCGGGCGCGGATATCATCGAAACCAACACCTTCAACTCCACGACCATTGCGATGGCGGATTACCAGATGGAATCCCTGTCGGCGGAAATCAACTTTGCGGCGGCGAAACTGGCGCGAGCTTGTGCTGACGAGTGGACCGCGCGCACg                                                                                                                                              >  2:495326/1‑150 (MQ=255)
          aaGCGGGCGCGGATATCATCGAAACCAACACCTTCAACTCCACGACCATTGCGATGGCGGATTACCAGATGGAATCCCTGTCGGCGGAAATCAACTTTGCGGCGGCGAAACTGGCGCGAGCTTGTGCTGACGAGTGGACCGCGCGCACGc                                                                                                                                             >  2:1997675/1‑150 (MQ=255)
               ggCGCGGATATCATCGAAACCAACACCTTCAACTCCACGACCATTGCGATGGCGGATTACCAGATGGAATCCCTGTCGGCGGAAATCAACTTTGCGGCGGCGAAACTGGCGCGAGCTTGTGCTGACGAGTGGACCGCGCGCACGCCagag                                                                                                                                        <  1:3862570/150‑1 (MQ=255)
                 cgcgGATATCATCGAAACCAACACCTTCAACTCCACGACCATTGCGATGGCGGATTACCAGATGGAATCCCTGTCGGCGGAAATCAACTTTGCGGCGGCGAAACTGGCGCGAGCTTGTGCTGACGAGTGGACCGCGCGCACGCCAGAGaa                                                                                                                                      >  1:4528673/1‑150 (MQ=255)
                   aggATATCATCGAAACCAACACCTTAAACTCCACGACCATTGCTATGGCGGATTACCAGAAGGAATCCCTGTCGGCGGAAATCAACTTTGCGGCGGCGAAACTGGCGCGAGCTTGTGCTGACGAGTGGACCGCGCGCACGCCAGAGacac                                                                                                                                    >  1:4096433/2‑150 (MQ=255)
                    ggATATCATCGAAACCAACACCTTCAACTCCACGACCATTGCGATGGCGGATTACCAGATGGAATCCCTGTCGGCGGAAATCAACTTTGCGGCGGCGAAACTGGCGCGAGCTTGTGCTGACGAGTGGACCGCGCGCACGCCAGAGAAAcc                                                                                                                                   <  1:542538/150‑1 (MQ=255)
                     gATATCATCGAAACCAACACCTTCAACTCCACGACCATTGCGATGGCGGATTACCAGATGGAATCCCTGTCGGCGGAAATCAACTTTGCGGCGGCGAAACTGGCGCGAGCTTGTGCTGACGAGTGGACCGCGCGCACGCCAGAGAAACcg                                                                                                                                  <  2:2512437/150‑1 (MQ=255)
                     gATATCATCGAAACCAACACCTTCAACTCCACGACCATTGCGATGGCGGATTACCAGATGGAATCCCTGTCGGCGGAAATCAACTTTGCGGCGGCGAAACTGGCGCGAGCTTGTGCTGACGAGTGGACCGCGCGCACGCCAGAGAAACcg                                                                                                                                  >  1:4637947/1‑150 (MQ=255)
                      atatCATCGAAACCAACACCTTCAACTCCACGACCATTGCGATGGCGGATTACCAGATGGAATCCCTGTCGGCGGAAATCAACTTTGCGGCGGCGAAACTGGCGCGAGCTTGTGCTGACGAGTGGACCGCGCGCACGCCAGAGAAACcgc                                                                                                                                 >  2:4174114/1‑150 (MQ=255)
                      atatCATCGAAACCAACACCATCAACTCCACGACCATTGCGATGGCGGATTACCAGATGGAATCCCTGTCGGCGGAAATCAACTTTGCGGCGGCGAAACTGGCGCGAGCTTGTGCTGACGAGTGGACCGCGCGCACGCCAGAGAAACcgc                                                                                                                                 >  1:3777752/1‑150 (MQ=255)
                       tatCATCGAAACCAACACCTTCAACTCCACGACCATTGCGATGGCGGATTACCAGATGGAATCCCTGTCGGCGGAAATCAACTTTGCGGCGGCGAAACTGGCGCGAGCTTGTGCTGACGAGTGGACCGCGCGCACGCCAGAGAAACcgcg                                                                                                                                >  2:2416727/1‑150 (MQ=255)
                       tatCATCGAAACCAACACCTTCAACTCCACGACCATTGCGATGGCGGATTACCAGATGGAATCCCTGTCGGCGGAAATCAACTTTGCGGCGGCGAAACTGGCGCGAGCTTGTGCTGACGAGTGGACCGCGCGCACGCCAGAGAAACcgcg                                                                                                                                >  1:1547893/1‑150 (MQ=255)
                        atcatcGAAACCAACACCTTCAACTCCACGACCATTGCGATGGCGGATTACCAGATGGAATCCCTGTCGGCGGAAATCAACTTTGCGGCGGCGAAACTGGCGCGAGCTTGTGCTGACGAGTGGACCGCGCGCACGCCAGAGAAACCgcgc                                                                                                                               >  2:3777698/1‑150 (MQ=255)
                         tcatcGAAACCAACACCTTCAACTCCACGACCATTGCGATGGCGGATTACCAGATGGAATCCCTGTCGGCGGAAATCAACTTTGCGGCGGCGAAACTGGCGCGAGCTTGTGCTGACGAGTGGACCGCGCGCACGCCAGAGAAACCGCGCt                                                                                                                              <  1:2194937/150‑1 (MQ=255)
                         tcatcGAAACCAACACCTTCAACTCCACGACCATTGCGATGGCGGATTACCAGATGGAATCCCTGTCGGCGGAAATCAACTTTGCGGCGGCGAAACTGGCGCGAGCTTGTGCTGACGAGTGGACCGCGCGCACGCCAGAGAAACCGCGCt                                                                                                                              <  2:1632059/150‑1 (MQ=255)
                           atcGAAACCAACACCTTCAACTCCACGACCATTGCGATGGCGGATTACCAGATGGAATCCCTGTCGGCGGAAATCAACTTTGCGGCGGCGAAACTGGCGCGAGCTTGTGCTGACGAGTGGACCGCGCGCACGCCAGAGAAACCGCGCTAc                                                                                                                            >  1:5361551/1‑150 (MQ=255)
                                aaCCAACACCTTCAACTCCACGACCATTGCGATGGCGGATTACCAGATGGAATCCCTGTCGGCGGAAATCAACTTTGCGGCGGCGAAACTGGCGCGAGCTTGTGCTGACGAGTGGACCGCGCGCACGCCAGAGAAACCGCGCTACGTTGc                                                                                                                       >  2:451062/1‑150 (MQ=255)
                                aaCCAACACCTTCAACTCCACGACCATTGCGATGGCGGATTACCAGATGGAATCCCTGTCGGCGGAAATCAACTTTGCGGCGGCGAAACTGGCGCGAGCTTGTGCTGACGAGTGGACCGCGCGCACGCCAGAGAAACCGCGCTACGTTGc                                                                                                                       >  1:521926/1‑150 (MQ=255)
                                aaCCAACACCTTCAACTCCACGACCATTGCGATGGCGGATTACCAGATGGAATCCCTGTCGGCGGAAATCAACTTTGCGGCGGCGAAACTGGCGCGAGCTTGTGCTGACGAGTGGACCGCGCGCACGCCAGAGAAACCGCGCTACGTTGc                                                                                                                       >  2:4577637/1‑150 (MQ=255)
                                 aCCAACACCTTCAACTCCACGACCATTGCGATGGCGGATTACCAGATGGAATCCCTGTCGGCGGAAATCAACTTTGCGGCGGCGAAACTGGCGCGAGCTTGTGCTGACGAGTGGACCGCGCGCACGCCAGAGAAACCGCGCTACGTTGcc                                                                                                                      >  1:69603/1‑150 (MQ=255)
                                 aCCAACACCTTCAACTCCACGACCATTGCGATGGCGGATTACCAGATGGAATCCCTGTCGGCGGAAATCAACTTTGCGGCGGCGAAACTGGCGCGAGCTTGTGCTGACGAGTGGACCGCGCGCACGCCAGAGAAACCGCGCTACGTTGcc                                                                                                                      <  1:1618155/150‑1 (MQ=255)
                                    aaCACCTTCAACTCCACGACCATTGCGATGGCGGATTACCAGATGGAATCCCTGTCGGCGGAAATCAACTTTGCGGCGGCGAAACTGGCGCGAGCTTGTGCTGACGAGTGGACCGCGCGCACGCCAGAGAAACCGCGCTACGTTGCCGgt                                                                                                                   <  1:20238/150‑1 (MQ=255)
                                          ttCAACTCCACGACCATTGCGATGGCGGATTACCAGATGGAATCCCTGTCGGCGGAAATCAACTTTGCGGCGGCGAAACTGGCGCGAGCTTGTGCTGACGAGTGGACCGCGCGCACGCCAGAGAAACCGCGCTACGTTGCCGGTGTtctc                                                                                                             <  2:1433102/150‑1 (MQ=255)
                                          ttCAACTCCACGACCATTGCGATGGCGGATTACCAGATGGAATCCCTGTCGGCGGAAATCAACTTTGCGGCGGCGAAACTGGCGCGAGCTTGTGCTGACGAGTGGACAGCGCGCACGCCAGAGAAACCGCGCTACGTTGCCGGTGTtctc                                                                                                             >  2:2960219/1‑150 (MQ=255)
                                             aaCTCCACGACCATTGCGATGGCGGATTACCAGATGGAATCCCTGTCGGCGGAAATCAACTTTGCGGCGGCGAAACTGGCGCGAGCTTGTGCTGACGAGTGGACCGCGCGCACGCCAGAGAAACCGCGCTACGTTGCCGGTGTTCTCGGc                                                                                                          >  1:5030593/1‑150 (MQ=255)
                                             aaCTCCACGACCATTGCGATGGCGGATTACCAGATGGAATCCCTGTCGGCGGAAATCAACTTTGCGGCGGCGAAACTGGCGCGAGCTTGTGCTGACGAGTGGACCGCGCGCACGCCAGAGAAACCGCGCTACGTTGCCGGTGTTCTCGGc                                                                                                          >  2:2652585/1‑150 (MQ=255)
                                              aCTCCACGACCATTGCGATGGCGGATTACCAGATGGAATCCCTGTCGGCGGAAATCAACTTTGCGGCGGCGAAACTGGCGCGAGCTTGTGCTGACGAGTGGACCGCGCGCACGCCAGAGAAACCGCGCTACGTTGCCGGTGTTCTCGGcc                                                                                                         <  1:3077672/150‑1 (MQ=255)
                                                tCCACGACCATTGCGATGGCGGATTACCAGATGGAATCCCTGTCGGCGGAAATCAACTTTGCGGCGGCGAAACTGGCGCGAGCTTGTGCTGACGAGTGGACCGCGCGCACGCCAGAGAAACCGCGCTACGTTGCCGGTGTTCTCGGCCCg                                                                                                       <  2:5109506/150‑1 (MQ=255)
                                                  cACGACCATTGCGATGGCGGATTACCAGATGGAATCCCTGTCGGCGGAAATCAACTTTGCGGCGGCGAAACTGGCGCGAGCTTGTGCTGACGAGTGGACCGCGCGCACGCCAGAGAAACCGCGCTACGTTGCCGGTGTTCTCGGCCCGAc                                                                                                     <  2:4231383/150‑1 (MQ=255)
                                                    cGACCATTGCGATGGCGGATTACCAGATGGAATCCCTGTCGGCGGAAATCAACTTTGCGGCGGCGAAACTGGCGCGAGCTTGTGCTGACGAGTGGACCGCGCGCACGCCAGAGAAACCGCGCTACGTTGCCGGTGTTCTCGGCCCGACCa                                                                                                   >  2:322582/1‑150 (MQ=255)
                                                       ccATTGCGATGGCGGATTACCAGATGGAATCCCTGTCGGCGGAAATCAACTTTGCGGCGGCGAAACTGGCGCGAGCTTGTGCTGACGAGTGGACAGCGCGCACGCCAGAGAAACCGCGCTACGTTGCCGGTGTTCTCGGCCCGACCAAcc                                                                                                <  2:2288379/150‑1 (MQ=255)
                                                         aTTGCGATGGCGGATTACCAGATGGAATCCCTGTCGGCGGAAATCAACTTTGCGGCGGCGAAACTGGCGCGAGCTTGTGCTGACGAGTGGACCGCGCGCACGCCAGAGAAACCGCGCTACGTTGCCGGTGTTCTCGGCCCGACCAACCGc                                                                                              >  1:1177364/1‑150 (MQ=255)
                                                          ttGCGATGGCGGATTACCAGATGGAATCCCTGTCGGCGGAAATCAACTTTGCGGCGGCGAAACTGGCGCGAGCTTGTGCTGACGAGTGGGCCGCGCGCACGCCAGAGAAACCGCGCTACGTTGCCGGTGTTCTCGGCCCGACCAACCGCa                                                                                             <  2:582702/150‑1 (MQ=255)
                                                           tGCGATGGCGGATTACCAGATGGAATCCCTGTCGGCGGAAATCAACTTTGCGGCGGCGAAACTGGCGCGAGCTTGTGCTGACGAGTGGACCGCGCGCACGCCAGAGAAACCGCGCTACGTTGCCGGTGTTCTCGGCCCGACCAACCGCAc                                                                                            <  2:3182678/150‑1 (MQ=255)
                                                              gATGGCGGATTACCAGATGGAATCCCTGTCGGCGGAAATCAACTTTGCGGCGGCGAAACTGGCGCGAGCTTGTGCTGACGAGTGGACCGCGCGCACGCCAGAGAAACCGCGCTACGTTGCCGGTGTTCTCGGCCCGACCAACCGCACGGc                                                                                         <  1:822254/150‑1 (MQ=255)
                                                                tGGCGGATTACCAGATGGAATCCCTGTCGGCGGAAATCAACTTTGCGGCGGCGAAACTGGCGCGAGCTTGTGCTGACGAGTGGACCGCGCGCACGCCAGAGAAACCGCGCTACGTTGCCGGTGTTCTCGGCCCGACCAACCGCACGGCGt                                                                                       <  2:3643268/150‑1 (MQ=255)
                                                                    ggATTACCAGATGGAATCCCTGTCGGCGGAAATCAACTTTGCGGCGGCGAAACTGGCGCGAGCTTGTGCTGACGAGTGGACCGCGCGCACGCCAGAGAAACCGCGCTCCGTTGCCGGTGTTCTCGGCCCGACCAACCGCACGGCGTCTAt                                                                                   >  1:4202400/1‑150 (MQ=255)
                                                                    ggATTACCAGATGGAATCCCTGTCGACGGAAATCAACTTTGCGGCGGCGAAACTGGCGCGAGCTTGTGCTGACGAGTGGACCGCGCGCACGCCAGAGAAACCGCGCTACGTTGCCGGTGTTCTCGGCCCGACCAACCGCACGGCGTCTAt                                                                                   <  1:2391835/150‑1 (MQ=255)
                                                                     gATTACCAGATGGAATCCCTGTCGGCGGAAATCAACTTTGCGGCGGCGAAACTGGCGCGAGCTTGTGCTGACGAGTGGACCGCGCGCACGCCAGAGAAACCGCGCTACGTTGCCGGTGTTCTCGGCCCGACCAACCGCACGGCGTCTAtt                                                                                  >  2:2045020/1‑150 (MQ=255)
                                                                      aTTACCAGATGGAATCCCTGTCGGCGGAAATCAACTTTGCGGCGGCGAAACTGGCGCGAGCTTGTGCTGACGAGTGGACCGCGCGCACGCCAGAGAAACCGCGCTACGTTGCCGGTGTTCTCGGCCCGACCAACCGCACGGCGTCTAttt                                                                                 >  1:2512288/1‑150 (MQ=255)
                                                                         aCCAGATGGAATCCCTGTCGGCGGAAATCAACTTTGCGGCGGCGAAACTGGCGCGAGCTTGTGCTGACGAGTGGACCGCGCGCACGCCAGAGAAACCGCGCTACGTTGCCGGTGTTCTCGGCCCGACCAACCGCACGGCGTCTATTtctc                                                                              <  1:4985164/150‑1 (MQ=255)
                                                                          ccAGATGGAATCCCTGTCGGCGGAAATCAACTTTGCGGCGGCGAAACTGGCGCGAGCTTGTGCTGACGAGTGGACCGCGCGCACGCCAGAGAAACCGCGCTACGTTGCCGGTGTTCTCGGCCCGACCAACCGCACGGCGTCTATTTCTcc                                                                             >  2:4887771/1‑150 (MQ=255)
                                                                          ccAGATGGAATCCCTGTCGGCGGAAATCAACTTTGCGGCGGCGAAACTGGCGCGAGCTTGTGCTGACGAGTGGACCGCGCGCACGCCAGAGAAACCGCGCTACGTTGCCGGTGTTCTCGGCCCGACCAACCGCACGGCGTCTATTTCTcc                                                                             <  1:5581211/150‑1 (MQ=255)
                                                                                ggAATCCCTGTCGGCGGAAATCAACTTTGCGGCGGCGAAACTGGCGCGAGCTTGTGCTGACGAGTGGACAGCGCGCACGCCAGAGAAACCGCGCTACGTTGCCGGTGTTCTCGGCCCGACCAACCGCACGGCGTCTATTTCTCCGGACGt                                                                       >  2:478887/1‑150 (MQ=255)
                                                                                  aaTCCCTGTCGGCGGAAATCAACTTTGCGGCGGCGAAACTGGCGCGAGCTTGTGCTGACGAGTGGACCGCGCGCACGCCAGAGAAACCGCGCTACGTTGCCGGTGTTCTCGGCCCGACCAACCGCACGGCGTCTATTTCTCCGGACGTCa                                                                     >  2:3317002/1‑150 (MQ=255)
                                                                                  aaTCCCTGTCGGCGGAAATCAACTTTGCGGCGGCGAAACTGGCGCGAGCTTGTGCTGACGAGTGGACCGCGCGCACGCCAGAGAAACCGCGCTACGTTGCCGGTGTTCTCGGCCCGACCAACCGCACGGCGTCTATTTCTCCGGACGTCa                                                                     <  1:4719734/150‑1 (MQ=255)
                                                                                  aaTCCCTGTCGGCGGAAATCAACTTTGCGGCGGCGAAACTGGCGCGAGCTTGTGCTGACGAGTGGACCGCGCGCACGCCAGAGAAACCGCGCTACGTTGCCGGTGTTCTCGGCCCGACCAACCGCACGGCGTCTATTTCTCCGGACGTCa                                                                     <  1:1666640/150‑1 (MQ=255)
                                                                                   aTCCCTGTCGGCGGAAATCAACTTTGCGGCGGCGAAACTGGCGCGAGCTTGTGCTGACGAGTGGACCGCGCGCACGCCAGAGAAACCGCGCTACGTTGCCGGTGTTCTCGGCCCGACCAACCGCACGGCGTCTATTTCTCCGGACGTCaa                                                                    <  1:1973892/150‑1 (MQ=255)
                                                                                    tCCCTGTCGGCGGAAATCAACTTTGCGGCGGCGAAACTGGCGCGAGCTTGTGCTGACGAGTGGACAGCGCGCACGCCAGAGAAACCGCGCTACGTTGCCGGTGTTCTCGGCCCGACCAACCGCACGGCGTCTATTTCTCCGGACGTCAAc                                                                   <  2:4373536/150‑1 (MQ=255)
                                                                                     cccTGTCGGCGGAAATCAACTTTGCGGCGGCGAAACTGGCGCGAGCTTGTGCTGACGAGTGGACCGCGCGCACGCCAGAGAAACCGCGCTACGTTGCCGGTGTTCTCGGCCCGACCAACCGCACGGCGTCTATTTCTCCGGACGTCAACg                                                                  >  2:5236512/1‑150 (MQ=255)
                                                                                         gTCGGCGGAAATCAACTTTGCGGCGGCGAAACTGGCGCGAGCTTGTGCTGACGAGTGGACCGCGCGCACGCCAGAGAAACCGCGCTACGTTGCCGGTGTTCTCGGCCCGACCAACCGCACGGCGTCTATTTCTCCGGACGTCAACGATcc                                                              <  2:2235850/150‑1 (MQ=255)
                                                                                         gTCGGCGGAAATCAACTTTGCGGCGGCGAAACTGGCGCGAGCTTGTGCTGACGAGTGGACCGCGCGCACGCCAGAGAAACCGCGCTACGTTGCCGGTGTTCTCGGCCCGACCAACCGCACGGCGTCTATTTCTCCGGACGTCAACGATcc                                                              <  1:3196294/150‑1 (MQ=255)
                                                                                         gTCGGCGGAAATCAACTTTGCGGCGGCGAAACTGGCGCGAGCTTGTGCTGACGAGTGGACCGCGCGCACGCCAGAGAAACCGCGCTACGTTGCCGGTGTTCTCGGCCCGACCAACCGCACGGCGTCTATTTCTCCGGACGTCAACGATcc                                                              <  2:2378618/150‑1 (MQ=255)
                                                                                          tCGGCTGAAATCAACTTTGCGGCGGCGAAACTGGCGCGAGCTTGTGCTGACGAGTGGACCGCGCGCACGCCAGAGAAACAGCGCTACGTTGCCGGTGTTCTCGGCCCGACCAACCGCTCGGCGTCTATTTCTCCGGACGTCAACGATCCg                                                             <  2:2942887/150‑1 (MQ=255)
                                                                                           cggcggAAATCAACTTTGCGGCGGCGAAACTGGCGCGAGCTTGTGCTGACGAGTGGACAGCGCGCACGCCAGAGAAACCGCGCTACGTTGCCGGTGTTCTCGGCCCGACCAACCGCACGGCGTCTATTTCTCCGGACGTCAACGATCCgg                                                            >  1:1027907/1‑150 (MQ=255)
                                                                                           cggcggAAATCAACTTTGCGGCGGCGAAACTGGCGCGAGCTTGTGCTGACGAGTGGACAGCGCGCACGCCAGAGAAACCGCGCTACGTTGCCGGTGTTCTCGGCCCGACCAACCGCACGGCGTCTATTTCTCCGGACGTCAACGATCCgg                                                            >  1:3873288/1‑150 (MQ=255)
                                                                                            ggcggAAATCAACTTTGCGGCGGCGAAACTGGCGCGAGCTTGTGCTGACGAGTGGACCGCGCGCACGCCAGAGAAACCGCGCTACGTTGCCGGTGTTCTCGGCCCGACCAACCGCACGGCGTCTATTTCTCCGGACGTCAACGATCCGGc                                                           >  2:3146071/1‑150 (MQ=255)
                                                                                             gcggAAATCAACTTTGCGGCGGCGAAACTGGCGCGAGCTTGTGCTGACGAGTGGACCGCGCGCACGCCAGAGAAACCGCGCTACGTTGCCGGTGTTCTCGGCCCGACCAACCGCACGGCGTCTATTTCTCCGGACGTCAACGATCCGGCa                                                          <  2:2152759/150‑1 (MQ=255)
                                                                                              cggAAATCAACTTTGCGGCGGCGAAACTGGCGCGAGCTTGTGCTGACGAGTGGACCGCGCGCACGCCAGAGAAACCGCGCTACGTTGCCGGTGTTCTCGGCCCGACCAACCGCACGGCGTCTATTTCTCCGGACGTCAACGATCCGGCAt                                                         >  1:2216365/1‑150 (MQ=255)
                                                                                              cggAAATCAACTTTGCGGCGGCGAAACTGGCGCGAGCTTGTGCTGACGAGTGGACCGCGCGCACGCCAGAGAAACCGCGCTACGTTGCCGGTGTTCTCGGCCCGACCAACCGCACGGCGTCTATTTCTCCGGACGTCAACGATCCGGCAt                                                         >  1:3243195/1‑150 (MQ=255)
                                                                                              cggAAATCAACTTTGCGGCGGCGAAACTGGCGCGAGCTTGTGCTGACGAGTGGACAGCGCGCACGCCAGAGAAACCGCGCTACGTTGCCGGTGTTCTCGGCCCGACCAACCGCACGGCGTCTATTTCTCCGGACGTCAACGATCCGGCAt                                                         <  1:3568803/150‑1 (MQ=255)
                                                                                               ggAAATCAACTTTGCGGCGGCGAAACTGGCGCGAGCTTGTGCTGACGAGTGGACCGCGCGCACGCCAGAGAAACCGCGCTACGTTGCCGGTGTTCTCGGCCCGACCAACCGCACGGCGTCTATTTCTCCGGACGTCAACGATCCGGCAtt                                                        >  2:3475649/1‑150 (MQ=255)
                                                                                               ggAAATCAACTGTGCGGCGGCGAAACTGGCGCGAGCTTGTGCTGACGAGTGGACCGCGCGCACGCCAGAGAAACCGCGCTACGTTGCCGGTGTTCTCGGCCCGACCAACCGCACGGCGTCTATTTCTCCGGACGTCAACGATCCGGCAtt                                                        <  1:3510723/150‑1 (MQ=255)
                                                                                                gAAATCAACTTTGCGGCGGCGAAACTGGCGCGAGCTTGTGCTGACGAGTGGACCGCGCGCACGCCAGAGAAACCGCGCTACGTTGCCGGTGTTCTCGGCCCGACCAACCGCACGGCGTCTATTTCTCCGGACGTCAACGATCCGGCAttt                                                       <  2:888494/150‑1 (MQ=255)
                                                                                                 aaaTCAACTTTGCGGCGGCGAAACTGGCGCGAGCTTGTGCTGACGAGTGGACCGCGCGCACGCCAGAGAAACCGCGCTACGGTGCCGGTGTTCTCGGCCCGACCAACCGCACGGCGTCTATTTCTCCGGACGTCAACGATCCGGCATTTc                                                      <  2:4417687/150‑1 (MQ=255)
                                                                                                  aaTCAACTTTGCGGCGGCGAAACTGGCGCGAGCTTGTGCTGACGAGTGGACCGCGCGCACGCCAGAGAAACCGCGCTACGTTGCCGGTGTTCTCGGCCCGACCAACCGCACGGCGTCTATTTCTCCGGACGTCAACGATCCGGCATTTCg                                                     >  2:1798581/1‑150 (MQ=255)
                                                                                                   aTCAACTTTGCGGCGGCGAAACTGGCGCGAGCTTGTGCTGACGAGTGGACAGCGCGCACGCCAGAGAAACCGCGCTACGTTGCCGGTGTTCTCGGCCCGACCAACCGCACGGCGTCTATTTCTCCGGACGTCAACGATCCGGCATTTCGt                                                    >  1:2101595/1‑150 (MQ=255)
                                                                                                     cAACTTTGCGGCGGCGAAACTGGCGCGAGCTTGTGCTGACGAGTGGACCGCGCGCACGCCAGAGAAACCGCGCTACGTTGCCGGTGTTCTCGGCCCGACCAACCGCACGGCGTCTATTTCTCCGGACGTCAACGATCCGGCATTTCGTaa                                                  <  2:5135062/150‑1 (MQ=255)
                                                                                                      aaCTTTGCGGCGGCGAAACTGGCGCGAGCTTGTGCTGACGAGTGGACCGCGCGCACGCCAGAGAAACCGCGCTACGTTGCCGGTGTTCTCGGCCCGACCAACCGCACGGCGTCTATTTCTCCGGACGTCAACGATCCGGCATTTCGTAat                                                 >  2:1524133/1‑150 (MQ=255)
                                                                                                          ttGCGGCGGCGAAACTGGCGCGAGCTTGTGCTGACGAGTGGACCGCGCGCACGCCAGAGAAACCGCGCTACGTTGCCGGTGTTCTCGGCCCGACCAACCGCACGGCGTCTATTTCTCCGGACGTCAACGATCCGGCATTTCGTAATATCa                                             <  1:2901742/150‑1 (MQ=255)
                                                                                                               gcggcgAAACTGGCGCGAGCTTGTGCTGACGAGTGGACCGCGCGCACGCCAGAGAAACCGCGCTACGTTGCCGGTGTTCTCGGCCCGACCAACCGCACGGCGTCTATTTCTCCGGACGTCAACGATCCGGCATTTCGTAATATCACtttt                                        <  2:673553/150‑1 (MQ=255)
                                                                                                               gcggcgAAACTGGCGCGAGCTTGTGCTGACGAGTGGACCGCGCGCACGCCAGAGAAACCGCGCTACGTTGCCGGTGTTCTCGGCCCGACCAACCGCACGGCGTCTATTTCTCCGGACGTCAACGATCCGGCATTTCGTAATATCACtttt                                        <  1:2604679/150‑1 (MQ=255)
                                                                                                               gcggcgAAACTGGCGCGAGCTTGTGCTGACGAGTGGACAGCGCGCACGCCAGAGAAACCGCGCTACGTTGCCGGTGTTCTCGGCCCGACCAACCGCACGGCGTCTATTTCTCCGGACGTCAACGATCCGGCATTTCGTAATATCACtttt                                        <  1:5145436/150‑1 (MQ=255)
                                                                                                                 ggcgAAACTGGCGCGAGCTTGTGCTGACGAGTGGACCGCGCGCACGCCAGAGAAACCGCGCTACGTTGCCGGTGTTCTCGGCCCGACCAACCGCACGGCGTCTATTTCTCCGGACGTCAACGATCCGGCATTTCGTAATATCACTTTTGa                                      >  1:2998482/1‑150 (MQ=255)
                                                                                                                  gcgAAACTGGCGCGAGCTTGTGCTGACGAGTGGACCGCGCGCACGCCAGAGAAACCGCGCTACGTTGCCGGTGTTCTCGGCCCGACCAACCGCACGGCGTCTATTTCTCCGGACGTCAACGATCCGGCATTTCGTAATATCACTTTTGAc                                     >  1:4300370/1‑150 (MQ=255)
                                                                                                                    gAAACTGGCGCGAGCTTGTGCTGACGAGTGGACCGCGCGCACGCCAGAGAAACCGCGCTACGTTGCCGGTGTTCTCGGCCCGACCAACCGCACGGCGTCTATTTCTCCGGACGTCAACGATCCGGCATTTCGTAATATCACTTTTGACgg                                   >  1:4508136/1‑150 (MQ=255)
                                                                                                                     aaaCTGGCGCGAGCTTGTGCTGACGAGTGGACCGCGCGCACGCCAGAGAAACCGCGCTACGTTGCCGGTGTTCTCGGCCCGACCAACCGCACGGCGTCTATTTCTCCGGACGTCAACGATCCGGCATTTCGTAATATCACTTTTGACggg                                  <  1:4779360/150‑1 (MQ=255)
                                                                                                                      aaCTGGCGCGAGCTTGTGCTGACGAGTGGACCGCGCGCACGCCAGAGAAACCGCGCTACGTTGCCGGTGTTCTCGGCCCGACCAACCGCACGGCGTCTATTTCTCCGGACGTCAACGATCCGGCATTTCGTAATATCACTTTTGACGGGc                                 <  2:31619/150‑1 (MQ=255)
                                                                                                                      aaCTGGCGCGAGCTTGTGCTGACGAGTGGACCGCGCGCACGCCAGAGAAACCGCGCTACGTTGCCGGTGTTCTCGGCCCGACCAACCGCACGGCGTCTATTTCTCCGGACGTCAACGATCCGGCATTTCGTAATATCACTTTTGACGGGc                                 >  1:1835018/1‑150 (MQ=255)
                                                                                                                        cTGGCGCGAGCTTGTGCTGACGAGTGGACCGCGCGCACGCCAGAGAAACCGCGCTACGTTGCCGGTGTTCTCGGCCCGACCAACCGCACGGCGTCTATTTCTCCGGACGTCAACGATCCGGCATTTCGTAATATCACTTTTGACGGGCtg                               <  1:4462256/150‑1 (MQ=255)
                                                                                                                          ggCGCGAGCTTGTGCTGACGAGTGGACCGCGCGCACGCCAGAGAAACCGCGCTACGTTGCCGGTGTTCTCGGCCCGACCAACCGCACGGCGTCTATTTCTCCGGACGTCAACGATCCGGCATTTCGTAATATCACTTTTGACGGGCtggt                             <  1:3174890/150‑1 (MQ=255)
                                                                                                                             gcgAGCTTGTGCTGACGAGTGGACCGCGCGCACGCCAGAGAAACCGCGCTACGTTGCCGGTGTTCTCGGCCCGACCAACCGCACGGCGTCTATTTCTCCGGACGTCAACGATCCGGCATTTCGTAATATCACTTTTGACGGGCTGGTggc                          <  2:126794/150‑1 (MQ=255)
                                                                                                                             gcgAGCTTGTGCTGACGAGTGGACCGCGCGCACGCCAGAGAAACCGCGCTACGTTGCCGGTGTTCTCGGCCCGACCAACCGCACGGCGTCTATTTCTCCGGACGTCAACGATCCGGCATTTCGTAATATCACTTTTGACGGGCTGGTggc                          <  2:4690248/150‑1 (MQ=255)
                                                                                                                                  cTTGTGCTGACGAGTGGACCGCGCGCACGCCAGAGAAACCGCGCTACGTTGCCGGTGTTCTCGGCCCGACCAACCGCACGGCGTCTATTTCTCCGGACGTCAACGATCCGGCATTTCGTAATATCACTTTTGACGGGCTGGTGGCGGCtt                     >  2:5383591/1‑150 (MQ=255)
                                                                                                                                  cTTGTGCTGACGAGTGGACCGCGCGCACGCCAGAGAAACCGCGCTACGTTGCCGGTGTTCTCGGCCCGACCAACCGCACGGCGTCTATTTCTCCGGACGTCAACGATCCGGCATTTCGTAATATCACTTTTGACGGGCTGGTGGCGGCtt                     <  1:3444820/150‑1 (MQ=255)
                                                                                                                                  cTTGTGCTGACGAGTGGACCGCGCGCACGCCAGAGAAACCGCGCTACGTTGCCGGTGTTCTCGGCCCGACCAACCGCACGGCGTCTATTTCTCCGGACGTCAACGATCCGGCATTTCGTAATATCACTTTTGACGGGCTGGTGGCGGCtt                     >  2:3906131/1‑150 (MQ=255)
                                                                                                                                  cTTGTGCTGACGAGTGGACCGCGCGCACGCCAGAGAAACCGCGCTACGTTGCCGGTGTTCTCGGCCCGACCAACCGCACGGCGTCTATTTCTCCGGACGTCAACGATCCGGCATTTCGTAATATCACTTTTGACGGGCTGGTGGCGGCtt                     <  1:655665/150‑1 (MQ=255)
                                                                                                                                    tgtgCTGACGAGTGGACCGCGCGCACGCCAGAGAAACCGCGCTACGTTGCCGGTGTTCTCGGCCCGACCAACCGCACGGCGTCTATTTCTCCGGACGTCAACGATCCGGCATTTCGTAATATCACTTTTGACGGGCTGGTGGCGGCTTAt                   <  2:3080178/150‑1 (MQ=255)
                                                                                                                                       gCTGACGAGTGGACCGCGCGCACGCCAGAGAAACCGCGCTACGTTGCCGGTGTTCTCGGCCCGACCAACCGCACGGCGTCTATTTCTCCGGACGTCAACGATCCGGCATTTCGTAATATCACTTTTGACGGGCTGGTGGCGGCTTATCga                <  1:1093087/150‑1 (MQ=255)
                                                                                                                                         tGACGAGTGGACCGCGCGCACGCCAGAGAAACCGCGCTACGTTGCCGGTGTTCTCGGCCCGACCAACCGCACGGCGTCTATTTCTCCGGACGTCAACGATCCGGCATTTCGTAATATCACTTTTGACGGGCTGGTGGCGGCTTATCgaga              >  1:1449757/1‑150 (MQ=255)
                                                                                                                                         tGACGAGTGGACCGCGCGCACGCCAGAGAAACCGCGCTACGTTGCCGGTGTTCTCGGCCCGACCAACCGCACGGCGTCTATTTCTCCGGACGTCAACGATCCGGCATTTCGTAATATCACTTTTGACGGGCTGGTGGCGGCTTATCgaga              >  2:2326783/1‑150 (MQ=255)
                                                                                                                                          gACGAGTGGACCGCGCGCACGCCAGAGAAACCGCGCTACGTTGCCGGTGTTCTCGGCCCGACCAACCGCACGGCGTCTATTTCTCCGGACGTCAACGATCCGGCATTTCGTAATATCACTTTTGACGGGCTGGTGGCGGCTTATCGagag             <  1:4174114/150‑1 (MQ=255)
                                                                                                                                           aCGAGTGGACCGCGCGCACGCCAGAGAAACCGCGCTACGTTGCCGGTGTTCTCGGCCCGACCAACCGCACGGCGTCTATTTCTCCGGACGTCAACGATCCGGCATTTCGTAATATCACTTTTGACGGGCTGGTGGCGGCTTATCGAGAGt            <  2:286575/150‑1 (MQ=255)
                                                                                                                                             gAGTGGACAGCGCGCACGCCAGAGAAACCGCGCTACGTTGCCGGTGTTCTCGGCCCGACCAACCGCACGGCGTCTATTTCTCCGGACGTCAACGATCCGGCATTTCGTAATATCACTTTTGACGGGCTGGTGGCGGCTTATCGAGAGTcc          <  1:867211/150‑1 (MQ=255)
                                                                                                                                              aGTGGACCGCGCGCACGCCAGAGAAACCGCGCTACGTTGCCGGTGTTCTCGGCCCGACCAACCGCACGGCGTCTATTTCTCCGGACGTCAACGATCCGGCATTTCGTAATATCACGTTTGACGGGCTGGTGGCGGCTTATCGAGAGTcca         >  2:3386632/1‑150 (MQ=255)
                                                                                                                                                 ggACCGCGCGCACGCCAGAGAAACCGCGCTACGTTGCCGGTGTTCTCGGCCCGACCAACCGCACGGCGTCTATTTCTCCGGACGTCAACGATCCGGCATTTCGTAATATCACTTTTGACGGGCTGGTGGCGGCTTATCGAGAGTccacca      <  1:486650/150‑1 (MQ=255)
                                                                                                                                                     cgcgcgCACGCCAGAGAAACCGCGCTACGTTGCCGGTGTTCTCGGCCCGACCAACCGCACGGCGTCTATTTCTCCGGACGTCAACGATCCGGCATTTCGTAATATCACTTTTGACGGGCTGGTGGCGGCTTATCGAGAGTCCACCAAAgc  <  2:1168861/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GCCTACTTTGAAGCGGGCGCGGATATCATCGAAACCAACACCTTCAACTCCACGACCATTGCGATGGCGGATTACCAGATGGAATCCCTGTCGGCGGAAATCAACTTTGCGGCGGCGAAACTGGCGCGAGCTTGTGCTGACGAGTGGACCGCGCGCACGCCAGAGAAACCGCGCTACGTTGCCGGTGTTCTCGGCCCGACCAACCGCACGGCGTCTATTTCTCCGGACGTCAACGATCCGGCATTTCGTAATATCACTTTTGACGGGCTGGTGGCGGCTTATCGAGAGTCCACCAAAGC  >  CP009273/4213963‑4214261

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: