Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,019,203 G→T 27.1% E27* (GAG→TAG)  guaD → guanine deaminase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,019,2030GT27.1% 595.0 / 191.8 267E27* (GAG→TAG) guaDguanine deaminase
Reads supporting (aligned to +/- strand):  ref base G (107/87);  new base T (33/40);  total (140/127)
Fisher's exact test for biased strand distribution p-value = 1.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 279 total aligned reads displayed.

TACCGACAGGAAAACGTTCTGCCTGGTATTACCTCAGCGGAAGAGATGGATTAACAGTAAAGGAGTCAATGATGTCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCTGACACCATTCGCGTGCGCGACTATCGCGGCAAACTGATAGTACCG  >  CP009273/3019057‑3019352
                                                                                                                                                  |                                                                                                                                                     
tACCGACAGGAAAACGTTCTGCCTGGTATTACCTCAGCGGAAGAGATGGATTAACAGTAAAGGAGTCAATGATGTCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAtaga                                                                                                                                                    >  1:4447080/1‑150 (MQ=255)
tACCGACAGGAAAACGTTCTGCCTGGTATTACCTCAGCGGAAGAGATGGATTAACAGTAAAGGAGTCAATGATGTCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAgaga                                                                                                                                                    >  2:3207563/1‑150 (MQ=255)
   cGACAGGAAAACGTTCTGCCTGGTATTACCTCAGCGGAAGAGATGGATTAACAGTAAAGGAGTCAATGATGTCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAATAGATTg                                                                                                                                                 >  2:1025979/1‑150 (MQ=255)
      cAGGAAAACGTTCTGCCTGGTATTACCTCAGCGGAAGAGATGGATTAACAGTAAAGGAGTCAATGATGTCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAATAGATTGCct                                                                                                                                              <  2:892285/150‑1 (MQ=255)
        ggAAAACGTTCTGCCTGGTATTACCTCAGCGGAAGAGATGGATTAACAGTAAAGGAGTCAATGATGTCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCctct                                                                                                                                            <  1:1292383/150‑1 (MQ=255)
        ggAAAACGTTCTGCCTGGTATTACCTCAGCGGAAGAGATGGATTAACAGTAAAGGAGTCAATGATGTCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCctct                                                                                                                                            <  2:2999037/150‑1 (MQ=255)
         gAAAACGTTCTGCCTGGTATTACCTCAGCGGAAGAGATGGATTAACAGTAAAGGAGTCAATGATGTCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTctg                                                                                                                                           <  2:2048558/150‑1 (MQ=255)
               gTTCTGCCTGGTATTACCTCAGCGGAAGAGATGGATTAACAGTAAAGGAGTCAATGATGTCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAATAGATTGCCTctgcgctgc                                                                                                                                     <  2:603834/150‑1 (MQ=255)
               gTTCTGCCTGGTATTACCTCAGCGGAAGAGATGGATTAACAGTAAAGGAGTCAATGATGTCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTctgcgctgc                                                                                                                                     >  2:3957214/1‑150 (MQ=255)
               gTTCTGCCTGGTATTACCTCAGCGGAAGAGATGGATTAACAGTAAAGGAGTCAATGATGTCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTctgcgctgc                                                                                                                                     >  1:2547056/1‑150 (MQ=255)
                  cTGCCTGGTATTACCTCAGCGGAAGAGATGGATTAACAGTAAAGGAGTCAATGATGTCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGt                                                                                                                                  >  1:4180631/1‑150 (MQ=255)
                  cTGCCTGGTATTACCTCAGCGGAAGAGATGGATTAACAGTAAAGGAGTCAATGATGTCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGt                                                                                                                                  <  2:1252919/150‑1 (MQ=255)
                     ccTGGTATTACCTCAGCGGAAGAGATGGATTAACAGTAAAGGAGTCAATGATGTCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTa                                                                                                                               <  1:900375/150‑1 (MQ=255)
                      cTGGTATTACCTCAGCGGAAGAGATGGATTAACAGTAAAGGAGTCAATGATGTCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTAt                                                                                                                              <  2:4241082/150‑1 (MQ=255)
                      cTGGTATTACCTCAGCGGAAGAGATGGATTAACAGTAAAGGAGTCAATGATGTCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTAt                                                                                                                              <  2:4484533/150‑1 (MQ=255)
                       tGGTATTACCTCAGCGGAAGAGATGGATTAACAGTAAAGGAGTCAATGATGTCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTAtt                                                                                                                             <  2:3669685/150‑1 (MQ=255)
                        ggTATTACCTCAGCGGACGAGATGGATTAACAGTAAAGGAGTCAATGATGTCAGGAGAACACACGTTAAAAGCGGTCCGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGAGTGCCGCTGCGCTGCgggggcttg                                                                                                                            >  2:3430515/1‑143 (MQ=255)
                         gTATTACCTCAGCGGAAGAGATGGATTAACAGTAAAGGAGTCAATGATGTCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAATAGATTGCCTCTGCGCTGCGGTTTATTGa                                                                                                                           >  2:3032687/1‑150 (MQ=255)
                           aTTACCTCAGCGGAAGAGATGGATTAACAGTAAAGGAGTCAATGATGTCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAgg                                                                                                                         >  1:2824626/1‑150 (MQ=255)
                            ttACCTCAGCGGAAGAGATGGATTAACAGTAAAGGAGTCAATGATGTCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGa                                                                                                                        >  1:1819752/1‑150 (MQ=255)
                                cTCAGCGGAAGAGATGGATTAACAGTAAAGGAGTCAATGATGTCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGt                                                                                                                    >  2:3158841/1‑150 (MQ=255)
                                 tCAGCGGACGAGATGGATTAACAGTAAAGGAGTCAATGATGTCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGtt                                                                                                                   >  1:1062395/1‑150 (MQ=255)
                                  cAGCGGAAGAGATGGATTAACAGTAAAGGAGTCAATGATGTCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAATAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGttt                                                                                                                  <  1:1220192/150‑1 (MQ=255)
                                  cAGCGGAAGAGATGGATTAACAGTAAAGGAGTCAATGATGTCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGttt                                                                                                                  >  1:4838745/1‑150 (MQ=255)
                                      ggAAGAGATGGATTAACAGTAAAGGAGTCAATGATGTCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTttatta                                                                                                              >  1:2103485/1‑150 (MQ=255)
                                       gAAGAGATGGATTAACAGTAAAGGAGTCAATGATGTCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAATAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTAc                                                                                                             >  2:4163916/1‑150 (MQ=255)
                                        aaGAGATGGATTAACAGTAAAGGAGTCAATGATGTCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACt                                                                                                            >  2:4335209/1‑150 (MQ=255)
                                          gagaTGGATTAACAGTAAAGGAGTCAATGATGTCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAATAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCa                                                                                                          <  1:39144/150‑1 (MQ=255)
                                          gagaTGGATTAACAGTAAAGGAGTCAATGATGTCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCa                                                                                                          >  2:4678251/1‑150 (MQ=255)
                                            gaTGGATTAACAGTAAAGGAGTCAATGATGTCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCAtt                                                                                                        <  2:2921668/150‑1 (MQ=255)
                                             aTGGATTAACAGTAAAGGAGTCAATGATGTCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTa                                                                                                       >  1:321817/1‑150 (MQ=255)
                                                   tAACAGTAAAGGAGTCAATGATGTCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAgg                                                                                                 <  2:256329/150‑1 (MQ=255)
                                                    aaCAGTAAAGGAGTCAATGATGTCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAATAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAggg                                                                                                >  1:1475583/1‑150 (MQ=255)
                                                    aaCAGTAAAGGAGTCAATGATGTCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAggg                                                                                                <  2:468288/150‑1 (MQ=255)
                                                    aaCAGTAAAGGAGTCAATGATGTCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAggg                                                                                                <  2:1855864/150‑1 (MQ=255)
                                                      cAGTAAAGGAGTCAATGATGTCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAATAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGaa                                                                                              <  2:5138184/150‑1 (MQ=255)
                                                      cAGTAAAGGAGTCAATGATGTCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAATAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGaa                                                                                              >  2:1726678/1‑150 (MQ=255)
                                                         tAAAGGAGTCAATGATGTCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAg                                                                                           >  2:1861757/1‑150 (MQ=255)
                                                               aGTCAATGATGTCAGGATAACACACGTTAAAAGCGGTACGAGGCATTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGTTTTATTACTCATTAAACAGGGAAAAGTGGAAt                                                                                     <  2:859147/150‑1 (MQ=255)
                                                                gTCAATGATGTCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCCTTAAACAGGGACAAGTGGAATg                                                                                    >  2:1147726/1‑150 (MQ=255)
                                                                gTCAATGATGTCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATg                                                                                    >  2:387343/1‑150 (MQ=255)
                                                                gTCAATGAGGTCAGGAGAACACACGTTACAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATg                                                                                    <  2:4538053/150‑1 (MQ=255)
                                                                 tCAATGATGTCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATgg                                                                                   >  2:3087571/1‑150 (MQ=255)
                                                                  cAATGATGTCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAATAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGt                                                                                  <  1:3521260/150‑1 (MQ=255)
                                                                  cAATGATGTCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGt                                                                                  >  2:2742476/1‑150 (MQ=255)
                                                                  cAATGATGTCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGt                                                                                  >  1:4896889/1‑150 (MQ=255)
                                                                   aaTGATGTCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAATAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGtt                                                                                 >  1:310340/1‑150 (MQ=255)
                                                                   aaTGATGTCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGtt                                                                                 >  1:3599208/1‑150 (MQ=255)
                                                                    atgatgTCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGttt                                                                                <  1:1412673/150‑1 (MQ=255)
                                                                     tgatgTCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAATAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTg                                                                               >  1:4721596/1‑150 (MQ=255)
                                                                     tgatgTCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAATAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTg                                                                               <  1:4538095/150‑1 (MQ=255)
                                                                       atgTCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAATAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGc                                                                             <  1:150007/150‑1 (MQ=255)
                                                                         gTCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGa                                                                           <  1:1701578/150‑1 (MQ=255)
                                                                            aGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATg                                                                        <  1:4948075/150‑1 (MQ=255)
                                                                            aGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATg                                                                        <  2:2403514/150‑1 (MQ=255)
                                                                               agaACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGa                                                                     <  1:3203864/150‑1 (MQ=255)
                                                                                 aaCACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGaaa                                                                   <  2:1603803/150‑1 (MQ=255)
                                                                                        gTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGaaa                                                            >  2:779141/1‑150 (MQ=255)
                                                                                          tAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAATAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGc                                                          <  2:331368/150‑1 (MQ=255)
                                                                                            aaaGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCAt                                                        >  1:1005359/1‑150 (MQ=255)
                                                                                            aaaGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCAt                                                        >  1:30567/1‑150 (MQ=255)
                                                                                             aaGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATc                                                       <  2:206527/150‑1 (MQ=255)
                                                                                              aGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCa                                                      <  2:1208352/150‑1 (MQ=255)
                                                                                                cGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAATAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCaaa                                                    <  2:3607611/150‑1 (MQ=255)
                                                                                                cGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCaaa                                                    <  2:3764634/150‑1 (MQ=255)
                                                                                                 ggTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAATAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGACAACGGAAAGCATCAAAt                                                   >  2:4842128/1‑150 (MQ=255)
                                                                                                     cGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAATAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCt                                               <  1:2021688/150‑1 (MQ=255)
                                                                                                      gAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCTg                                              >  1:4360931/1‑150 (MQ=255)
                                                                                                        ggCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCTGac                                            >  1:380605/1‑150 (MQ=255)
                                                                                                           aGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCTGACAcc                                         <  1:2170579/150‑1 (MQ=255)
                                                                                                             ttttATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCTGACACCAt                                       <  2:2958370/150‑1 (MQ=255)
                                                                                                               ttATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCTGACACCATTc                                     <  1:3461737/150‑1 (MQ=255)
                                                                                                                 aTTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCTGACACCATTcgc                                   >  1:4061747/1‑150 (MQ=255)
                                                                                                                 aTTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCTGACACCATTcgc                                   >  2:1114925/1‑150 (MQ=255)
                                                                                                                  ttGATGTCACCCGTACGATCGATAACCCGGAATAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCTGACACCATTcgcg                                  <  2:793162/150‑1 (MQ=255)
                                                                                                                   tGATGTCACCCGTACGATCGATAACCCGGAATAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCTGACACCATTCGCGt                                 >  1:719556/1‑150 (MQ=255)
                                                                                                                    gATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCTGACACCATTCGCGTg                                <  1:2260529/150‑1 (MQ=255)
                                                                                                                     aTGTCACCCGTACGATCGATAACCCGGAATAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCTGACACCATTCGCGTgc                               >  2:2523252/1‑150 (MQ=255)
                                                                                                                      tGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCTGACACCATTCGCGTgcg                              <  1:2784609/150‑1 (MQ=255)
                                                                                                                         cACCCGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCTGACACCATTCGCGTGCGCGa                           >  1:502353/1‑150 (MQ=255)
                                                                                                                               tACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCTGACACCATTCGCGTGCGCGACTATcg                     >  2:3738806/1‑150 (MQ=255)
                                                                                                                                    tcgatAACCCGGAATAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCTGACACCATTCGCGTGCGCGACTATCGCGGCa                <  2:2121621/150‑1 (MQ=255)
                                                                                                                                    tcgatAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCTGACACCATTCGCGTGCGCGACTATCGCGGCa                >  2:1532568/1‑150 (MQ=255)
                                                                                                                                     cgatAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCTGACACCATTCGCGTGCGCGACTATCGCGGCaa               <  2:2188000/150‑1 (MQ=255)
                                                                                                                                     cgatAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCTGACACCATTCGCGTGCGCGACTATCGCGGCaa               <  1:877538/150‑1 (MQ=255)
                                                                                                                                        tAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCTGACACCATTCGCGTGCGCGACTATCGCGGCAAACt            >  2:5067802/1‑150 (MQ=255)
                                                                                                                                        tAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCTGACACCATTCGCGTGCGCGACTATCGCGGCAAACt            <  1:1169972/150‑1 (MQ=255)
                                                                                                                                         aaCCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCTGACACCATTCGCGTGCGCGACTATCGCGGCAAACTg           >  2:3978521/1‑150 (MQ=255)
                                                                                                                                         aaCCCGGAAGAGAGTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCTGACACCATTCGCGTGCGCGACTATCGCGGCAAACTg           <  1:192419/150‑1 (MQ=255)
                                                                                                                                          aCCCGGAATAGAGTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCTGACACCATTCGCGTGCGCGACTATCGCGGCAAACTGa          <  1:4347089/150‑1 (MQ=255)
                                                                                                                                          aCCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCTGACACCATTCGCGTGCGCGACTATCGCGGCAAACTGa          <  1:4678251/150‑1 (MQ=255)
                                                                                                                                           cccGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCTGACACCATTCGCGTGCGCGACTATCGCGGCAAACTGAt         <  2:1819752/150‑1 (MQ=255)
                                                                                                                                            ccGGAATAGAGTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCTGACACCATTCGCGTGCGCGACTATCGCGGCAAACTGATa        <  2:4423814/150‑1 (MQ=255)
                                                                                                                                             cGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCTGACACCATTCGCGTGCGCGACTATCGCGGCAAACTGATAg       <  2:4835677/150‑1 (MQ=255)
                                                                                                                                             cGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCTGACACCATTCGCGTGCGCGACTATCGCGGCAAACTGATAg       <  2:2152755/150‑1 (MQ=255)
                                                                                                                                              ggAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCTGACACCATTCGCATGCGCGACTATCGCGGCAAACTGATAGt      <  1:240905/150‑1 (MQ=255)
                                                                                                                                               gAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCTGACACCATTCGCGTGCGCGACTATCGCGGCAAACTGATAGTa     <  2:4751068/150‑1 (MQ=255)
                                                                                                                                                 agagATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCTGACACCATTCGCGTGCGCGACTATCGCGGCAAACTGATAGTAcc   <  1:4568385/150‑1 (MQ=255)
                                                                                                                                                 agagATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCTGACACCATTCGCGTGCGCGACTATCGCGGCAAACTGATAGTAcc   <  2:3540956/150‑1 (MQ=255)
                                                                                                                                                 agagATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCTGACACCATTCGCGTGCGCGACTATCGCGGCAAACTGATAGTAcc   <  2:2322828/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
TACCGACAGGAAAACGTTCTGCCTGGTATTACCTCAGCGGAAGAGATGGATTAACAGTAAAGGAGTCAATGATGTCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCTGACACCATTCGCGTGCGCGACTATCGCGGCAAACTGATAGTACCG  >  CP009273/3019057‑3019352

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: