Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,838,587 T→A 12.9% S422T (TCC→ACC)  ade → cryptic adenine deaminase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,838,5870TA12.9% 942.7 / 109.6 310S422T (TCC→ACC) adecryptic adenine deaminase
Reads supporting (aligned to +/- strand):  ref base T (112/158);  new base A (15/25);  total (127/183)
Fisher's exact test for biased strand distribution p-value = 7.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 331 total aligned reads displayed.

AAGAGTCGGCGAGACTGGCACAATCCGCTCCGCCATATGGCAACACCATTGCCCGCCAGCCAGTTTCCGCCAGCGACTTTGCCCTGCAATTTACGCCCGGAAAACGCTATCGGGTCATTGACGTCATCCATAACGAATTGATTACGCACTCCCACTCCAGCGTCTACAGCGAAAATGGTTTTGATCGCGATGATGTGAGCTTTATTGCCGTACTTGAGCGTTACGGGCAACGGCTGGCTCCGGCTTGTGGTTTGCTTGGCGGCTTTGGACTGAATGAAGGTGCGCTGGCTGCGACG  >  CP009273/3838438‑3838733
                                                                                                                                                     |                                                                                                                                                  
aaGAGTCGGCGAGACTGGCACAATCCGCTCCGCCATATGGCAACACCATTGCCCGCCAGCCAGTTTCCGCCAGCGACTTTGCCCTGCAATTTACGCCCGGAAAACGCTATCGGGTCATTGACGTCATCCATAACGAATTGATTACGcact                                                                                                                                                    >  2:547664/1‑150 (MQ=255)
 agagTCGGCGAGACTGGCACAATCCGCTCCGCCATATGGCAACACCATTGCCCGCCAGCCAGTTTCCGCCAGCGACTTTGCCCTGCAATTTACGCCCGGAAAACGCTATCGGGTCATTGACGTCATCCATAACGAATTGATTACGcactc                                                                                                                                                   >  2:139270/1‑150 (MQ=255)
  gagTCGGCGAGACTGGCACAATCCGCTCCGCCATATGGCAACACCATTGCCCGCCAGCCAGTTTCCGCCAGCGACTTTGCCCTGCAATTTACGCCCGGAAAACGCTATCGGGTCATTGACGTCATCCATAACGAATTGATTACGcactcc                                                                                                                                                  <  1:1200952/150‑1 (MQ=255)
          gagaCTGGCACAATCCGCTCCGCCATATGGCAACACCATTGCCCGCCAGCCAGTTTCCGCCAGCGACTTTGCCCTGCAATTTACGCCCGGAAAACGCTATCGGGTCATTGACGTCATCCATAACGAATTGATTACGCACTCCCACTCCAg                                                                                                                                          <  2:4247241/150‑1 (MQ=255)
          gagaCTGGCACAATCCGCTCCGCCATATGGCAACACCATTGCCCGCCAGCCAGTTTCCGCCAGCGACTTTGCCCTGCAATTTACGCCCGGAAAACGCTATCGGGTCATTGACGTCATCCATAACGAATTGATTACGCACTCCCACTCCAg                                                                                                                                          <  2:3274916/150‑1 (MQ=255)
          gagaCTGGCACAATCCGCTCCGCCATATGGCAACACCATTGCCCGCCAGCCAGTTTCCGCCAGCGACTTTGCCCTGCAATTTACGCCCGGAAAACGCTATCGGGTCATTGACGTCATCCATAACGAATTGATTACGCACTCCCACTCCAg                                                                                                                                          <  2:3126283/150‑1 (MQ=255)
          gagaCTGGCACAATCCGCTCCGCCATATGGCAACACCATTGCCCGCCAGCCAGTTTCCGCCAGCGACTTTGCCCTGCAATTTACGCCCGGAAAACGCTATCGGGTCATTGACGTCATCCATAACGAATTGATTACGCACTCCCACTCCAg                                                                                                                                          >  2:2572102/1‑150 (MQ=255)
          gagaCTGGCACAATCCGCTCCGCCATATGGCAACACCATTGCCCGCCAGCCAGTTTCCGCCAGCGACTTTGCCCTGCAATTTACGCCCGGAAAACGCTATCGGGTCATTGACGTCATCCATAACGAATTGATTACGCACTCCCACTCCAg                                                                                                                                          >  1:2074605/1‑150 (MQ=255)
           agaCTGGCACAATCCGCTCCGCCATATGGCAACACCATTGCCCGCCAGCCAGTTTCCGCCAGCGACTTTGCCCTGCAATTTACGCCCGGAAAACGCTATCGGGTCATTGACGTCATCCATAACGAATTGATTACGCACTCCCACTCCAGc                                                                                                                                         <  2:1179955/150‑1 (MQ=255)
           agaCTGGCACAATCCGCTCCGCCATATGGCAACACCATTGCCCGCCAGCCAGTTTCCGCCAGCGACTTTGCCCTGCAATTTACGCCCGGAAAACGCTATCGGGTCATTGACGTCATCCATAACGAATTGATTACGCACTCCCACTCCAGc                                                                                                                                         <  1:4238369/150‑1 (MQ=255)
            gaCTGGCACAATCCGCTCCGCCATATGGCAACACCATTGCCCGCCAGCCAGTTTCCGCCAGCGACTTTGCCCTGCAATTTACGCCCGGAAAACGCTATCGGGTCATTGACGTCATCCATAACGAATTGATTACGCACTCCCACTCCAGCg                                                                                                                                        >  2:236001/1‑150 (MQ=255)
             aCTGGCACAATCCGCTCCGCCATATGGCAACACCATTGCCCGCCAGCCAGTTTCCGCCAGCGACTTTGCCCTGCAATTTACGCCCGGAAAACGCTATCGGGTCATTGACGTCATCCATAACGAATTGATTACGCACTCCCACTCCAGCGt                                                                                                                                       <  2:1957963/150‑1 (MQ=255)
               tGGCACAATCCGCTCCGCCATATGGCAACACCATTGCCCGCCAGCCAGTTTCCGCCAGCGACTTTGCCCTGCAATTTACGCCCGGAAAACGCTATCGGGTCATTGACGTCATCCATAACGAATTGATTACGCACTCCCACTCCAGCGTCt                                                                                                                                     >  2:4596103/1‑150 (MQ=255)
                ggCACAATCCGCTCCGCCATATGGCAACACCATTGCCCGCCAGCCAGTTTCCGCCAGCGACTTTGCCCTGCAATTTACGCCCGGAAAACGCTATCGGGTCATTGACGTCATCCATAACGAATTGATTACGCACTCCCACTCCAGCGTCTa                                                                                                                                    >  2:907219/1‑150 (MQ=255)
                ggCACAATCCGCTCCGCCATATGGCAACACCATTGCCCGCCAGCCAGTTTCCGCCAGCGACTTTGCCCTGCAATTTACGCCCGGAAAACGCTATCGGGTCATTGACGTCATCCATAACGAATTGATTACGCACTCCCACTCCAGCGTCTa                                                                                                                                    <  1:3899809/150‑1 (MQ=255)
                 gCACAATCCGCTCCGCCATATGGCAACACCATTGCCCGCCAGCCAGTTTCCGCCAGCGACTTTGCCCTGCAATTTACGCCCGGAAAACGCTATCGGGTCATTGACGTCATCCATAACGAATTGATTACGCACTCCCACTCCAGCGTCTAc                                                                                                                                   <  2:3769093/150‑1 (MQ=255)
                  cacaATCCGCTCCGCCATATGGCAACACCATTGCCCGCCAGCCAGTTTCCGCCAGCGACTTTGCCCTGCAATTTACGCCCGGAAAACGCTATCGGGTCATTGACGTCATCCATAACGAATTGATTACGCACTCCCACTCCAGCGTCTACa                                                                                                                                  <  2:248322/150‑1 (MQ=255)
                     aaTCCGCTCCGCCATATGGCAACACCATTGCCCGCCAGCCAGTTTCCGCCAGCGACTTTGCCCTGCAATTTACGCCCGGAAAACGCTATCGGGTCATTGACGTCATCCATAACGAAGTGATTACGCACTCCCACTCCAGCGTCTACAGCg                                                                                                                               >  1:2542332/1‑150 (MQ=255)
                       tccgctccgcCATATGGCAACACCATTGCCCGCCAGCCAGTTTCCGCCAGCGACTTTGCCCTGCAATTTACGCCCGGAAAACGCTATCGGGTCATTGACGTCATCCATAACGAATTGATTACGCACTCCCACTCCAGCGTCTACAGCGaa                                                                                                                             <  2:3590359/150‑1 (MQ=255)
                        ccgctccgcCATATGGCAACACCATTGCCCGCCAGCCAGTTTCCGCCAGCGACTTTGCCCTGCAATTTACGCCCGGAAAACGCTATCGGGTCATTGACGTCATCCATAACGAATTGATTACGCACTCCCACTCCAGCGTCTACAGCGaaa                                                                                                                            <  1:1250503/150‑1 (MQ=255)
                           ctccgcCATATGGCAACACCATTGCCCGCCAGCCAGTTTCCGCCAGCGACTTTGCCCTGCAATTTACGCCCGGAAAACGCTATCGGGTCATTGACGTCATCCATAACGAATTGATTACGCACTCCCACTCCAGCGTCTACAGCGAAAATg                                                                                                                         <  1:1362635/150‑1 (MQ=255)
                            tccgcCATATGGCAACACCATTGCCCGCCAGCCAGTTTCCGCCAGCGACTTTGCCCTGCAATTTACGCCCGGAAAACGCTATCGGGTCATTGACGTCATCCATAACGAATTGATTACGCACTCCCACTCCAGCGTCTACAGCGAAAATgg                                                                                                                        >  2:3910988/1‑150 (MQ=255)
                             ccgcCATATGGCAACACCATTGCCCGCCAGCCAGTTTCCGCCAGCGACTTTGCCCTGCAATTTACGCCCGGAAAACGCTATCGGGTCATTGACGTCATCCATAACGAATTGATTACGCACTCCCACTCCAGCGTCTACAGCGAAAATGGt                                                                                                                       <  2:1217505/150‑1 (MQ=255)
                             ccgcCATATGGCAACACCATTGCCCGCCAGCCAGTTTCCGCCAGCGACTTTGCCCTGCAATTTACGCCCGGAAAACGCTATCGGGTCATTGACGTCATCCATAACGAATTGATTACGCACTCCCACTCCAGCGTCTACAGCGAAAATGGt                                                                                                                       <  2:2577677/150‑1 (MQ=255)
                              cgcCATATGGCAACACCATTGCCCGCCAGCCAGTTTCCGCCAGCGACTTTGCCCTGCAATTTACGCCCGGAAAACGCTATCGGGTCATTGACGTCATCCATAACGAATTGATTACGCACTCCCACTCCAGCGTCTACAGCGAAAATGGtt                                                                                                                      >  2:1765426/1‑150 (MQ=255)
                              cgcCATATGGCAACACCATTGCCCGCCAGCCAGTTTCCGCCAGCGACTTTGCCCTGCAATTTACGCCCGGAAAACGCTATCGGGTCATTGACGTCATCCATAACGAATTGATTACGCACTCCCACTCCAGCGTCTACAGCGAAAATGGtt                                                                                                                      >  2:2210405/1‑150 (MQ=255)
                                ccATATGGCAACACCATTGCCCGCCAGCCAGTTTCCGCCAGCGACTTTGCCCTGCAATTTACGCCCGGAAAACGCTATCGGGTCATTGACGTCATCCATAACGAATTGATTACGCACTCCCACTCCAGCGTCTACAGCGAAAATGGtttt                                                                                                                    <  1:907926/150‑1 (MQ=255)
                                 cATATGGCAACACCATTGCCCGCCAGCCAGTTTCCGCCAGCGACTTTGCCCTGCAATTTACGCCCGGAAAACGCTATCGGGTCATTGACGTCATCCATAACGAATTGATTACGCACTCCCACTCCAGCGTCTACAGCGAAAATGGTTTTg                                                                                                                   <  2:744962/150‑1 (MQ=255)
                                 cATATGGCAACACCATTGCCCGCCAGCCAGTTTCCGCCAGCGACTTTGCCCTGCAATTTACGCCCGGAAAACGCTATCGGGTCATTGACGTCATCCATAACGAATTGATTACGCACTCCCACTCCAGCGTCTACAGCGAAAATGGTTTTg                                                                                                                   >  1:3396822/1‑150 (MQ=255)
                                  atatGGCAACACCATTGCCCGCCAGCCAGTTTCCGCCAGCGACTTTGCCCTGCAATTTACGCCCGGAAAACGCTATCGGGTCATTGACGTCATCCATAACGAATTGATTACGCACTCCCACTCCAGCGTCTACAGCGAAAATGGTTTTGa                                                                                                                  <  1:198136/150‑1 (MQ=255)
                                  atatGGCAACACCATTGCCCGCCAGCCAGTTTCCGCCAGCGACTTTGCCCTGCAATTTACGCCCGGAAAACGCTATCGGGTCATTGACGTCATCCATAACGAATTGATTACGCACACCCACTCCAGCGTCTACAGCGAAAATGGTTTTGa                                                                                                                  <  1:1297867/150‑1 (MQ=255)
                                  atatGGCAACACCATTGCCCGCCAGCCAGTTTCCGCCAGCGACTTTGCCCTGCAATTTACGCCCGGAAAACGCTATCGGGTCATTGACGTCATCCATAACGAATTGATTACGCACACCCACTCCAGCGTCTACAGCGAAAATGGTTTTGa                                                                                                                  <  2:2819666/150‑1 (MQ=255)
                                   tatGGCAACACCATTGCCCGCCAGCCAGTTTCCGCCAGCGACTTTGCCCTGCAATTTACGCCCGGAAAACGCTATCGGGTCATTGACGTCATCCATAACGAATTGATTACGCACTCCCACTCCAGCGTCTACAGCGAAAATGGTTTTGAt                                                                                                                 >  1:2930814/1‑150 (MQ=255)
                                    atGGCAACACCATTGCCCGCCAGCCAGTTTCCGCCAGCGACTTTGCCCTGCAATTTACGCCCGGAAAACGCTATCGGGTCATTGACGTCATCCATAACGAATTGATTACGCACTCCCACTCCAGCGTCTACAGCGAAAATGGTTTTGATc                                                                                                                <  2:465735/150‑1 (MQ=255)
                                    atGGCAACACCATTGCCCGCCAGCCAGTTTCCGCCAGCGACTTTGCCCTGCAATTTACGCCCGGAAAACGCTATCGGGTCATTGACGTCATCCATAACGAATTGATTACGCACTCCCACTCCAGCGTCTACAGCGAAAATGGTTTTGATc                                                                                                                >  1:729427/1‑150 (MQ=255)
                                     tGGCAACACCATTGCCCGCCAGCCAGTTTCCGCCAGCGACTTTGCCCTGCAATTTACGCCCGGAAAACGCTATCGGGTCATTGACGTCATCCATAACGAATTGATTACGCACTCCCACTCCAGCGTCTACAGCGAAAATGGTTTTGATcg                                                                                                               <  2:2220011/150‑1 (MQ=255)
                                     tGGCAACACCATTGCCCGCCAGCCAGTTTCCGCCAGCGACTTTGCCCTGCAATTTACGCCCGGAAAACGCTATCGGGTCATTGACGTCATCCATAACGAATTGATTACGCACTCCCACTCCAGCGTCTACAGCGAAAATGGTTTTGATcg                                                                                                               <  1:746778/150‑1 (MQ=255)
                                      ggCAACACCATTGCCCGCCAGCCAGTTTCCGCCAGCGACTTTGCCCTGCAATTTACGCCCGGAAAACGCTATCGGGTCATTGACGTCATCCATAACGAATTGATTACGCACTCCCACTCCAGCGTCTACAGCGAAAATGGTTTTGATcgc                                                                                                              >  1:3212249/1‑150 (MQ=255)
                                      ggCAACACCATTGCCCGCCAGCCAGTTTCCGCCAGCGACTTTGCCCTGCAATTTACGCCCGGAAAACGCTATCGGGTCATTGACGTCATCCATAACGAATTGATTACGCACACCCACTCCAGCGTCTACAGCGAAAATGGGTTTGATcgc                                                                                                              >  2:1581780/1‑150 (MQ=255)
                                       gCAACACCATTGCCCGCCAGCCAGTTTCCGCCAGCGACTTTGCCCTGCAATTTACGCCCGGAAAACGCTATCGGGTCATTGACGTCATCCATAACGAATTGATTACGCACTCCCACTCCAGCGTCTACAGCGAAAATGGTTTTGATcgcg                                                                                                             <  1:2717990/150‑1 (MQ=255)
                                          acacCATTGCCCGCCAGCCAGTTTCCGCCAGCGACTTTGCCCTGCAATTTACGCCCGGAAAACGCTATCGGGTCATTGACGTCATCCATAACGAATTGATTACGCACTCCCACTCCAGCGTCTACAGCGAAAATGGTTTTGATCGCgatg                                                                                                          <  2:3517124/150‑1 (MQ=255)
                                          acacCATTGCCCGCCAGCCAGTTTCCGCCAGCGACTTTGCCCTGCAATTTACGCCCGGAAAACGCTATCGGGTCATTGACGTCATCCATAACGAATTGATTACGCACTCCCACTCCAGCGTCTACAGCGAAAATGGTTTTGATCGCgatg                                                                                                          <  2:1896480/150‑1 (MQ=255)
                                           cacCATTGCCCGCCAGCCAGTTTCCGCCAGCGACTTTGCCCTGCAATTTACGCCCGGAAAACGCTATCGGGTCATTGACGTCATCCATAACGAATTGATTACGCACTCCCACTCCAGCGTCTACAGCGAAAATGGTTTTGATCGCgatga                                                                                                         >  1:2595686/1‑150 (MQ=255)
                                             ccATTGCCCGCCAGCCAGTTTCCGCCAGCGACTTTGCCCTGCAATTTACGCCCGGAAAACGCTATCGGGTCATTGACGTCATCCATAACGAATTGATTACGCACTCCCACTCCAGCGTCTACAGCGAAAATGGTTTTGATCGCGatgatg                                                                                                       <  2:1844796/150‑1 (MQ=255)
                                              cATTGCCCGCCAGCCAGTTTCCGCCAGCGACTTTGCCCTGCAATTTACGCCCGGAAAACGCTATCGGGTCATTGACGTCATCCATAACGAATTGATTACGCACTCCCACTCCAGCGTCTACAGCGAAAATGGTTTTGATCGCGATGAtgt                                                                                                      >  1:34623/1‑150 (MQ=255)
                                                   cccGCCAGCCAGTTTCCGCCAGCGACTTTGCCCTGCAATTTACGCCCGGAAAACGCTATCGGGTCATTGACGTCATCCATAACGAATTGATTACGCACTCCCACTCCAGCGTCTACAGCGAAAATGGTGTTGATCGCGATGATGTGAGCt                                                                                                 >  1:485266/1‑150 (MQ=255)
                                                   cccGCCAGCCAGTTTCCGCCAGCGACTTTGCCCTGCAATTTACGCCCGGAAAACGCTATCGGGTCATTGACGTCATCCATAACGAATTGATTACGCACACCCACTCCAGCGTCTACAGCGAAAATGGTTTTGATCGCGATGATGTGAGCt                                                                                                 <  1:1579472/150‑1 (MQ=255)
                                                     cGCCAGCCAGTTTCCGCCAGCGACTTTGCCCTGCAATTTACGCCCGGAAAACGCTATCGGGTCATTGACGTCATCCATAACGAATTGATTACGCACTCCCACTCCAGCGTCTACAGCGAAAATGGTTTTGATCGCGATGATGTGAGCttt                                                                                               <  2:2198917/150‑1 (MQ=255)
                                                     cGCCAGCCAGTTTCCGCCAGCGACTTTGCCCTGCAATTTACGCCCGGAAAACGCTATCGGGTCATTGACGTCATCCATAACGAATTGATTACGCACTCCCACTCCAGCGTCTACAGCGAAAATGGTTTTGATCGCGATGATGTGAGCttt                                                                                               <  2:3832917/150‑1 (MQ=255)
                                                     cGCCAGCCAGTTTCCGCCAGCGACTTTGCCCTGCAATTTACGCCCGGAAAACGCTATCGGGTCATTGACGTCATCCATAACGAATTGATTACGCACACCCACTCCAGCGTCTACAGCGAAAATGGTTTTGATCGCGATGATGTGAGCttt                                                                                               >  1:1044652/1‑150 (MQ=255)
                                                       ccagccagTTTCCGCCAGCGACTTTGCCCTGCAATTTACGCCCGGAAAACGCTATCGGGTCATTGACGTCATCCATAACGAATTGATTACGCACTCCCACTCCAGCGTCTACAGCGAAAATGGTTTTGATCGCGATGATGTGAGCTTTAt                                                                                             <  1:2865169/150‑1 (MQ=255)
                                                         agccagTTTCCGCCAGCGACTTTGCCCTGCAATTTACGCCCGGAAAACGCTATCGGGTCATTGACGTCATCCATAACGAATTGATTACGCACTCCCACTCCAGCGTCTACAGCGAAAATGGTTTTGATCGCGATGATGTGAGCTTTATTg                                                                                           >  1:1933779/1‑150 (MQ=255)
                                                            cagTTTCCGCCAGCGACTTTGCCCTGCAATTTACGCCCGGAAAACGCTATCGGGTCATTGACGTCATCCATAACGAATTGATTACGCACTCCCACTCCAGCGTCTACAGCGAAAATGGTTTTGATCGCGATGATGTGAGCTTTATTGCCg                                                                                        <  1:3291449/150‑1 (MQ=255)
                                                            cagTTTCCGCCAGCGACTTTGCCCTGCAATTTACGCCCGGAAAACGCTATCGGGTCATTGACGTCATCCATAACGAATTGATTACGCACTCCCACTCCAGCGTCTACAGCGAAAATGGTTTTGATCGCGATGATGTGAGCTTTATTGCCg                                                                                        >  1:1257957/1‑150 (MQ=255)
                                                            cagTTTCCGCCAGCGACTTTGCCCTGCAATTTACGCCCGGAAAACGCTATCGGGTCATTGACGTCATCCATAACGAATTGATTACGCACTCCCACTCCAGCGTCTACAGCGAAAATGGTTTTGATCGCGATGATGTGAGCTTTATTGCCg                                                                                        >  2:2037527/1‑150 (MQ=255)
                                                             agTTTCCGCCAGCGACTTTGCCCTGCAATTTACGCCCGGAAAACGCTATCGGGTCATTGACGTCATCCATAACGAATTGATTACGCACACCCACTCCAGCGTCTACAGCGAAAATGGTTTTGATCGCGATGATGTGAGCTTTATTGCCGt                                                                                       <  1:2112778/150‑1 (MQ=255)
                                                               tttCCGCCAGCGACTTTGCCCTGCAATTTACGCCCGGAAAACGCTATCGGGTCATTGACGTCATCCATAACGAATTGATTACGCACTCCCACTCCAGCGTCTACAGCGAAAATGGTTTTGATCGCGATGATGTGAGCTTTATTGCCGTAc                                                                                     >  2:4360446/1‑150 (MQ=255)
                                                                ttCCGCCAGCGACTTTGCCCTGCAATTTACGCCCGGAAAACGCTATCGGGTCATTGACGTCATCCATAACGAATTGATTACGCACACCCACTCCAGCGTCTACAGCGAAAATGGTTTTGATCGCGATGATGTGAGCTTTATTGCCGTACt                                                                                    <  1:502771/150‑1 (MQ=255)
                                                                 tCCGCCAGCGACTTTGCCCTGCAATTTACGCCCGGAAAACGCTATCGGGTCATTGACGTCATCCATAACGAATTGATTACGCACTCCCACTCCAGCGTCTACAGCGAAAATGGTTTTGATCGCGATGATGTGAGCTTTATTGCCGTACtt                                                                                   <  1:1708675/150‑1 (MQ=255)
                                                                 tCCGCCAGCGACTTTGCCCTGCAATTTACGCCCGGAAAACGCTATCGGGTCATTGACGTCATCCATAACGAATTGATTACGCACTCCCACTCCAGCGTCTACAGCGAAAATGGTTTTGATCGCGATGATGTGAGCTTTATTGCCGTACtt                                                                                   <  2:2791974/150‑1 (MQ=255)
                                                                   cGCCAGCGACTTTGCCCTGCAATTTACGCCCGGAAAACGCTATCGGGTCATTGACGTCATCCATAACGAATTGATTACGCACTCCCACTCCAGCGTCTACAGCGAAAATGGTTTTGATCGCGATGATGTGAGCTTTATTGCCGTACTTGa                                                                                 <  2:681528/150‑1 (MQ=255)
                                                                     ccAGCGACTTTGCCCTGCAATTTACGCCCGGAAAACGCTATCGGGTCATTGACGTCATCCATAACGAATTGATTACGCACTCCCACTCCAGCGTCTACAGCGAAAATGGTTTTGATCGCGATGATGTGAGCTTTATTGCCGTACTTGAGc                                                                               <  2:4295695/150‑1 (MQ=255)
                                                                     ccAGCGACTTTGCCCTGCAATTTACGCCCGGAAAACGCTATCGGGTCATTGACGTCATCCATAACGAATTGATTACGCACTCCCACTCCAGCGTCTACAGCGAAAATGGTTTTGATCGCGATGATGTGAGCTTTATTGCCGTACTTGAGc                                                                               <  1:2357768/150‑1 (MQ=255)
                                                                      cAGCGACTTTGCCCTGCAATTTACGCCCGGAAAACGCTATCGGGTCATTGACGTCATCCATAACGAATTGATTACGCACTCCCACTCCAGCGTCTACAGCGAAAATGGTTTTGATCGCGATGATGTGAGCTTTATTGCCGTACTTGAGCg                                                                              <  1:3878397/150‑1 (MQ=255)
                                                                          gACTTTGCCCTGCAATTTACGCCCGGAAAACGCTATCGGGTCATTGACGTCATCCATAACGAATTGATTACGCACTCCCACTCCAGCGTCTACAGCGAAAATGGTTTTGATCGCGATGATGTGAGCTTTATTGCCGTACTTGAGCGTTAc                                                                          >  1:2737976/1‑150 (MQ=255)
                                                                          gACTTTGCCCTGCAATTTACGCCCGGAAAACGCTATCGGGTCATTGACGTCATCCATAACGAATTGATTACGCACTCCCACTCCAGCGTCTACAGCGAAAATGGTTTTGATCGCGATGATGTGAGCTTTATTGCCGTACTTGAGCGTTAc                                                                          >  1:761599/1‑150 (MQ=255)
                                                                           aCTTTGCCCTGCAATTTACGCCCGGAAAACGCTATCGGGTCATTGACGTCATCCATAACGAATTGATTACGCACTCCCACTCCAGCGTCTACAGCGAAAATGGTTTTGATCGCGATGATGTGAGCTTTATTGCCGTACTTGAGCGTTACg                                                                         <  1:4182939/150‑1 (MQ=255)
                                                                           aCTTTGCCCTGCAATTTACGCCCGGAAAACGCTATCGGGTCATTGACGTCATCCATAACGAATTGATTACGCACTCCCACTCCAGCGTCTACAGCGAAAATGGTTTTGATCGCGATGATGTGAGCTTTATTGCCGTACTTGAGCGTTACg                                                                         <  2:2415871/150‑1 (MQ=255)
                                                                            cTTTGCCCTGCAATTTACGCCCGGAAAACGCTATCGGGTCATTGACGTCATCCATAACGAATTGATTACGCACTCCCACTCCAGCGTCTACAGCGAAAATGGTTTTGATCGCGATGATGTGAGCTTTATTGCCGTACTTGAGCGTTACgg                                                                        <  2:1110687/150‑1 (MQ=255)
                                                                                gCCCTGCAATTTACGCCCGGAAAACGCTATCGGGTCATTGACGTCATCCATAACGAATTGATTACGCACTCCCACTCCAGCGTCTACAGCGAACATGGTTTTGATCGCGATGATGTGAGCTTTATTGCCGTACTTGAGCGTTACGGGCaa                                                                    >  1:4541411/1‑150 (MQ=255)
                                                                                 cccTGCAATTTACGCCCGGAAAACGCTATCGGGTCATTGACGTCATCCATAACGAATTGATTACGCACTCCCACTCCAGCGTCTACAGCGAAAATGGTTTTGATCGCGATGATGTGAGCTTTATTGCCGTACTTGAGCGTTACGGGCAAc                                                                   >  1:3994840/1‑150 (MQ=255)
                                                                                 cccTGCAATTTACGCCCGGAAAACGCTATCGGGTCATTGACGTCATCCATAACGAATTGATTACGCACTCCCACTCCAGCGTCTACAGCGAAAATGGTTTTGATCGCGATGATGTGAGCTTTATTGCCGTACTTGAGCGTTACGGGCAAc                                                                   >  2:2897815/1‑150 (MQ=255)
                                                                                      cAATTTACGCCCGGAAAACGCTATCGGGTCATTGACGTCATCCATAACGAATTGATTACGCACTCCCACTCCAGCGTCTACAGCGAAAATGGTTTTGATCGCGATGATGTGAGCTTTATTGCCGTACTTGAGCGTTACGGGCAACggctg                                                              <  2:876044/150‑1 (MQ=255)
                                                                                            aCGCCCGGAAAACGCTATCGGGTCATTGACGTCATCCATAACGAATTGATTACGCACTCCCACTCCAGCGTCTACAGCGAAAATGGTTTTGATCGCGATGATGTGAGCTTTATTGCCGTACTTGAGCGTTACGGGCAACGGCTGGCTCCg                                                        <  1:3734891/150‑1 (MQ=255)
                                                                                            aCGCCCGGAAAACGCTATCGGGTCATTGACGTCATCCATAACGAATTGATTACGCACTCCCACTCCAGCGTCTACAGCGAAAATGGTTTTGATCGCGATGATGTGAGCTTTATTGCCGTACTTGAGCGTTACGGGCAACGGCTGGCTCCg                                                        <  2:3040052/150‑1 (MQ=255)
                                                                                            aCGCCCGGAAAACGCTATCGGGTCATTGACGTCATCCATAACGAATTGATTACGCACTCCCACTCCAGCGTCTACAGCGAAAATGGTTTTGATCGCGATGATGTGAGCTTTATTGCCGTACTTGAGCGTTACGGGCAACGGCTGGCTCCg                                                        <  1:2463853/150‑1 (MQ=255)
                                                                                             cGCCCGGAAAACGCTATCGGGTCATTGACGTCATCCATAACGAATTGATTACGCACTCCCACTCCAGCGTCTACAGCGAAAATGGTTTTGATCGCGATGATGTGAGCTTTATTGCCGTACTTGAGCGTTACGGGCAACGGCTGGCTCCgg                                                       <  1:3152292/150‑1 (MQ=255)
                                                                                               cccGGAAAACGCTATCGGGTCATTGACGTCATCCATAACGAATTGATTACGCACTCCCACTCCAGCGTCTACAGCGAAAATGGTTTTGATCGCGATGATGTGAGCTTTATTGCCGTACTTGAGCGTTACGGGCAACGGCTGGCTCCGGCt                                                     >  2:4109588/1‑150 (MQ=255)
                                                                                                  ggAAAACGCTATCGGGTCATTGACGTCATCCATAACGAATTGATTACGCACACCCACTCCAGCGTCTACAGCGAAAATGGTTTTGATCGCGATGATGTGAGCTTTATTGCCGTACTTGAGCGTTACGGGCAACGGCTGGCTCCGGCTtgt                                                  <  2:3327859/150‑1 (MQ=255)
                                                                                                      aaCGCTATCGGGTCATTGACGTCATCCATAACGAATTGATTACGCACTCCCACTCCAGCGTCTACAGCGAAAATGGTTTTGATCGCGATGATGTGAGCTTTATTGCCGTACTTGAGCGTTACGGGCAACGGCTGGCTCCGGCTTGTGGtt                                              >  1:1244981/1‑150 (MQ=255)
                                                                                                       aCGCTATCGGGTCATTGACGTCATCCATAACGAATTGATTACGCACTCCCACTCCAGCGTCTACAGCGAAAATGGTTTTGATCGCGATGATGTGAGCTTTATTGCCGTACTTGAGCGTTACGGGCAACGGCTGGCTCCGGCTTGTGGttt                                             <  2:4187218/150‑1 (MQ=255)
                                                                                                                   cATTGACGTCATCCATAACGAATTGATTACGCACTCCCACTCCAGCGTCTACAGCGAAAATGGTTTTGATCGCGATGATGTGAGCTTTATTGCCGTACTTGAGCGTTACGGGCAACGGCTGGCTCCGGCTTGTGGTTTGCTTGGCGGCtt                                 <  2:2074605/150‑1 (MQ=255)
                                                                                                                    aTTGACGTCATCCATAACGAATTGATTACTCACTCCCACTCCAGCGTCAACAGCGAAAATGGTTAAGATCGCGAAGATGTGAGCTTTATTGCCGTACTAGAGCGTTACGGGCAACGGCTGGCTACGGCTTGTGGTTTGCTTGGCGGCttt                                >  2:811686/1‑150 (MQ=255)
                                                                                                                    aTTGACGTCATCCATAACGAATTGATTACGCACTCCCACTCCAGCGTCTACAGCGAAAATGGTTTTGATCGCGATGATGTGAGCTTTATTGCCGTACTTGAGCGTTACGGGCAACGGCTGGCTCCGGCTTGTGGTTTGCTTGGCGGCttt                                <  1:2341267/150‑1 (MQ=255)
                                                                                                                    aTTGACGTCATCCATAACGAATTGATTACGCACTCCCACTCCAGCGTCTACAGCGAAAATGGTTTTGATCGCGATGATGTGAGCTTTATTGCCGTACTTGAGCGTTACGGGCAACGGCTGGCTCCGGCTTGTGGTTTGCTTGGCGGCttt                                >  2:1528732/1‑150 (MQ=255)
                                                                                                                    aTTGACGTCATCCATAACGAATTGATTACGCACTCCCACTCCAGCGTCTACAGCGAAAATGGTTTTGATCGCGATGATGTGAGCTTTATTGCCGTACTTGAGCGTTACGGGCAACGGCTGGCTCCGGCTTGTGGTTTGCTTGGCGGCttt                                <  1:619591/150‑1 (MQ=255)
                                                                                                                     ttGACGTCATCCATAACGAATTGATTACGCACTCCCACTCCAGCGTCTACAGCGAAAATGGTTTTGATCGCGATGATGTGAGCTTTATTGCCGTACTTGAGCGTTACGGGCAACGGCTGGCTCCGGCTTGTGGTTTGCTTGGCGGCTTTg                               >  2:1012423/1‑150 (MQ=255)
                                                                                                                          gTCATCCATAACGAATTGATTACGCACTCCCACTCCAGCGTCTACAGCGAAAATGGTTTTGATCGCGATGATGTGAGCTTTATTGCCGTACTTGAGCGTTACGGGCAACGGCTGGCTCCGGCTTGTGGTTTGCTTGGCGGCTTTGGACtg                          >  2:4439561/1‑150 (MQ=255)
                                                                                                                          gTCATCCATAACGAATTGATTACGCACTCCCACTCCAGCGTCTACAGCGAAAATGGTTTTGATCGCGATGATGTGAGCTTTATTGCCGTACTTGAGCGTTACGGGCAACGGCTGGCTCCGGCTTGTGGTTTGCTTGGCGGCTTTGGACtg                          >  2:634957/1‑150 (MQ=255)
                                                                                                                                 aTAACGAATTGATTACGCACTCCCACTCCAGCGTCTACAGCGAAAATGGTTTTGATCGCGATGATGTGAGCTTTATTGCCGTACTTGAGCGTTACGGGCAACGGCTGGCTCCGGCTTGTGGTTTGCTTGGCGGCTTTGGACTGAATGAAg                   <  1:1900013/150‑1 (MQ=255)
                                                                                                                                 aTAACGAATTGATTACGCACTCCCACTCCAGCGTCTACAGCGAAAATGGTTTTGATCGCGATGATGTGAGCTTTATTGCCGTACTTGAGCGTTACGGGCAACGGCTGGCTCCGGCTTGTGGTTTGCTTGGCGGCTTTGGACTGAATGAAg                   <  1:2839112/150‑1 (MQ=255)
                                                                                                                                 aTAACGAATTGATTACGCACACCCACTCCAGCGTCTACAGCGAAAATGGTTTTGATCGCGATGATGTGAGCTTTATTGCCGTACTTGAGCGTTACGGGCAACGGCTGGCTCCGGCTTGTGGTTTGCTTGGCGGCTTTGGACTGAATGAAg                   >  1:4332090/1‑150 (MQ=255)
                                                                                                                                    aCGAATTGATTACGCACTCCCACTCCAGCGTCTACAGCGAAAATGGTTTTGATCGCGATGATGTGAGCTTTATTGCCGTACTTGAGCGTTACGGGCAACGGCTGGCTCCGGCTTGTGGTTTGCTTGGCGGCTTTGGACTGAATGAAGGTg                <  2:2910507/150‑1 (MQ=255)
                                                                                                                                      gAATTGATTACGCACTCCCACTCCAGCGTCTACAGCGAAAATGGTTTTGATCGCGATGATGTGAGCTTTATTGCCGTACTTGAGCGTTACGGGCAACGGCTGGCTCCGGCTTGTGGTTTGCTTGGCGGCTTTGGACTGAATGAAGGTgcg              <  1:3062377/150‑1 (MQ=255)
                                                                                                                                        aTTGATTACGCACTCCCACTCCAGCGTCTACAGCGAAAATGGTTTTGATCGCGATGATGTGAGCGTTATTGCCGTACTTGAGCGTTACGGGCAACGGCTGGCTCCGGCTTGTGGTTTGCTTGGCGGCTTTGGACTGAATGAAGGTGCgct            >  2:4363164/1‑150 (MQ=255)
                                                                                                                                         ttGATTACGCACTCCCACTCCAGCGTCTACAGCGAAAATGGTTTTGATCGCGATGATGTGAGCTTTATTGCCGTACTTGAGCGTTACGGGCAACGGCTGGCTCCGGCTTGTGGTTTGCTTGGCGGCTTTGGACTGAATGAAGGTGCgctg           <  2:3670362/150‑1 (MQ=255)
                                                                                                                                                cGCACTCCCACTCCAGCGTCTACAGCGAAAATGGTTTTGATCGCGATGATGTGAGCTTTATTGCCGTACTTGAGCGTTACGGGCAACGGCTGGCTCCGGCTTGTGGTTTGCTTGGCGGCTTTGGACTGAATGAAGGTGCGCTGGCTGCGa    <  2:1864478/150‑1 (MQ=255)
                                                                                                                                                cGCACACCCACTCCAGCGTCTACAGCGAAAATGGTTTTGATCGCGATGATGTGAGCTTTATTGCCGTACTTGAGCGTTACGGGCAACGGCTGGCTCCGGCTTGTGGTTTGCTTGGCGGCTTTGGACTGAATGAAGGTGCGCTGGCTGCGa    <  1:1636210/150‑1 (MQ=255)
                                                                                                                                                 gCACTCCCACTCCAGCGTCTACAGCGAAAATGGTTTTGATCGCGATGATGTGAGCTTTATTGCCGTACTTGAGCGTTACGGGCAACGGCTGGCTCCGGCTTGTGGTTTGCTTGGCGGCTTTGGACTGAATGAAGGTGCGCTGGCTGCGAc   <  1:2206825/150‑1 (MQ=255)
                                                                                                                                                  cactcccactccAGCGTCTACAGCGAAAATGGTTTTGATCGCGATGATGTGAGCTTTATTGCCGTACTTGAGCGTTACGGGCAACGGCTGGCTCCGGCTTGTGGTTTGCTTGGCGGCTTTGGACTGAATGAAGGTGCGCTGGCTGCGACg  <  2:1262880/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
AAGAGTCGGCGAGACTGGCACAATCCGCTCCGCCATATGGCAACACCATTGCCCGCCAGCCAGTTTCCGCCAGCGACTTTGCCCTGCAATTTACGCCCGGAAAACGCTATCGGGTCATTGACGTCATCCATAACGAATTGATTACGCACTCCCACTCCAGCGTCTACAGCGAAAATGGTTTTGATCGCGATGATGTGAGCTTTATTGCCGTACTTGAGCGTTACGGGCAACGGCTGGCTCCGGCTTGTGGTTTGCTTGGCGGCTTTGGACTGAATGAAGGTGCGCTGGCTGCGACG  >  CP009273/3838438‑3838733

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: