Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,839,034 C→A 8.1% 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,0340CA 8.1% 963.3 / 37.4 270Q571K (CAG→AAG) adecryptic adenine deaminase
Reads supporting (aligned to +/- strand):  ref base C (121/126);  new base A (11/12);  total (132/138)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 272 total aligned reads displayed.

GATGAGCACCGACACGGCGCAGTCGCTGGCGGAACAAATTGACGCCTTGAAAGCCGCCGCCCGTGAATGCGGTCCGTTACCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTT  >  CP009273/3838886‑3839180
                                                                                                                                                    |                                                                                                                                                  
gatgaGCACCGACACGGCGCAGTCGCTGGCGGAACAAATTGACGCCTTGAAAGCCGCCGCCCGTGAATGCGGTCCGTTACCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTc                                                                                                                                                    <  1:3150118/149‑1 (MQ=255)
 atgaGCACCGACACGGCGCAGTCGCTGGCGGAACAAATTGACGCCTTGAAAGCCGCCGCCCGTGAATGCGGTCCGTTACCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCa                                                                                                                                                   <  1:1927524/149‑1 (MQ=255)
   gaGCACCGACACGGCGCAGTCGCTGGCGGAACAAATTGACGCCTTGAAAGCCGCCGCCCGTGAATGCGGTCCGTTACCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAgg                                                                                                                                                 <  2:169949/149‑1 (MQ=255)
    aGCACCGACACGGCGCAGTCGCTGGCGGAACAAATTGACGCCTTGAAAGCCGCCGCCCGTGAATGCGGTCCGTTACCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAggg                                                                                                                                                <  1:1068065/149‑1 (MQ=255)
    aGCACCGACACGGCGCAGTCGCTGGCGGAACAAATTGACGCCTTGAAAGCCGCCGCCCGTGAATGCGGTCCGTTACCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAggg                                                                                                                                                <  2:2791406/149‑1 (MQ=255)
     gCACCGACACGGCGCAGTCGCTGGCGGAACAAATTGACGCCTTGAAAGCCGCCGCCCGTGAATGCGGTCCGTTACCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAgggg                                                                                                                                               <  1:3853582/149‑1 (MQ=255)
       aCCGACACGGCGCAGTCGCTGGCGGAACAAATTGACGCCTTGAAAGCCGCCGCCCGTGAATGCGGTCCGTTACCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCt                                                                                                                                             >  1:3897108/1‑149 (MQ=255)
          gACACGGCGCAGTCGCTGGCGGAACAAATTGACGCCTTGAAAGCCGCCGCCCGTGAATGCGGTCCGTGACCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTAtt                                                                                                                                          <  2:3336314/149‑1 (MQ=255)
             acGGCGCAGTCGCTGGCGGAACAAATTGACGCCTTGAAAGCCGCCGCCCGTGAATGCGGTCCGTTACCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGa                                                                                                                                       <  2:2422051/149‑1 (MQ=255)
                    aGTCGCTGGCGGAACAAATTGACGCCTTGAAAGCCGCCGCCCGTGAATGCGGTCCGTTACCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTTCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCgag                                                                                                                                >  1:4441074/1‑149 (MQ=255)
                    aGTCGCTGGCGGAACAAATTGACGCCTTGAAAGCCGCCGCCCGTGAATGCGGTCCGTTACCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCgag                                                                                                                                >  1:3434032/1‑149 (MQ=255)
                    aGTCGCTGGCGGAACAAATTGACGCCTTGAAAGCCGCCGCCCGTGAATGCGGTCCGTTACCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCgag                                                                                                                                >  1:4231176/1‑149 (MQ=255)
                    aGTCGCTGGCGGAACAAATTGACGCCTTGAAAGCCGCCGCCCGTGAATGCGGTCCGTTACCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCgag                                                                                                                                >  1:2049753/1‑149 (MQ=255)
                       cGCTGGCGGAACAAATTGACGCCTTGAAAGCCGCCGCCCGTGAATGCGGTCCGTTACCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAg                                                                                                                             <  2:767661/149‑1 (MQ=255)
                          tGGCGGAACAAATTGACGCCTTGAAAGCCGCCGCCCGTGAATGCGGTCCGTTACCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGttt                                                                                                                          <  2:2852335/149‑1 (MQ=255)
                            gCGGAACAAATTGACGCCTTGAAAGCCGCCGCCCGTGAATGCGGTCCGTTACCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGc                                                                                                                        >  2:1468023/1‑149 (MQ=255)
                            gCGGAACAAATTGACGCCTTGAAAGCCGCCGCCCGTGAATGCGGTCCGTTACCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTAAGGGGCTATTTGATGGCGAGAAGTTTGc                                                                                                                        <  1:3670453/149‑1 (MQ=255)
                            gCGGAACAAATTGACGCCTTGAAAGCCGCCGCCCGTGAATGCGGTCCGTTACCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTAAGGGGCTATTTGATGGCGAGAAGTTTGc                                                                                                                        <  1:3670397/149‑1 (MQ=255)
                              ggAACAAATTGACGCCTTGAAAGCCGCCGCCCGTGAATGCGGTCCGTTACCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCt                                                                                                                      <  2:1194277/149‑1 (MQ=255)
                              ggAACAAATTGACGCCTTGAAAGCCGCCGCCCGTGAATGCGGTCCGTTACCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCt                                                                                                                      <  1:2259101/149‑1 (MQ=255)
                               gAACAAATTGACGCCTTGAAAGCCGCCGCCCGTGAATGCGGTCCGTTACCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCtt                                                                                                                     >  2:1226203/1‑149 (MQ=255)
                               gAACAAATTGACGCCTTGAAAGCCGCCGCCCGTGAATGCGGTCCGTTACCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCtt                                                                                                                     >  1:73125/1‑149 (MQ=255)
                               gAACAAATTGACGCCTTGAAAGCCGCCGCCCGTGAATGCGGTCCGTTACCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCtt                                                                                                                     >  1:2734848/1‑149 (MQ=255)
                                    aaTTGCCGCCTTGAAAGCCGCCGCCCGTGAATGCGGTCCGTTACCCGACGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCGAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTa                                                                                                                <  1:804196/149‑1 (MQ=255)
                                     aTTGACGCCTTGAAAGCCGCCGCCCGTGAATGCGGTCCGTTACCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTAc                                                                                                               >  1:2622296/1‑149 (MQ=255)
                                     aTTGACGCCTTGAAAGCCGCCGCCCGTGAATGCGGTCCGTTACCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTAc                                                                                                               >  1:4310718/1‑149 (MQ=255)
                                        gACGCCTTGAAAGCCGCCGCCCGTGAATGCGGTCCGTTACCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCt                                                                                                            >  2:899222/1‑149 (MQ=255)
                                        gACGCCTTGAAAGCCGCCGCCCGTGAATGCGGTCCGTTACCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCt                                                                                                            >  2:322126/1‑149 (MQ=255)
                                        gACGCCTTGAAAGCCGCCGCCCGTGAATGCGGTCCGTTACCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCt                                                                                                            >  2:2314990/1‑149 (MQ=255)
                                         aCGCCTTGAAAGCCGCCGCCCGTGAATGCGGTCCGTTACCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTAAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTg                                                                                                           >  2:3730377/1‑149 (MQ=255)
                                          cGCCTTGAAAGCCGCCGCCCGTGAATGCGGTCCGTTACCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTgg                                                                                                          <  1:775009/149‑1 (MQ=255)
                                          cGCCTTGAAAGCCGCCGCCCGTGAATGCGGTCCGTTACCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTgg                                                                                                          <  1:3365875/149‑1 (MQ=255)
                                           gCCTTGAAAGCCGCCGCCCGTGAATGCGGTCCGTTACCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGa                                                                                                         <  1:131835/149‑1 (MQ=255)
                                           gCCTTGAAAGCCGCCGCCCGTGAATGCGGTCCGTTACCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGa                                                                                                         >  1:3805113/1‑149 (MQ=255)
                                                gAAAGCCGCCGCCCGTGAATGCGGTCCGTTACCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCa                                                                                                    <  2:1814136/149‑1 (MQ=255)
                                                 aaaGCCGCCGCCCGTGAATGCGGTCCGTTACCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCAc                                                                                                   <  2:47046/149‑1 (MQ=255)
                                                      cgccgccCGTGAATGCGGTCCGTTACCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGaat                                                                                              >  1:675433/1‑149 (MQ=255)
                                                        ccgccCGTGAATGCGGTCCGTTACCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGaataa                                                                                            <  1:1497687/149‑1 (MQ=255)
                                                        ccgccCGTGAATGCGGTCCGTTACCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGaataa                                                                                            <  1:4704193/149‑1 (MQ=255)
                                                        ccgccCGTGAATGCGGTCCGTTACCCGATGAGCCGTTTATTAAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGATATTTGATGGCGAGAATTTTGCCTTCAATAAGCTGGAAGTCACGGaataa                                                                                            >  1:4362647/1‑149 (MQ=255)
                                                         cgccCGTGAATGCGGTCCGTTACCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGaataat                                                                                           <  2:702804/149‑1 (MQ=255)
                                                                gAATGCGGTCCGTTACCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATaaaaaaa                                                                                    >  2:1679286/1‑149 (MQ=255)
                                                                 aaTGCGGTCCGTTACCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAg                                                                                   <  1:1253257/149‑1 (MQ=255)
                                                                  aTGCGGTCCGTTACCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGc                                                                                  >  1:2355940/1‑149 (MQ=255)
                                                                   tGCGGTCCGTTACCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGcc                                                                                 <  1:1958202/149‑1 (MQ=255)
                                                                      ggTCCGTTACCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCgg                                                                              <  1:2540519/149‑1 (MQ=255)
                                                                        tCCGTTACCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCg                                                                            >  2:1798553/1‑149 (MQ=255)
                                                                          cGTTACCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTc                                                                          <  1:2913042/149‑1 (MQ=255)
                                                                          cGTTACCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTc                                                                          <  2:2988201/149‑1 (MQ=255)
                                                                             tACCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATg                                                                       >  2:4499876/1‑149 (MQ=255)
                                                                             tACCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTAAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATg                                                                       >  2:3687236/1‑149 (MQ=255)
                                                                              aCCCGATGAGCCGTTTAGTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGc                                                                      <  2:4668536/149‑1 (MQ=255)
                                                                               cccGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCAAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGcc                                                                     <  2:1988690/149‑1 (MQ=255)
                                                                                    tGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCa                                                                >  2:691776/1‑149 (MQ=255)
                                                                                    tGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCa                                                                >  2:2687728/1‑149 (MQ=255)
                                                                                    tGAGCCGTTTATTCAGATGGCGTTTCTTTATCTGCCAGTGATCCCCGCGCTAAAAAAAACCAGTCAGGGGCTATTTTATGGCGAGAAGTTTTCCATCACTACGCAGGAAGTCACGGAAAAATAAAAAAAGCCCCGCGTCATGCCGGGCa                                                                >  1:3737030/1‑149 (MQ=255)
                                                                                       gCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCaaaa                                                             <  2:4096147/149‑1 (MQ=255)
                                                                                        ccGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATATGATGGCGAGAAGTTTGCCTACACTACGCTGGAAGACACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAg                                                            >  2:2193011/1‑149 (MQ=255)
                                                                                          gTTTATTCAGATGGCGTTTCTTTCTCTGCCTGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTc                                                          <  1:2871798/149‑1 (MQ=255)
                                                                                          gTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTc                                                          >  1:962086/1‑149 (MQ=255)
                                                                                           tttATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCa                                                         <  1:2840845/149‑1 (MQ=255)
                                                                                           tttATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCa                                                         >  1:3979020/1‑149 (MQ=255)
                                                                                            ttATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCAc                                                        >  1:1648331/1‑149 (MQ=255)
                                                                                            ttATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCAc                                                        >  1:3318226/1‑149 (MQ=255)
                                                                                              aTTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCa                                                      <  2:4317950/149‑1 (MQ=255)
                                                                                                  aGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTa                                                  >  1:2494902/1‑149 (MQ=255)
                                                                                                  aGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTa                                                  <  2:1590105/149‑1 (MQ=255)
                                                                                                   gATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTAc                                                 >  2:560630/1‑149 (MQ=255)
                                                                                                   gATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTAc                                                 >  2:3376247/1‑149 (MQ=255)
                                                                                                   gATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTAc                                                 >  1:1089295/1‑149 (MQ=255)
                                                                                                    aTGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTAAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACg                                                <  1:3478537/149‑1 (MQ=255)
                                                                                                        cGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTa                                            <  1:4061407/149‑1 (MQ=255)
                                                                                                         gTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTAt                                           >  2:3151738/1‑149 (MQ=255)
                                                                                                         gTTTCTTTCTCTGCCAGTGATCCCCGCGATAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTAt                                           <  2:3994493/149‑1 (MQ=255)
                                                                                                            tctttctCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGcc                                        <  2:4655415/149‑1 (MQ=255)
                                                                                                            tctttctCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTAAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGcc                                        >  2:1977856/1‑149 (MQ=255)
                                                                                                              tttctCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCAc                                      <  1:2985452/149‑1 (MQ=255)
                                                                                                                 ctctGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGt                                   >  2:2226521/1‑149 (MQ=255)
                                                                                                                  tctGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTc                                  <  2:3897108/149‑1 (MQ=255)
                                                                                                                  tctGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGGTTATGCCACTGTc                                  >  1:3414128/1‑149 (MQ=255)
                                                                                                                   ctGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCa                                 <  1:527565/149‑1 (MQ=255)
                                                                                                                    tGCCAGTGATCCCCGCGCTAAAACTAACCAGTAAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCaa                                >  2:424807/1‑149 (MQ=255)
                                                                                                                       cAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGACCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAActg                             >  1:1023064/1‑149 (MQ=255)
                                                                                                                        aGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAActgc                            <  2:3583691/149‑1 (MQ=255)
                                                                                                                        aGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAActgc                            <  1:3710253/149‑1 (MQ=255)
                                                                                                                         gTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAActgct                           >  1:4778127/1‑149 (MQ=255)
                                                                                                                            aTCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCCCTGTCAACTGCTGaa                        >  2:3510300/1‑149 (MQ=255)
                                                                                                                            aTCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGaa                        >  2:1701637/1‑149 (MQ=255)
                                                                                                                              ccccGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAAtt                      >  2:4488059/1‑149 (MQ=255)
                                                                                                                                ccGCGCTAAAACTAACCAGTAAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAAtttt                    >  2:4672535/1‑149 (MQ=255)
                                                                                                                                    gcTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTtc                >  2:2825642/1‑149 (MQ=255)
                                                                                                                                    gcTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTtc                >  1:4331562/1‑149 (MQ=255)
                                                                                                                                     cTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTtct               >  2:2572976/1‑149 (MQ=255)
                                                                                                                                      tAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTtctc              >  2:1123858/1‑149 (MQ=255)
                                                                                                                                        aaaCTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTcgc            >  1:268460/1‑149 (MQ=255)
                                                                                                                                        aaaCTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTcgc            >  2:1266859/1‑149 (MQ=255)
                                                                                                                                           cTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCgcggc         <  1:2120530/149‑1 (MQ=255)
                                                                                                                                           cTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCgcggc         <  2:1260066/149‑1 (MQ=255)
                                                                                                                                           cTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCgcggc         >  1:745793/1‑149 (MQ=255)
                                                                                                                                               ccAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGAt     <  1:710361/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
GATGAGCACCGACACGGCGCAGTCGCTGGCGGAACAAATTGACGCCTTGAAAGCCGCCGCCCGTGAATGCGGTCCGTTACCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTT  >  CP009273/3838886‑3839180

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: