Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 557,887 A→T 6.4% S32S (TCA→TCT)  sfmH → FimA homolog

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273557,8870AT 6.4% 813.3 / 34.7 236S32S (TCA→TCT) sfmHFimA homolog
Reads supporting (aligned to +/- strand):  ref base A (118/103);  new base T (9/6);  total (127/109)
Fisher's exact test for biased strand distribution p-value = 7.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.07e-01

Only 100 of 248 total aligned reads displayed.

CACACATCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCAGATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTAATAGCACCAATAACCAGCCAGGGCAGATTGTTGTTTTACCGGAAAAATCCGGCTGGGTAGGTGTCTCAGCAATTTGTCCACCCGGTACGCTGGTGAATTATACATACCGTAGTTATGTCACCAACTTTATTGTTCAGGA  >  CP009273/557738‑558036
                                                                                                                                                     |                                                                                                                                                     
cacacaTCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCAGATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCa                                                                                                                                                       >  2:520801/1‑150 (MQ=255)
 acacaTCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCAGATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCaa                                                                                                                                                      >  2:447723/1‑150 (MQ=255)
      tCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCAGATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATc                                                                                                                                                 >  2:2149233/1‑150 (MQ=255)
        cTGGCTCATAAAGGATATTATCAATAAGATTATCTGCAGATTATTATTGGCGATGGCATGTTTGTGTCTGGCTAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGTACACTATGATCTCTCAAATATCtt                                                                                                                                               <  2:2968239/150‑1 (MQ=255)
        cTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCAGATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCtt                                                                                                                                               <  1:2642272/150‑1 (MQ=255)
         tGGCTCATAAAGGAAATTATCAATAAGATAATCTGCAGATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATTCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCttt                                                                                                                                              >  1:4123857/1‑150 (MQ=255)
         tGGCTCATAAAGGAAATTATCAATAAGATAATCTGCAGATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCttt                                                                                                                                              <  1:1014954/150‑1 (MQ=255)
           gCTCATAAAGGAAATTATCAATAAGATACTCTGCAGATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTaa                                                                                                                                            >  2:480851/1‑150 (MQ=255)
           gCTCATAAAGGAAATTATCAATAAGATAATCTGCAGATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTaa                                                                                                                                            >  1:4187399/1‑150 (MQ=255)
            cTCATAAAGGAAATTATCAATAAGATAATCTGCAGATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATTCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTAAt                                                                                                                                           >  2:833167/1‑150 (MQ=255)
            cTCATAAAGGAAATTATCAATAAGATAATCTGCAGATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTAAt                                                                                                                                           <  2:1773376/150‑1 (MQ=255)
             tCATAAAGGAAATTATCAATAAGATAATCTGCAGATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATTCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTAATa                                                                                                                                          <  2:3797905/150‑1 (MQ=255)
              cATAAAGGAAATTATCAATAAGATAATCTGCAGATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTAATAg                                                                                                                                         >  1:4028033/1‑150 (MQ=255)
               aTACAGGAAATTATCAATACGATAATCTGCAGATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTAATAGc                                                                                                                                        <  1:4865065/150‑1 (MQ=255)
                 aaaGTAAATTATCAATAAGATAATCTGCAGATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTAATAGCAc                                                                                                                                      <  2:4831555/150‑1 (MQ=255)
                 aaaGGAAATTATCAATAAGATAATCTGCAGATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTAATAGCAc                                                                                                                                      >  2:4174111/1‑150 (MQ=255)
                     gAAATTATCAATAAGATAATCTGCAGATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTAATAGCACCAAt                                                                                                                                  <  1:667044/150‑1 (MQ=255)
                       aaTTATCAATAAGATAATCTGCAGATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTAATAGCACCAATaa                                                                                                                                <  1:736749/150‑1 (MQ=255)
                           aTCAATAAGATACTCTGCAGATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTAATAGCACCAATAAccag                                                                                                                            >  2:3720073/1‑150 (MQ=255)
                             cAATAAGATAATCTGCATATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTAATAGCACCAATAAccagcc                                                                                                                          <  2:1612566/150‑1 (MQ=255)
                               aTAAGATAATCTGCAGATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTAATAGCACCAATAAccagccag                                                                                                                        >  2:2574138/1‑150 (MQ=255)
                               aTAAGATAATCTGCAGATTATTATAGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTAATAGCACCAATAAccagacag                                                                                                                        >  1:1475586/1‑150 (MQ=255)
                                 aaGATAATCTGCAGATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTAATAGCACCAATAACCAGCCAggg                                                                                                                      <  1:4224371/150‑1 (MQ=255)
                                 aaGATAATCTGCAGATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTAATAGCACCAATAACCAGCCAggg                                                                                                                      <  2:2838710/150‑1 (MQ=255)
                                    aTAATCTGCAGATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTAATAGCACCAATAACCAGCCAGGGCAg                                                                                                                   <  1:579680/150‑1 (MQ=255)
                                     tAATCTGCAGATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCTAATATCTTTAATAGCACCAATAACCAGCCAGGGCAGa                                                                                                                  >  1:1004894/1‑150 (MQ=255)
                                      aaTCTGCAGATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTAATAGCACCAATAACCAGCCAGGGCAGAt                                                                                                                 >  2:4608196/1‑150 (MQ=255)
                                         cTGCAGATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTAATAGCACCAATAACCAGCCAGGGCAGAttgt                                                                                                              >  2:4182283/1‑150 (MQ=255)
                                          tGCAGATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTACTAGCACCAATAACCAGCCAGGGCAGAttgtt                                                                                                             <  1:4606511/150‑1 (MQ=255)
                                             aGATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTAATAGCACCAATAACCAGCCAGGGCAGATTgttgtt                                                                                                          >  1:2074393/1‑150 (MQ=255)
                                             aGATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTAATAGCACCAATAACCAGCCAGGGCAGATTgttgtt                                                                                                          >  1:985180/1‑150 (MQ=255)
                                                  attattGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTAATAGCACCAATAACCAGCCAGGGCAGATTGTTGTTTTAcc                                                                                                     >  1:600266/1‑150 (MQ=255)
                                                   ttattGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTAATAGCACCAATAACCAGCCAGGGCAGATTGTTGTTTTACCg                                                                                                    >  2:3950778/1‑150 (MQ=255)
                                                     attGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCTAATATCTTTAATAGCACCAATAACCAGCCAGGGCAGATTGTTGTTTTACCGGa                                                                                                  >  2:2129258/1‑150 (MQ=255)
                                                           gATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTAATAGCACCAATAACCAGCCAGGGCAGATTGTTGTTTTACCGGAAAAATc                                                                                            >  1:3902978/1‑150 (MQ=255)
                                                              ggCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTAATAGCACCAATAACCAGCCAGGGCAGATTGTTGTTTTACCGGAAAAATCCgg                                                                                         >  1:1287033/1‑150 (MQ=255)
                                                                 aTGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTAATAGCACCAATAACCAGCCAGGGCAGATTGTTGTTTTACCGGAAAAATCCGGCTg                                                                                      <  1:1729935/150‑1 (MQ=255)
                                                                 aTGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTAATAGCACCAATAACCAGCCAGGGCAGATTGTTGTTTTACCGGAAAAATCCGGCTg                                                                                      >  2:3158578/1‑150 (MQ=255)
                                                                 aTGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTAATAGCACCAATAACCAGCCAGGGCAGATTGTTGTTTTACCGGAAAAATCCGGCTg                                                                                      <  2:1791575/150‑1 (MQ=255)
                                                                  tGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACGATGATCTCTCAAATATCTTTAATAGCACCAATAACCAGCCAGGGCAGATTGTTGTTTTACCGGAAAAATCCGGCTgg                                                                                     <  1:3537705/150‑1 (MQ=255)
                                                                   gTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTAATAGCACCAATAACCAGCCAGGGCAGATTGTTGTTTTACCGGAAAAATCCGGCTggg                                                                                    <  2:3174010/150‑1 (MQ=255)
                                                                    tttGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTAATAGCACCAATAACCAGCCAGGGCAGATTGTTGTTTTACCGGAAAAATCCGGCTGGGt                                                                                   <  1:4704843/150‑1 (MQ=255)
                                                                       gtgtCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTAATAGCACCAATAACCAGCCAGGGCAGATTGTTGTTTTACCGGAAAAATCCGGCTGGGTAgg                                                                                >  1:1322616/1‑150 (MQ=255)
                                                                        tgtCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTAATAGCACCAATAACCAGCCAGGGCAGATTGTTGTTTTACCGGAAAAATCCGGCTGGGTAGgt                                                                               >  2:1583600/1‑150 (MQ=255)
                                                                            tGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTAATAGCACCAATAACCAGCCAGGGCAGATTGTTGTTTTACCGGAAAAATCCGGCTGGGTAGGTGtct                                                                           <  1:3206459/150‑1 (MQ=255)
                                                                            tGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTAATAGCACCAATAACCAGCCAGGGCAGATTGTTGTTTTACCGGAAAAATCCGGCTGGGTAGGTGtct                                                                           >  2:3620092/1‑150 (MQ=255)
                                                                             ggCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTAATAGCACCAATAACCAGCCAGGGCAGATTGTTGTTTTACCGGAAAAATCCGGCTGGGTAGGTGtctc                                                                          <  1:4594640/150‑1 (MQ=255)
                                                                             ggCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTAATAGCACCAATAACCAGCCAGGGCAGATTGTTGTTTTACCGGAAAAATCCGGCTGGGTAGGTGtctc                                                                          <  1:2596437/150‑1 (MQ=255)
                                                                               cAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCTAATATCTTTAATAGCACCAATAACCAGCCAGGGCAGATTGTTGTTTTACCGGAAAAATCCGGCTGGGTAGGTGTCTCAg                                                                        >  1:2368256/1‑150 (MQ=255)
                                                                               cAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTAATAGCACCAATAACCAGCCAGGGCAGATTGTTGTTTTACCGGAAAAATCCGGCTGGGTAGGTGTCTCAg                                                                        <  1:3334300/150‑1 (MQ=255)
                                                                                aaaCATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTAATAGCACCAATAACCAGCCAGGGCAGATTGTTGTTTTACCGGAAAAATCCGGCTGGGTAGGTGTCTCAGc                                                                       >  1:1023637/1‑150 (MQ=255)
                                                                                  aCATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCTAATATCTTTAATAGCACCAATAACCAGCCAGGGCAGATTGTTGTTTTACCGGAAAAATCCGGCTGGGTAGGTGTCTCAGCaa                                                                     >  2:4840270/1‑150 (MQ=255)
                                                                                   cATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTAATAGCACCAATAACCAGCCAGGGCAGATTGTTGTTTTACCGGAAAAATCCGGCTGGGTAGGTGTCTCAGCAAt                                                                    <  1:2335740/150‑1 (MQ=255)
                                                                                      atCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTAATAGCACCAATAACCAGCCAGGGCAGATTGTTGTTTTACCGGAAAAATCCGGCTGGGTAGGTGTCTCAGCAATTTg                                                                 <  2:1423143/150‑1 (MQ=255)
                                                                                         cTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTAATAGCACCAATAACCAGCCAGGGCAGATTGTTGTTTTACCGGAAAAATCCGGCTGGGTAGGTGTCTCAGCAATTTGTcc                                                              >  1:2473649/1‑150 (MQ=255)
                                                                                          tGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTAATAGCACCAATAACCAGCCAGGGCAGATTGTTGTTTTACCGGAAAAATCCGGCTGGGTAGGTGTCTCAGCAATTTGTCCa                                                             >  2:3414518/1‑150 (MQ=255)
                                                                                          tGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTAATAGCACCAATAACCAGCCAGGGCAGATTGTTGTTTTACCGGAAAAATCCGGCTGGGTAGGTGTCTCAGCAATTTGTCCa                                                             <  1:1165566/150‑1 (MQ=255)
                                                                                               tACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTAATAGCACCAATAACCAGCCAGGGCAGATTGTTGTTTTACCGGAAAAATCCGGCTGGGTAGGTGTCTCAGCAATTTGTCCACCCgg                                                        >  2:275437/1‑150 (MQ=255)
                                                                                                 cTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTAATAGCACCAATAACCAGCCAGGGCAGATTGTTGTTTTACCGGAAAAATCCGGCTGGGTAGGTGTCTCAGCAATTTGTCCACCCGGTa                                                      >  2:3299753/1‑150 (MQ=255)
                                                                                                  tGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTAATAGCACCAATAACCAGCCAGGGCAGATTGTTGTTTTACCGGAAAAATCCGGCTGGGTAGGTGTCTCAGCAATTTGTCCACCCGGTAc                                                     >  1:2984165/1‑150 (MQ=255)
                                                                                                   gTTTTTGCATTTAGTTCTGGCGTAGCAGAAGATTAACACTATGATCTCTCATATATCTTTAATAGCACCAATACCCAGCCAGGGCAGATTGTTGTTTTACCGGAAAAATCCGGCTGGTTAGTTGTCTCAGCAATTTGTCCACCCGGTTCg                                                    <  2:680495/150‑1 (MQ=255)
                                                                                                     ttGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAACTATCTTTAATAGCACCAATAACCAGCCAGGGCAGATTGTTGTTTTACCGGAAAAATCCGGCTGGGTAGGTGTCTCAGCAATTTGTCCACCCGGTACGCt                                                  <  1:3824778/150‑1 (MQ=255)
                                                                                                      tgtgCAAATAGTTCTGGCGTAGCAGAAGATGAACACTATGATCTCTCTAATATCTTTAATAGCACCAATAACCAGCCAGGGCAGATTGTTGTTTTACCGGAAAAATCCGGCTGGGTAGGTGTCTCAGCAATTTGTCCACCCGGTACGCTg                                                 <  2:3009235/150‑1 (MQ=255)
                                                                                                       gtgCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTAATAGCACCAATAACCAGCCAGGGCAGATTGTTGTTTTACCGGAAAAATCCGGCTGGGTAGGTGTCTCAGCAATTTGTCCACCCGGTACGCTgg                                                <  1:4045933/150‑1 (MQ=255)
                                                                                                       gtgCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTAATAGCACCAATAACCAGCCAGGGCAGATTGTTGTTTTACCGGAAAAATCCGGCTGGGTAGGTGTCTCAGCAATTTGTCCACCCGGTACGCTgg                                                <  1:3687613/150‑1 (MQ=255)
                                                                                                          cAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTAATAGCACCAATAACCAGCCAGGGCAGATTGTTGTTTTACCGGAAAAATCCGGCTGGGTAGGTGTCTCAGCAATTTGTCCACCCGGTACGCTGGTGa                                             <  1:5013679/150‑1 (MQ=255)
                                                                                                          cAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTAATAGCACCAATAACCAGCCAGGGCAGATTGTTGTTTTACCGGAAAAATCCGGCTGGGTAGGTGTCTCAGCAATTTGTCCACCCGGTACGCTGGTGa                                             <  1:834198/150‑1 (MQ=255)
                                                                                                           aaaTAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTAATAGCACCAATAACCAGCCAGGGCAGATTGTTGTTTTACCGGAAAAATCCGGCTGGGTAGGTGTCTCAGCAATTTGTCCACCCGGTACGCTGGTGaa                                            >  2:2461358/1‑150 (MQ=255)
                                                                                                            aaTAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTAATAGCACCAATAACCAGCCAGGGCAGATTGTTGTTTTACCGGAAAAATCCGGCTGGGTAGGTGTCTCAGCAATTTGTCCACCCGGTACGCTGGTGAAt                                           >  2:2416186/1‑150 (MQ=255)
                                                                                                            aaTAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTAATAGCACCAATAACCAGCCAGGGCAGATTGTTGTTTTACCGGAAAAATCCGGCTGGGTAGGTGTCTCAGCAATTTGTCCACCCGGTACGCTGGTGAAt                                           <  2:788945/150‑1 (MQ=255)
                                                                                                                gTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTAATAGCACCAATAACCAGCCAGGGCAGATTGTTGTTTTACCGGAAAAATCCGGCTGGGTAGGTGTCTCAGCAATTTGTCCACCCGGTACGCTGGTGAATtata                                       >  2:4343294/1‑150 (MQ=255)
                                                                                                                      gCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTAATAGCACCAATAACCAGCCAGGGCAGATTGTTGTTTTACCGGAAAAATCCGGCTGGGTAGGTGTCTCAGCAATTTGTCCACCCGGTACGCTGGTGAATTATACATAcc                                 >  1:2753965/1‑150 (MQ=255)
                                                                                                                       cGTAGCAGAAGATGAACACTATGATCTCTCTAATATCTTTAATAGCACCAATAACCAGCCAGGGCAGATTGTTGTTTTACCGGAAAAATCCGGCTGGGTAGGTGTCTCAGCAATTTGTCCACCCGGTACGCTGGTGAATTATACATACCg                                >  1:1969532/1‑150 (MQ=255)
                                                                                                                           gCAGAAGATGAACACTATGATCTCTCAAATATCTTTAATAGCACCAATAACCAGCCAGGGCAGATTGTTGTTTTACCGGAAAAATCCGGCTGGGTAGGTGTCTCAGCAATTTGTCCACCCGGTACGCTGGTGAATTATACATACCGTAGt                            <  2:1023568/150‑1 (MQ=255)
                                                                                                                            cAGAAGATGAACACTATGATCTCTCAAATATCTTTAATAGCACCAATAACCAGCCAGGGCAGATTGTTGTTTTACCGGAAAAATCCGGCTGGGTAGGTGTCTCAGCAATTTGTCCACCCGGTACGCTGGTGAATTATACATACCGTAGtt                           >  2:3489703/1‑150 (MQ=255)
                                                                                                                              gaagaTGAACACTATGATCTCTCTAATATCTTTAATAGCACCAATAACCAGCCAGGGCAGATTGTTGTTTTACCGGAAAAATCCGGCTGGGTAGGTGTCTCAGCAATTTGTCCACCCGGTACGCTGGTGAATTATACATACCGTAGTTAt                         <  2:4442826/150‑1 (MQ=255)
                                                                                                                              gaagaTGAACACTATGATCTCTCAAATATCTTTAATAGCACCAATAACCAGCCAGGGCAGATTGTTGTTTTACCGGAAAAATCCGGCTGGGTAGGTGTCTCAGCAATTTGTCCACCCGGTACGCTGGTGAATTATACATACCGTAGTTAt                         >  1:1442270/1‑150 (MQ=255)
                                                                                                                              gaagaTGAACACTATGATCTCTCAAATATCTTTAATAGCACCAATAACCAGCCAGGGCAGATTGTTGTTTTACCGGAAAAATCCGGCTGGGTAGGTGTCTCAGCAATTTGTCCACCCGGTACGCTGGTGAATTATACATACCGTAGTTAt                         >  1:884197/1‑150 (MQ=255)
                                                                                                                               aagaTGAACACTATGATCTCTCAAATATCTTTAATAGCACCAATAACCAGCCAGGGCAGATTGTTGTTTTACCGGAAAAATCCGGCTGGGTAGGTGTCTCAGCAATTTGTCCACCCGGTACGCTGGTGAATTATACATACCGTAGTTATg                        <  1:3152698/150‑1 (MQ=255)
                                                                                                                               aagaTGAACACTATGATCTCTCAAATATCTTTAATAGCACCAATAACCAGCCAGGGCAGATTGTTGTTTTACCGGAAAAATCCGGCTGGGTAGGTGTCTCAGCAATTTGTCCACCCGGTACGCTGGTGAATTATACATACCGTAGTTATg                        <  1:566891/150‑1 (MQ=255)
                                                                                                                               aagaTGAACACTATGATCTCTCAAATATCTTTAATAGCACCAATAACCAGCCAGGGCAGATTGTTGTTTTACCGGAAAAATCCGGCTGGGTAGGTGTCTCAGCAATTTGTCCACCCGGTACGCTGGTGAATTATACATACCGTAGTTATg                        >  2:4510701/1‑150 (MQ=255)
                                                                                                                                 gaTGAACACTATGATCTCTCTAATATCTTTAATAGCACCAATAACCAGCCAGGGCAGATTGTTGTTTTACCGGAAAAATCCGGCTGGGTAGGTGTCTCAGCAATTTGTCCACCCGGTACGCTGGTGAATTATACATACCGTAGTTATGTc                      <  2:1969532/150‑1 (MQ=255)
                                                                                                                                  aTGAACACTATGATCTCTCAAATATCTTTAATAGCACCAATAACCAGCCAGGGCAGATTGTTGTTTTACCGGAAAAATCCGGCTGGGTAGGTGTCTCAGCAATTTGTCCACCCGGTACGCTGGTGAATTATACATACCGTAGTTATGTCa                     <  2:3878050/150‑1 (MQ=255)
                                                                                                                                    gAACACTATGATCTCTCAAATATCTTTAATAGCACCAATAACCAGCCAGGGCAGATTGTTGTTTTACCGGAAAAATCCGGCTGGGTAGGTGTCTCAGCAATTTGTCCACCCGGTACGCTGGTGAATTATACATACCGTAGTTATGTCAcc                   >  1:2297919/1‑150 (MQ=255)
                                                                                                                                    gAACACTATGATCTCTCAAATATCTTTAATAGCACCAATAACCAGCCAGGGCAGATTGTTGTTTTACCGGAAAAATCCGGCTGGGTAGGTGTCTCAGCAATTTGTCCACCCGGTACGCTGGTGAATTATACATACCGTAGTTATGTCAcc                   >  2:5118072/1‑150 (MQ=255)
                                                                                                                                      acacTATGATCTCTCAAATATCTTTAATAGCACCAATAACCAGCCAGGGCAGATTGTTGTTTTACCGGAAAAATCCGGCTGGGTAGGTGTCTCAGCAATTTGTCCACCCGGTACGCTGGTGAATTATACATACCGTAGTTATGTCACCaa                 >  2:2604809/1‑150 (MQ=255)
                                                                                                                                       cacTAAGATCTCTCAAATATCTTTAATAGCACCAATAACCAGCCAGGGCAGATTGTTGTTTTACCGGAAAAATCCGGCTGGGTAGGTGTCTCAGCAATTTGTCCACCCGGTACGCTGGTGACTTATACATACCGTAGTTATGTCACCAAc                >  1:2247629/1‑150 (MQ=255)
                                                                                                                                        acTATGATCTCTCAAATATCTTTAATAGCACCAATAACCAGCCAGGGCAGATTGTTGTTTTACCGGAAAAATCCGGCTGGGTAGGTGTCTCAGCAATTTGTCCACCCGGTACGCTGGTGAATTATACATACCGTAGTTATGTCACCAACt               <  2:1676035/150‑1 (MQ=255)
                                                                                                                                           aTGATCTCTCAAATATCTTTAATAGCACCAATAACCAGCCAGGGCAGATTGTTGTTTTACCGGAAAAATCCGGCTGGGTAGGTGTCTCAGCAATTTGTCCACCCGGTACGCTGGTGAATTATACATACCGTAGTTATGTCACCAACTTTa            <  2:3348160/150‑1 (MQ=255)
                                                                                                                                             gATCTCTCAAATATCTTTAATAGCACCAATAACCAGCCAGGGCAGATTGTTGTTTTACCGGAAAAATCCGGCTGGGTAGGTGTCTCAGCAATTTGTCCACCCGGTACGCTGGTGAATTATACATACCGTAGTTATGTCACCAACTTTAtt          >  2:4825194/1‑150 (MQ=255)
                                                                                                                                              aTCTCTCAAATATCTTTAATAGCACCAATAACCAGCCAGGGCAGATTGTTGTTTTACCGGAAAAATCCGGCTGGGTAGGTGTCTCAGCAATTTGTCCACCCGGTACGCTGGTGAATTATACATACCGTAGTTATGTCACCAACTTTATTg         >  1:1390685/1‑150 (MQ=255)
                                                                                                                                               tctctcAAATATCTTTAATAGCACCAATAACCAGCCAGGGCAGATTGTTGTTTTACCGGAAAAATCCGGCTGGGTAGGTGTCTCAGCAATTTGTCCACCCGGTACGCTGGTGAATTATACATACCGTAGTTATGTCACCAACTTTATTGt        >  2:65726/1‑150 (MQ=255)
                                                                                                                                                ctctcAAATATCTTTAATAGCACCAATAACCAGCCAGGGCAGATTGTTGTTTTACCGGAAAAATCCGGCTGGGTAGGTGTCTCAGCAATTTGTCCACCCGGTACGCTGGTGAATTATACATACCGTAGTTATGTCACCAACTTTATTGtt       <  1:67576/150‑1 (MQ=255)
                                                                                                                                                 tctcAAATATCTTTAATAGCACCAATAACCAGCCAGGGCAGATTGTTGTTTTACCGGAAAAATCCGGCTGGGTAGGTGTCTCAGCAATTTGTCCACCCGGTACGCTGGTGAATTATACATACCGTAGTTATGTCACCAACTTTATTGTTc      <  1:2824874/150‑1 (MQ=255)
                                                                                                                                                 tctcAAATATCTTTAATAGCACCAATAACCAGCCAGGGCAGATTGTTGTTTTACCGGAAAAATCCGGCTGGGTAGGTGTCTCAGCAATTTGTCCACCCGGTACGCTGGTGAATTATACATACCGTAGTTATGTCACCAACTTTATTGTTc      >  1:330707/1‑150 (MQ=255)
                                                                                                                                                  ctcAAATATCTTTAATAGCACCAATAACCAGCCAGGGCAGATTGTTGTTTTACCGGAAAAATCCGGCTGGGTAGGTGTCTCAGCAATTTGTCCACCCGGTACGCTGGTGAATTATACATACCGTAGTTATGTCACCAACTTTATTGTTCa     >  1:1756851/1‑150 (MQ=255)
                                                                                                                                                  ctcAAATATCTTTAATAGCACCAATAACCAGCCAGGGCAGATTGTTGTTTTACCGGAAAAATCCGGCTGGGTAGGTGTCTCAGCAATTTGTCCACCCGGTACGCTGGTGAATTATACATACCGTAGTTATGTCACCAACTTTATTGTTCa     <  1:1890284/150‑1 (MQ=255)
                                                                                                                                                   tcAAATATCTTTAATAGCACCAATAACCAGCCAGGGCAGATTGTTGTTTTACCGGAAAAATCCGGCTGGGTAGGTGTCTCAGCAATTTGTCCACCCGGTACGCTGGTGAATTATACATACCGTAGTTATGTCACCAACTTTATTGTTCAg    >  2:1290812/1‑150 (MQ=255)
                                                                                                                                                     aaaTATCTTTAATAGCACCAATAACCAGCCAGGGCAGATTGTTGTTTTACCGGAAAAATCCGGCTGGGTAGGTGTCTCAGCAATTTGTCCACCCGGTACGCTGGTGAATTATACATACCGTAGTTATGTCACCAACTTTATTGTTCAGGa  >  2:4640312/1‑150 (MQ=255)
                                                                                                                                                     aaaTATCTTTAATAGCACCAATAACCAGCCAGGGCAGATTGTTGTTTTACCGGAAAAATCCGGCTGGGTAGGTGTCTCAGCAATTTGTCCACCCGGTACGCTGGTGAATTATACATACCGTAGTTATGTCACCAACTTTATTGTTCAGGa  <  1:4280157/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CACACATCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCAGATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTAATAGCACCAATAACCAGCCAGGGCAGATTGTTGTTTTACCGGAAAAATCCGGCTGGGTAGGTGTCTCAGCAATTTGTCCACCCGGTACGCTGGTGAATTATACATACCGTAGTTATGTCACCAACTTTATTGTTCAGGA  >  CP009273/557738‑558036

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