Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 554,715 G→T 5.0% G96C (GGT→TGT)  sfmC → putative periplasmic pilus chaperone

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273554,7150GT 5.0% 1080.1 / 21.8 277G96C (GGT→TGT) sfmCputative periplasmic pilus chaperone
Reads supporting (aligned to +/- strand):  ref base G (139/124);  new base T (6/8);  total (145/132)
Fisher's exact test for biased strand distribution p-value = 5.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.69e-01

Only 100 of 280 total aligned reads displayed.

GATTAGAAATAGCCATACCAATGAGCGCTTTCTGGTCAATTCGTGGATTGAAAACAGCAGCGGTGTAAAAGAAAAGTCATTCATCATTACACCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTGGCAGCAGATCGTGAGTCTCTGTTCTGGATGAATGTTAAGACGATCCCTTCGGTAGATAAAAATGCATTGAACGGCAGGAATGTTTTGCAACTGGCGATTTTATCGCGCATGAAATTATTTCTCCGTCCAATTCAA  >  CP009273/554567‑554861
                                                                                                                                                    |                                                                                                                                                  
 aTTAGAAATAGCCATACCAATGAGCGCTTTCTGGTCAATTCGTGGATTGAAAACAGCAGCGGTGTAAAAGAAAAGTCATTCATCATTACACCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGt                                                                                                                                                  >  1:4628419/1‑150 (MQ=255)
   tAGAAATAGCCATACCAATGAGCGCTTTCTGGTCAATTCGTGGATTGAAAACAGCAGCGGTGTAAAAGAAAAGTCATTCATCATTACACCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTcc                                                                                                                                                <  2:609625/150‑1 (MQ=255)
     gAAATAGCCATACCAATGAGCGCTTTCTGGTCAATTCGTGGATTGAAAACAGCAGCGGTGTAAAAGAAAAGTCATTCATCATTACACCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCTGTCCAc                                                                                                                                              <  2:1761644/150‑1 (MQ=255)
      aaaTAGCCATACCAATGAGCGCTTTCTGGTCAATTCGTGGATTGAAAACAGCAGCGGTGTAAAAGAAAAGTCATTCATCATTACACCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCAcc                                                                                                                                             <  2:2841180/150‑1 (MQ=255)
       aaTAGCCATACCAAGGAGCGCTTTCTGGTCAATTCGTGGATTGAAAACAGCAGCGGTGTAAAAGAAAAGTCATTCATCATTACACCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCg                                                                                                                                            >  1:892017/1‑150 (MQ=255)
         tAGCCATACCAATGAGCGCTTTCTGGTCAATTCGTGGATTGAAAACAGCAGCGGTGTACAAGAAAAGTCATTCATCATTACACCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCt                                                                                                                                          >  2:3670427/1‑150 (MQ=255)
            ccATACCAATGAGCGCTTTCTGGTCAATTCGTGGATTGAAAACAGCAGCGGTGTAAAAGAAAAGTCATTCATCATTACACCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTGgc                                                                                                                                       <  1:2095587/150‑1 (MQ=255)
             cATACCAANGAGCGCTNTCTGGTCAATTCGTGGATTGAAAACAGCAGCGGTGTAAAAGAAAAGTCATTCATCATTACACCGCCACTGTTTGTTAGTGAACCCAAAAGCNAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTGgca                                                                                                                                      >  2:2716158/1‑150 (MQ=255)
              aTACCAATGAGCGCTTTCTGGTCAATTCGTGGATTGAAAACAGCAGCGGTGTAAAAGAAAAGTCATTCATCATTACACCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTGgcag                                                                                                                                     <  1:3084287/150‑1 (MQ=255)
                 ccAATGAGCGCTTTCTGGTCAATTCGTGGATTGAAAACAGCAGCGGTGTAAAAGAAAAGTCATTCATCATTACACCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTGGcagcag                                                                                                                                  >  1:3243103/1‑150 (MQ=255)
                  cAATGAGCGCTTTCTGGTCAATTCGTGGATTGAAAACAGCAGCGGTGTAAAAGAAAAGTCATTCATCATTACACCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTGGCAGCAGa                                                                                                                                 >  1:3960941/1‑150 (MQ=255)
                   aaTGAGCGCTTTCTGGTCAATTCGTGGATTGAAAACAGCAGCGGTGTAAAAGAAAAGTCATTCATCATTACACCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTGGCAGCAGAt                                                                                                                                >  1:2025005/1‑150 (MQ=255)
                    aTGAGCGCTTTCTGGTCAATTCGTGGATTGAAAACAGCAGCGGTGTAAAAGAAAAGTCATTCATCATTACACCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTGGCAGCAGATc                                                                                                                               <  2:3474201/150‑1 (MQ=255)
                      gAGCGCTTTCTGGTCAATTCGTGGATTGAAAACAGCAGCGGTGTAAAAGAAAAGTCATTCATCATTACACCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTGGCAGCAGATCGt                                                                                                                             <  1:2462067/150‑1 (MQ=255)
                      gAGCGCTTTCTGGTCAATTCGTGGATTGAAAACAGCAGCGGTGTAAAAGAAAAGTCATTCATCATTACACCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTGGCAGCAGATCGt                                                                                                                             >  1:278683/1‑150 (MQ=255)
                        gcgcTTTCTGGTCAATTCGTGGATTGAAAACAGCAGCGGTGTAAAAGAAAAGTCATTCATCATTACACCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTGGCAGCAGATCGTGa                                                                                                                           <  1:2892429/150‑1 (MQ=255)
                            tttCTGGTCAATTCGTGGATTGAAAACAGCAGCGGTGTAAAAGAAAAGTCATTCATCATTACACCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTGGCAGCAGATCGTGAGtct                                                                                                                       >  2:1352984/1‑150 (MQ=255)
                            tttCTGGTCAATTCGTGGATTGAAAACAGCAGCGGTGTAAAAGAAAAGTCATTCATCATTACACCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTGGCAGCAGATCGTGAGtct                                                                                                                       >  2:2177794/1‑150 (MQ=255)
                             ttCTGGTCAATTCGTGGATTGAAAACAGCAGCGGTGTAAAAGAAAAGTCATTCATCATTACACCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTGGCAGCAGATCGTGAGtctc                                                                                                                      <  2:199745/150‑1 (MQ=255)
                              tCTGGTCAATTCGTGGATTGAATACAGCAGCGGTGTAAAAGAAAAGTCATTCATCATTACACCGCCAATGTTTGTTAGTGAACCCAAAAGCGAAAATACTTGGCGTATTATTTACACCGGTCCACCGCTGGCAGCAGATCGTGAGTctct                                                                                                                     <  2:198108/150‑1 (MQ=255)
                               cTGGTCAATTCGTGGATTGAAAACAGCAGCGGTGTAAAAGAAAAGTCATTCATCATTACACCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTGGCAGCAGATCGTGAGTCTCTg                                                                                                                    >  2:3110372/1‑150 (MQ=255)
                                tGGTCAATTCGTGGATTGAAAACAGCAGCGGTGTAAAAGAAAAGTCATTCATCATTACACCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCTGTCCACCGCTGGCAGCAGATCGTGAGTCTCTGt                                                                                                                   >  2:520228/1‑150 (MQ=255)
                                 ggTCAATTCGTGGATTGAAAACAGCAGCGGTGTAAAAGAAAAGTCATTCATCATTACACCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTGGCAGCAGATCGTGAGTCTCTGtt                                                                                                                  <  2:3564337/150‑1 (MQ=255)
                                 ggTCAATTCGTGGATTGAAAACAGCAGCGGTGTAAAAGAAAAGTCATTCATCATTACACCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTGGCAGCAGATCGTGAGTCTCTGtt                                                                                                                  <  2:2547562/150‑1 (MQ=255)
                                    cAATTCGTGGATTGAAAACAGCAGCGGTGTAAAATAAAAGTCATTCATCATTACACCGCCACTGTTTGTTAGTGAAACCAAAAGCGAAAAAACTTTGCGTATTATTTACACCGTTCCACCGCAGCCAGCAGATCGTTAGTCTCTGTTCTg                                                                                                               >  2:2369359/1‑150 (MQ=255)
                                      aTTCGTGGATTGAAAACAGCAGCGGTGTAAAAGAAAAGTAATTCATCATTACACCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTGGCAGCAGATCGTGAGTCTCTGTTCTGGa                                                                                                             >  1:2628511/1‑150 (MQ=255)
                                         cGTGGATTGAAAACAGCAGCGGTGTAAAAGAAAAGTCATTCATCATTACACCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTGGCAGCAGATCGTGAGTCTCTGTTCTGGATGa                                                                                                          >  2:297977/1‑150 (MQ=255)
                                           tGGATTGAAAACAGCAGCGGTGTAAAAGAAAAGTCATTCATCATTACACCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTGGCAGCAGATCGTGAGTCTCTGTTCTGGATGAAt                                                                                                        >  2:4534096/1‑150 (MQ=255)
                                           tGGATTGAAAACAGCAGCGGTGTAAAAGAAAAGTCATTCATCATTACACCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTGGCAGCAGATCGTGAGTCTCTGTTCTGGATGAAt                                                                                                        <  2:1665324/150‑1 (MQ=255)
                                            ctaTTGAAAACAGCAGCGGTTTAAATTAAAAGTCATTCATCATTACACCGCAACTGTTTGTTAGTGAACACAAAAGATAAAATACTTTGCTTATTATTTACACCGGTCCACAGCTGGAAGCAGATCGTGAGTCTCTGTTCTGTAGGAATg                                                                                                       >  2:1768837/3‑150 (MQ=255)
                                                tGAAAACAGCAGCGGTGTAAAAGAAAAGTCATTCATCATTACACCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTGGCAGCAGATCGTGAGTCTCTGTTCTGGATGAATGTTaa                                                                                                   <  1:1933145/150‑1 (MQ=255)
                                                  aaaaCAGCAGCGGTGTAAAAGAAAAGTCATTCATCATTACACCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTGGCAGCAGATCGTGAGTCTCTGTTCTGGATGAATGTTAAGa                                                                                                 <  1:3159865/150‑1 (MQ=255)
                                                     aCAGCAGCGGTGTAAAAGAAAAGTCATTCATCATTACACCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTGGCAGCAGATCGTGAGTCTCTGTTCTGGATGAATGTTAAGACGa                                                                                              >  2:779788/1‑150 (MQ=255)
                                                          agcGGTGTAAAAGAAAAGTCATTCATCATTACACCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTGGCAGCAGATCGTGAGTCTCTGTTCTGGATGAATGTTAAGACGATCCCt                                                                                         <  2:1308080/150‑1 (MQ=255)
                                                          agcGGTGTAAAAGAAAAGTCATTCATCATTACACCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTGGCAGCAGATCGTGAGTCTCTGTTCTGGATGAATGTTAAGACGATCCCt                                                                                         <  2:3051669/150‑1 (MQ=255)
                                                           gcGGTGTAAAAGAAAAGTCATTCATCATTACACCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTGGCAGCAGATCGTGAGTCTCTGTTCTGGATGAATGTTAAGACGATCCCtt                                                                                        <  2:401932/150‑1 (MQ=255)
                                                            cGGTGTAAAAGAAAAGTCATTCATCATTACACCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTGGCAGCAGATCGTGAGTCTCTGTTCTGGATGAATGTTAAGACGATCCCTTc                                                                                       <  2:123696/150‑1 (MQ=255)
                                                             ggTGTAAAAGAAAAGTCATTCATCATTACACCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTGGCAGCAGATCGTGAGTCTCTGTTCTGGATGAATGTTAAGACGATCCCTTCg                                                                                      >  1:696713/1‑150 (MQ=255)
                                                              gtgtAAAAGAAAAGTCATTCATCATTACACCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTGGCAGCAGATCGTGAGTCTCTGTTCTGGATGAATGTTAAGACGATCCCTTCgg                                                                                     >  2:4479807/1‑150 (MQ=255)
                                                                gtAAAAGAAAAGTCATTCATCATTACACCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTGGCAGCAGATCGTGAGTCTCTGTTCTGGATGAATGTTAAGACGATCCCTTCGGTa                                                                                   >  2:64632/1‑150 (MQ=255)
                                                                   aaagaaaagTCATTCATCATTACACCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTGGCAGCAGATCGTGAGTCTCTGTTCTGGATGAATGTTAAGACGATCCCTTCGGTAGAt                                                                                <  1:3952759/150‑1 (MQ=255)
                                                                    aagaaaagTCATTCATCATTACACCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTGGCAGCAGATCGTGAGTCTCTGTTCTGGATGAATGTTAAGACGATCCCTTCGGTAGATa                                                                               >  1:630520/1‑150 (MQ=255)
                                                                    aagaaaagTCATTCATCATTACACCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTGGCAGCAGATCGTGAGTCTCTGTTCTGGATGAATGTTAAGACGATCCCTTCGGTAGATa                                                                               >  1:3964746/1‑150 (MQ=255)
                                                                              attcatCATTACACCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTGGCAGCAGATCGTGAGTCTCTGTTCTGGATGAATGTTAAGACGATCCCTTCGGTAGATAAAAATGCAtt                                                                     <  1:2212941/150‑1 (MQ=255)
                                                                                 catcatTACACCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTGGCAGCAGATCGTGAGTCTCTGTTCTGGATGAATGTTAAGACGATCCCTTCGGTAGATAAAAATGCATTGaa                                                                  >  2:154956/1‑150 (MQ=255)
                                                                                    catTACACCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTGGCAGCAGATCGTGAGTCTCTGTTCTGGATGAATGTTAAGACGATCCCTTCGGTAGATAAAAATGCATTGAACgg                                                               >  2:1437846/1‑150 (MQ=255)
                                                                                      ttACACCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTGGCAGCAGATCGTGAGTCTCTGTTCTGGATGAATGTTAAGACGATCCCTTCGGTAGATAAAAATGCATTGAACGGCa                                                             >  2:704761/1‑150 (MQ=255)
                                                                                       tACACCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTGGCAGCAGATCGTGAGTCTCTGTTCTGGATGAATGTTAAGACGATCCCTTCGGTAGATAAAAATGCATTGAACGGCAg                                                            <  1:3010495/150‑1 (MQ=255)
                                                                                        acacCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTGGCAGCAGATCGTGAGTCTCTGTTCTGGATGAATGTTAAGACGATCCCTTCGGTAGATAAAAATGCATTGAACGGCAgg                                                           <  1:2152525/150‑1 (MQ=255)
                                                                                        acacCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTGGCAGCAGATCGTGAGTCTCTGTTCTGGATGAATGTTAAGACGATCCCTTCGGTAGATAAAAATGCATTGAACGGCAgg                                                           <  2:4178109/150‑1 (MQ=255)
                                                                                        acacCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTGGCAGCAGATCGTGAGTCTCTGTTCTGGATGAATGTTAAGACGATCCCTTCGGTAGATAAAAATGCATTGAACGGCAgg                                                           >  2:3361099/1‑150 (MQ=255)
                                                                                         cacCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTGGCAGCAGATCGTGAGTCTCTGTTCTGGATGACTGTTAAGACGATCCCTTCGGTAGATAAAAATGCATTGAACGGCAGGa                                                          >  1:3631929/1‑150 (MQ=255)
                                                                                         cacCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTGGCAGCAGATCGTGAGTCTCTGTTCTGGATGAATGTTAAGACGATCCCTTCGGTAGATAAAAATGCATTGAACGGCAGGa                                                          >  1:1064557/1‑150 (MQ=255)
                                                                                            cGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTGGCAGCAGATCGTGAGTCTCTGTTCTGGATGAATGTTAAGACGATCCCTTCGGTAGATAAAAATGCATTGAACGGCAGGAATg                                                       >  1:1097546/1‑150 (MQ=255)
                                                                                              ccACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTGGCAGCAGATCGTGAGTCTCTGTTCTGGATGAATGTTAAGACGATCCCTTCGGTAGATAAAAATGCATTGAACGGCAGGAATGtt                                                     <  1:1527147/150‑1 (MQ=255)
                                                                                               cACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTGGCAGCAGATCGTGAGTCTCTGTTCTGGATGAATGTTAAGACGATCCCTTCGGTAGATAAAAATGCATTGAACGGCAGGAATGttt                                                    <  2:1273955/150‑1 (MQ=255)
                                                                                               cACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTGGCAGCAGATCGTGAGTCTCTGTTCTGGATGAATGTTAAGACGATCCCTTCGGTAGATAAAAATGCATTGAACGGCAGGAATGttt                                                    <  1:3103580/150‑1 (MQ=255)
                                                                                                 cTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTGGCAGCAGATCGTGAGTCTCTGTTCTGGATGAATGTTAAGACGATCCCTTCGGTAGATAAAAATGCATTGAACGGCAGGAATGTTTTg                                                  >  1:2361398/1‑150 (MQ=255)
                                                                                                  tgtttgttAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTGGCAGCAGATCGTGAGTCTCTGTTCTGGATGAATGTTAAGACGATCCCTTCGGTAGATAAAAATGCATTGAACGGCAGGAATGTTTTGc                                                 <  1:2693712/150‑1 (MQ=255)
                                                                                                   gtttgttAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTGGCAGCAGATCGTGAGTCTCTGTTCTGGATGAATGTTAAGACGATCCCTTCGGTAGATAAAAATGCATTGAACGGCAGGAATGTTTTGCa                                                <  1:1014391/150‑1 (MQ=255)
                                                                                                        ttAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTGGCAGCAGATCGTGAGTCTCTGTTCTGGATGAATGTTAAGACGATCCCTTCGGTAGATAAAAATGCATTGAACGGCAGGAATGTTTTGCAACTgg                                           >  1:1537603/1‑150 (MQ=255)
                                                                                                         tAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTGGCAGCAGATCGTGAGTCTCTGTTCTGGATGAATGTTAAGACGATCCCTTCGGTAGATAAAAATGCATTGAACGGCAGGAATGTTTTGCAACTGGc                                          >  1:1874950/1‑150 (MQ=255)
                                                                                                          aGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTGGCAGCAGATCGTGAGTCTCTGTTCTGGATGAATGTTAAGACGATCCCTTCGGTAGATAAAAATGCATTGAACGGCAGGAATGTTTTGCAACTGGCg                                         >  2:3231646/1‑150 (MQ=255)
                                                                                                            tGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTGGCAGCAGATCGTGAGTCTCTGTTCTGGATGAATGTTAAGACGATCCCTTCGGTAGATAAAAATGCATTGAACGGCAGGAATGTTTTGCAACTGGCGAt                                       <  2:1187546/150‑1 (MQ=255)
                                                                                                             gAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTGGCAGCAGATCGTGAGTCTCTGTTCTGGATGAATGTTAAGACGATCCCTTCGGTAGATAAAAATGCATTGAACGGCAGGAATGTTTTGCAACTGGCGAtt                                      >  1:408338/1‑150 (MQ=255)
                                                                                                              aaCCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTGGCAGCAGATCGTGAGTCTCTGTTCTGGATGAATGTTAAGACGATCCCTTCGGTAGATAAAAATGCATTGAACGGCAGGAATGTTTTGCAACTGGCGAttt                                     >  1:4160755/1‑150 (MQ=255)
                                                                                                                cccAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTGGCAGCAGATCGTGAGTCTCTGTTCTGGATGAATGTTAAGACGATCCCTTCGGTAGATAAAAATGCATTGAACGGCAGGAATGTTTTGCAACTGGCGATTTTa                                   <  2:938259/150‑1 (MQ=255)
                                                                                                                cccAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTGGCAGCAGATCGTGAGTCTCTGTTCTGGATGAATGTTAAGACGATCCCTTCGGTAGATAAAAATGCATTGAACGGCAGGAATGTTTTGCAACTGGCGATTTTa                                   <  2:3963715/150‑1 (MQ=255)
                                                                                                                 ccAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTGGCAGCAGATCGTGAGTCTCTGTTCTGGATGAATGTTAAGACGATCCCTTCGGTAGATAAAAATGCATTGAACGGCAGGAATGTTTTGCAACTGGCGATTTTAt                                  >  2:2203140/1‑150 (MQ=255)
                                                                                                                  cAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTGGCAGCAGATCGTGAGTCTCTGTTCTGGATGAATGTTAAGACGATCCCTTCGGTAGATAAAAATGCATTGAACGGCAGGAATGTTTTGCAACTGGCGATTTTATc                                 <  1:4422066/150‑1 (MQ=255)
                                                                                                                  cAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTGGCAGCAGATCGTGAGTCTCTGTTCTGGATGAATGTTAAGACGATCCCTTCGGTAGATAAAAATGCATTGAACGGCAGGAATGTTTTGCAACTGGCGATTTTATc                                 <  1:3411867/150‑1 (MQ=255)
                                                                                                                    aaaGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTGGCAGCAGATCGTGAGTCTCTGTTCTGGATGAATGTTAAGACGATCCCTTCGGTAGATAAAAATGCATTGAACGGCAGGAATGTTTTGCAACTGGCGATTTTATcgc                               <  1:1345934/150‑1 (MQ=255)
                                                                                                                    aaaGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTGGCAGCAGATCGTGAGTCTCTGTTCTGGATGAATGTTAAGACGATCCCTTCGGTAGATAAAAATGCATTGAACGGCAGGAATGTTTTGCAACTGGCGATTTTATcgc                               >  1:3429463/1‑150 (MQ=255)
                                                                                                                     aaGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTGGCAGCAGATCGTGAGTCTCTGTTCTGGATGAATGTTAAGACGATCCCTTCGGTAGATAAAAATGCATTGAACGGCAGGAATGTTTTGCAACTGGCGATTTTATcgcg                              >  1:3520883/1‑150 (MQ=255)
                                                                                                                     aaGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTGGCAGCAGATCGTGAGTCTCTGTTCTGGATGAATGTTAAGACGATCCCTTCGGTAGATAAAAATGCATTGAACGGCAGGAATGTTTTGCAACTGGCGATTTTATcgcg                              <  1:1754301/150‑1 (MQ=255)
                                                                                                                     aaGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTGGCAGCAGATCGTGAGTCTCTGTTCTGGATGAATGTTAAGACGATCCCTTCGGTAGATAAAAATGCATTGAACGGCAGGAATGTTTTGCAACTGGCGATTTTATcgcg                              >  2:4111357/1‑150 (MQ=255)
                                                                                                                       gCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTGGCAGCAGATCGTGAGTCTCTGTTCTGGATGAATGTTAAGACGATCCCTTCGGTAGATAAAAATGCATTGAACGGCAGGAATGTTTTGCAACTGGCGATTTTATCGCGCa                            <  2:3067189/150‑1 (MQ=255)
                                                                                                                        cGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTGGCAGCAGATCGTGAGTCTCTGTTCTGGATGAATGTTAAGACGATCCCTTCGGTAGATAAAAATGCATTGAACGGCAGGAATGTTTTGCAACTGGCGATTTTATCGCGCAt                           >  2:152983/1‑150 (MQ=255)
                                                                                                                        cGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTGGCAGCAGATCGTGAGTCTCTGTTCTGGATGAATGTTAAGACGATCCCTTCGGTAGATAAAAATGCATTGAACGGCAGGAATGTTTTGCAACTGGCGATTTTATCGCGCAt                           <  2:2658443/150‑1 (MQ=255)
                                                                                                                              tACTTTGCGTATTATTTACACCTGTCCACCGCTGGCAGCAGATCGTGAGTCTCTGTTCTGGATGAATGTTAAGACGATCCCTTCGGTAGATAAAAATGCATTGCACGGCAGGAATGTTTTGCAACTGGCGATTTTATCGCGCATGAAatt                     <  1:2839119/150‑1 (MQ=255)
                                                                                                                               aCTTTGCGTATTATTTACACCGGTCCACCGCTGGCAGCAGATCGTGAGTCTCTGTTCTGGATGAATGTTAAGACGATCCCTTCGGTAGATAAAAATGCATTGAACGGCAGGAATGTTTTGCAACTGGCGATTTTATCGCGCATGAAatta                    >  2:4307252/1‑150 (MQ=255)
                                                                                                                                cTTTGCGTATTATTTACACCGGTCCACCGCTGGCAGCAGATCGTGAGTCTCTGTTCTGGATGAATGTTAAGACGATCCCTTCGGTAGATAAAAATGCATTGAACGGCAGGAATGTTTTGCAACTGGCGATTTTATCGCGCATGAAattat                   <  2:2253312/150‑1 (MQ=255)
                                                                                                                                 tttGCGTATTATTTACACCGGTCCACCGCTGGCAGCAGATCGTGAGTCTCTGTTCTGGATGAATGTTAAGACGATCCCTTCGGTATATAAAAATGCATTGAACGGCAGGAATGTTTTGCAACTGGCGATTTTATCGCGCATGAAattatt                  <  2:1059230/150‑1 (MQ=255)
                                                                                                                                 tttGCGTATTATTTACACCGGTCCACCGCTGGCAGCAGATCGTGAGTCTCTGTTCTGGATGAATGTTAAGACGATCCCTTCGGTAGATAAAAATGCATTGAACGGCAGGAATGTTTTGCAACTGGCGATTTTATCGCGCATGAAattatt                  <  1:4591973/150‑1 (MQ=255)
                                                                                                                                  ttGCGTATTATTTACACCTGTCCACCGCTGGCAGCAGATCGTGAGTCTCTGTTCTGGATGAATGTTAAGACGATCCCTTCGGTAGATAAAAATGCATTGAACGGCAGGAATGTTTTGCAACTGGCGATTTTATCGCGCATGAAATTAttt                 <  1:619435/150‑1 (MQ=255)
                                                                                                                                   tGCGTATTATTTACACCGGTCCACCGCTGGCAGCAGATCGTGAGTCTCTGTTCTGGATGAATGTTAAGACGATCCCTTCGGTAGATAAAAATGCATTGAACGGCAGGAATGTTTTGCAACTGGCGATTTTATCGCGCATGAAATTATTtc                <  2:3519357/150‑1 (MQ=255)
                                                                                                                                   tGCGTATTATTTACACCGGTCCACCGCTGGCAGCAGATCGTGAGTCTCTGTTCTGGATGAATGTTAAGACGATCCCTTCGGTAGATAAAAATGCATTGAACGGCAGGAATGTTTTGCAACTGGCGATTTTATCGCGCATGAAATTATTtc                >  1:3561981/1‑150 (MQ=255)
                                                                                                                                     cGTATTATTTACACCGGTCCACCGCTGGCAGCAGATCGTGAGTCTCTGTTCTGGATGAATGTTAAGACGATCCCTTCGGTAGATAAAAATGCATTGAACGGCAGGAATGTTTTGCAACTGGCGATTTTATCGCGCATGAAATTATTtctc              <  2:3599215/150‑1 (MQ=255)
                                                                                                                                     cGTATTATTTACACCGGTCCACCGCTGGCAGCAGATCGTGAGTCTCTGTTCTGGATGAATGTTAAGACGATCCCTTCGGTAGATAAAAATGCATTGAACGGCAGGAATGTTTTGCAACTGGCGATTTTATCGCGCATGAAATTATTtctc              >  2:3100120/1‑150 (MQ=255)
                                                                                                                                       tattatTTACACCGGTCCACCGCTGGCAGCAGATCGTGAGTCTCTGTTCTGGATGAATGTTAAGACGATCCCTTCGGTAGATAAAAATGCATTGAACGGCAGGAATGTTTTGCAACTGGCGATTTTATCGCGCATGAAATTATTTCTCCg            >  2:3715810/1‑150 (MQ=255)
                                                                                                                                        attattTACACCGGTCCACCGCTGGCAGCAGATCGTGAGTCTCTGTTCTGGATGAATGTTAAGACGATCCCTTCGGTAGATAAAAATGCATTGAACGGCAGGAATGTTTTGCAACTGGCGATTTTATCGCGCATGAAATTATTTCTCCGt           >  2:1119289/1‑150 (MQ=255)
                                                                                                                                         ttattTACACCGGTCCACCGCTGGCAGCAGATCGTGAGTCTCTGTTCTGGATGAATGTTAAGACGATCCCTTCGGTAGATAAAAATGCATTGAACGGCAGGAATGTTTTGCAACTGGCGATTTTATCGCGCATGAAATTATTTCTCCGTc          <  1:2421754/150‑1 (MQ=255)
                                                                                                                                         ttattTACACCGGTCCACCGCTGGCAGCAGATCGTGAGTCTCTGTTCTGGATGAATGTTAAGACGATCCCTTCGGTAGATAAAAATGCATTGAACGGCAGGAATGTTTTGCAACTGGCGATTTTATCGCGCATGAAATTATTTCTCCGTc          <  1:2021247/150‑1 (MQ=255)
                                                                                                                                          tattTACACCGGTCCACCGCTGGCAGCAGATCGTGAGTCTCTGTTCTGGATGAATGTTAAGACGATCCCTTCGGTAGATAAAAATGCATTGAACGGCAGGAATGTTTTGCAACTGGCGATTTTATCGCGCATGAAATTATTTCTCCGTcc         <  1:534703/150‑1 (MQ=255)
                                                                                                                                           attTACACCGGTCCACCGCTGGCAGCAGATCGTGAGTCTCTGTTCTGGATGAATGTTAAGACGATCCCTTCGGTAGATAAAAATGCATTGAACGGCAGGAATGTTTTGCAACTGGCGATTTTATCGCGCATGAAATTATTTCTCCGTCca        <  1:1595430/150‑1 (MQ=255)
                                                                                                                                           attTACACCGGTCCACCGCTGGCAGCAGATCGTGAGTCTCTGTTCTGGATGAATGTTAAGACGATCCCTTCGGTAGATAAAAATGCATTGAACGGCAGGAATGTTTTGCAACTGGCGATTTTATCGCGCATGAAATTATTTCTCCGTCca        <  2:2430717/150‑1 (MQ=255)
                                                                                                                                              tACACCGGTCCACCGCTGGCAGCAGATCGTGAGTCTCTGTTCTGGATGAATGTTAAGACGATCCCTTCGGTAGATAAAAATGCATTGAACGGCAGGAATGTTTTGCAACTGGCGATTTTATCGCGCATGAAATTATTTCTCCGTCcaatt     >  1:3796691/1‑150 (MQ=255)
                                                                                                                                              tACACCGGTCCACCGCTGGCAGCAGATCGTGAGTCTCTGTTCTGGATGAATGTTAAGACGATCCCTTCGGTAGATAAAAATGCATTGAACGGCAGGAATGTTTTGCAACTGGCGATTTTATCGCGCATGAAATTATTTCTCCGTCcaatt     >  1:4402262/1‑150 (MQ=255)
                                                                                                                                                cacCGGTCCACCGCTGGCAGCAGATCGTGAGTCTCTGTTCTGGATGAATGTTAAGACGATCCCTTCGGTAGATAAAAATGCATTGAACGGCAGGAATGTTTTGCAACTGGCGATTTTATCGCGCATGAAATTATTTCTCCGTCcaattca   >  1:3646021/1‑150 (MQ=255)
                                                                                                                                                 acCGGTCCACCGCTGGCAGCAGATCGTGAGTCTCTGTTCTGGATGAATGTTAAGACGATCCCTTCGGTAGATAAAAATGCATTGAACGGCAGGAATGTTTTGCAACTGGCGATTTTATCGCGCATGAAATTATTTCTCCGTCcaattcaa  <  1:64797/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
GATTAGAAATAGCCATACCAATGAGCGCTTTCTGGTCAATTCGTGGATTGAAAACAGCAGCGGTGTAAAAGAAAAGTCATTCATCATTACACCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTGGCAGCAGATCGTGAGTCTCTGTTCTGGATGAATGTTAAGACGATCCCTTCGGTAGATAAAAATGCATTGAACGGCAGGAATGTTTTGCAACTGGCGATTTTATCGCGCATGAAATTATTTCTCCGTCCAATTCAA  >  CP009273/554567‑554861

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: