Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,839,034 C→A 21.5% Q571K (CAG→AAG)  ade → cryptic adenine deaminase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,839,0340CA21.5% 716.3 / 139.3 264Q571K (CAG→AAG) adecryptic adenine deaminase
Reads supporting (aligned to +/- strand):  ref base C (104/103);  new base A (21/36);  total (125/139)
Fisher's exact test for biased strand distribution p-value = 9.89e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

Only 100 of 272 total aligned reads displayed.

TGATGAGCACCGACACGGCGCAGTCGCTGGCGGAACAAATTGACGCCTTGAAAGCCGCCGCCCGTGAATGCGGTCCGTTACCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTTGC  >  CP009273/3838885‑3839182
                                                                                                                                                     |                                                                                                                                                    
tgatgaGCACCGACACGGCGCAGTCGCTGGCGGAACAAATTGACGCCTTGAAAGCCGCCGCCCGTGAATGCGGTCCGTTACCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTc                                                                                                                                                      >  2:2937564/1‑150 (MQ=255)
tgatgaGCACCGACACGGCGCAGTCGCTGGCGGAACAAATTGACGCCTTGAAAGCCGCCGCCCGTGAATGCGGTCCGTTACCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTc                                                                                                                                                      >  2:1558172/1‑150 (MQ=255)
  atgaGCACCGACACGGCGCAGTCGCTGGCGGAACAAATTGACGCCTTGAAAGCCGCCGCCCGTGAATGCGGTCCGTTACCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAg                                                                                                                                                    <  1:3182322/150‑1 (MQ=255)
   tgaGCACCGACACGGCGCAGTCGCTGGCGGAACAAATTGACGCCTTGAAAGCCGCCGCCCGTGAATGCGGTCCGTTACCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAgg                                                                                                                                                   <  1:2345296/150‑1 (MQ=255)
   tgaGCACCGACACGGCGCAGTCGCTGGCGGAACAAATTGACGCCTTGAAAGCCGCCGCCCGTGAATGCGGTCCGTTACCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAAACAGTCAgg                                                                                                                                                   <  2:3173304/150‑1 (MQ=255)
    gaGCACCGACACGGCGCAGTCGCTGGCGGAACAAATTGACGCCTTGAAAGCCGCCGCCCGTGAATGCGGTCCGTTACCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAggg                                                                                                                                                  <  1:3686018/150‑1 (MQ=255)
    gaGCACCGACACGGCGCAGTCGCTGGCGGAACAAATTGACGCCTTGAAAGCCGCCGCCCGTGAATGCGGTCCGTTACCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTAAggg                                                                                                                                                  >  1:1695843/1‑150 (MQ=255)
      gCACCGACACGGCGCAGTCGCTGGCGGAACAAATTGACGCCTTGAAAGCCGCCGCCCGTGAATGCGGTCCGTTACCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGc                                                                                                                                                <  2:3475316/150‑1 (MQ=255)
       cACCGACACGGCGCAGTCGCTGGCGGAACAAATTGACGCCTTGAAAGCCGCCGCCCGTGAATGCGGTCCGTTACCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCt                                                                                                                                               <  2:3736510/150‑1 (MQ=255)
         ccGACACGGCGCAGTCGCTGGCGGAACAAATTGACGCCTTGAAAGCCGCCGCCCGTGAATGCGGTCCGTTACCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTAAGGGGCTAt                                                                                                                                             <  2:3113518/150‑1 (MQ=255)
          cGACACGGCGCAGTCGCTGGCGGAACAAATTGACGCCTTGAAAGCCGCCGCCCGTGAATGCGGTCCGTTACCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTAtt                                                                                                                                            <  2:9192/150‑1 (MQ=255)
           gACACGGCGCAGTCGCTGGCGGAACAATTTGACGCCTTGAAAGCCGCCGCCCGTGAATGCGGTCCGTTACCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAACACTAACCAGTCAGGGGCTAttt                                                                                                                                           <  1:353808/150‑1 (MQ=255)
           gACACGGCGCAGTCGCTGGCGGAACAAATTGACGCCTTGAAAGCCGCCGCCCGTGAATGCGGTCCGTTACCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTAttt                                                                                                                                           <  1:1847658/150‑1 (MQ=255)
           gACACGGCGCAGTCGCTGGCGGAACAAATTGACGCCTTGAAAGCCGCCGCCCGTGAATGCGGTCCGTTACCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTAttt                                                                                                                                           >  1:40038/1‑150 (MQ=255)
           gACACGGCGCAGTCGCTGGCGGAACAAATTGACGCCTTGAAAGCCGCCGCCCGTGAATGCGGTCCGTTACCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTAttt                                                                                                                                           <  2:3728044/150‑1 (MQ=255)
           gACACGGCGCAGTCGCTGGCGGAACAAATTGACGCCTTGAAAGCCGCCGCCCGTGAATGCGGTCCGTTACCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTAAGGGGCTAttt                                                                                                                                           >  2:333708/1‑150 (MQ=255)
                 gcgcAGTCGCTGGCGGAACAAATTGACGCCTTGAAAGCCGCCGCCCGTGAATGCGGTCCGTTACCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTAAGGGGCTATTTGATGGc                                                                                                                                     <  1:702284/150‑1 (MQ=255)
                    cAGTCGCTGGCGGAACAAATTGACGCCTTGAAAGCCGCCGCCCGTGAATGCGGTCCGTTACCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCgag                                                                                                                                  >  1:4032026/1‑150 (MQ=255)
                     aGTCGCTGGCGGAACAAATTGACGCCTTGAAAGCCGCCGCCCGTGAATGCGGTCCGTTACCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCgaga                                                                                                                                 <  1:2218105/150‑1 (MQ=255)
                      gTCGCTGGCGGAACAAATTGACGCCTTGAAAGCCGCCGCCCGTGAATGCGGTCCGTTACCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGaa                                                                                                                                >  1:2590383/1‑150 (MQ=255)
                            ggCGGAACAAATTGACGCCTTGAAAGCCGCCGCCCGTGAATGCGGTCCGTTACCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTAAGGGGCTATTTGATGGCGAGAAGTTTGc                                                                                                                          <  2:4203705/150‑1 (MQ=255)
                             gCGGAACAAATTGACGCCTTGACAGCCGCCGCCCGTGAATGCGGTCCGTTACCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGcc                                                                                                                         <  2:597858/150‑1 (MQ=255)
                               ggAACAAATTGACGCCTTGAAAGCCGCCGCCCGTGAATGCGGTCCGTTACCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTAAGGGGCTATTTGATGGCGAGAAGTTTGCCtt                                                                                                                       >  2:1832206/1‑150 (MQ=255)
                               ggAACAAATTGACGCCTTGAAAGCCGCCGCCCGTGAATGCGGTCCGTTACCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTATCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCtt                                                                                                                       >  1:501981/1‑150 (MQ=255)
                                  aCAAATTGACGCCTTGAAAGCCGCCGCCCGTGAATGCGGTCCGTTACCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTAAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCAc                                                                                                                    <  1:253653/150‑1 (MQ=255)
                                      aTTGACGCCTTGAAAGCCGCCGCCCGTGAATGCGGTCCGTTACCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACg                                                                                                                <  2:4503310/150‑1 (MQ=255)
                                      aTTGACGCCTTGAAAGCCGCCGCCCGTGAATGCGGTCCGTTACCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACg                                                                                                                <  2:2995500/150‑1 (MQ=255)
                                       ttGACGCCTTGAAAGCCGCCGCCCGTGAATGCGGTCCGTTACCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTAAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGc                                                                                                               >  2:748641/1‑150 (MQ=255)
                                        tGACGCCTTGAAAGCCGCCGCCCGTGAATGCGGTCCGTTACCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTAAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCt                                                                                                              <  2:3501540/150‑1 (MQ=255)
                                         gACGCCTTGAAAGCCGCCGCCCGTGAATGCGGTCCGTTACCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTg                                                                                                             <  1:4392569/150‑1 (MQ=255)
                                         gACGCCTTGAAAGCCGCCGCCCGTGAATGCGGTCCGTTACCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCACGCGCTAAAACTAACCAGTAAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTg                                                                                                             >  1:2695156/1‑150 (MQ=255)
                                           cGCCTTGAAAGCCGCCGCCCGTGAATGCGGTCCGTTACCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTAAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGa                                                                                                           <  2:2007424/150‑1 (MQ=255)
                                              cTTGAAAGCCGCCGCCCGTGAATGCGGTCCGTTACCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGt                                                                                                        <  2:2221626/150‑1 (MQ=255)
                                               ttGAAAGCCGCCGCCCGTGAATGCGGTCCGTTACCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTAAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTc                                                                                                       <  1:820287/150‑1 (MQ=255)
                                                 gAAAGCCGCCGCCCGTGAATGCGGTCCGTTACCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCAc                                                                                                     <  1:3069509/150‑1 (MQ=255)
                                                 gAAAGCCGCCGCCCGTGAATGCGGTCCGTTACCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCAc                                                                                                     <  1:912419/150‑1 (MQ=255)
                                                   aaGCCGCCGCCCGTGAATGCGGTCCGTTACCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTAAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACgg                                                                                                   <  1:1140223/150‑1 (MQ=255)
                                                    aGCCGCCGCCCGTGAATGCGGTCCGTTACCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGa                                                                                                  <  2:4076430/150‑1 (MQ=255)
                                                            cccGTGACTGCGGTCCGTTACCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATaaa                                                                                          <  1:867830/150‑1 (MQ=255)
                                                               gTGAATGCGGTCCGTTACCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATaaaaaa                                                                                       >  2:2364899/1‑150 (MQ=255)
                                                               gTGAATGCGGTCCGTTACCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATaaaaaa                                                                                       >  1:4064491/1‑150 (MQ=255)
                                                               gTGAATGCGGTCAGTTACCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGGCAGTGATCACCGCGCTAAAACTAACCAGACAGGGGCTATTTGATGGCGAGAAGTTTGCCAGCACTACGCTGGAAGTCACGGAATAaaaaaaaa                                                                                       >  1:2115292/1‑150 (MQ=255)
                                                                tGAATGCGGTCCGTTACCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATaaaaaaa                                                                                      >  1:2155180/1‑150 (MQ=255)
                                                                   aTGCGGTCCGTTACCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGcc                                                                                   >  1:2491112/1‑150 (MQ=255)
                                                                    tGCGGTCCGTTACCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGccc                                                                                  >  1:3131403/1‑150 (MQ=255)
                                                                       ggTCCGTTACCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCTAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGc                                                                               <  2:1664499/150‑1 (MQ=255)
                                                                       ggTCCGTTACCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGc                                                                               >  1:3936682/1‑150 (MQ=255)
                                                                       ggTCCGTTACCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGc                                                                               >  1:4545858/1‑150 (MQ=255)
                                                                       ggTCCGTTACCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGc                                                                               >  1:3793765/1‑150 (MQ=255)
                                                                        gTCCGTTACCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTAAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCg                                                                              <  1:3850881/150‑1 (MQ=255)
                                                                         tCCGTTACCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGt                                                                             >  1:4229119/1‑150 (MQ=255)
                                                                         tCCGTTACACGTTGAGCCGTTTTTTCATATGGAGTTTCTTTCTCTGCCAGTGATCACCGCTCTAAAAATAACCAGTCATGGGCTATTTGATGGCGATATGTTTGCCTTCAATACGCTGGAAGTCACGGAATAATTATAAAAGCCCGGCGt                                                                             <  2:2767454/150‑1 (MQ=255)
                                                                              tACCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTAAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGc                                                                        <  2:1120057/150‑1 (MQ=255)
                                                                                 ccGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCgg                                                                     <  1:2651399/150‑1 (MQ=255)
                                                                                 ccGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCgg                                                                     >  1:3322019/1‑150 (MQ=255)
                                                                                  cGATGAGCCGTTTATTCAGATGGCGTTTCATTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTAAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCggg                                                                    >  2:3669097/1‑150 (MQ=255)
                                                                                        gCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAg                                                              >  1:1474299/1‑150 (MQ=255)
                                                                                        gCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAg                                                              >  1:1514953/1‑150 (MQ=255)
                                                                                           gTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCa                                                           <  1:1511806/150‑1 (MQ=255)
                                                                                              tATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCa                                                        >  1:3476466/1‑150 (MQ=255)
                                                                                                 tCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAATTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGtt                                                     <  2:1170556/150‑1 (MQ=255)
                                                                                                 tCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGtt                                                     <  1:2965213/150‑1 (MQ=255)
                                                                                                  cAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTAAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTa                                                    <  1:3860628/150‑1 (MQ=255)
                                                                                                      tGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTAAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGtt                                                <  1:1241947/150‑1 (MQ=255)
                                                                                                       ggCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGttt                                               >  2:676627/1‑150 (MQ=255)
                                                                                                         cGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTAt                                             <  1:4426294/150‑1 (MQ=255)
                                                                                                            ttctttctCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGcc                                          >  1:1960305/1‑150 (MQ=255)
                                                                                                            ttctttctCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTAAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGcc                                          >  2:3403977/1‑150 (MQ=255)
                                                                                                             tctttctCTGCCAGTGATCCCCGCGCTAAACCTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCa                                         <  1:1815182/150‑1 (MQ=255)
                                                                                                               tttctCTGCCAGTGATCCCCGCGCTAAAACTACCCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACt                                       <  2:3950266/150‑1 (MQ=255)
                                                                                                                   tctGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCa                                   <  2:248328/150‑1 (MQ=255)
                                                                                                                    ctGCCAGTGATCCCCGCGCTAAAACTAACCAGTAAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCaa                                  <  2:1765422/150‑1 (MQ=255)
                                                                                                                     tGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAAc                                 <  2:942622/150‑1 (MQ=255)
                                                                                                                     tGCCAGTGATCCCCGCGCTAAAACTAACCAGTAAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAAc                                 >  1:2976798/1‑150 (MQ=255)
                                                                                                                      gCCAGTGATCCCTGCGCTAAAACTAACCAGTAAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAAct                                <  1:1189208/150‑1 (MQ=255)
                                                                                                                      gCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAAct                                <  1:2881841/150‑1 (MQ=255)
                                                                                                                      gCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAAct                                <  1:4495402/150‑1 (MQ=255)
                                                                                                                          gTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAActgctg                            <  1:1341046/150‑1 (MQ=255)
                                                                                                                          gTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAActgctg                            <  2:1293228/150‑1 (MQ=255)
                                                                                                                              tCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAAtt                        <  1:3345894/150‑1 (MQ=255)
                                                                                                                              tCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAAtt                        <  1:3345974/150‑1 (MQ=255)
                                                                                                                              tCCCCGCGCTAAAACTAACCAGTAAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAAtt                        <  2:4222934/150‑1 (MQ=255)
                                                                                                                                cccGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAAtttt                      >  1:1072607/1‑150 (MQ=255)
                                                                                                                                 ccGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAAttttt                     >  2:1450707/1‑150 (MQ=255)
                                                                                                                                     gcTAAAACTAACCAGTAAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTtct                 <  2:2045381/150‑1 (MQ=255)
                                                                                                                                      cTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTtctc                >  2:197690/1‑150 (MQ=255)
                                                                                                                                       tAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTcg               <  1:731870/150‑1 (MQ=255)
                                                                                                                                        aaaaCTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTcgc              >  2:4577562/1‑150 (MQ=255)
                                                                                                                                         aaaCTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTcgcg             >  1:1580337/1‑150 (MQ=255)
                                                                                                                                         aaaCTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTcgcg             >  1:3148295/1‑150 (MQ=255)
                                                                                                                                         aaaCTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTcgcg             <  2:3030558/150‑1 (MQ=255)
                                                                                                                                         aaaCTAACCAGTCAGGGGATATTTGATGGCGATAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGACCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGAAACTGTCAAATGCAGAATTTTTTTCTcgcg             >  1:953933/1‑150 (MQ=255)
                                                                                                                                            cTAACCAGTAAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCgcggcg          <  2:2695156/150‑1 (MQ=255)
                                                                                                                                              aaCCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGa        >  1:1103285/1‑150 (MQ=255)
                                                                                                                                               aCCAGTAAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGAt       <  1:27598/150‑1 (MQ=255)
                                                                                                                                                ccAGTAAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGAtt      >  2:370327/1‑150 (MQ=255)
                                                                                                                                                 cAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGAttt     <  2:4521366/150‑1 (MQ=255)
                                                                                                                                                 cAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGAttt     >  1:666834/1‑150 (MQ=255)
                                                                                                                                                   gTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTTg   <  2:2513999/150‑1 (MQ=255)
                                                                                                                                                    tCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTTGc  >  1:3427361/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
TGATGAGCACCGACACGGCGCAGTCGCTGGCGGAACAAATTGACGCCTTGAAAGCCGCCGCCCGTGAATGCGGTCCGTTACCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTTGC  >  CP009273/3838885‑3839182

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: