Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,697,123 C→A 10.7% R212S (CGT→AGT)  add → adenosine deaminase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,697,1230CA10.7% 1009.8 / 61.1 278R212S (CGT→AGT) addadenosine deaminase
Reads supporting (aligned to +/- strand):  ref base C (99/149);  new base A (15/15);  total (114/164)
Fisher's exact test for biased strand distribution p-value = 3.28e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.72e-01

Only 100 of 286 total aligned reads displayed.

AGATTACCGCACTTGATTTAGCCGGTGATGAACTTGGTTTCCCGGGAAGTCTGTTCCTTTCTCACTTCAACCGCGCGCGTGATGCGGGCTGGCATATTACCGTCCATGCAGGCGAAGCTGCCGGGCCGGAAAGCATCTGGCAGGCGATTCGTGAACTGGGTGCGGAGCGTATTGGACATGGCGTAAAAGCCATTGAAGATCGGGCGCTGATGGATTTTCTCGCCGAGCAACAAATTGGTATTGAATCCTGTCTGACCTCCAATATTCAGACCAGCACCGTAGCAGAGCTGGCTGCAC  >  CP009273/1696974‑1697270
                                                                                                                                                     |                                                                                                                                                   
aGATTACCGCACTTGATTTAGCCGGTGATGAACTTGGTTTCCCGGGAAGTCTGTTCCTTTCTCACTTCAACCGCGCGCGTGATGCGGGCTGGCATATTACCGTCCATGCAGGCGAAGCTGCCGGGCCGGAAAGCATCTGGCAGGCGATTc                                                                                                                                                     <  2:588015/150‑1 (MQ=255)
aGATTACCGCACTTGATTTAGCCGGTGATGAACTTGGTTTCCCGGGAAGTCTGTTCCTTTCTCACTTCAACCGCGCGCGTGATGCGGGCTGGCATATTACCGTCCATGCAGGCGAAGCTGCCGGGCCGGAAAGCATCTGGCAGGCGATTc                                                                                                                                                     <  2:3581252/150‑1 (MQ=255)
aGATTACCGCACTTGATTTAGCCGGTGATGAACTTGGTTTCCCGGGAAGTCTGTTCCTTTCTCACTTCAACCGCGCGCGTGATGCGGGCTGGCATATTACCGTCCATGCAGGCGAAGCTGCCGGGCCGGAAAGCATCTGGCAGGCGATTc                                                                                                                                                     <  2:2382663/150‑1 (MQ=255)
  aTTACCGCACTTGATTTAGCCGGTGATGAACTTGGTTTCCCGGGAAGTCTGTTCCTTTCTCACTTCAACCGCGCGCGTGATGCGGGCTGGCATATTACCGTCCATGCAGGCGAAGCTGCCGGGCCGGAAAGCATCTGGCAGGCGATTAGt                                                                                                                                                   <  2:1612772/150‑1 (MQ=255)
     aCCGCACTTGATTTAGCCGGTGATGAACTTGGTTTCCCGGGAAGTCTGTTCCTTTCTCACTTCAACCGCGCGCGTGATGCGGGCTGGCATATTACCGTCCATGCAGGCGAAGCTGCCGGGCCGGAAAGCATCTGGCAGGCGATTCGTGaa                                                                                                                                                >  1:2741489/1‑150 (MQ=255)
      ccGCACTTGATTTAGCCGGTGATGAACTTGGTTTCCCGGGAAGTCTGTTCCTTTCTCACTTCAACCGCGCGCGTGATGCGGGCTGGCATATTACCGTCCATGCAGGCGAAGCTGCCGGGCCGGAAAGCATCTGGCAGGCGATTCGTGAAc                                                                                                                                               >  1:2834691/1‑150 (MQ=255)
           cTTGATTTAGCCGGTGATGAACTTGGTTTCCCGGGAAGTCTGTTCCTTTCTCACTTCAACCGCGCGCGTGATGCGGGCTGGCATATTACCGTCCATGCAGGCGAAGCTGCCGGGCCGGAAAGCATCTGGCAGGCGATTCGTGAACTGGGt                                                                                                                                          >  2:960529/1‑150 (MQ=255)
            ttGATTTAGCCGGTGATGAACTTGGTTTCCCGGGAAGTCTGTTCCTTTCTCACTTCAACCGCGCGCGTGATGCGGGCTGGCATATTACCGTCCATGCAGGCGAAGCTGCCGGGCCGGAAAGCATCTGGCAGGCGATTCGTGAACTGGGTg                                                                                                                                         <  1:4354423/150‑1 (MQ=255)
             tGATTTAGCCGGTGATGAACTTGGTTTCCCGGGAAGTCTGTTCCTTTCTCACTTCAACCGCGCGCGTGATGCGGGCTGGCATATTACCGTCCATGCAGGCGAAGCTGCCGGGCCGGAAAGCATCTGGCAGGCGATTCGTGAACTGGGTGc                                                                                                                                        <  2:340731/150‑1 (MQ=255)
             tGATTTAGCCGGTGATGAACTTGGTTTCCCGGGAAGTCTGTTCCTTTCTCACTTCAACCGCGCGCGTGATGCGGGCTGGCATATTACCGTCCATGCAGGCGAAGCTGCCGGGCCGGAAAGCATCTGGCAGGCGATTCGTGAACTGGGTGc                                                                                                                                        <  2:2850330/150‑1 (MQ=255)
              gATTTAGCCGGTGATGAACTTGGTTTCCCGGGAAGTCTGTTCCTTTCTCACTTCAACCGCGCGCGTGATGCGGGCTGGCATATTACCGTCCATGCAGGCGAAGCTGCCGGGCCGGAAAGCATCTGGCAGGCGATTCGTGAACTGGGTGCg                                                                                                                                       <  1:2690345/150‑1 (MQ=255)
                     ccGGTGATGAACTTGGTTTCCCGGGAAGTCTGTTCCTTTCTCACTTCAACCGCGCGCGTGATGCGGGCTGGCATATTACCGTCCATGCAGGCGAAGCTGCCGGGCCGGAAAGCATCTGGCAGGCGATTCGTGAACTGGGTGCGGAGCGTa                                                                                                                                <  1:3628407/150‑1 (MQ=255)
                      cGGTGATGAACTTGGTTTCCCGGGAAGTCTGTTCCTTTCTCACTTCAACCGCGCGCGTGATGCGGGCTGGCATATTACCGTCCATGCAGGCGAAGCTGCCGGGCCGGAAAGCATCTGGCAGGCGATTCGTGAACTGGGTGCGGAGCGTAt                                                                                                                               >  2:4234532/1‑150 (MQ=255)
                        gTGATGAACTTGGTTTCCCGGGAAGTCTGTTCCTTTCTCACTTCAACCGCGCGCGTGATGCGGGCTGGCATATTACCGTCCATGCAGGCGAAGCTGCCGGGCCGGAAAGCATCTGGCAGGCGATTCGTGAACTGGGTGCGGAGCGTATTg                                                                                                                             <  2:1208130/150‑1 (MQ=255)
                          gatgaACTTGGTTTCCCGGGAAGTCTGTTCCTTTCTCACTTCAACCGCGCGCGTGATGCGGGCTGGCATATTACCGTCCATGCAGGCGAAGCTGCCGGGCCGGAAAGCATCTGGCAGGCGATTCGTGAACTGGGTGCGGAGCGTATTGGa                                                                                                                           <  2:204499/150‑1 (MQ=255)
                           atgaACTTGGTTTCCCGGGAAGTCTGTTCCTTTCTCACTTCAACCGCGCGCGTGATGCGGGCTGGTATATTACCGTCCATGCAGGCGAAGCTGCCGGGCCGGAAAGCATCTGGCAGGCGATTCGTGAACTGGGTGCGGAGCGTATTGGAc                                                                                                                          <  1:1460430/150‑1 (MQ=255)
                           atgaACTTGGTTTCCCGGGAAGTCTGTTCCTTTCTCACTTCAACCGCGCGCGTGATGCGGGCTGGCCTATTACCGTCCATGCAGGCGAAGCTGCCGGGCCGGAACGCATCTGGCAGGCGATTCGTGAACTGGGTGCGGAGCGTATTGGAc                                                                                                                          >  2:3105765/1‑150 (MQ=255)
                              aaCTTGGTTTCCCGGGAAGTCTGTTCCTTTCTCACTTCAACCGCGCGCGTGATGCGGGCTGGCATATTACCGTCCATGCAGGCGAAGCTGCCGGGCCGGAAAGCATCTGGCAGGCGATTCGTGAACTGGGTGCGGAGCGTATTGGACATg                                                                                                                       >  2:2235139/1‑150 (MQ=255)
                                cTTGGTTTCCCGGGAAGTCTGTTCCTTTCTCACTTCAACCGCGCGCGTGATGCGGGCTGGCATATTACCGTCCATGCAGGCGAAGCTGCCGGGCCGGAAAGCATCTGGCAGGCGATTCGTGAACTGGGTGCGGAGCGTATTGGACATGGc                                                                                                                     <  1:2128535/150‑1 (MQ=255)
                                cTTGGTTTCCCGGGAAGTCTGTTCCTTTCTCACTTCAACCGCGCGCGTGATGCGGGCTGGCATATTACCGTCCATGCAGGCGAAGCTGCCGGGCCGGAAAGCATCTGGCAGGCGATTCGTGAACTGGGTGCGGAGCGTATTGGACATGGc                                                                                                                     <  1:3281749/150‑1 (MQ=255)
                                  tGGTTTCCCGGGAAGTCTGTTCCTTTCTCACTTCAACCGCGCGCGTGATGCGGGCTGGCATATTACCGTCCATGCAGGCGAAGCTGCCGGGCCGGAAAGCATCTGGCAGGCGATTCGTGAACTGGGTGCGGAGCGTATTGGACATGGCGt                                                                                                                   <  1:974429/150‑1 (MQ=255)
                                   ggTTTCCCGGGAAGTCTGTTCCTTTCTCACTTCAACCGCGCGCGTGATGCGGGCTGGCATATTACCGTCCATGCAGGCGAAGCTGCCGGGCCGGAAAGCATCTGGCAGGCGATTCGTGAACTGGGTGCGGAGCGTATTGGACATGGCGTa                                                                                                                  <  2:3766584/150‑1 (MQ=255)
                                      ttCCCGGGAAGTCTGTTCCTTTCTCACTTCAACCGCGCGCGTGATGCGGGCTGGCATATTACCGTCCATGCAGGCGAAGCTGCCGGGCCGGAAAGCATCTGGCAGGCGATTCGTGAACTGGGTGCGGAGCGTATTGGACATGGCGTaaaa                                                                                                               <  2:2302174/150‑1 (MQ=255)
                                         ccGGGAAGTCTGTTCCTTTCTCACTTCAACCGCGCGCGTGATGCGGGCTGGCATATTACCGTCCATGCAGGCGAAGCTGCCGGGCCGGAAAGCATCTGGCAGGCGATTCGTGAACTGGGTGCGGAGCGTATTGGACATGGCGTAAAAGcc                                                                                                            <  2:3267617/150‑1 (MQ=255)
                                          cGGGAAGTCTGTTCCTTTCTCACTTCAACCGCGCGCGTGATGCGGGCTGGCATATTACCGTCCATGCAGGCGAAGCTGCCGGGCCGGAAAGCATCTGGCAGGCGATTCGTGAACTGGGTGCGGAGCGTATTGGACATGGCGTAAAAGCCa                                                                                                           >  1:1432836/1‑150 (MQ=255)
                                           gggAAGTCTGTTCCTTTCTCACTTCAACCGCGCGCGTGATGCGGGCTGGCATATTACCGTCCATGCAGGCGAAGCTGCCGGGCCGGAAAGCATCTGGCAGGCGATTCGTGAACTGGGTGCGGAGCGTATTGGACATGGCGTAAAAGCCAt                                                                                                          >  1:69066/1‑150 (MQ=255)
                                                gTCTGTTCCTTTCTCACTTCAACCGCGCGCGTGATGCGGGCTGGCATATTACCGTCCATGCAGGCGAAGCTGCCGGGCCGGAAAGCATCTGGCAGGCGATTCGTGAACTGGGTGCGGAGCGTATTGGACATGGCGTAAAAGCCATTGAAg                                                                                                     <  2:1356459/150‑1 (MQ=255)
                                                     ttCCTTTCTCACTTCAACCGCGCGCGTGATGCGGGCTGGCATATTACCGTCCATGCAGGCGAAGCTGCCGGGCCGGAAAGCATCTGGCAGGCGATTCGTGAACTGGGTGCGGAGCGTATTGGACATGGCGTAAAAGCCATTGAAGATCgg                                                                                                <  1:1178722/150‑1 (MQ=255)
                                                      tCCTTTCTCACTTCAACCGCGCGCGTGATGCGGGCTGGCATATTACCGTCCATGCAGGCGAAGCTGCCGGGCCGGAAAGCATCTGGCAGGCGATTAGTGAACTGGGTGCGGAGCGTATTGGACATGGCGTAAAAGCCATTGAAGATCggg                                                                                               <  2:3293470/150‑1 (MQ=255)
                                                        cTTTCTCACTTCAACCGCGCGCGTGATGCGGGCTGGCATATTACCGTCCATGCAGGCGAAGCTGCCGGGCCGGAAAGCATCTGGCAGGCGATTCGTGAACTGGGTGCGGAGCGTATTGGACATGGCGTAAAAGCCATTGAAGATCGGgcg                                                                                             >  2:1543145/1‑150 (MQ=255)
                                                         tttCTCACTTCAACCGCGCGCGTGATGCGGGCTGGCATATTACCGTCCATGCAGGCGAAGCTGCCGGGCCGGAAAGCATCTGGCAGGCGATTCGTGAACTGGGTGCGGAGCGTATTGGACATGGCGTAAAAGCCATTGAAGATCGGgcgc                                                                                            <  1:1068924/150‑1 (MQ=255)
                                                         tttCTCACTTCAACCGCGCGCGTGATGCGGGCTGGCATATTACCGTCCATGCAGGCGAAGCTGCCGGGCCGGAAAGCATCTGGCAGGCGATTCGTGAACTGGGTGCGGAGCGTATTGGACATGGCGTAAAAGCCATTGAAGATCGGgcgc                                                                                            >  1:1037120/1‑150 (MQ=255)
                                                          ttCTCACTTCAACCGCGCGCGTGATGCGGGCTGGCATATTACCGTCCATGCAGGCGAAGCTGCCGGGCCGGAAAGCATCTGGCAGGCGATTCGTGAACTGGGTGCGGAGCGTATTGGACATGGCGTAAAAGCCATTGAAGATCGGGCGCt                                                                                           >  1:1393556/1‑150 (MQ=255)
                                                            ctcACTTCAACCGCGCGCGTGATGCGGGCTGGCATATTACCGTCCATGCAGGCGAAGCTGCCGGGCCGGAAAGCATCTGGCAGGCGATTCGTGAACTGGGTGCGGAGCGTATTGGACATGGCGTAAAAGCCATTGAAGATCGGGCGCTGa                                                                                         <  2:32584/150‑1 (MQ=255)
                                                               aCTTCAACCGCGCGCGTGATGCGGGCTGGCATATTACCGTCCATGCAGGCGAAGCTGCCGGGCCGGAAAGCATCTGGCAGGCGATTCGTGAACTGGGTGCGGAGCGTATTGGACATGGCGTAAAAGCCATTGAAGATCGGGCGCTGATgg                                                                                      <  2:3198207/150‑1 (MQ=255)
                                                                cTTCAACCGCGCGCGTGATGCGGGCTGGCATATTACCGTCCATGCAGGCGAAGCTGCCGGGCCGGAAAGCATCTGGCAGGCGATTCGTGAACTGGGTGCGGAGCGTATTGGACATGGCGTAAAAGCCATTGAAGATCGGGCGCTGATGGa                                                                                     >  2:4606590/1‑150 (MQ=255)
                                                                 ttCAACCGCGCGCGTGATGCGGGCTGGCATATTACCGTCCATGCAGGCGAAGCTGCCGGGCCGGAAAGCATCTGGCAGGCGATTCGTGAACTGGGTGCGGAGCGTATTGGACATGGCGTAAAAGCCATTGAAGATCGGGCGCTGATGGAt                                                                                    >  1:2565942/1‑150 (MQ=255)
                                                                   cAACCGCGCGCGTGATGCGGGCTGGCATATTACCGTCCATGCAGGCGAAGCTGCCGGGCCGGAAAGCATCTGGCAGGCGATTCGTGAACTGGGTGCGGAGCGTATTGGACATGGCGTAAAAGCCATTGAAGATCGGGCGCTGATGGAttt                                                                                  >  2:480363/1‑150 (MQ=255)
                                                                   cAACCGCGCGCGTGATGCGGGCTGGCATATTACCGTCCATGCAGGCGAAGCTGCCGGGCCGGAAAGCATCTGGCAGGCGATTCGTGAACTGGGTGCGGAGCGTATTGGACATGGCGTAAAAGCCATTGAAGATCGGGCGCTGATGGAttt                                                                                  <  1:188911/150‑1 (MQ=255)
                                                                   cAACCGCGCGCGTGATGCGGGCTGGCATATTACCGTCCATGCAGGCGAAGCTGCCGGGCCGGAAAGCATCTGGCAGGCGATTCGTGAACTGGGTGCGGAGCGTATTGGACATGGCGTAAAAGCCATTGAAGATCGGGCGCTGATGGAttt                                                                                  >  1:1916628/1‑150 (MQ=255)
                                                                    aaCCGCGCGCGTGATGCGGGCTGGCATATTACCGTCCATGCAGGCGAAGCTGCCGGGCCGGAAAGCATCTGGCAGGCGATTCGTGAACTGGGTGCGGAGCGTATTGGACATGGCGTAAAAGCCATTGAAGATCGGGCGCTGATGGAtttt                                                                                 <  2:3010048/150‑1 (MQ=255)
                                                                     aCCGCGCGCGTGATGCGGGCTGGCATATTACCGTCCATGCAGGCGAAGCTGCCGGGCCGGAAAGCATCTGGCAGGCGATTCGTGAACTGGGTGCGGAGCGTATTGGACATGGCGTAAAAGCCATTGAAGATCGGGCGCTGATGGATTTtc                                                                                <  2:1225865/150‑1 (MQ=255)
                                                                     aCCGCGCGCGTGATGCGGGCTGGCATATTACCGTCCATGCAGGCGAAGCTGCCGGGCCGGAAAGCATCTGGCAGGCGATTCGTGAACTGGGTGCGGAGCGTATTGGACATGGCGTAAAAGCCATTGAAGATCGGGCGCTGATGGATTTtc                                                                                <  1:161421/150‑1 (MQ=255)
                                                                      ccGCGCGCGTGATGCGGGCTGGCATATTACCGTCCATGCAGGCGAAGCTGCCGGGCCGGAAAGCATCTGGCAGGCGATTAGTGAACTGGGTGCGGAGCGTATTGGACATGGCGTAAAAGCCATTGAAGATCGGGCGCTGATGGATTTtct                                                                               >  1:217770/1‑150 (MQ=255)
                                                                       cgcgcgcgTGATGCGGGCTGGCATATTACCGTCCATGCAGGCGAAGCTGCCGGGCCGGAAAGCATCTGGCAGGCGATTAGTGAACTGGGTGCGGAGCGTATTGGACATGGCGTAAAAGCCATTGAAGATCGGGCGCTGATGGATTTtctc                                                                              <  2:3346606/150‑1 (MQ=255)
                                                                       cgcgcgcgTGATGCGGGCTGGCATATTACCGTCCATGCAGGCGAAGCTGACGGGCCGGAAAGCATCTGGCAGGCGATTCGTGAACTGGGTGCGGAGCGTATTGGACATGGCGTAAAAGCCATTGAAGATCGGGCGCTGATGGATTTtctc                                                                              >  2:2716004/1‑150 (MQ=255)
                                                                          gcgcgTGATGCGGGCTGGCATATTACCGTCCATGCAGGCGAAGCTGCCGGGCCGGAAAGCATCTGGCAGGCGATTCGTGAACTGGGTGCGGAGCGTATTGGACATGGCGTAAAAGCCATTGAAGATCGGGCGCTGATGGATTTTCTCGcc                                                                           <  1:43869/150‑1 (MQ=255)
                                                                          gcgcgTGATGCGGGCTGGCATATTACCGTCCATGCAGGCGAAGCTGCCGGGCCGGAAAGCATCTGGCAGGCGATTCGTGAACTGGGTGCGGAGCGTATTGGACATGGCGTAAAAGCCATTGAAGATCGGGCGCTGATGGATTTTCTCGcc                                                                           <  1:1705172/150‑1 (MQ=255)
                                                                               tGATGCGGGCTGGCATATTACCGTCCATGCAGGCGAAGCTGCCGGGCCGGAAAGCATCTGGCAGGCGATTCGTGAACTGGGTGCGGAGCGTATTGGACATGGCGTAAAAGCCATTGAAGATCGGGCGCTGATGGATTTTCTCGCCGAGca                                                                      >  1:1109743/1‑150 (MQ=255)
                                                                                gATGCGGGCTGGCATATTACCGTCCATGCAGGCGAAGCTGCCGGGCCGGAAAGCATCTGGCAGGCGATTCGTGAACTGGGTGCGGAGCGTATTGGACATGGCGTAAAAGCCATTGAAGATCGGGCGCTGATGGATTTTCTCGCCGAGcaa                                                                     >  1:884592/1‑150 (MQ=255)
                                                                                  tGCGGGCTGGCATATTACCGTCCATGCAGGCGAAGCTGCCGGGCCGGAAAGCATCTGGCAGGCGATTAGTGAACTGGGTGCGGAGCGTATTGGACATGGCGTAAAAGCCATTGAAGATCGGGCGCTGATGGATTTTCTCGCCGAGcaaca                                                                   >  2:3839647/1‑150 (MQ=255)
                                                                                      ggCTGGCATATTACCGTCCATGCAGGCGAAGCTGCCGGGCCGGAAAGCATCTGGCAGGCGATTCGTGAACTGGGTGCGGAGCGTATTGGACATGGCGTAAAAGCCATTGAAGATCGGGCGCTGATGGATTTTCTCGCCGAGCAACAAAtt                                                               <  2:119021/150‑1 (MQ=255)
                                                                                        cTGGCATATTACCGTCCATGCAGGCGAAGCTGCCGGGCCGGAAAGCATCTGGCAGGCGATTCGTGAACTGGGTGCGGAGCGTATTGGACATGGCGTAAAAGCCATTGAAGATCGGGCGCTGATGGATTTTCTCGCCGAGCAACAAATTgg                                                             <  1:3793002/150‑1 (MQ=255)
                                                                                        cTGGCATATTACCGTCCATGCAGGCGAAGCTGCCGGGCCGGAAAGCATCTGGCAGGCGATTCGTGAACTGGGTGCGGAGCGTATTGGACATGGCGTAAAAGCCATTGAAGATCGGGCGCTGATGGATTTTCTCGCCGAGCAACAAATTgg                                                             >  2:839153/1‑150 (MQ=255)
                                                                                           gCATATTACCGTCCATGCAGGCGAAGCTGCCGGGCCGGAAAGCATCTGGCAGGCGATTCGTGAACTGGGTGCGGAGCGTATTGGACATGGCGTAAAAGCCATTGAAGATCGGGCGCTGATGGATTTTCTCGCCGAGCAACAAATTGGTAt                                                          <  2:705603/150‑1 (MQ=255)
                                                                                             atatTACCGTCCATGCAGGCGAAGCTGCCGGGCCGGAAAGCATCTGGCAGGCGATTCGTGAACTGGGTGCGGAGCGTATTGGACATGGCGTAAAAGCCATTGAAGATCGGGCGCTGATGGATTTTCTCGCCGAGCAACAAATTGGTATTg                                                        >  1:767444/1‑150 (MQ=255)
                                                                                             atatTACCGTCCATGCAGGCGAAGCTGCCGGGCCGGAAAGCATCTGGCAGGCGATTAGTGAACTGGGTGCGGAGCGTATTGGACATGGCGTAAAAGCCATTGAAGATCGGGCGCTGATGGATTTTCTCGCCGAGCAACAAATTGGTATTg                                                        <  1:2257475/150‑1 (MQ=255)
                                                                                             atatTACCGTCCATGCAGGCGAAGCTGCCGGGACGGAAAGCATTTGGCAGGCGATTCGTGAACTGGGTGCGGAGCGTATTGGACATGGCGAAAAAGCCATTGAAGATCGGGCGCTGATGGATTTTCTCGCCGAGCAACAAATTGGTATTg                                                        >  2:3652480/1‑150 (MQ=255)
                                                                                              tatTACCGTCCATGCAGGCGAAGCTGCCGGGCCGGAAAGCATCTGGCAGGCGATTCGTGAACTGGGTGCGGAGCGTATTGGACATGGCGTAAAAGCCATTGAAGATCGGGCGCTGATGGATTTTCTCGCCGAGCAACAAATTGGTATTGa                                                       <  1:386328/150‑1 (MQ=255)
                                                                                                 tACCGTCCATGCAGGCGAAGCTGCCGGGCCGGAAAGCATCTGGCAGGCGATTCGTGAACTGGGTGCGGAGCGTATTGGACATGGCGTAAAAGCCATTGAAGATCGGGCGCTGATGGATTTTCTCGCCGAGCAACAAATTGGTATTGAATc                                                    <  2:392541/150‑1 (MQ=255)
                                                                                                 tACCGTCCATGCAGGCGAAGCTGCCGGGCCGGAAAGCATCTGGCAGGCGATTCGTGAACTGGGTGCGGAGCGTATTGGACATGGCGTAAAAGCCATTGAAGATCGGGCGCTGATGGATTTTCTCGCCGAGCAACAAATTGGTATTGAATc                                                    <  2:1435491/150‑1 (MQ=255)
                                                                                                 tACCGTCCATGCAGGCGAAGCTGCCGGGCCGGAAAGCATCTGGCAGGCGATTCGTGAACTGGGTGCGGAGCGTATTGGACATGGCGTAAAAGCCATTGAAGATCGGGCGCTGATGGATTTTCTCGCCGAGCAACAAATTGGTATTGAATc                                                    <  2:644186/150‑1 (MQ=255)
                                                                                                   ccGTCCATGCAGGCGAAGCTGCCGGGCCGGAAAGCATCTGGCAGGCGATTCGTGAACTGGGTGCGGAGCGTATTGGACATGGCGTAAAAGTCATTGAAGATCGGGCGCTGATGGATTTTCTCGCCGAGCAACAAATTGGTATTGAATCCt                                                  <  1:2639236/150‑1 (MQ=255)
                                                                                                    cGTCCATGCAGGCGAAGCTGCCGGGCCGGAAAGCATCTGGCAGGCGATTCGTGAACTGGGTGCGGAGCGTATTGGACATGGCGTAAAAGCCATTGAAGATCGGGCGCTGATGGATTTTCTCGCCGAGCAACAAATTGGTATTGAATCCTg                                                 >  1:2996135/1‑150 (MQ=255)
                                                                                                    cGTCCATGCAGGCGAAGCTGCCGGGCCGGAAAGCATCTGGCAGGCGATTCGTGAACTGGGTGCGGAGCGTATTGGACATGGCGTAAAAGCCATTGAAGATCGGGCGCTGATGGATTTTCTCGCCGAGCAACAAATTGGTATTGAATCCTg                                                 <  2:4064945/150‑1 (MQ=255)
                                                                                                    cGTCCATGCAGGCGAAGCTGCCGGGCCGGAAAGCATCTGGCAGGCGATTCGTGAACTGGGTGCGGAGCGTATTGGACATGGCGTAAAAGCCATTGAAGATCGGGCGCTGATGGATTTTCTCGCCGAGCAACAAATTGGTATTGAATCCTg                                                 >  2:2574627/1‑150 (MQ=255)
                                                                                                     gTCCATGCAGGCGAAGCTGCCGGGCCGGAAAGCATCTGGCAGGCGATTCGTGAACTGGGTGCGGAGCGTATTGGACATGGCGTAAAAGCCATTGAAGATCGGGCGCTGATGGATTTTCTCGCCGAGCAACAAATTGGTATTGAATCCTGt                                                >  2:1079630/1‑150 (MQ=255)
                                                                                                            cAGGCGAAGCTGCCGGGCCGGAAAGCATCTGGCAGGCGATTCGTGAACTGGGTGCGGAGCGTATTGGACATGGCGTAAAAGCCATTGAAGATCGGGCGCTGATGGATTTTCTCGCCGAGCAACAAATTGGTATTGAATCCTGTCTGACCt                                         <  2:1787968/150‑1 (MQ=255)
                                                                                                            cAGGCGAAGCTGCCGGGCCGGAAAGCATCTGGCAGGCGATTCGTGAACTGGGTGCGGAGCGTATTGGACATGGCGTAAAAGCCATTGAAGATCGGGCGCTGATGGATTTTCTCGCCGAGCAACAAATTGGTATTGAATCCTGTCTGACCt                                         <  1:3355991/150‑1 (MQ=255)
                                                                                                            cAGGCGAAGCTGCCGGGCCGGAAAGCATCTGGCAGGCGATTCGTGAACTGGGTGCGGAGCGTATTGGACATGGCGTAAAAGCCATTGAAGATCGGGCGCTGATGGATTTTCTCGCCGAGCAACAAATTGGTATTGAATCCTGTCTGACCt                                         <  2:4514091/150‑1 (MQ=255)
                                                                                                            cAGGCGAAGCTGCCGGGCCGGAAAGCATCTGGCAGGCGATTAGTGAACTGGGTGCGGAGCGTATTGGACATGGCGTAAAAGCCATTGAAGATCGGGCGCTGATGGATTTTCTCGCCGAGCAACAAATTGGTATTGAATCCTGTCTGACCt                                         >  1:2364886/1‑150 (MQ=255)
                                                                                                            cAGGCGAAGCTGCCGGGCCGGAAAGCATCGGGCAGGCGATTCGTGAACTGGGTGCGGAGCGTATTGGACATGGCGTAAAAGCCATTGAAGATCGGGCGCTGATGGATTTTCTCGCCGAGCAACAAATTGGTATTGAATCCTGTCTGACCt                                         <  1:3478851/150‑1 (MQ=255)
                                                                                                              ggCGAAGCTGCCGGGCCGGAAAGCATCTGGCAGGCGATTCGTGAACTGGGTGCGGAGCGTATTGGACATGGCGTAAAAGCCATTGAAGATCGGGCGCTGATGGATTTTCTCGCCGAGCAACAAATTGGTATTGAATCCTGTCTGACCTcc                                       <  1:3506204/150‑1 (MQ=255)
                                                                                                                cGAAGCTGCCGGGCCGGAAAGCATCTGGCAGGCGATTCGTGAACTGGGTGCGGAGCGTATTGGACATGGCGTAAAAGCCATTGAAGATCGGGCGCTGATGGATTTTCTCGCCGAGCAACAAATTGGTATTGAATCCTGTCTGACCTCCaa                                     >  1:1748634/1‑150 (MQ=255)
                                                                                                                 gAAGCTGCCGGGCCGGAAAGCATCTGGCAGGCGATTCGTGAACTGGGTGCGGAGCGTATTGGACATGGCGTAAAAGCCATTGAAGATCGGGCGCTGATGGATTTTCTCGCCGAGCAACAAATTGGTATTGAATCCTGTCTGACCTCCAat                                    <  2:2635374/150‑1 (MQ=255)
                                                                                                                 gAAGCTGCCGGGCCGGAAAGCATCTGGCAGGCGATTCGTGAACTGGGTGCGGAGCGTATTGGACATGGCGTAAAAGCCATTGAAGATCGGGCGCTGATGGATTTTCTCGCCGAGCAACAAATTGGTATTGAATCCTGTCTGACCTCCAat                                    >  2:160450/1‑150 (MQ=255)
                                                                                                                 gAAGCTGCCGGGCCGGAAAGCATCTGGCAGGCGATTCGTGAACTGGGTGCGGAGCGTATTGGACATGGCGTAAAAGCCATTGAAGATCGGGCGCTGATGGATTTTCTCGCCGAGCAACAAATTGGTATTGAATCCTGTCTGACCTCCAat                                    <  2:2892020/150‑1 (MQ=255)
                                                                                                                  aaGCTGCCGGGCCGGAAAGCATCTGGCAGGCGATTCGTGAACTGGGTGCGGAGCGTATTGGACATGGCGTAAAAGCCATTGAAGATCGGGCGCTGATGGATTTTCTCGCCGAGCAACAAATTGGTATTGAATCCTGTCTGACCTCCAata                                   <  2:317726/150‑1 (MQ=255)
                                                                                                                  aaGCTGCCGGGCCGGAAAGCATCTGGCAGGCGATTAGTGAACTGGGTGCGGAGCGTATTGGACATGGCGTAAAAGCCATTGAAGATCGGGCGCTGATGGATTTTCTCGCCGAGCAACAAATTGGTATTGAATCCTGTCTGACCTCCAata                                   >  1:1228503/1‑150 (MQ=255)
                                                                                                                    gCTGCCGGGCCGGAAAGCATCTGGCAGGCGATTCGTGAACTGGGTGCGGAGCGTATTGGACATGGCGTAAAAGCCATTGAAGATCGGGCGCTGATGGATTTTCTCGCCGAGCAACAAATTGGTATTGAATCCTGTCTGACCTCCAATAtt                                 <  2:1000045/150‑1 (MQ=255)
                                                                                                                     cTGCCGGGCCGGAAAGCATCTGGCAGGCGATTCGTGAACTGGGTGCGGAGCGTATTGGACATGGCGTAAAAGCCATTGAAGATCGGGCGCTGATGGATTTTCTCGCCGAGCAACAAATTGGTATTGAATCCTGTCTGACCTCCAATATTc                                <  1:2823343/150‑1 (MQ=255)
                                                                                                                     cTGCCGGGCCGGAAAGCATCTGGCAGGCGATTCGTGAACTGGGTGCGGAGCGTATTGGACATGGCGTAAAAGCCATTGAAGATCGGGCGCTGATGGATTTTCTCGCCGAGCAACAAATTGGTATTGAATCCTGTCTGACCTCCAATATTc                                <  1:1318686/150‑1 (MQ=255)
                                                                                                                       gccgggccggAAAGCATCTGGCAGGCGATTCGTGAACTGGGTGCGGAGCGTATTGGACATGGCGTAAAAGCCATTGAAGATCGGGCGCTGATGGATTTTCTCGCCGAGCAACAAATTGGTATTGAATCCTGTCTGACCTCCAATATTCAg                              <  2:2433301/150‑1 (MQ=255)
                                                                                                                          gggccggAAAGCATCTGGCAGGCGATTCGTGAACTGGGTGCGGAGCGTATTGGACATGGCGTAAAAGCCATTGAAGATCGGGCGCTGATGGATTTTCTCGCCGAGCAACAAATTGGTATTGAATCCTGTCTGACCTCCAATATTCAGAcc                           >  2:3027600/1‑150 (MQ=255)
                                                                                                                          gggccggAAAGCATCTGGCAGGCGATTCGTGAACTGGGTGCGGAGCGTATTGGACATGGCGTAAAAGCCATTGAAGATCGGGCGCTGATGGATTTTCTCGCCGAGCAACAAATTGGTATTGAATCCTGTCTGACCTCCAATATTCAGAcc                           >  2:4138958/1‑150 (MQ=255)
                                                                                                                               ggAAAGCATCTGGCAGGCGATTCGTGAACTGGGTGCGGAGCGTATTGGACATGGCGTAAAAGCCATTGAAGATCGGGCGCTGATGGATTTTCTCGCCGAGCAACAAATTGGTATTGAATCCTGTCTGACCTCCAATATTCAGACCAGCAc                      >  1:498208/1‑150 (MQ=255)
                                                                                                                                  aaGCATCTGGCAGGCGATTCGTGAACTGGGTGCGGAGCGTATTGGACATGGCGTAAAAGCCATTGAAGATCGGGCGCTGATGGATTTTCTCGCCGAGCAACAAATTGGTATTGAATCCTGTCTGACCTCCAATATTCAGACCAGCACCGt                   >  2:3503327/1‑150 (MQ=255)
                                                                                                                                   aGCATCTGGCAGGCGATTCGTGAACTGGGTGCGGAGCGTATTGGACATGGCGTAAAAGCCATTGAAGATCGGGCGCTGATGGATTTTCTCGCCGAGCAACAAATTGGTATTGAATCCTGTCTGACCTCCAATATTCAGACCAGCACCGTa                  <  1:3278659/150‑1 (MQ=255)
                                                                                                                                    gCATCTGGCAGGCGATTCGTGAACTGGGTGCGGAGCGTATTGGACATGGCGTAAAAGCCATTGAAGATCGGGCGCTGATGGATTTTCTCGCCGAGCAACAAATTGGTATTGAATCCTGTCTGACCTCCAATATTCAGACCAGCACCGTAg                 <  1:2102622/150‑1 (MQ=255)
                                                                                                                                     cATCTGGCAGGCGATTCGTGAACTGGGTGCGGAGCGTATTGGACATGGCGTAAAAGCCATTGAAGATCGGGCGCTGATGGATTTTCTCGCCGAGCAACAAATTGGTATTGAATCCTGTCTGACCTCCAATATTCAGACCAGCACCGTAGc                >  2:1673247/1‑150 (MQ=255)
                                                                                                                                      aTCTGGCAGGAGATTCGTGAACTGGGTGCGGAGCGTATTGGACATGGCGTAAAAGCCATTGAAGATCGGGCGCTGATGGATTTTCTCGCCGAGCAACAAATTGGTATTGAATCCTGTCTGACCTCCAATATTCAGACCAGCACCGTAGCa               <  2:293107/150‑1 (MQ=255)
                                                                                                                                           gCAGGCGATTCGTGAACTGGGTGCGGAGCGTATTGGACATGGCGTAAAAGCCATTGAAGATCGGGCGCTGATGGATTTTCTCGCCGAGCAACAAATTGGTATTGAATCCTGTCTGACCTCCAATATTCAGACCAGCACCGTAGCAGAgct          <  1:4021332/150‑1 (MQ=255)
                                                                                                                                           gCAGGCGATTAGTGAACTGGGTGCGGAGCGTATTGGACATGGCGTAAAAGCCATTGAAGATCGGGCGCTGATGGATTTTCTCGCCGAGCAACAAATTGGTATTGAATCCTGTCTGACCTCCAATATTCAGACCAGCACCGTAGCAGAgct          >  1:677945/1‑150 (MQ=255)
                                                                                                                                             aGGCGATTCGTGAACTGGGTGCGGAGCGTATTGGACATGGCGTAAAAGCCATTGAAGATCGGGCGCTGATGGATTTTCTCGCCGAGCAACAAATTGGTATTGAATCCTGTCTGACCTCCAATATTCAGACCAGCACCGTAGCAGAgctgg        <  2:1408282/150‑1 (MQ=255)
                                                                                                                                             aGGCGATTCGTGAACTGGGTGCGGAGCGTATTGGACATGGCGTAAAAGCCATTGAAGATCGGGCGCTGATGGATTTTCTCGCCGAGCAACAAATTGGTATTGAATCCTGTCTGACCTCCAATATTCAGACCAGCACCGTAGCAGAgctgg        >  1:3185994/1‑150 (MQ=255)
                                                                                                                                               gCGATTCGTGAACTGGGTGCGGAGCGTATTGGACATGGCGTAAAAGCCATTGAAGATCGGGCGCTGATGGATTTTCTCGCCGAGCAACAAATTGGTATTGAATCCTGTCTGACCTCCAATATTCAGACCAGCACCGTAGCAGAgctggct      <  2:2542253/150‑1 (MQ=255)
                                                                                                                                               gCGATTCGTGAACTGGGTGCGGAGCGTATTGGACATGGCGTAAAAGCCATTGAAGATCGGGCGCTGATGGATTTTCTCGCCGAGCAACAAATTGGTATTGAATCCTGTCTGACCTCCAATATTCAGACCAGCACCGTAGCAGAgctggct      <  1:468004/150‑1 (MQ=255)
                                                                                                                                               gCGATTCGTGAACTGGGTGCGGAGCGTATTGGACATGGCGTAAAAGCCATTGAAGATCGGGCGCTGATGGATTTTCTCGCCGAGCAACAAATTGGTATTGAATCCTGTCTGACCTCCAATATTCAGACCAGCACCGTAGCAGAgctggct      <  1:1058131/150‑1 (MQ=255)
                                                                                                                                                cGATTCGTGAACTGGGTGCGGAGCGTATTGGACATGGCGTAAAAGCCATTGAAGATCGGGCGCTGATGGATTTTCTCGCCGAGCAACAAATTGGTATTGAATCCTGTCTGACCTCCAATATTCAGACCAGCACCGTAGCAGAgctggctg     <  1:4150634/150‑1 (MQ=255)
                                                                                                                                                 gATTCGTGAACTGGGTGCGGAGCGTATTGGACATGGCGTAAAAGCCATTGAAGATCGGGCGCTGATGGATTTTCTCGCCGAGCAACAAATTGGTATTGAATCCTGTCTGACCTCCAATATTCAGACCAGCACCGTAGCAGAGCTGGCTGc    >  2:2107414/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
AGATTACCGCACTTGATTTAGCCGGTGATGAACTTGGTTTCCCGGGAAGTCTGTTCCTTTCTCACTTCAACCGCGCGCGTGATGCGGGCTGGCATATTACCGTCCATGCAGGCGAAGCTGCCGGGCCGGAAAGCATCTGGCAGGCGATTCGTGAACTGGGTGCGGAGCGTATTGGACATGGCGTAAAAGCCATTGAAGATCGGGCGCTGATGGATTTTCTCGCCGAGCAACAAATTGGTATTGAATCCTGTCTGACCTCCAATATTCAGACCAGCACCGTAGCAGAGCTGGCTGCAC  >  CP009273/1696974‑1697270

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: ‑