Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,838,587 T→A 19.4% 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,5870TA19.4% 726.6 / 164.5 293S422T (TCC→ACC) adecryptic adenine deaminase
Reads supporting (aligned to +/- strand):  ref base T (95/141);  new base A (30/27);  total (125/168)
Fisher's exact test for biased strand distribution p-value = 1.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 309 total aligned reads displayed.

AAGAGTCGGCGAGACTGGCACAATCCGCTCCGCCATATGGCAACACCATTGCCCGCCAGCCAGTTTCCGCCAGCGACTTTGCCCTGCAATTTACGCCCGGAAAACGCTATCGGGTCATTGACGTCATCCATAACGAATTGATTACGCACTCCCACTCCAGCGTCTACAGCGAAAATGGTTTTGATCGCGATGATGTGAGCTTTATTGCCGTACTTGAGCGTTACGGGCAACGGCTGGCTCCGGCTTGTGGTTTGCTTGGCGGCTTTGGACTGAATGAAGGTGCGCTGGCTGCGACGGTC  >  CP009273/3838438‑3838736
                                                                                                                                                     |                                                                                                                                                     
 agagTCGGCGAGACTGGCACAATCCGCTCCGCCATATGGCAACACCATTGCCCGCCAGCCAGTTTCCGCCAGCGACTTTGCCCTGCAATTTACGCCCGGAAAACGCTATCGGGTCATTGACGTCATCCATAACGAATTGATTACGcactc                                                                                                                                                      >  1:2849766/1‑150 (MQ=255)
       ggCGAGACTGGCACAATCCGCTCCGCCATATGGCAACACCATTGCCCGCCAGCCAGTTTCCGCCAGCGACTTTGCCCTGCAATTTACGCCCGGAAAACGCTATCGGGTCATTGACGTCATCCATAACGAATTGATTACGcactcccactc                                                                                                                                                <  2:3130313/150‑1 (MQ=255)
         cGAGACTGGCACAATCCGCTCCGCCATATGGCAACACCATTGCCCGCCAGCCAGTTTCCGCCAGCGACTTTGCCCTGCAATTTACGCCCGGAAAACGCTATCGGGTCATTGACGTCATCCATAACGAATTGATTACGCACTCCCACTCCa                                                                                                                                              >  1:4626529/1‑150 (MQ=255)
         cGAGACTGGCACAATCCGCTCCGCCATATGGCAACACCATTGCCCGCCAGCCAGTTTCCGCCAGCGACTTTGCCCTGCAATTTACGCCCGGAAAACGCTATCGGGTCATTGACGTCATCCATAACGAATTGATTACGCACACCCACTCCa                                                                                                                                              >  2:3498043/1‑150 (MQ=255)
         cGAGACTGGCACAATCCGCTCCGCCATATGGCAACACCATTGCCCGCCAGCCAGTTTCCGCCAGCGACTTTGCCCTGCAATTTACGCCCGGAAAACGCTATCGGGTCATTGACGTCATCCATAACGAATTGATTACGCACACCCACTCCa                                                                                                                                              <  1:3509527/150‑1 (MQ=255)
          gagaCTGGCACAATCCGCTCCGCCATATGGCAACACCATTGCCCGCCAGCCAGTTTCCGCCAGCGACTTTGCCCTGCAATTTACGCCCGGAAAACGCTATCGGGTCATTGACGTCATCCATAACGAATTGATTACGCACTCCCACTCCAg                                                                                                                                             <  2:2564332/150‑1 (MQ=255)
           agaCTGGCACAATCCGCTCCGCCATATGGCAACACCATTGCCCGCCAGCCAGTTTCCGCCAGCGACTTTGCCCTGCAATTTACGCCCGGAAAACGCTATCGGGTCATTGACGTCATCCATAACGAATTGATTACGCACTCCCACTCCAGc                                                                                                                                            <  2:821244/150‑1 (MQ=255)
           agaCTGGCACAATCCGCTCCGCCATATGGCAACACCATTGCCCGCCAGCCAGTTTCCGCCAGCGACTTTGCCCTGCAATTTACGCCCGGAAAACGCTATCGGGTCATTGACGTCATCCATAACGAATTGATTACGCACTCCCACTCCAGc                                                                                                                                            >  1:4340549/1‑150 (MQ=255)
            gaCTGGCACAATCCGCTCCGCCATATGGCAACACCATTGCCCGCCAGCCAGTTTCCGCCAGCGACTTTGCCCTGCAATTTACGCCCGGAAAACGCTATCGGGTCATTGACGTCATCCATAACGAATTGATTACGCACTCCCACTCCAGCg                                                                                                                                           <  1:3339621/150‑1 (MQ=255)
            gaCTGGCACAATCCGCTCCGCCATATGGCAACACCATTGCCCGCCAGCCAGTTTCCGCCAGCGACTTTGCCCTGCAATTTACGCCCGGAAAACGCTATCGGGTCATTGACGTCATCCATAACGAATTGATTACGCACTCCCACTCCAGCg                                                                                                                                           <  1:1485374/150‑1 (MQ=255)
            gaCTGGCACAATCCGCTCCGCCATATGGCAACACCATTGCCCGCCAGCCAGTTTCCGCCAGCGACTTTGCCCTGCAATTTACGCCCGGAAAACGCTATCGGGTCATTGACGTCATCCATAACGAATTGATTACGCACTCCCACTCCAGCg                                                                                                                                           <  1:1733567/150‑1 (MQ=255)
              cTGGCACAATCCGCTCCGCCATATGGCAACACCATTGCCCGCCAGCCAGTTTCCGCCAGCGACTTTGCCCTGCAATTTACGCCCGGAAAACGCTATCGGGTCATTGACGTCATCCATAACGAATTGATTACGCACTCCCACTCCAGCGTc                                                                                                                                         >  2:3008772/1‑150 (MQ=255)
               tGGCACAATCCGCTCCGCCATATGGCAACACCATTGCCCGCCAGCCAGTTTCCGCCAGCGACTTTGCCCTGCAATTTACGCCCGGAAAACGCTATCGGGTCATTGACGTCATCCATAACGAATTGATTACGCACTCCCACTCCAGCGTCt                                                                                                                                        >  1:3949082/1‑150 (MQ=255)
                ggCACAATCCGCTCCGCCATATGGCAACACCATTGCCCGCCAGCCAGTTTCCGCCAGCGACTTTGCCCTGCAATTTACGCCCGGAAAACGCTATCGGGTCATTGACGTCATCCATAACGAATTGATTACGCACTCCCACTCCAGCGTCTa                                                                                                                                       <  2:2198832/150‑1 (MQ=255)
                 gCACAATCCGCTCCGCCATATGGCAACACCATTGCCCGCCAGCCAGTTTCCGCCAGCGACTTTGCCCTGCAATTTACGCCCGGAAAACGCTATCGGGTCATTGACGTCATCCATAACGAATTGATTACGCACTCCCACTCCAGCGTCTAc                                                                                                                                      <  1:1548181/150‑1 (MQ=255)
                   acaATCCGCTCCGCCATATGGCAACACCATTGCCCGCCAGCCAGTTTCCGCCAGCGACTTTGCCCTGCAATTTACGCCCGGAAAACGCTATCGGGTCATTGACGTCATCCATAACGAATTGATTACGCACTCCCACTCCAGCGTCTACAg                                                                                                                                    <  2:1199839/150‑1 (MQ=255)
                      aTCCGCTCCGCCATATGGCAACACCATTGCCCGCCAGCCAGTTTCCGCCAGCGACTTTGCCCTGCAATTTACGCCCGGAAAACGCTATCGGGTCATTGACGTCATCCATAACGAATTGATTACGCACTCCCACTCCAGCGTCTACAGCGa                                                                                                                                 >  2:2074820/1‑150 (MQ=255)
                       tccgctccgcCATATGGCAACACCATTGCCCGCCAGCCAGTTTCCGCCAGCGACTTTGCCCTGCAATTTACGCCCGGAAAACGCTATCGGGTCATTGACGTCATCCATAACGAATTGATTACGCACTCCCACTCCAGCGTCTACAGCGaa                                                                                                                                >  1:3653213/1‑150 (MQ=255)
                         cgctccgcCATATGGCAACACCATTGCCCGCCAGCCAGTTTCCGCCAGCGACTTTGCCCTGCAATTTACGCCCGGAAAACGCTATCGGGTCATTGACGTCATCCATAACGAATTGATTACGCACTCCCACTCCAGCGTCTACAGCGaaaa                                                                                                                              <  2:1210341/150‑1 (MQ=255)
                          gctccgcCATATGGCAACACCATTGCCCGCCAGCCAGTTTCCGCCAGCGACTTTGCCCTGCAATTTACGCCCGGAAAACGCTATCGGGTCATTGACGTCATCCATAACGAATTGATTACGCACTCCCACTCCAGCGTCTACAGCGAAAAt                                                                                                                             <  1:3126449/150‑1 (MQ=255)
                           ctccgcCATATGGCAACACCATTGCCCGCCAGCCAGTTTCCGCCAGCGACTTTGCCCTGCAATTTACGCCCGGAAAACGCTATCGGGTCATTGACGTCATCCATAACGAATTGATTACGCACTCCCACTCCAGCGTCTACAGCGAAAATg                                                                                                                            <  1:1254902/150‑1 (MQ=255)
                           ctccgcCATATGGCAACACCATTGCCCGCCAGCCAGTTTCCGCCAGCGACTTTGCCCTGCAATTTACGCCCGGAAAACGCTATCGGGTCATTGACGTCATCCATAACGAATTGATTACGCACACCCACTCCAGCGTCTACAGCGAACATg                                                                                                                            >  2:522763/1‑150 (MQ=255)
                           ctccgcCATATGGCAACACCATTGCCCGCCAGCCAGTTTCCGCCAGCGACTTTGCCCTGCAATTTACGCCCGGAAAACGCTATCGGGTCATTGACGTCATCCATAACGAATTGATTACGCACACCCACTCCAGCGTCTACAGCGAAAATg                                                                                                                            <  1:2668758/150‑1 (MQ=255)
                            tccgcCATATGGCAACACCATTGCCCGCCAGCCAGTTTCCGCCAGCGACTTTGCCCTGCAATTTACGCCCGGAAAACGCTATCGGGTCATTGACGTCATCCATAACGAATTGATTACGCACTCCCACTCCAGCGTCTACAGCGAAAATgg                                                                                                                           >  1:1040739/1‑150 (MQ=255)
                            tccgcCATATGGCAACACCATTGCCCGCCAGCCAGTTTCCGCCAGCGACTTTGCCCTGCAATTTACGCCCGGAAAACGCTATCGGGTCATTGACGTCATCCATAACGAATTGATTACGCACTCCCACTCCAGCGTCTACAGCGAAAATgg                                                                                                                           <  1:2414994/150‑1 (MQ=255)
                            tccgcCATATGGCAACACCATTGCCCGCCAGCCAGTTTCCGCCAGCGACTTTGCCCTGCAATTTACGCCCGGAAAACGCTATCGGGTCATTGACGTCATCCATAACGAATTGATTACGCACTCCCACTCCAGCGTCTACAGCGAAAATgg                                                                                                                           <  2:583431/150‑1 (MQ=255)
                              cgcCATATGGCAACACCATTGCCCGCCAGCCAGTTTCCGCCAGCGACTTTGCCCTGCAATTTACGCCCGGAAAACGCTATCGGGTCATTGACGTCATCCATAACGAATTGATTACGCACTCCCACTCCAGCGTCTACAGCGAAAATGGtt                                                                                                                         <  1:3297716/150‑1 (MQ=255)
                                  atatGGCAACACCATTGCCCGCCAGCCAGTTTCCGCCAGCGACTTTGCCCTGCAATTTACGCCCGGAAAACGCTATCGGGTCATTGACGTCATCCATAACGAATTGATTACGCACTCCCACTCCAGCGTCTACAGCGAAAATGGTTTTGa                                                                                                                     >  1:1073076/1‑150 (MQ=255)
                                  atatGGCAACACCATTGCCCGCCAGCCAGTTTCCGCCAGCGACTTTGCCCTGCAATTTACGCCCGGAAAACGCTATCGGGTCATTGACGTCATCCATAACGAATTGATTACGCACACCCACTCCAGCGTCTACAGCGAAAATGGTTTTGa                                                                                                                     <  2:1879849/150‑1 (MQ=255)
                                   tatGGCAACACCATTGCCCGCCAGCCAGTTTCCGCCAGCGACTTTGCCCTGCAATTTACGCCCGGAAAACGCTATCGGGTCATTGACGTCATCCATAACGAATTGATTACGCACTCCCACTCCAGCGTCTACAGCGAAAATGGTTTTGAt                                                                                                                    <  1:223403/150‑1 (MQ=255)
                                   tatGGCAACACCATTGCCCGCCAGCCAGTTTCCGCCAGCGACTTTGCCCTGCAATTTACGCCCGGAAAACGCTATCGGGTCATTGACGTCATCCATAACGAATTGATTACGCACACCCACTCCAGCGTCTACAGCGAAAATGGTTTTGAt                                                                                                                    >  1:3584216/1‑150 (MQ=255)
                                    atGGCAACACCATTGCCCGCCAGCCAGTTTCCGCCAGCGACTTTGCCCTGCAATTTACGCCCGGAAAACGCTATCGGGTCATTGACGTCATCCATAACGAATTGATTACGCACTCCCACTCCAGCGTCTACAGCGAAAATGGTTTTGATc                                                                                                                   >  1:1129119/1‑150 (MQ=255)
                                     tGGCAACACCATTGCCCGCCAGCCAGTTTCCGCCAGCGACTTTGCCCTGCAATTTACGCCCGGAAAACGCTATCGGGTCATTGACGTCATCCATAACGAATTGATTACGCACTCCCACTCCAGCGTCTACAGCGAAAATGGTTTTGATcg                                                                                                                  <  2:3777662/150‑1 (MQ=255)
                                     tGGCAACACCATTGCCCGCCAGCCAGTTTCCGCCAGCGACTTTGCCCTGCAATTTACGCCCGGAAAACGCTATCGGGTCATTGACGTCATCCATAACGAATTGATTACGCACTCCCACTCCAGCGTCTACAGCGAAAATGGTTTTGATcg                                                                                                                  >  2:1539158/1‑150 (MQ=255)
                                        cAACACCATTGCCCGCCAGCCAGTTTCCGCCAGCGACTTTGCCCTGCAATTTACGCCCGGAAAACGCTATCGGGTCATTGACGTCATCCATAACGAATTGATTACGCACACCCACTCCAGCGTCTACAGCGAAAATGGTTTTGATCGCga                                                                                                               >  1:918104/1‑150 (MQ=255)
                                         aaCACCATTGCCCGCCAGCCAGTTTCCGCCAGCGACTTTGCCCTGCAATTTACGCCCGGAAAACGCTATCGGGTCATTGACGTCATCCATAACGAATTGATTACGCACTCCCACTCCAGCGTCTACAGCGAAAATGGTTTTGATCGCgat                                                                                                              <  2:49916/150‑1 (MQ=255)
                                         aaCACCATTGCCCGCCAGCCAGTTTCCGCCAGCGACTTTGCCCTGCAATTTACGCCCGGAAAACGCTATCGGGTCATTGACGTCATCCATAACGAATTGATTACGCACACCCACTCCAGCGTCTACAGCGAAAATGGTTTTGATCGCgat                                                                                                              >  2:2447967/1‑150 (MQ=255)
                                             ccATTGCCCGCCAGCCAGTGTCCGCCAGCGACTTTGCCCTGCAATTTACGCCCGGAAAACGCTATCGGGTCATTGACGTCATCCATAACGAATTGATTACGCACTCCCACTCCAGCGTCTACAGCGAAAATGGTTTTGATCGCGatgatg                                                                                                          <  1:569386/150‑1 (MQ=255)
                                              cATTGCCCGCCAGCCAGTTTCCGCCAGCGACTTTGCCCTGCAATTTACGCCCGGAAAACGCTATCGGGTCATTGACGTCATCCATAACGAATTGATTACGCACTCCCACTCCAGCGTCTACAGCGAAAATGGTTTTGATCGCGATGAtgt                                                                                                         >  2:4448190/1‑150 (MQ=255)
                                              cATTGCCCGCCAGCCAGTTTCCGCCAGCGACTTTGCCCTGCAATTTACGCCCGGAAAACGCTATCGGGTCATTGACGTCATCCATAACGAATTGATTACGCACTCCCACTCCAGCGTCTACAGCGAAAATGGTTTTGATCGCGATGAtgt                                                                                                         >  2:2005603/1‑150 (MQ=255)
                                                ttGCCCGCCAGCCAGTTTCCGCCAGCGACTTTGCCCTGCAATTTACGCCCGGAAAACGCTATCGGGTCATTGACGTCATCCATAACGAATTGATTACGCACTCCCACTCCAGCGTCTACAGCGAAAATGGTTTTGATCGCGATGATGTGa                                                                                                       >  1:1825295/1‑150 (MQ=255)
                                                ttGCCCGCCAGCCAGTTTCCGCCAGCGACTTTGCCCTGCAATTTACGCCCGGAAAACGCTATCGGGTCATTGACGTCATCCATAACGAATTGATTACGCACACCCACTCCAGCGTCTACAGCGAAAATGGTTTTGATCGCGATGATGTGa                                                                                                       >  2:325889/1‑150 (MQ=255)
                                                 tGCCCGCCAGCCAGTTTCCGCCAGCGACTTTGCCCTGCAATTTACGCCCGGAAAACGCTATCGGGTCATTGACGTCATCCATAACGAATTGATTACGCACTCCCACTCCAGCGTCTACAGCGAAAATGGTTTTGATCGCGATGATGTGAg                                                                                                      <  2:701714/150‑1 (MQ=255)
                                                   cccGCCAGCCAGTTTCCGCCAGCGACTTTGCCCTGCAATTTACGCCCGGAAAACGCTATCGGGTCATTGACGTCATCTATAACGAATTGATTACGCACTCCCACTCCAGCGTCTACAGCGAAAATGGTTTTGATCGCGATGATGTGAGCt                                                                                                    <  1:3635710/150‑1 (MQ=255)
                                                   cccGCCAGCCAGTTTCCGCCAGCGACTTTGCCCTGCAATTTACGCCCGGAAAACGCTATCGGGTCATTGACGTCATCCATAACGAATTGATTACGCACTCCCACTCCAGCGTCTACAGCGAAAATGGTTTTGATCGCGATGATGTGAGCt                                                                                                    <  2:1117371/150‑1 (MQ=255)
                                                   cccGCCAGCCAGTTTCCGCCAGCGACTTTGCCCTGCAATTTACGCCCGGAAAACGCTATCGGGTCATTGACGTCATCCATAACGAATTGATTACGCACTCCCACTCCAGCGTCTACAGCGAAAATGGTTTTGATCGCGATGATGTGAGCt                                                                                                    <  1:2435915/150‑1 (MQ=255)
                                                   cccGCCAGCCAGTTTCCGCCAGCGACTTTGCCCTGCAATTTACGCCCGGAAAACGCTATCGGGTCATTGACGTCATCCATAACGAATTGATTACGCACTCCCACTCCAGCGTCTACAGCGAAAATGGTTTTGATCGCGATGATGTGAGCt                                                                                                    <  1:130724/150‑1 (MQ=255)
                                                    ccGCCAGCCAGTTTCCGCCAGCGACTTTGCCCTGCAATTTACGCCCGGAAAACGCTATCGGGTCATTGACGTCATCCATAACGAATTGATTACGCACTCCCACTCCAGCGTCTACAGCGAAAATGGTTTTGATCGCGATGATGTGAGCtt                                                                                                   <  2:2293163/150‑1 (MQ=255)
                                                    ccGCCAGCCAGTTTCCGCCAGCGACTTTGCCCTGCAATTTACGCCCGGAAAACGCTATCGGGTCATTGACGTCATCCATAACGAATTGATTACGCACTCCCACTCCAGCGTCTACAGCGAAAATGGTTTTGATCGCGATGATGTGAGCtt                                                                                                   <  1:3081744/150‑1 (MQ=255)
                                                       ccagccagTTTCCGCCAGCGACTTTGCCCTGCAATTTACGCCCGGAAAACGCTATCGGGTCATTGACGTCATCCATAACGAATTGATTACGCACACCCACTCCAGCGTCTACAGCGAAAATGGTTGTGATCGCGAGGATGTGAGCTTTAt                                                                                                >  2:4189045/1‑150 (MQ=255)
                                                        cagccagTTTCCGCCAGCGACTTTGCCCTGCAATTTACGCCCGGAAAACGCTATCGGGTCATTGACGTCATCCATAACGAATTGATTACGCACTCCCACTCCAGCGTCTACAGCGAAAATGGTTTTGATCGCGATGATGTGAGCTTTAtt                                                                                               <  1:3857914/150‑1 (MQ=255)
                                                        cagccagTTTCCGCCAGCGACTTTGCCCTGCAATTTACGCCCGGAAAACGCTATCGGGTCATTGACGTCATCCATAACGAATTGATTACGCACTCCCACTCCAGCGTCTACAGCGAAAATGGTGTTGATCGCGATGATGTGAGCTTTAtt                                                                                               >  1:3225592/1‑150 (MQ=255)
                                                           ccagTTTCCGCCAGCGACTTTGCCCTGCAATTTACGCCCGGAAAACGCTATCGGGTCATTGACGTCATCCATAACGAATTGATTACGCACTCCCACTCCAGCGTCTACAGCGAAAATGGTTTTGATCGCGATGATGTGAGCTTTATTGcc                                                                                            <  2:1025072/150‑1 (MQ=255)
                                                            cagTTTCCGCCAGCGACTTTGCCCTGCAATTTACGCCCGGAAAACGCTATCGGGTCATTGACGTCATCCATAACGAATTGATTACGCACTCCCACTCCAGCGTCTACAGCGAAAATGGTTTTGATCGCGATGATGTGAGCTTTATTGCCg                                                                                           <  1:1495609/150‑1 (MQ=255)
                                                              gTTTCCGCCAGCGACTTTGCCCTGCAATTTACGCCCGGAAAACGCTATCGGGTCATTGACGTCATCCATAACGAATTGATTACGCACTCCCACTCCAGCGTCTACAGCGAAAATGGTTTTGATCGCGATGATGTGAGCTTTATTGCCGTa                                                                                         >  1:1941474/1‑150 (MQ=255)
                                                              gTTTCCGCCAGCGACTTTGCCCTGCAATTTACGCCCGGAAAACGCTATCGGGTCATTGACGTCATCCATAACGAATTGATTACGCACACCCACTCCAGCGTCTACAGCGAAAATGGTTTTGATCGCGATGATGTGAGCTTTATTGCCGTa                                                                                         <  2:3313191/150‑1 (MQ=255)
                                                               tttCCGCCAGCGACTTTGCCCTGCAATTTACGCCCGGAAAACGCTATCGGGTCATTGACGTCATCCATAACGAATTGATTACGCACACCCACTCCAGCGTCTACAGCGAACATGGTTTTGATCGCGATGATGTGAGCTTTATTGCCGTAc                                                                                        >  2:3912649/1‑150 (MQ=255)
                                                                ttCCGCCAGCGACTTTGCCCTGCAATTTACGCCCGGAAAACGCTATCGGGTCATTGACGTCATCCATAACGAATTGATTACGCACACCCACTCCAGCGTCTACAGCGAAAATGGTTTTGATCGCGATGATGTGAGCTTTATTGCCGTACt                                                                                       <  2:3016428/150‑1 (MQ=255)
                                                                  ccGCCAGCGACTTTGCCCTGCAATTTACGCCCGGAAAACGCTATCGGGTCATTGACGTCATCCATAACGAATTGATTACGCACTCCCACTCCAGCGTCTACAGCGAAAATGGTTTTGATCGCGATGATGTGAGCTTTATTGCCGTACTTg                                                                                     <  2:1385624/150‑1 (MQ=255)
                                                                       aGCGACTTTGCCCTGCAATTTACGCCCGGAAAACGCTATCGGGTCATTGACGTCATCCATAACGAATTGATTACGCACTCCCACTCCAGCGTCTACAGCGAAAATGGTTTTGATCGCGATGATGTGAGCTTTATTGCCGTACTTGAGCGt                                                                                >  1:4073198/1‑150 (MQ=255)
                                                                          gACTTTGCCCTGCAATTTACGCCCGGAAAACGCTATCGGGTCATTGACGTCATCCATAACGAATTGATTACGCACTCCCACTCCAGCGTCTACAGCGAAAATGGTTTTGATCGCGATGATGTGAGCTTTATTGCCGTACTTGAGCGTTAc                                                                             >  2:735008/1‑150 (MQ=255)
                                                                           aCTTTGCCCTGCAATTTACGCCCGGAAAACGCTATCGGGTCATTGACGTCATCCATAACGAATTGATTACGCACTCCCACTCCAGCGTCTACAGCGAAAATGGTTTTGATCGCGATGATGTGAGCTTTATTGCCGTACTTGAGCGTTACg                                                                            <  1:2498725/150‑1 (MQ=255)
                                                                                   cTGCAATTTACGCCCGGAAAACGCTATCGGGTCATTGACGTCATCCATAACGAATTGATTACGCACTCCCACTCCAGCGTCTACAGCGAAAATGGTTTTGATCGCGATGATGTGAGCTTTATTGCCGTACTTGAGCGTTACGGGCAACgg                                                                    >  1:2061736/1‑150 (MQ=255)
                                                                                      cAATTTACGCCCGGAAAACGCTATCGGGTCATTGACGTCATCCATAACGAATTGATTACGCACTCCCACTCCAGCGTCTACAGCGAAAATGGTTTTGATCGCGATGATGTGAGCTTTATTGCCGTACTTGAGCGTTACGGGCAACggctg                                                                 >  2:950589/1‑150 (MQ=255)
                                                                                       aaTTTACGCCCGGAAAACGCTATCGGGTCATTGACGTCATCCATAACGAATTGATTACGCACTCCCACTCCAGCGTCTACAGCGAAAATGGTTTTGATCGCGATGATGTGAGCTTTATTGCCGTACTTGAGCGTTACGGGCAACggctgg                                                                >  2:2354999/1‑150 (MQ=255)
                                                                                             cGCCCGGAAAACGCTATCGGGTCATTGACGTCATCCATAACGAAGTGATTACGCACTCCCACTCCAGCGTCTACAGCGAAAATGGTTTTGATCGCGATGATGTGAGCTTTATTGCCGTACTTGAGCGTTACGGGCAACGGCTGGCTCCgg                                                          <  2:2838681/150‑1 (MQ=255)
                                                                                              gCCCGGAAAACGCTATCGGGTCATTGACGTCATCCATAACGAATTGATTACGCACTCCCACTCCAGCGTCTACAGCGAAAATGGTTTTGATCGCGATGATGTGAGCTTTATTGCCGTACTTGAGCGTTACGGGCAACGGCTGGCTCCGGc                                                         <  1:3532795/150‑1 (MQ=255)
                                                                                                ccGGAAAACGCTATCGGGTCATTGACGTCATCCATAACGAATTGATTACGCACTCCCACTCCAGCGTCTACAGCGAAAATGGTTTTGATCGCGATGATGTGAGCTTTATTGCCGTACTTGAGCGTTACGGGCAACGGCTGGCTCCGGCtt                                                       <  1:667854/150‑1 (MQ=255)
                                                                                                    aaaaCGCTATCGGGTCATTGACGTCATCCATAACGAATTGATTACGCACACCCACTCCAGCGTCTACAGCGAAAATGGTTTTGATCGCGATGATGTTAGCTTTATTGCAGTACTTGAGCGTTACGGGCAACGGCTGGCTCCGGCTTGTgg                                                   >  2:799368/1‑150 (MQ=255)
                                                                                                      aaCGCTATCGGGTCATTGACGTCATCCATAACGAATTGATTACGCACACCCACTCCCGCGTCTACAGCGAAAATGGTTTTGATCGCGATGATGTGAGCTTTATTGCCGTACTTGAGCGTTACGGGCAACGGCTGGCTCCGGCTTGTGGtt                                                 <  2:55054/150‑1 (MQ=255)
                                                                                                       aCGCTATCGGGTCATTGACGTCATCCATAACGAATTGAGTACGCACTCCCACTCCAGCGTCTACAGCGAAAATGGTTTTGATCGCGATGATGTGAGCTTTATTGCCGTACTTGAGCGTTACGGGCAACGGCTGGCTCCGGCTTGTGGttt                                                <  2:2621355/150‑1 (MQ=255)
                                                                                                         gCTATCGGGTCATTGACGTCATCCATAACGAATTGATTACGCACTCCCACTCCAGCGTCTACAGCGAAAATGGTTTTGATCGCGATGATGTGAGCTTTATTGCCGTACTTGAGCGTTACGGGCAACGGCTGGCTCCGGCTTGTGGTTTGc                                              <  1:2975118/150‑1 (MQ=255)
                                                                                                          cTATCGGGTCATTGACGTCATCCATAACGAATTGATTACGCACTCCCACTCCAGCGTCTACAGCGAAAATGGTTTTGATCGCGATGATGTGAGCTTTATTGCCGTACTTGAGCGTTACGGGCAACGGCTGGCTCCGGCTTGTGGTTTGCt                                             <  2:167990/150‑1 (MQ=255)
                                                                                                             tCGGGTCATTGACGTCATCCATAACGAATTGATTACGCACTCCCACTCCAGCGTCTACAGCGAAAATGGTTTTGATCGCGATGATGTGAGCTTTATTGCCGTACTTGAGCGTTACGGGCAACGGCTGGCTCCGGCTTGTGGTTTGCTTgg                                          >  1:3490880/1‑150 (MQ=255)
                                                                                                              cGGGTCATTGACGTCATCCATAACGAATTGATTACGCACTCCCACTCCAGCGTCTACAGCGAAAATGGTTTTGATCGCGATGATGTGAGCTTTATTGCCGTACTTGAGCGTTACGGGCAACGGCTGGCTCCGGCTTGTGGTTTGCTTggc                                         <  2:2379199/150‑1 (MQ=255)
                                                                                                                  tCATTGACGTCATCCATAACGAATTGATTACGCACTCCCACTCCAGCGTCTACAGCGAAAATGGTTTTGATCGCGATGATGTGAGCTTTATTGCCGTACTTGAGCGTTACGGGCAACGGCTGGCTCCGGCTTGTGGTTTGCTTGGCGGCt                                     <  2:66730/150‑1 (MQ=255)
                                                                                                                   cATTGACGTCATCCATAACGAATTGATTACGCACTCCCACTCCAGCGTCTACAGCGAAAATGGTTTTGATCGCGATGATGTGAGCTTTATTGCCGTACTTGAGCGTTACGGGCAACGGCTGGCTCCGGCTTGTGGTTTGCTTGGCGGCtt                                    <  1:1021817/150‑1 (MQ=255)
                                                                                                                   cATTGACGTCATCCATAACGAATTGATTACGCACACCCACTCCAGCGTCTACAGCGAAAATGGTTTTGATCGCGATGATGTGAGCTTTATTGCCGTACTTGAGCGTTACGGGCAACGGCTGGCTCCGGCTTGTGGTTTGCTTGGCGGCtt                                    >  1:1458541/1‑150 (MQ=255)
                                                                                                                    aTTGACGTCATCCATAACGAATTGATTACGCACTCCCACTCCAGCGTCTACAGCGAAAATGGTTTTGATCGCGATGATGTGAGCTTTATTGCCGTACTTGAGCGTTACGGGCAACGGCTGGCTCCGGCTTGTGGTTTGCTTGGCGGCttt                                   <  2:35703/150‑1 (MQ=255)
                                                                                                                    aTTGACGTCATCCATAACGAATTGATTACGCACTCCCACTCCAGCGTCTACAGCGAAAATGGTTTTGATCGCGATGATGTGAGCTTTATTGCCGTACTTGAGCGTTACGGGCAACGGCTGGCTCCGGCTTGTGGTTTGCTTGGCGGCttt                                   >  1:4558463/1‑150 (MQ=255)
                                                                                                                      tGACGTCATCCATAACGAATTGATTACGCACTCCCACTCCAGCGTCTACAGCGAAAATGGTTTTGATCGCGATGATGTGAGCTTTATTGCCGTACTTGAGCGTTACGGGCAACGGCTGGCTCCGGCTTGTGGTTTGCTTGGCGGCTTTgg                                 >  2:431009/1‑150 (MQ=255)
                                                                                                                       gACGTCATCCATAACGAATTGATTACGCACACCCACTCCAGCGTCTACAGCGAAAATGGGTTTGATCGCGATGATGTGAGCTTTATTGCCGTACTTGAGCGTTACGGGCAACGGCTGGCTCCGGCTTGTGGTTTGCTTGGCGGCTTTGGa                                >  2:1157888/1‑150 (MQ=255)
                                                                                                                             aTCCATAACGAATTGATTACGCACTCCCACTCCAGCGTCTACAGCGAAAATGGTTTTGATCGCGATGATGTGAGCTTTATTGCCGTACTTGAGCGTTACGGGCAACGGCTGGCTCCGGCTTGTGGTTTGCTTGGCGGCTTTGGACtgaat                          <  2:4011628/150‑1 (MQ=255)
                                                                                                                               ccATAACGAATTGATTACGCACTCCCACTCCAGCGTCTACAGCGAAAATGGTTTTGATCGCGATGATGTGAGCTTTATTGCCGTACTTGAGCGTTACGGGCAACGGCTGGCTCCGGCTTGTGGTTTGCTTGGCGGCTTTGGACtgaatga                        >  2:2585961/1‑150 (MQ=255)
                                                                                                                               ccATAACGAATTGATTACGCACTCCCACTCCAGCGTCTACAGCGAAAATGGTTTTGATCGCGATGATGTGAGCTTTATTGCCGTACTTGAGCGTTACGGGCAACGGCTGGCTCCGGCTTGTGGTTTGCTTGGCGGCTTTGGACtgaatga                        <  2:3545421/150‑1 (MQ=255)
                                                                                                                                 aTAACGAATTGATTACGCACTCCCACTCCAGCGTCTACAGCGAAAATGGTTTTGATCGCGATGATGTGAGCTTTATTGCCGTACTTGAGCGTTACGGGCAACGGCTGGCTCCGGCTTGTGGTTTGCTTGGCGGCTTTGGACTGAATGAAg                      >  1:2977528/1‑150 (MQ=255)
                                                                                                                                 aTAACGAATTGATTACGCACTCCCACTCCAGCGTCTACAGCGAAAATGGTTTTGATCGCGATGATGTGAGCTTTATTGCCGTACTTGAGCGTTACGGGCAACGGCTGGCTCCGGCTTGTGGTTTGCTTGGCGGCTTTGGACTGAATGAAg                      >  1:2366931/1‑150 (MQ=255)
                                                                                                                                  tAACGAATTGATTACGCACTCCCACTCCAGCGTCTACAGCGAAAATGGTTTTGATCGCGATGATGTGAGCTTTATTGCCGTACTTGAGCGTTACGGGCAACGGCTGGCTCCGGCTTGTGGTTTGCTTGGCGGCTTTGGACTGAATGAAgg                     >  2:1007697/1‑150 (MQ=255)
                                                                                                                                  tAACGAATTGATTACGCACTCCCACTCCAGCGTCTACAGCGAAAATGGTTTTGATCGCGATGATGTGAGCTTTATTGCCGTACTTGAGCGTTACGGGCAACGGCTGGCTCCGGCTTGTGGTTTGCTTGGCGGCTTTGGACTGAATGAAgg                     >  2:1645728/1‑150 (MQ=255)
                                                                                                                                    aCGAATTGATTACGCACTCCCACTCCAGCGTCTACAGCGAAAATGGTTTTGATCGCGATGATGTGAGCTTTATTGCCGTACTTGAGCGTTACGGGCAACGGCTGGCTCCGGCTTGTGGTTTGCTTGGCGGCTTTGGACTGAATGAAGGTg                   >  1:4391272/1‑150 (MQ=255)
                                                                                                                                     cGAATTGATTACGCACTCCCACTCCAGCGTCTACAGCGAAAATGGTTTTGATCGCGATGATGTGAGCTTTATTGCCGTACTTGAGCGTTACGGGCAACGGCTGGCTCCGGCTTGTGGTTTGCTTGGCGGCTTTGGACTGAATGAAGGTgc                  <  2:918169/150‑1 (MQ=255)
                                                                                                                                      gAATTGATTACGCACTCCCACTCCAGCGTCTACAGCGAAAATGGTTTTGATCGCGATGATGTGAGCTTTATTGCCGTACTTGAGCGTTACGGGCAACGGCTGGCTCCGGCTTGTGGTTTGCTTGGCGGCTTTGGACTGAATGAAGGTgcg                 >  1:2298777/1‑150 (MQ=255)
                                                                                                                                       aaTTGATTACGCACTCCCACTCCAGCGTCTACAGCGAAAATGGTTTTGATCGCGATGATGTGAGCTTTATTGCCGTACTTGAGCGTTACGGGCAACGGCTGGCTCCGGCTTGTGGTTTGCTTGGCGGCTTTGGACTGAATGAAGGTgcgc                <  1:93000/150‑1 (MQ=255)
                                                                                                                                         ttGATTACGCACTCCCACTCCAGCGTCTACAGCGAAAATGGTTTTGATCGCGATGATGTGAGCTTTATTGCCGTACTTGAGCGTTACGGGCAACGGCTGGCTCCGGCGTGTGGTTTGCTTGGCGGCTTTGGACTGAAGGAAGGTGCgctg              >  2:2255912/1‑150 (MQ=255)
                                                                                                                                            aTTACGCACTCCCACTCCAGCGTCTACAGCGAAAATGGTTTTGATCGCGATGATGTGAGCTTTATTGCCGTACTTGAGCGTTACGGGCAACGGCTGGCTCCGGCTTGTGGTTTGCTTGGCGGCTTTGGACTGAATGAAGGTGCgctggct           <  1:2121783/150‑1 (MQ=255)
                                                                                                                                            aTTACGCACACCCACTCCAGCGTCTACAGCGAAAATGGTTTTGATCGCGATGATGTGAGCTTTATTGCCGTACTTGAGCGTTACGGGCAACGGCTGGCTCCGGCTTGTGGTTTGCTTGGCGGCTTTGGACTGAATGAAGGTGCgctggct           <  2:4030972/150‑1 (MQ=255)
                                                                                                                                                 gCACTCCCACTCCAGCGTCTACAGCGAAAATGGTTTTGATCGCGATGATGTGAGCTTTATTGCCGTACTTGAGCGTTACGGGCAACGGCTGGCTCCGGCTTGTGGTTTGCTTGGCGGCTTTGGACTGAATGAAGGTGCGCTGGCTGCGAc      <  2:4540762/150‑1 (MQ=255)
                                                                                                                                                  cactcccactccAGCGTCTACAGCGAAAATGGTTTTGATCGCGATGATGTGAGCTTTATTGCCGTACTTGAGCGTTACGGGCAACGGCTGGCTCCGGCTTGTGGTTTGCTTGGCGGCTTTGGACTGAATGAAGGTGCGCTGGCTGCGACg     >  2:1297725/1‑150 (MQ=255)
                                                                                                                                                   actcccactccAGCGTCTACAGCGAAAATGGTTTTGATCGCGATGATGTGAGCTTTATTGCCGTACTTGAGCGTTACGGGCAACGGCTGGCTCCGGCTTGTGGTTTGCTTGGCGGCTTTGGACTGAATGAAGGTGCGCTGGCTGCGACgg    >  1:200129/1‑150 (MQ=255)
                                                                                                                                                    ctcccactccAGCGTCTACAGCGAAAATGGTTTTGATCGCGATGATGTGAGCTTTATTGCCGTACTTGAGCGTTACGGGCAACGGCTGGCTCCGGCTTGTGGTTTGCTTGGCGGCTTTGGACTGAATGAAGGTGCGCTGGCTGCGACGGt   >  1:4257107/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AAGAGTCGGCGAGACTGGCACAATCCGCTCCGCCATATGGCAACACCATTGCCCGCCAGCCAGTTTCCGCCAGCGACTTTGCCCTGCAATTTACGCCCGGAAAACGCTATCGGGTCATTGACGTCATCCATAACGAATTGATTACGCACTCCCACTCCAGCGTCTACAGCGAAAATGGTTTTGATCGCGATGATGTGAGCTTTATTGCCGTACTTGAGCGTTACGGGCAACGGCTGGCTCCGGCTTGTGGTTTGCTTGGCGGCTTTGGACTGAATGAAGGTGCGCTGGCTGCGACGGTC  >  CP009273/3838438‑3838736

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: