Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,867,537 C→A 11.3% R414S (CGT→AGT)  yeaJ → putative diguanylate cyclase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,867,5370CA11.3% 1014.8 / 71.6 289R414S (CGT→AGT) yeaJputative diguanylate cyclase
Reads supporting (aligned to +/- strand):  ref base C (137/119);  new base A (21/12);  total (158/131)
Fisher's exact test for biased strand distribution p-value = 3.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 297 total aligned reads displayed.

AGCGGTTGCAGAAACTGGTGCAATCCGGTTCTTCGGTGATGTTTATTGCTATTGACATGGACAAGTTAAAGCAAATAAATGACACCCTCGGTCATCAGGAGGGGGATTTAGCGATTACGTTATTAGCTCAGGCGATTAAACAGTCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAACGTCTGCAACATATCGCGCCGCAGAAAGAGATCGGCTTCTCTTCC  >  CP009273/1867388‑1867686
                                                                                                                                                     |                                                                                                                                                     
agcGGTTGCAGAAACTGGTGCAATCCGGTTCTTCGGTGATGTTTATTGCTATTGACATGGACAAGTTAAAGCAAATAAATGACACCCTCGGTCATCAGGAGGGGGATTTAGCGATTACGTTATTAGCTCAGGCGATTAAACAGTCGATTc                                                                                                                                                       >  2:3026819/1‑150 (MQ=255)
agcGGTTGCAGAAACTGGTGCAATCCGGTTCTTCGGTGATGTTTATTGCTATTGACATGGACAAGTTAAAGCAAATAAATGACACCCTCGGTCATCAGGAGGGGGATTTAGCGATTACGTTATTAGCTCAGGCGATTAAACAGTCGATTc                                                                                                                                                       >  1:2471002/1‑150 (MQ=255)
 gcGGTTGCAGAAACTGGTGCAATCCGGTTCTTCGGTGATGTTTATTGCTATTGACATGGACAAGTTAAAGCAAATAAATGACACCCTCGGTCATCAGGAGGGGGATTTAGCGATTACGTTATTAGCTCAGGCGATTAAACAGTCGATTCg                                                                                                                                                      >  2:364713/1‑150 (MQ=255)
  cGGTTGCAGAAACTGGTGCAATCCGGTTCTTCGGTGATGTTTATTGCTATTGACATGGACAAGTTAAAGCAAATAAATGACACCCTCGGTCATCAGGAGGGGGATTTAGCGATTACGTTATTAGCTCAGGCGATTAAACAGTCGATTCGt                                                                                                                                                     >  2:1449816/1‑150 (MQ=255)
     ttGCAGAAACTGGTGCAATCCGGTTCTTCGGTGATGTTTATTGCTATTGACATGGACAAGTTAAAGCAAATAAATGACACCCTCGGTCATCAGGAGGGGGATTTAGCGATTACGTTATTAGCTCAGGCGATTAAACAGTCGATTCGTaaa                                                                                                                                                  <  1:1320818/150‑1 (MQ=255)
     ttGCAGAAACTGGTGCAATCCGGTTCTTCGGTGATGTTTATTGCTATTGACATGGACAAGTTAAAGCAAATAAATGACACCCTCGGTCATCAGGAGGGGGATTTAGCGATTACGTTATTAGCTCAGGCGATTAAACAGTCGATTCGTaaa                                                                                                                                                  >  2:4022158/1‑150 (MQ=255)
      tGCAGAAACTGGTGCAATCCGGTTCTTCGGTGATGTTTATTGCTATTGACATGGACAAGTTAAAGCAAATAAATGACACCCTCGGTCATCAGGAGGGGGATTTAGCGATTACGTTATTAGCTCAGGCGATTAAACAGTCGATTCGTaaaa                                                                                                                                                 <  1:5061108/150‑1 (MQ=255)
         aGAAACTGGTGCAATCCGGTTCTTCGGTGATGTTTATTGCTATTGACATGGACAAGTTAAAGCAAATAAATGACACCCTCGGTCATCAGGAGGGGGATTTAGCGATTACGTTATTAGCTCAGGCGATTAAACAGTCGATTAGTAAAAGTg                                                                                                                                              >  1:3113994/1‑150 (MQ=255)
                  tGCAATCCGGTTCTTCGGTGATGTTTATTGCTATTGACATGGACAAGTTAAAGCAAATAAATGACACCCTCGGTCATCAGGAGGGGGATTTAGCGATTACGTTATTAGCTCAGGCGATTAAACAGTCGATTCGTAAAAGTGATTATGCCa                                                                                                                                     >  2:1531973/1‑150 (MQ=255)
                     aaTCCGGTTCTTCGGTGATGTTTATTGCTATTGACATGGACAAGTTAAAGCAAATAAATGACACCCTCGGTCATCAGGAGGGGGATTTAGCGATTACGTTATTAGCTCAGGCGATTAAACAGTCGATTCGTAAAAGTGATTATGCCATcc                                                                                                                                  <  1:3965968/150‑1 (MQ=255)
                     aaTCCGGTTCTTCGGTGATGTTTATTGCTATTGACATGGACAAGTTAAAGCAAATAAATGACACCCTCGGTCATCAGGAGGGGGATTTAGCGATTACGTTATTAGCTCAGGCGATTAAACAGTCGATTCGTAAAAGTGATTATGCCATcc                                                                                                                                  <  1:4990929/150‑1 (MQ=255)
                      aTCCGGTTCTTCGGTGATGTTTATTGCTATTGACATGGACAAGTTAAAGCAAATAAATGACACCCTCGGTCATCAGGAGGGGGATTTAGCGATTACGTTATTAGCTCAGGCGATTAAACAGTCGATTCGTAAAAGTGATTATGCCATCCg                                                                                                                                 >  1:3454673/1‑150 (MQ=255)
                       tCCGGTTCTTCGGTGATGTTTATTGCTATTGACATGGACAAGTTAAAGCAAATAAATGACACCCTCGGTCATCAGGAGGGGGATTTAGCGATTACGTTATTAGCTCAGGCGATTAAACAGTCGATTCGTAAAAGTGATTATGCCATCCGa                                                                                                                                >  2:5047110/1‑150 (MQ=255)
                          ggTTCTTCGGTGATGTTTATTGCTATTGACATGGACAAGTTAAAGCAAATAAATGACACCCTCGGTCATCAGGAGGGGGATTTAGCGATTACGTTATTAGCTCAGGCGATTAAACAGTCGATTCGTAAAAGTGATTATGCCATCCGACTc                                                                                                                             <  1:1690798/150‑1 (MQ=255)
                              cttcGGTGATGTTTATTGCTATTGACATGGACAAGTTAAAGCAAATAAATGACACCCTCGGTCATCAGGAGGGGGATTTAGCGATTACGTTATTAGCTCAGGCGATTAAACAGTCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTg                                                                                                                         <  1:4529078/150‑1 (MQ=255)
                                tcGGTGATGTTTATTGCTATTGACATGGACAAGTTAAAGCAAATAAATGACACCCTCGGTCATCAGGAGGGGGATTTAGCGATTACGTTATTAGCTCAGGCGATTAAACAGTCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGc                                                                                                                       >  1:4585316/1‑150 (MQ=255)
                                   gTGATGTTTATTGCTATTGACATGGACAAGTTAAAGCAAATAAATGACACCCTCGGTCATCAGGAGGGGGATTTAGCGATTACGTTATTAGCTCAGGCGATTAAACAGTCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGAt                                                                                                                    >  1:892937/1‑150 (MQ=255)
                                     gATGTTTATTGCTATTGACATGGACAAGTTAAAGCAAATAAATGACACCCTCGGTCATCAGGAGGGGGATTTAGCGATTACGTTATTAGCTCAGGCGATTAAACAGTCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGa                                                                                                                  >  1:2658451/1‑150 (MQ=255)
                                      aTGTTTATTGCTATTGACATGGACAAGTTAAAGCAAATAAATGACACCCTCGGTCATCAGGAGGGGGATTTAGCGATTACGTTATTAGCTCAGGCGATTAAACAGTCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGaa                                                                                                                 >  1:1998776/1‑150 (MQ=255)
                                      aTGTTTATTGCTATTGACATGGACAAGTTAAAGCAAATAAATGACACCCTCGGTCATCAGGAGGGGGATTTAGCGATTACGTTATTAGCTCAGGCGATTAAACAGTCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGaa                                                                                                                 >  2:1942593/1‑150 (MQ=255)
                                      aTGTTTATTGCTATTGACATGGACAAGTTAAAGCAAATAAATGACACCCTCGGTCATCAGGAGGGGGATTTAGCGATTACGTTATTAGCTCAGGCGATTAAACAGTCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGaa                                                                                                                 <  2:376952/150‑1 (MQ=255)
                                      aTGTTTATTGCTATTGACATGGACAAGTTAAAGCAAATAAATGACACCCTCGGTCATCAGGAGGGGGATTTAGCGATTACGTTATTAGCTCAGGCGATTAAACAGTCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGaa                                                                                                                 <  2:750281/150‑1 (MQ=255)
                                        gTTTATTGCTATTGACATGGACAAGTTAAAGCAAATAAATGACACCCTCGGTCATCAGGAGGGGGATTTAGCGATTACGTTATTAGCTCAGGCGATTAAACAGTCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAAtt                                                                                                               <  2:4879782/150‑1 (MQ=255)
                                          ttATTGCTATTGACATGGACAAGTTAAAGCAAATAAATGACACCCTCGGTCATCAGGAGGGGGATTTAGCGATTACGTTATTAGCTCAGGCGATTAAACAGTCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCt                                                                                                             <  1:1083608/150‑1 (MQ=255)
                                           tATTGCTATTGACATGGACAAGTTAAAGCAAATAAATGACACCCTCGGTCATAAGGAGGGGGATTTAGCGATTACGTTATTAGCTCAGGCGATTAAACAGTCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTg                                                                                                            >  2:3605457/1‑150 (MQ=255)
                                            aTTGCTATTGACATGGACAAGTTAAAGCAAATAAATGACACCCTCGGTCATCAGGAGGGGGATTTAGCGATTACGTTATTAGCTCAGGCGATTAAACAGTCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGc                                                                                                           >  1:1512321/1‑150 (MQ=255)
                                            aTTGCTATTGACATGGACAAGTTAAAGCAAATAAATGACACCCTCGGTCATAAGGAGGGGGATTTAGCGATTACGTTATTAGCTCAGGCGATTAAACAGTCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGc                                                                                                           >  2:4215871/1‑150 (MQ=255)
                                              tGCTATTGACATGGACAAGTTAAAGCAAATAAATGACACCCTCGGTCATCAGGAGGGGGATTTAGCGATTACGTTATTAGCTCAGGCGATTAAACAGTCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGcat                                                                                                         >  1:2704920/1‑150 (MQ=255)
                                                cTATTGACATGGACAAGTTAAAGCAAATAAATGACACCCTCGGTCATCAGGAGGGGGATTTAGCGATTACGTTATTAGCTCAGGCGATTAAACAGTCGATTAGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGcatca                                                                                                       >  2:3065875/1‑150 (MQ=255)
                                                 tATTGACATGGACAAGTTAAAGCAAATAAATGACACCCTCGGTCATCAGGAGGGGGATTTAGCGATTACGTTATTAGCTCAGGCGATTAAACAGTCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGcatcat                                                                                                      >  2:793052/1‑150 (MQ=255)
                                                    tGACATGGACAAGTTAAAGCAAATAAATGACACCCTCGGTCATCAGGAGGGGGATTTAGCGATTACGTTATTAGCTCAGGCGATTAAACAGTCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCt                                                                                                   <  2:1009819/150‑1 (MQ=255)
                                                    tGACATGGACAAGTTAAAGCAAATAAATGACACCCTCGGTCATCAGGAGGGGGATTTAGCGATTACGTTATTAGCTCAGGCGATTAAACAGTCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCt                                                                                                   <  2:2374926/150‑1 (MQ=255)
                                                      aCATGGACAAGTTAAAGCAAATAAATGACACCCTCGGTCATCAGGAGGGGGATTTAGCGATTACGTTATTAGCTCAGGCGATTAAACAGTCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTg                                                                                                 <  1:1364881/150‑1 (MQ=255)
                                                       cATGGACAAGTTAAAGCAAATAAATGACACCCTCGGTCATCAGGAGGGGGATTTAGCGATTACGTTATTAGCTCAGGCGATTAAACAGTCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGt                                                                                                >  1:3197998/1‑150 (MQ=255)
                                                         tGGACAAGTTAAAGCAAATAAATGACACCCTCGGTCATCAGGAGGGGGATTTAGCGATTACGTTATTAGCTCAGGCGATTAAACAGTCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCg                                                                                              >  1:2603990/1‑150 (MQ=255)
                                                            aCAAGTTAAAGCAAATAAATGACACCCTCGGTCATCAGGAGGGGGATTTAGCGATTACGTTATTAGCTCAGGCGATTAAACAGTCGATTAGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGAtt                                                                                           <  1:1889398/150‑1 (MQ=255)
                                                             cAAGTTAAAGCAAATAAATGACACCCTCGGTCATCAGGAGGGGGATTTAGCGATTACGTTATTAGCTCAGGCGATTAAACAGTCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTc                                                                                          <  2:1641095/150‑1 (MQ=255)
                                                             cAAGTTAAAGCAAATAAATGACACCCTCGGTAATCAGGAGGGGGATTTAGCGATTACGTTATTAGCTCAGGCGATTAAACAGTCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTc                                                                                          >  2:2470939/1‑150 (MQ=255)
                                                              aaGTTAAAGCAAATAAATGACACCCTCGGTCATCAGGAGGGGGATTTAGCGATTACGTTATTAGCTCAGGCGATTAAACAGTCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCg                                                                                         >  2:626188/1‑150 (MQ=255)
                                                               aGTTAAAGCAAATAAATGACACCCTCGGTCATCAGGAGGGGGATTTAGCGATTACGTTATTAGCTCAGGCGATTAAACAGTCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGa                                                                                        >  2:2670571/1‑150 (MQ=255)
                                                                  tAAAGCAAATAAATGACACCCTCGGTCATCAGGAGGGGGATTTAGCGATTACGTTATTAGCTCAGGCGATTAAACAGTCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGAcgc                                                                                     >  1:103258/1‑150 (MQ=255)
                                                                   aaaGCAAATAAATGACACCCTCGGTCATCAGGAGGGGGATTTAGCGATTACGTTATTAGCTCAGGCGATTAAACAGTCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGAcgcc                                                                                    >  1:328295/1‑150 (MQ=255)
                                                                    aaGCAAATAAATGACACCCTCGGTCATCAGGAGGGGGATTTAGCGATTACGTTATTAGCTCAGGCGATTAAACAGTCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGAcgccg                                                                                   >  1:3657823/1‑150 (MQ=255)
                                                                    aaGCAAATAAATGACACCCTCGGTCATCAGGAGGGGGATTTAGCGATTACGTTATTAGCTCAGGCGATTAAACAGTCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGAcgccg                                                                                   <  2:4315704/150‑1 (MQ=255)
                                                                      gCAAATAAATGACACCCTCGGTCATCAGGAGGGGGATTTAGCGATTACGTTATTAGCTCAGGCGATTAAACAGTCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCa                                                                                 >  2:4717524/1‑150 (MQ=255)
                                                                        aaataaatGACACCCTCGGTCATCAGGAGGGGGATTTAGCGATTACGTTATTAGCTCAGGCGATTAAACAGTCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCaaa                                                                               >  1:2294726/1‑150 (MQ=255)
                                                                           taaatGACACCCTCGGTCATCAGGAGGGGGATTTAGCGATTACGTTATTAGCTCAGGCGATTAAACAGTCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTg                                                                            <  2:3959261/150‑1 (MQ=255)
                                                                           taaatGACACCCTCGGTCATCAGGAGGGGGATTTAGCGATTACGTTATTAGCTCAGGCGATTAAACAGTCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTg                                                                            >  1:62062/1‑150 (MQ=255)
                                                                           taaatGACACCCTCGGTCATCAGGAGGGGGATTTAGCGATTACGTTATTAGCTCAGGCGATTAAACAGTCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTg                                                                            >  2:3372292/1‑150 (MQ=255)
                                                                               tGACACCCTCGGTCATAAGGAGGGGGATTTAGCGATTACGTTATTAGCTCAGGCGATTAAACAGTCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTgcagc                                                                        <  1:2886120/150‑1 (MQ=255)
                                                                                gACACCCTCGGTCATCAGGAGGGGGATTTAGCGATTACGTTATTAGCTCAGGCGATTAAACAGTCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTgcagca                                                                       >  1:2173516/1‑150 (MQ=255)
                                                                                  cacCCTCGGTCATCAGGAGGGGGATTTAGCGATTACGTTATTAGCTCAGGCGATTAAACAGTCGATTAGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGcaca                                                                     <  2:3293712/150‑1 (MQ=255)
                                                                                     ccTCGGTCATCAGGAGGGGGATTTAGCGATTACGTTATTAGCTCAGGCGATTAAACAGTCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACt                                                                  <  2:1881522/150‑1 (MQ=255)
                                                                                      cTCGGTCATAAGGAGGGGGATTTAGCGATTACGTTATTAGCTCAGGCGATTAAACAGTCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTg                                                                 >  1:1133974/1‑150 (MQ=255)
                                                                                         ggTCATCAGGAGGGGGATTTAGCGATTACGTTATTAGCTCAGGCGATTAAACAGTCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCt                                                              >  1:4376822/1‑150 (MQ=255)
                                                                                         ggTCATCAGGAGGGGGATTTAGCGATTACGTTATTAGCTCAGGCGATTAAACAGTCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCt                                                              >  2:4832179/1‑150 (MQ=255)
                                                                                           tcatcaGGAGGGGGATTTAGCGATTACGTTATTAGCTCAGGCGATTAAACAGTCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGa                                                            <  1:2131243/150‑1 (MQ=255)
                                                                                            catcaGGAGGGGGATTTAGCGATTACGTTATTAGCTCAGGCGATTAAACAGTCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGaa                                                           >  2:1318652/1‑150 (MQ=255)
                                                                                            catcaGGAGGGGGATTTAGCGATTACGTTATTAGCTCAGGCGATTAAACAGTCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGaa                                                           <  1:4141468/150‑1 (MQ=255)
                                                                                            catcaGGAGGGGGATTTAGCGATTACGTTATTAGCTCAGGCGATTAAACAGTCGATTAGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGaa                                                           <  2:3093457/150‑1 (MQ=255)
                                                                                               caGGAGGGGGATTTAGCGATTACGTTATTAGCTCAGGCGATTAAACAGTCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGt                                                        >  2:2161197/1‑150 (MQ=255)
                                                                                                aggaggGGGATTTAGCGATTACGTTATTAGCTCAGGCGATTAAACAGTCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTa                                                       >  2:4918524/1‑150 (MQ=255)
                                                                                                aggaggGGGATTTAGCGATTACGTTATTAGCTCAGGCGATTAAACAGTCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTa                                                       >  2:2942225/1‑150 (MQ=255)
                                                                                                 ggaggGGGAGTTAGCGATTACGTTATTAGCTCAGGCGATTAAACAGTCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTAt                                                      <  2:1695384/150‑1 (MQ=255)
                                                                                                     ggggATTTAGCGATTACGTTATTAGCTCAGGCGATTAAACAGTCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGaa                                                  >  1:1827652/1‑150 (MQ=255)
                                                                                                      gggATTTAGCGATTACGTTATTAGCTCAGGCGATTAAACAGTCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGaaa                                                 >  1:1358000/1‑150 (MQ=255)
                                                                                                      gggATTTAGCGATTACGTTATTAGCTCAGGCGATTAAACAGTCGATTAGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATGATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGaaa                                                 <  1:3797083/150‑1 (MQ=255)
                                                                                                                gATTACGTTATTAGCTCAGGCGATTAAACAGTCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAACGTCTGCa                                       >  2:4381850/1‑150 (MQ=255)
                                                                                                                  ttACGTTATTAGCTCAGGCGATTAAACAGTCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAACGTCTGCAAc                                     <  2:916867/150‑1 (MQ=255)
                                                                                                                  ttACGTTATTAGCTCAGGCGATTAAACAGTCGATTAGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAACGTCTGCAAc                                     <  1:1631414/150‑1 (MQ=255)
                                                                                                                     cGTTATTAGCTCAGGCGATTAAACAGTCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAACGTCTGCAACata                                  <  2:555886/150‑1 (MQ=255)
                                                                                                                     cGTTATTAGCTCAGGCGATTAAACAGTCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAACGTCTGCAACata                                  >  2:2853978/1‑150 (MQ=255)
                                                                                                                      gTTATTAGCTCAGGCGATTAAACAGTCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAACGTCTGCAACatat                                 <  1:2906164/150‑1 (MQ=255)
                                                                                                                       ttattaGCTCAGGCGATTAAACAGTCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAACGTCTGCAACATATc                                <  2:843672/150‑1 (MQ=255)
                                                                                                                       ttattaGCTCAGGCGATTAAACAGTCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAACGTCTGCAACATATc                                <  1:913802/150‑1 (MQ=255)
                                                                                                                        tattaGCTCAGGCGATTAAACAGTCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAACGTCTGCAACATATcg                               <  1:2400916/150‑1 (MQ=255)
                                                                                                                        tattaGCTCAGGCGATTAAACAGTCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAACGTCTGCAACATATcg                               >  2:1109310/1‑150 (MQ=255)
                                                                                                                          ttaGCTCAGGCGATTAAACAGTCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCGGCATCATTCTTGTCGATTCGACGCCGCTAATTGCAGCACAACTGCCTGAACGTATCGAAAAACGTCTGCAACATATcgcg                             <  1:2562934/150‑1 (MQ=255)
                                                                                                                           taGCTCAGGCGATTAAACAGTCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAACGTCTGCAACATATCgcgc                            >  1:4151883/1‑150 (MQ=255)
                                                                                                                              cTCAGGCGATTAAACAGTCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAACGTCTGCAACATATCGcgccgc                         >  1:1767158/1‑150 (MQ=255)
                                                                                                                               tCAGGCGATTAAACAGTCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAACGTCTGCAACATATCGCGCCGCa                        >  2:3471661/1‑150 (MQ=255)
                                                                                                                                cAGGCGATTAAACAGTCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAACGTCTGCAACATATCGCGCCGCAg                       <  2:4502278/150‑1 (MQ=255)
                                                                                                                                cAGGCGATTAAACAGTCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAACGTCTGCAACATATCGCGCCGCAg                       >  2:2074997/1‑150 (MQ=255)
                                                                                                                                    cGATTAAACAGTCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAACGTCTGCAACATATCGCGCCGCAGAAag                   >  1:3817414/1‑150 (MQ=255)
                                                                                                                                    cGATTAAACAGTCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAACGTCTGCAACATATCGCGCCGCAGAAag                   >  1:3564106/1‑150 (MQ=255)
                                                                                                                                     gATTAAACAGTCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAACGTCTGCAACATATCGCGCCGCAGAAaga                  >  1:528232/1‑150 (MQ=255)
                                                                                                                                     gATTAAACAGTCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAACGTCTGCAACATATCGCGCCGCAGAAaga                  >  2:4606133/1‑150 (MQ=255)
                                                                                                                                       ttAAACAGTCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAACGTCTGCAACATATCGCGCCGCAGAAAgaga                >  2:826017/1‑150 (MQ=255)
                                                                                                                                        tAAACAGTCGATTCGTAAAAGTGAGTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAACGTCTGCAACATATCGCGCCGCAGAAAGAGAt               <  1:3711238/150‑1 (MQ=255)
                                                                                                                                         aaaCAGTCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAACGTCTGCAACATATCGCGCCGCAGAAAGAGATc              <  2:3454673/150‑1 (MQ=255)
                                                                                                                                              gTCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAACGTCTGCAACATATCGCGCCGCAGAAAGAGATCGGCtt         <  2:2704920/150‑1 (MQ=255)
                                                                                                                                              gTCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAACGTCTGCAACATATCGCGCCGCAGAAAGAGATCGGCtt         <  2:2578955/150‑1 (MQ=255)
                                                                                                                                              gTCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAACGTCTGCAACATATCGCGCCGCAGAAAGAGATCGGCtt         <  1:1449816/150‑1 (MQ=255)
                                                                                                                                              gTCGATTAGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAACGTCTGCAACATATCGCGCCGCAGAAAGAGATCGGCtt         >  1:3834968/1‑150 (MQ=255)
                                                                                                                                               tCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAACGTCTGCAACATATCGCGCCGCAGAAAGAGATCGGCTtc        >  1:4717955/1‑150 (MQ=255)
                                                                                                                                                cGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAACGTCTGCAACATATCGCGCCGCAGAAAGAGATCGGCTtct       >  1:2628664/1‑150 (MQ=255)
                                                                                                                                                cGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAACGTCTGCAACATATCGCGCCGCAGAAAGAGATCGGCTtct       >  1:400604/1‑150 (MQ=255)
                                                                                                                                                 gATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAACGTCTGCAACATATCGCGCCGCAGAAAGAGATCGGCTtctc      >  2:4976888/1‑150 (MQ=255)
                                                                                                                                                   ttCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAACGTCTGCAACATATCGCGCCGCAGAAAGAGATCGGCTTCTCtt    <  1:2797128/150‑1 (MQ=255)
                                                                                                                                                     cGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAACGTCTGCAACATATCGCGCCGCAGAAAGAGATCGGCTTCTCTTcc  <  1:2375674/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AGCGGTTGCAGAAACTGGTGCAATCCGGTTCTTCGGTGATGTTTATTGCTATTGACATGGACAAGTTAAAGCAAATAAATGACACCCTCGGTCATCAGGAGGGGGATTTAGCGATTACGTTATTAGCTCAGGCGATTAAACAGTCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAACGTCTGCAACATATCGCGCCGCAGAAAGAGATCGGCTTCTCTTCC  >  CP009273/1867388‑1867686

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: