Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,644,550 C→A 43.2% pseudogene (216/765 nt) ydfE → Qin prophage; pseudogene;Phage or Prophage Related

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,644,5500CA43.2% 203.6 / 251.6 222pseudogene (216/765 nt)ydfEQin prophage; pseudogene;Phage or Prophage Related
Reads supporting (aligned to +/- strand):  ref base C (58/68);  new base A (37/59);  total (95/127)
Fisher's exact test for biased strand distribution p-value = 2.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 222 total aligned reads displayed.

GTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTCGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGATGAGGAGTCTGCTTCTTTTGTCCCGGTTGAGACTGAATCCGATCCGATGAACGTCACTTTTGACAAGCTGGCCCCTGAAGTACAGAACGCTGTCATGGTTAAGTTC  >  CP009273/1644404‑1644697
                                                                                                                                                  |                                                                                                                                                   
    ccACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTCGATTCa                                                                                                                                               >  1:2930022/1‑149 (MQ=255)
    ccACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTCGATTCa                                                                                                                                               <  1:3608977/149‑1 (MQ=255)
     cacaGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTCGATTCAt                                                                                                                                              >  2:24535/1‑149 (MQ=255)
     cacaGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTAGATTCAt                                                                                                                                              >  1:654518/1‑149 (MQ=255)
        aGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTAGATTCATggg                                                                                                                                           >  2:3866570/1‑149 (MQ=255)
         gCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTCGATTCATGGGa                                                                                                                                          <  1:808678/149‑1 (MQ=255)
                 aGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAATTTTATTGTCTGAGAGGATAAACCTGGCATACCCCGCCCTGCCCTCGATTCATGGGATGCTGAat                                                                                                                                  >  1:92004/1‑149 (MQ=255)
                 aGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTCGATTCATGGGATGCTGAat                                                                                                                                  >  1:593587/1‑149 (MQ=255)
                 aGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGAGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTCGATTCATGGGATGCTGAat                                                                                                                                  <  1:814385/149‑1 (MQ=255)
                  gCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCACTGCCCTAGATTCATGGGATGCTGAata                                                                                                                                 <  1:1505076/149‑1 (MQ=255)
                   cTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCTGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTAGATTCATGGGATGCTGAatat                                                                                                                                <  1:4544/149‑1 (MQ=255)
                     ggCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTAGATTCATGGGATGCTGAatatat                                                                                                                              <  1:1018585/149‑1 (MQ=255)
                       cgacgaACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTAGATTCATGGGATGCTGAATATATGc                                                                                                                            >  1:414010/1‑149 (MQ=255)
                        gacgaACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTAGATTCATGGGATGCTGAATATATGCa                                                                                                                           >  1:3182019/1‑149 (MQ=255)
                        gacgaACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTAGATTCATGGGATGCTGAATATATGCa                                                                                                                           <  2:872125/149‑1 (MQ=255)
                         acgaACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTCGATTCATGGGATGCTGAATATATGCAg                                                                                                                          >  1:154813/1‑149 (MQ=255)
                         acgaACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTCGATTCATGGGATGCTGAATATATGCAg                                                                                                                          >  1:3664345/1‑149 (MQ=255)
                           gaACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTCGATTCATGGGATGCTGAATATATGCAGGa                                                                                                                        >  1:2203182/1‑149 (MQ=255)
                            aaCGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTAGATTCATGGGATGCTGAATATATGCAGGaa                                                                                                                       >  2:1844533/1‑149 (MQ=255)
                              cGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTAGATTCATGGGATGCTGAATATATGCAGGaaaa                                                                                                                     <  1:2738519/149‑1 (MQ=255)
                               gCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTCGATTCATGGGATGCTGAATATATGCAGGAAAAc                                                                                                                    <  1:122644/149‑1 (MQ=255)
                               gCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTCGATTCATGGGATGCTGAATATATGCAGGAAAAc                                                                                                                    <  1:3856388/149‑1 (MQ=255)
                                ccGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTCGATTCATGGGATGCTGAATATATGCAGGAAAAcc                                                                                                                   <  1:2321238/149‑1 (MQ=255)
                                 cGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTCGATTCATGGGATGCTGAATATATGCAGGAAAACCg                                                                                                                  >  1:3976749/1‑149 (MQ=255)
                                  gTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTAGATTCATGGGATGCTGAATATATGCAGGAAAACCGc                                                                                                                 <  2:1301954/149‑1 (MQ=255)
                                       cGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTCGATTCATGGGATGCTGACTATATGCAGGAAAACCGCTGGGa                                                                                                            >  2:2497767/1‑149 (MQ=255)
                                       cGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTCGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGa                                                                                                            <  2:3078811/149‑1 (MQ=255)
                                       cGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTCGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGa                                                                                                            >  2:489908/1‑149 (MQ=255)
                                        gAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTAGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGAt                                                                                                           >  2:1596005/1‑149 (MQ=255)
                                        gAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTAGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGAt                                                                                                           >  1:1008477/1‑149 (MQ=255)
                                        gAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTAGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGAt                                                                                                           <  2:3093725/149‑1 (MQ=255)
                                            aaaaTTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTCGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGATgagg                                                                                                       >  2:2060743/1‑149 (MQ=255)
                                              caTTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTCGAGTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGATgaggag                                                                                                     <  2:4082758/148‑1 (MQ=255)
                                              aaTTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTAGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGATgaggag                                                                                                     >  2:3757936/1‑149 (MQ=255)
                                               aTTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTCGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGATGAGGAGt                                                                                                    >  1:1274505/1‑149 (MQ=255)
                                                tttCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTAGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGATGAGGAGTc                                                                                                   <  2:1661757/149‑1 (MQ=255)
                                                   cACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTCGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGATGAGGAGTCTGc                                                                                                >  1:499610/1‑149 (MQ=255)
                                                     cTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTAGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGATGAGGAGTCTGctt                                                                                              <  1:3601935/149‑1 (MQ=255)
                                                       ggCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTCGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGATGAGGAGTCTGcttct                                                                                            <  2:3210136/149‑1 (MQ=255)
                                                       ggCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTAGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGATGAGGAGTCTGcttct                                                                                            <  1:3321220/149‑1 (MQ=255)
                                                       ggCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTAGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGATGAGGAGTCTGcttct                                                                                            <  1:96969/149‑1 (MQ=255)
                                                            ttCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTCGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGATGAGGAGTCTGCTTCTTTTGt                                                                                       >  2:3902408/1‑149 (MQ=255)
                                                            ttCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTAGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGATGAGGAGTCTGCTTCTTTTGt                                                                                       >  2:3700790/1‑149 (MQ=255)
                                                            ttCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTAGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGATGAGGAGTCTGCTTCTTTTGt                                                                                       >  2:1105965/1‑149 (MQ=255)
                                                              ccTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTCGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGATGAGGAGTCTGCTTCTTTTGTcc                                                                                     <  2:3192843/149‑1 (MQ=255)
                                                              ccTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTAGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGATGAGGAGTCTGCTTCTTTTGTcc                                                                                     <  2:3481539/149‑1 (MQ=255)
                                                               cTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTCGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGATGAGGAGTCTGCTTCTTTTGTccc                                                                                    <  1:3177922/149‑1 (MQ=255)
                                                               cTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTCGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGATGAGGAGTCTGCTTCTTTTGTccc                                                                                    <  1:2043194/149‑1 (MQ=255)
                                                               cTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTCGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGATGAGGAGTCTGCTTCTTTTGTccc                                                                                    <  2:486968/149‑1 (MQ=255)
                                                                 ggAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTCGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGATGAGGAGTCTGCTTCTTTTGTCCCgg                                                                                  <  2:2843271/149‑1 (MQ=255)
                                                                  gAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTCGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGATGAGGAGTCTGCTTCTTTTGTCCCGGt                                                                                 <  2:1186325/149‑1 (MQ=255)
                                                                  gAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTAGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGATGAGGAGTCTGCTTCTTTTGTCCCGGt                                                                                 <  1:1934416/149‑1 (MQ=255)
                                                                   aaCATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTAGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGATGAGGAGTCTGCTTCTTTTGTCCCGGtt                                                                                <  2:601305/149‑1 (MQ=255)
                                                                   aaCATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTAGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGATGAGGAGTCTGCTTCTTTTGTCCCGGtt                                                                                <  2:601754/149‑1 (MQ=255)
                                                                   aaCATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTAGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGATGAGGAGTCTGCTTCTTTTGTCCCGGtt                                                                                <  2:1996589/149‑1 (MQ=255)
                                                                    aCATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTAGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGATGAGGAGTCTGCTTCTTTTGTCCCGGTTg                                                                               >  1:1250809/1‑149 (MQ=255)
                                                                     cATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTCGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGATGAGGAGTCTGCTTCTTTTGTCCCGGTTga                                                                              <  2:2538304/149‑1 (MQ=255)
                                                                     cATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTCGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGATGAGGAGTCTGCTTCTTTTGTCCCGGTTga                                                                              <  2:614361/149‑1 (MQ=255)
                                                                        tATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTCGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGATGAGGAGTCTGCTTCTTTTGTCCCGGTTGAGAc                                                                           <  2:909953/149‑1 (MQ=255)
                                                                        tATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTAGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGATGAGGAGTCTGCTTCTTTTGTCCCGGTTGAGAc                                                                           <  2:1033532/149‑1 (MQ=255)
                                                                         aTCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTAGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGATGAGGAGTCTGCTTCTTTTGTCCCGGTTGAGACt                                                                          <  1:2158939/149‑1 (MQ=255)
                                                                          tCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTCGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGATGAGGAGTCTGCTTCTTTTGTCCCGGTTGAGACTg                                                                         >  2:3843572/1‑149 (MQ=255)
                                                                          tCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTAGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGATGAGGAGTCTGCTTCTTTTGTCCCGGTTGAGACTg                                                                         >  1:3084576/1‑149 (MQ=255)
                                                                              gctgctCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTAGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGATGAGGAGTCTGCTTCTTTTGTCCCGGTTGAGACTGAatc                                                                     >  2:828073/1‑149 (MQ=255)
                                                                                      ggACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTCGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGATGAGGAGTCTGCTTCTTTTGTCCCGGTTGAGACTGAAtccgatccga                                                             <  1:72131/149‑1 (MQ=255)
                                                                                      ggACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTAGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGATGAGGAGTCTGCTTCTTTTGTCCCGGTTGAGACTGAAtccgatccga                                                             <  2:3275557/149‑1 (MQ=255)
                                                                                      ggACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTAGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGATGAGGAGTCTGCTTCTTTTGTCCCGGTTGAGACTGAAtccgatccga                                                             <  1:3692314/149‑1 (MQ=255)
                                                                                        aCTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTCGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGATGAGGAGTCTGCTTCTTTTGTCCCGGTTGAGACTGAATCCGATCCGATg                                                           <  1:1201723/149‑1 (MQ=255)
                                                                                         cTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTCGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGATGAGGAGTCTGCTTCTTTTGTCCCGGTTGAGACTGAATCCGATCCGATGa                                                          <  2:2394025/149‑1 (MQ=255)
                                                                                         cTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTCGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGATGAGGAGTCTGCTTCTTTTGTCCCGGTTGAGACTGAATCCGATCCGATGa                                                          <  1:1738431/149‑1 (MQ=255)
                                                                                             gcTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTAGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGATGAGGAGTCTGCTTCTTTTGTCCCGGTTGAGACTGAATCCGATCCGATGAACGt                                                      >  2:1129805/1‑149 (MQ=255)
                                                                                               ttATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTCGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGATGAGGAGTCTGCTTCTTTTGTCCCGGTTGAGACTGAATCCGATCCGATGAACGTCa                                                    >  1:2427517/1‑149 (MQ=255)
                                                                                               ttATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTAGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGATGAGGAGTCTGCTTCTTTTGTCCCGGTTGAGACTGAATCCGATCCGATGAACGTCa                                                    <  2:442771/149‑1 (MQ=255)
                                                                                                tataAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTCGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGATGAGGAGTCTGCTTCTTTTGTCCCGGTTGAGACTGAATCCGATCCGATGAACGTCAc                                                   >  1:1058693/1‑149 (MQ=255)
                                                                                                   aaaTTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTAGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGATGAGGAGTCTGCTTCTTTTGTCCCGGTTGAGACTGAATCCGATCCGATGAACGTCACttt                                                <  1:2797953/149‑1 (MQ=255)
                                                                                                     aTTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTCGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGATGAGGAGTCTGCTTCTTTTGTCCCGGTTGAGACTGAATCCGATCCGATGAACGTCACTTTTg                                              >  2:2087072/1‑149 (MQ=255)
                                                                                                        tATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTAGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGATGAGGAGTCTGCTTCTTTTGTCCCGGTTGAGACTGAATCCGATCCGATGAACGTCACTTTTGACa                                           <  1:2562259/149‑1 (MQ=255)
                                                                                                         aTTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTCGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGATGAGGAGTCTGCTTCTTTTGTCCCGGTTGAGACTGAATCCGATCCGATGAACGTCACTTTTGACaa                                          >  2:1390879/1‑149 (MQ=255)
                                                                                                            gTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTAGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGATGAGGAGTCTGCTTCTTTTGTCCCGGTTGAGACTGAATCCGATCCGATGAACGTCACTTTTGACAAGCt                                       <  1:2947798/149‑1 (MQ=255)
                                                                                                             tCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTCGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGATGAGGAGTCTGCTTCTTTTGTCCCGGTTGAGACTGAATCCGATCCGATGAACGTCACTTTTGACAAGCTg                                      >  2:757179/1‑149 (MQ=255)
                                                                                                                    ggATAAACCTGGCATACCCCGCCCTGCCCTCGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGATGAGGAGTCTGCTTCTTTTGTCCCGGTTGAGACTGAATCCGATCCGATGAACGTCACTTTTGACAAGCTGGCCCCTg                               >  1:1909361/1‑149 (MQ=255)
                                                                                                                      aTAAACCTGGCATACCCCGCCCTGCCCTAGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGATGAGGAGTCTGCTTCTTTTGTCCCGGTTGAGACTGAATCCGATCCGATGAACGTCACTTTTGACAAGCTGGCCCCTGaa                             >  1:2710879/1‑149 (MQ=255)
                                                                                                                       tAAACCTGGCATACCCCGCCCTGCCCTAGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGATGAGGAGTCTGCTTCTTTTGTCCCGGTTGAGACTGAATCCGATCCGATGAACGTCACTTTTGACAAGCTGGCCCCTGAAg                            <  2:1254992/149‑1 (MQ=255)
                                                                                                                        aaaCCTGGCATACCCCGCCCTGCCCTCGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGATGAGGAGTCTGCTTCTTTTGTCCCGGTTGAGACTGAATCCGATCCGATGAACGTCACTTTTGACAAGCTGGCCCCTGAAGt                           >  2:3534013/1‑149 (MQ=255)
                                                                                                                         aaCCTGGCATACCCCGCCCTGCCCTCGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGATGAGGAGTCTGCTTCTTTTGTCCCGGTTGAGACTGAATCCGATCCGATGAACGTCACTTTTGACAAGCTGGCCCCTGAAGTa                          <  1:3902408/149‑1 (MQ=255)
                                                                                                                         aaCCTGGCATACCCCGCCCTGCCCTAGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGATGAGGAGTCTGCTTCTTTTGTCCCGGTTGAGACTGAATCCGATCCGATGAACGTCACTTTTGACAAGCTGGCCCCTGAAGTa                          >  2:3709017/1‑149 (MQ=255)
                                                                                                                             tGGCATACCCCGCCCTGCCCTCGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGATGAGGAGTCTGCTTCTTTTGTCCCGGTTGAGACTGAATCCGATCCGATGAACGTCACTTTTGACAAGCTGGCCCCTGAAGTACAGa                      >  2:4044263/1‑149 (MQ=255)
                                                                                                                              ggCATACCGCGCCCTGCCCTCGTTTCATGGGATGCTGAATATAGGCTGGATAACCGGTGGGATGAGGAGTCTGCTTCTTTTGTCCCGGTTGAGACTGTATCCGTTCCTATGAACGTCACTTTTGACAAGCTGGGCCCTGTAGTACAGaa                     <  2:3333146/149‑1 (MQ=255)
                                                                                                                              ggCATACCCCGCCCTGCTCTCGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGATGAGGAGTCTGCTTCTTTTGTCCCGGTTGAGACTGAATCCGATCCGATGAACGTCACTTTTGACAAGCTGGCCCCTGAAGTACAGaa                     <  2:1282281/149‑1 (MQ=255)
                                                                                                                              ggCATACCCCGCCCTGCCCTCGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGATGAGGAGTCTGCTTCTTTTGTCCCGGTTGAGACTGAATCCGATCCGATGAACGTCACTTTTGACAAGCTGGCCCCTGAAGTACAGaa                     <  2:578104/149‑1 (MQ=255)
                                                                                                                              ggCATACCCCGCCCTGCCCTCGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGATGAGGAGTCTGCTTCTTTTGTCCCGGTTGAGACTGAATCCGATCCGATGAACGTCACTTTTGACAAGCTGGCCCCTGAAGTACAGaa                     <  2:2921577/149‑1 (MQ=255)
                                                                                                                              ggCATACCCCGCCCTGCCCTAGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGATGAGGAGTCTGCTTCTTTTGTCCCGGTTGAGACTGAATCCGATCCGATGAACGTCACTTTTGACAAGCTGGCCCCTGAAGTACAGaa                     >  1:198624/1‑149 (MQ=255)
                                                                                                                                 aTACCCCGCCCTGCCCTCGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGATGAGGAGTCTGCTTCTTTTGTCCCGGTTGAGACTGAATCCGATCCGATGAACGTCACTTTTGACAAGCTGGCCCCTGAAGTACAGAACGc                  >  2:633362/1‑149 (MQ=255)
                                                                                                                                 aTACCCCGCCCTGCCCTCGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGATGAGGAGTCTGCTTCTTTTGTCCCGGTTGAGACTGAATCCGATCCGATGAACGTCACTTTTGACAAGCTGGCCCCTGAAGTACAGAACGc                  <  2:2429197/149‑1 (MQ=255)
                                                                                                                                 aTACCCCGCCCTGCCCTCGAGTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGATGAGGAGTCTGCTTCTTTTGTCCCGGTTGAGACTGAATCCGATCCGATGAACGTCACTTTTGACAAGCTGGCCCCTGAAGTACAGAACGc                  <  2:2428734/149‑1 (MQ=255)
                                                                                                                                         ccctgccctAGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGATGAGGAGTCTGCTTCTTTTGTCCCGGTTGAGACTGAATCCGATCCGATGAACGTCACTTTTGACAAGCTGGCCCCTGAAGTACAGAACGCTGTCATgg          <  1:99907/149‑1 (MQ=255)
                                                                                                                                           ctgccctCGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGATGAGGAGTCTGCTTCTTTTGTCCCGGTTGAGACTGAATCCGATCCGATGAACGTCACTTTTGACAAGCTGGCCCCTGAAGTACAGAACGCTGTCATGGtt        >  1:1191783/1‑149 (MQ=255)
                                                                                                                                               cctCGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGATGAGGAGTCTGCTTCTTTTGTCCCGGTTGAGACTGAATCCGATCCGATGAACGTCACTTTTGACAAGCTGGCCCCTGAAGTACAGAACGCTGTCATGGTTAAGt    <  1:2041183/149‑1 (MQ=255)
                                                                                                                                               cctCGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGATGAGGAGTCTGCTTCTTTTGTCCCGGTTGAGACTGAATCCGATCCGATGAACGTCACTTTTGACAAGCTGGCCCCTGAAGTACAGAACGCTGTCATGGTTAAGt    >  2:1556753/1‑149 (MQ=255)
                                                                                                                                                ctCGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGATGAGGAGTCTGCTTCTTTTGTCCCGGTTGAGACTGAATCCGATCCGATGAACGTCACTTTTGACAAGCTGGCCCCTGAAGTACAGAACGCTGTCATGGTTAAGtt   <  2:4016302/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
GTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTCGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGATGAGGAGTCTGCTTCTTTTGTCCCGGTTGAGACTGAATCCGATCCGATGAACGTCACTTTTGACAAGCTGGCCCCTGAAGTACAGAACGCTGTCATGGTTAAGTTC  >  CP009273/1644404‑1644697

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: