Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 557,784 G→T 58.0% intergenic (+28/‑8) sfmD → / → sfmH putative outer membrane export usher protein/FimA homolog

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273557,7840GT58.0% 20.6 / inf 236intergenic (+28/‑8)sfmD/sfmHputative outer membrane export usher protein/FimA homolog
Reads supporting (aligned to +/- strand):  ref base G (64/35);  new base T (62/75);  total (126/110)
Fisher's exact test for biased strand distribution p-value = 3.67e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.45e-01

Only 100 of 236 total aligned reads displayed.

CTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGCTTACAGCAAGCCGTCACTGTTATTTCGGCAGTTTGCACACATCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCAGATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTAATAGCACCAATAACCAGCCAGGGCAGATTGTTGTT  >  CP009273/557637‑557932
                                                                                                                                                   |                                                                                                                                                    
  ggTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGCTTACAGCAAGCCGTCACTGTTATTTCGGCAGTTTGCACACATCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCATatta                                                                                                                                                  <  1:1406620/150‑1 (MQ=255)
     tCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGCTTACAGCAAGCCGTCACTGTTATTTCGGCAGTTTGCACACATCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCATattatta                                                                                                                                               >  2:3825079/1‑150 (MQ=255)
     tCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGAAAAGCTTACAGCAAGCCGTCACTGTTATTTCGGCAGTTTGCACACATCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCAGattatta                                                                                                                                               <  1:269309/150‑1 (MQ=255)
         tGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGCTTACAGCAAGCCGTCACTGTTATTTCGGCAGTTTGCACACATCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCATATTATTATTgg                                                                                                                                           <  1:948049/150‑1 (MQ=255)
          ggggAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGCTTACAGCAAGCCGTCACCGTTATTTCGGCAGTTTGCACACATCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCATATTATTATTGGc                                                                                                                                          >  1:2882768/1‑150 (MQ=255)
           gggAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGCTTACAGCAAGCCGTCACTGTTATTTCGGCAGTTTGCACACATCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCAGATTATTATTGGCg                                                                                                                                         >  2:563528/1‑150 (MQ=255)
              agaCGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGCTTACAGCAAGCCGTCACTGTTATTTCGGCAGTTTGCACACATCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCAGATTATTATTGGCGATg                                                                                                                                      >  2:3612706/1‑150 (MQ=255)
              agaCGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGCTTACAGCAAGCCGTCACTGTTATTTCGGCAGTTTGCACACATCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCAGATTATTATTGGCGATg                                                                                                                                      >  2:4550950/1‑150 (MQ=255)
               gaCGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGCTTACAGCAAGCCGTCACTGTTATTTCGGCAGTTTGCACACATCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCATATTATTATTGGCGATgg                                                                                                                                     <  2:751019/150‑1 (MQ=255)
                       gAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGCTTACAGCAAGCCGTCACTGTTATTTCGGCAGTTTGCACACATCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCATATTATTATTGGCGATGGCATGTTtg                                                                                                                             >  1:2661037/1‑150 (MQ=255)
                        aaCTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGCTTACAGCAAGCCGTCACTGTTATTTCGGCAGTTTGCACACATCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCAGATTATTATTGGCGATGGCATGTTtgt                                                                                                                            >  2:1599256/1‑150 (MQ=255)
                          cTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGCTTACAGCAAGCCGTCACTGTTATTTCGGCAGTTTGCACACATCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCATATTATTATTGGCGATGGCATGTTTgtgt                                                                                                                          <  2:3094728/150‑1 (MQ=255)
                           tCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGCTTACAGCAAGCCGTCACTGTTATTTCGGCAGTTTGCACACATCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCAGATTATTATTGGCGATGGCATGTTTGTGTc                                                                                                                         <  1:4085016/150‑1 (MQ=255)
                           tCACGCTGCATTGCGCACTATGTAGTGCCGAAGCAAAGCTTACAGCAAGCCGTCACTGTTATTTCGGCAGTTTGCACACATCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCAGATTATTATTGGCGATGGCATGTTTGTGTc                                                                                                                         <  1:714790/150‑1 (MQ=255)
                             aCGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGCTTACAGCAAGCCGTCACTGTTATTTCGGCAGTTTGCACACATCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCAGATTATTATTGGCGATGGCATGTTTGTGTCTg                                                                                                                       <  1:1041133/150‑1 (MQ=255)
                             aCGCTGCATTGCGCACTATGTAGTGCCGAAGCAAAGCTTACAGCAAGCCGTCACTGTTATTTCGGCAGTTTGCACACATCCTGGCTCATACAGGAAATTATCAATAAGATAATCTGCAGAGTATTATTGGCGATGGCATGTTTGTGTCTg                                                                                                                       <  2:4590171/150‑1 (MQ=255)
                                cTGCATTGCGCACTATGTATTGCCGAAGCAAAGCTTACAGCAAGCCGTCACTGTTATTTCGGCAGTTTGCACACATCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCAGATTATTATTGGCGATGGCATGTTTGTGTCTGGCa                                                                                                                    <  1:2203592/150‑1 (MQ=255)
                                 tGCATTGCGCACTATGTATTGCCGAAGCAAAGCTTACAGCAAGCCGTCACTGTTATTTCGGCAGTTTGCACACATCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCATATTATTATTGGCGATGGCATGTTTGTGTCTGGCaa                                                                                                                   <  1:1325803/150‑1 (MQ=255)
                                   cATTGCGCACTATGTATTGCCGAAGCAAAGCTTACAGCAAGCCGTCACTGTTATTTCGGCAGTTTGCACACATCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCAGATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAAc                                                                                                                 >  2:3030263/1‑150 (MQ=255)
                                    aTTGCGCACTATGTATTGCCGAAGCAAAGCTTACAGCAAGCCGTCACTGTTATTTCGGCAGTTTGCACACATCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCATATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACa                                                                                                                <  2:1627758/150‑1 (MQ=255)
                                    aTTGCGCACTATGTATTGCCGAAGCAAAGCTTACAGCAAGCCGTCACTGTTATTTCGGCAGTTTGCACACATCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCAGATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACa                                                                                                                >  2:3355744/1‑150 (MQ=255)
                                        cgcACTATGTATTGCCGAAGCAAAGCTTACAGCAAGCCGTCACTGTTATTTCGGCAGTTTGCACACATCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCAGATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATc                                                                                                            >  2:3916560/1‑150 (MQ=255)
                                        cgcACTATGTATTGCCGAAGCAAAGCTTACAGCAAGCAGTCACTGTTATTTCGGCAGTTTGCACACATCCTGGATCATAAAGGAAATTATCAATAAGATAATCTGCAGATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATc                                                                                                            <  1:3100898/150‑1 (MQ=255)
                                             tATGTATTGCCGAAGCAAAGCTTACAGCAAGCCGTCACTGTTATTTCGGCAGTTTGCACACATCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCATATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTggg                                                                                                       >  2:3841584/1‑150 (MQ=255)
                                               tGTATTGCCGAAGCAAAGCTTACAGCAAGCCGTCACTGTTATTTCGGCAGTTTGCACACATCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCAGATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCt                                                                                                     >  1:49802/1‑150 (MQ=255)
                                                gTATTGCCGAAGCAAAGCTTACAGCAAGCCGTCACTGTTATTTCGGCAGTTTGCACACATCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCAGATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTa                                                                                                    >  1:996205/1‑150 (MQ=255)
                                                      ccGAAGCAAAGCTTACAGCAAGCCGTCACTGTTATTTCGGCAGTTTGCACACATCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCATATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGttt                                                                                              >  2:2041656/1‑150 (MQ=255)
                                                       cGAAGCAAAGCTTACAGCAAGCCGTCACTGTTATTTCGGCAGTTTGCACACATCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCAGATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTtg                                                                                             >  1:2470942/1‑150 (MQ=255)
                                                       cGAAGCAAAGCTTACAGCAAGCCGTCACTGTTATTTCGGCAGTTTGCACACATCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCAGATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTtg                                                                                             >  1:1891927/1‑150 (MQ=255)
                                                       cGAAGCAAAGCTTACAGCAAGCCGTCACTGTTATTTCGGCAGTTTGCACACATCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCAGATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTtg                                                                                             >  2:1846026/1‑150 (MQ=255)
                                                       cGAAGCAAAGCTTACAGCAAGCCGTCACTGTTATTTCGGCAGTTTGCACACATCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCAGATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTtg                                                                                             >  2:3243999/1‑150 (MQ=255)
                                                        gAAGCAAAGCTTACAGCAAGCCGTCACTGTTATTTCGGCAGTTTGCACACATCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCATATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTtgt                                                                                            <  2:4434839/150‑1 (MQ=255)
                                                        gAAGCAAAGCTTACAGCAAGCCGTCACTGTTATTTCGGCAGTTTGCACACATCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCAGATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTtgt                                                                                            >  1:3554530/1‑150 (MQ=255)
                                                         aaGCAAAGCTTACAGCAAGCCGTCACTGTTATTTCGGCAGTTTGCACACATCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCATATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTtgtg                                                                                           >  2:926135/1‑150 (MQ=255)
                                                         aaGCAAAGCTTACAGCAAGCCGTCACTGTTATTTCGGCAGTTTGCACACATCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCAGATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTtgtg                                                                                           <  1:970463/150‑1 (MQ=255)
                                                           gCAAAGCTTACAGCAAGCCGTCACTGTTATTTCGGCAGTTTGCACACATCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCATATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCa                                                                                         <  2:241085/150‑1 (MQ=255)
                                                              aaGCTTACAGCAAGCCGTCACTGTTATTTCGGCAGTTTGCACACATCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCAGATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAAt                                                                                      >  2:2468414/1‑150 (MQ=255)
                                                               aGCTTACAGCAAGCCGTCACTGTTATTTCGGCAGTTTGCACACATCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCATATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATa                                                                                     <  1:2365560/150‑1 (MQ=255)
                                                                   tACAGCAAGCCGTCACTGTTATTTCGGCAGTTTGCACACATCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCATATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTAc                                                                                 >  1:4170920/1‑150 (MQ=255)
                                                                    aCAGCAAGCCGTCACTGTTATTTCGGCAGTTTGCACACAACCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCAGATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACAGTTTGTGCAAATAGTACt                                                                                >  1:3171989/1‑150 (MQ=255)
                                                                      aGCAAGCCGTCACTGTTATTTCGGCAGTTTGCACACATCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCAGATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTgg                                                                              >  2:12648/1‑150 (MQ=255)
                                                                       gCAAGCCGTCACTGTTATTTCGGCAGTTTGCACACATCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCAGATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGc                                                                             <  1:2802073/150‑1 (MQ=255)
                                                                             cGTCACTGTTATTTCGGCAGTTTGCACACATCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCATATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCa                                                                       <  2:1269431/150‑1 (MQ=255)
                                                                             cGTCACTGTTATTTCGGCAGTTTGCACACATCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCAGATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCa                                                                       >  1:4296520/1‑150 (MQ=255)
                                                                               tCACTGTTATTTCGGCAGTTTGCACACATCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCAGATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCaga                                                                     >  1:4341412/1‑150 (MQ=255)
                                                                                  cTGTTATTTCGGCAGTTTGCACACATCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCAGATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCagaaga                                                                  >  2:587431/1‑150 (MQ=255)
                                                                                  cTGTTATTTCGGCAGTTTGCACACATCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCAGATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCagaaga                                                                  >  2:2280889/1‑150 (MQ=255)
                                                                                   tGTTATTTCGGCAGTTTGCACACATCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCATATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGAt                                                                 <  2:1229282/150‑1 (MQ=255)
                                                                                     ttATTTCGGCAGTTTGCACACATCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCATATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGa                                                               >  2:3917311/1‑150 (MQ=255)
                                                                                       aTTTCGGCAGTTTGCACACATCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCATATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAac                                                             >  1:1095047/1‑150 (MQ=255)
                                                                                       aTTTCGGCAGTTTGCACACATCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCATATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAac                                                             >  1:1999180/1‑150 (MQ=255)
                                                                                        tttCGGCAGTTTGCACACATCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCATATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAaca                                                            <  1:1502897/150‑1 (MQ=255)
                                                                                         ttCGGCAGTTTGCACACATCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCATATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAacac                                                           >  1:3070831/1‑150 (MQ=255)
                                                                                            ggCAGTTTGCACACATCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCATATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTAt                                                        <  1:2550675/150‑1 (MQ=255)
                                                                                            ggCAGTTTGCACACATCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCAGATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTAt                                                        >  1:1384306/1‑150 (MQ=255)
                                                                                             gCAGTTTGCACACATCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCAGATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATg                                                       <  2:1483519/150‑1 (MQ=255)
                                                                                                  ttGCACACATCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCATATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGAtctc                                                  <  1:3999524/150‑1 (MQ=255)
                                                                                                   tGCACACATCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCAGATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGAtctct                                                 <  2:4621155/150‑1 (MQ=255)
                                                                                                     cacacaTCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCAGATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCa                                               >  1:134618/1‑150 (MQ=255)
                                                                                                      acacaTCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCAGATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCaa                                              <  2:2663593/150‑1 (MQ=255)
                                                                                                       cacaTCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCATATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCaaa                                             <  1:168104/150‑1 (MQ=255)
                                                                                                        acaTCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCATATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAat                                            <  1:2758183/150‑1 (MQ=255)
                                                                                                        acaTCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCATATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAat                                            <  1:1976469/150‑1 (MQ=255)
                                                                                                          aTCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCATATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAatat                                          <  2:3150355/150‑1 (MQ=255)
                                                                                                          aTCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCATATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAatat                                          >  1:3813733/1‑150 (MQ=255)
                                                                                                          aTCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCATATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAatat                                          >  1:3348319/1‑150 (MQ=255)
                                                                                                           tCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCATATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATc                                         >  1:1012084/1‑150 (MQ=255)
                                                                                                            ccTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCATATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCt                                        >  1:2912009/1‑150 (MQ=255)
                                                                                                            ccTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCAGATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCt                                        >  2:3520804/1‑150 (MQ=255)
                                                                                                            ccTGGCTCATAAACGAAATTATCAATAAGATAATCTGCAGATTATTATTGGCGATTGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGAGTAGCAGAAGATGAACACTATGATCTCTCAAATATCt                                        <  1:4627123/150‑1 (MQ=255)
                                                                                                               ggCTCATAAAGGAAATTATCAATAAGATAATCTGCATATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTa                                     <  1:2843575/150‑1 (MQ=255)
                                                                                                                gCTCATAAAGGAAATTATCAATAAGATAATCTGCATATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTaa                                    <  1:1042089/150‑1 (MQ=255)
                                                                                                                 cTCATAAAGGAAATTATCAATAAGATAATCTGCATATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTAAt                                   >  2:2760902/1‑150 (MQ=255)
                                                                                                                 cTCATAAAGGAAATTATCAATAAGATAATCTGCAGATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTAAt                                   >  2:2799415/1‑150 (MQ=255)
                                                                                                                  tCATAAAGGAAATTATCAATAAGATAATCTGCATATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTAATa                                  >  1:3861025/1‑150 (MQ=255)
                                                                                                                  tCATAAAGGAAATTATCAATAAGATAATCTGCATATTATTATTGGCGATGGCATGTGTGTGTCTGGCAAACATAACCTGGGCTACTGTTTGTGCAAATAGTACTGGCTTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTAATa                                  >  1:4590753/1‑150 (MQ=255)
                                                                                                                     tAAAGGAAATTATCAATAAGATAATCTGCAGATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTAATAGCa                               >  1:2062356/1‑150 (MQ=255)
                                                                                                                      aaaGGAAATTATCAATAAGATAATCTGCAGATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTAATAGCAc                              >  1:650325/1‑150 (MQ=255)
                                                                                                                       aaGGAAATTATCAATAAGATAATCTGCAGATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTAATAGCAcc                             <  1:1116837/150‑1 (MQ=255)
                                                                                                                         ggAAATTATCAATAAGATAATCTGCATATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTAATAGCACCaa                           >  1:1846729/1‑150 (MQ=255)
                                                                                                                         ggAAATTATCAATAAGATAATCTGCATATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTAATAGCACCaa                           <  1:3180564/150‑1 (MQ=255)
                                                                                                                         ggAAATTATCAATAAGATAATCTGCATATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTAATAGCACCaa                           >  2:4366094/1‑150 (MQ=255)
                                                                                                                         ggAAATTATCAATAAGATAATCTGCAGATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTAATAGCACCaa                           >  2:498436/1‑150 (MQ=255)
                                                                                                                             aTTATCAATAAGATAATCTGCATATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTAATAGCACCAATAAc                       <  1:1072046/150‑1 (MQ=255)
                                                                                                                              ttATCAATAAGATAATCTGCAGATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTAATAGCACCAATAAcc                      <  2:4049052/150‑1 (MQ=255)
                                                                                                                               tATCAATAAGATAATCTGCAGATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTAATAGCACCAATAAcca                     <  2:960690/150‑1 (MQ=255)
                                                                                                                                aTCAATAAGATAATCTGCATATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTAATAGCACCAATAAccag                    <  2:3421275/150‑1 (MQ=255)
                                                                                                                                aTCAATAAGATAATCTGCAGATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTAATAGCACCAATAAccag                    <  1:563528/150‑1 (MQ=255)
                                                                                                                                 tCAATAAGATAATCTGCAGATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTAATAGCACCAATAAccagc                   >  2:665408/1‑150 (MQ=255)
                                                                                                                                   aaTAAGATAATCTGCATATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTAATAGCACCAATAAccagcca                 <  1:2041679/150‑1 (MQ=255)
                                                                                                                                    aTAAGATAATCTGCATATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTAATAGCACCAATAAccagccag                >  1:3332453/1‑150 (MQ=255)
                                                                                                                                    aTAAGATAATCTGCAGATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTAATAGCACCAATAAccagccag                >  1:791214/1‑150 (MQ=255)
                                                                                                                                    aTAAGATAATCTGCAGATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTAATAGCACCAATAAccagccag                >  2:83908/1‑150 (MQ=255)
                                                                                                                                      aaGATAATCTGCATATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTAATAGCACCAATAACCAGCCAggg              >  1:2174577/1‑150 (MQ=255)
                                                                                                                                       aGATAATCTGCATATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTAATAGCACCAATAACCAGCCAGGGc             >  2:2810207/1‑150 (MQ=255)
                                                                                                                                       aGATAATCTGCAGATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTAATAGCACCAATAACCAGCCAGGGc             <  2:2057232/150‑1 (MQ=255)
                                                                                                                                        gATAATCTGCATATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTAATAGCACCAATAACCAGCCAGGGCa            <  1:4331779/150‑1 (MQ=255)
                                                                                                                                        gATAATCTGCATATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTAATAGCACCAATAACCAGCCAGGGCa            <  1:63473/150‑1 (MQ=255)
                                                                                                                                           aaTCTGCATATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTAATAGCACCAATAACCAGCCAGGGCAGAt         >  2:3256454/1‑150 (MQ=255)
                                                                                                                                                  aGATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTAATAGCACCAATAACCAGCCAGGGCAGATTgttgtt  >  2:1708412/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
CTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGCTTACAGCAAGCCGTCACTGTTATTTCGGCAGTTTGCACACATCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCAGATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTAATAGCACCAATAACCAGCCAGGGCAGATTGTTGTT  >  CP009273/557637‑557932

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: