Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 25,546 G→T 89.7% D114Y (GAC→TAC)  lspA → prolipoprotein signal peptidase (signal peptidase II)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP00927325,5460GT89.7% 866.2 / 94.6 311D114Y (GAC→TAC) lspAprolipoprotein signal peptidase (signal peptidase II)
Reads supporting (aligned to +/- strand):  ref base G (15/17);  new base T (136/143);  total (151/160)
Fisher's exact test for biased strand distribution p-value = 8.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.15e-01

Only 100 of 311 total aligned reads displayed.

GGCGGCTGGCAGCGTTGGTTCTTTGCCGGTATTGCGATTGGTATTAGCGTGATCCTGGCAGTGATGATGTATCGCTCGAAGGCCACGCAGAAGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCGACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAGGTTTTTTGCCTTCTAGAGCGAAA  >  CP009273/25399‑25692
                                                                                                                                                   |                                                                                                                                                  
ggcggcTGGCAGCGTTGGTTCTTTGCCGGTATTGCGATTGGTATTAGCGTGATCCTGGCAGTGATGATGTATCGCTCGAAGGCCACGCAGAAGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCTAc                                                                                                                                                  <  1:3319573/150‑1 (MQ=255)
   ggcTGGCAGCGTTGGTTCTTTGCCGGTATTGCGATTGGTATTAGCGTGATCCTGGCAGTGATGATGTATCGCTCGAAGGCCACGCAGAAGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCTACCGc                                                                                                                                               >  2:2888327/1‑150 (MQ=255)
         cAGCGTTGGTTCTTTGCCGGTATTGCGATTGGTATTAGCGTGATCCTGGCAGTGATGATGTATCGCTCGAAGGCCACGCAGAAGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTgg                                                                                                                                         >  1:88573/1‑150 (MQ=255)
         cAGCGTTGGTTCTTTGCCGGTATTGCGATTGGTATTAGCGTGATCCTGGCAGTGATGATGTATCGCTCGAAGGCCACGCAGAAGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTgg                                                                                                                                         <  1:3273812/150‑1 (MQ=255)
          aGCGTTGGTTCTTTGCCGGTATTGCGATTGGTATTAGCGTGATCCTGGCAGTGATGATGTATCGCTCGAAGGCCACGCAGAAGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGc                                                                                                                                        >  1:334696/1‑150 (MQ=255)
             gTTGGTTCTTTGCCGGTATTGCGATTGGTATTAGCGTGATCCTGGCAGTGATGATGTATCGCTCGAAGGCCACGCAGAAGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGCACg                                                                                                                                     >  2:2529979/1‑150 (MQ=255)
               tGGTTCTTTGCCGGTATTGCGATTGGTATTAGCGTGATCCTGGCAGTGATGATGTATCGCTCGAAGGCCACGCAGAAGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGCACGGc                                                                                                                                   <  1:4430921/150‑1 (MQ=255)
               tGGTTCTTTGCCGGTATTGCGATTGGTATTAGCGTGATCCTGGCAGTGATGATGTATCGCTCGAAGGCCACGCAGAAGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGCACGGc                                                                                                                                   >  1:2783761/1‑150 (MQ=255)
                ggTTCTTTGCCGGTATTGCGATTGGTATTAGCGTGATCCTGGCAGTGATGATGTATCGCTCGAAGGCCACGCAGAAGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGCACGGCt                                                                                                                                  <  1:3618751/150‑1 (MQ=255)
                ggTTCTTTGCCGGTATTGCGATTGGTATTAGCGTGATCCTGGCAGTGATGATGTATCGCTCGAAGGCCACGCAGAAGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGCACGGCt                                                                                                                                  <  2:2696559/150‑1 (MQ=255)
                  ttCTTTGCCGGTATTGCGATTGGTATTAGCGTGATCCTGGCAGTGATGATGTATCGCTCGAAGGCCACGCAGAAGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTc                                                                                                                                <  2:4518714/150‑1 (MQ=255)
                   tCTTTGCCGGTATTGCGATTGGTATTAGCGTGATCCTGGCAGTGATGATGTATCGCTCGAAGGCCACGCAGAAGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCg                                                                                                                               <  1:4176251/150‑1 (MQ=255)
                     tttGCCGGTATTGCGATTGGTATTAGCGTGATCCTGGCAGTGATGATGTATCGCTCGAAGGCCACGCAGAAGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCGtt                                                                                                                             <  2:313025/150‑1 (MQ=255)
                     tttGCCGGTATTGCGATTGGTATTAGCGTGATCCTGGCAGTGATGATGTATCGCTCGAAGGCCACGCAGAAGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCGtt                                                                                                                             <  2:4715749/150‑1 (MQ=255)
                          cGGTATTGCGATTGGTATTAGCGTGATCCTGGCAGTGATGATGTATCGCTCGAAGGCCACGCAGAAGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGa                                                                                                                        >  2:2546336/1‑150 (MQ=255)
                          cGGTATTGCGATTGGTATTAGCGTGATCCTGGCAGTGATGATGTATCGCTCGAAGGCCACGCAGAAGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGa                                                                                                                        >  2:1841671/1‑150 (MQ=255)
                           ggTATTGCGATTGGTATTAGCGTGATCCTGGCAGTGATGATGTATCGCTCGAAGGCCACGCAGAAGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGat                                                                                                                       <  1:3973035/150‑1 (MQ=255)
                           ggTATTGCGATTGGTATTAGCGTGATCCTGGCAGTGATGATGTATCGCTCGAAGGCCACGCAGAAGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGat                                                                                                                       <  1:2428258/150‑1 (MQ=255)
                           ggTATTGCGATTGGTATTAGCGTGATCCTGGCAGTGATGATGTATCGCTCGAAGGCCACGCAGAAGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCGACCGCCTGTGGCACGGCTTCGTTGTCGat                                                                                                                       >  1:4485999/1‑150 (MQ=255)
                              aTTGCGATTGGTATTAGCGTGATCCTGGCAGTGATGATGTATCGCTCGAAGGCCACGCAGAAGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATg                                                                                                                    >  2:837042/1‑150 (MQ=255)
                              aTTGCGATTGGTATTAGCGTGATCCTGGCAGTGATGATGTATCGCTCGAAGGCCACGCAGAAGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATg                                                                                                                    >  2:368696/1‑150 (MQ=255)
                                   gATTGGTATTAGCGTGATCCTGGCAGTGATGATGTATCGCTCGAAGGCCACGCAGAAGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGa                                                                                                               >  2:1381680/1‑150 (MQ=255)
                                      tGGTATTAGCGTGATCCTGGCAGTGATGATGTATCGCTCGAAGGCCACGCAGAAGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACtt                                                                                                            <  1:3776738/150‑1 (MQ=255)
                                         tATTAGCGTGATCCTGGCAGTGATGATGTATCGCTCGAAGGCCACGCAGAAGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTa                                                                                                         >  1:2893093/1‑150 (MQ=255)
                                         tATTAGCGTGATCCTGGCAGTGATGATGTATCGCTCGAAGGCCACGCAGAAGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTa                                                                                                         <  2:2750421/150‑1 (MQ=255)
                                         tATTAGCGTGATCCTGGCAGTGATGATGTATCGCTCGAAGGCCACGCAGAAGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTa                                                                                                         <  1:1376061/150‑1 (MQ=255)
                                           ttAGCGTGATCCTGGCAGTGATGATGTATCGCTCGAAGGCCACGCAGAAGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACg                                                                                                       <  2:4662878/150‑1 (MQ=255)
                                           ttAGCGTGATCCTGGCAGTGATGATGTATCGCTCGAAGGCCACGCAGAAGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACg                                                                                                       >  2:3760338/1‑150 (MQ=255)
                                            tAGCGTGATCCTGGCAGTGATGATGTATCGCTCGAAGGCCACGCAGAAGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGt                                                                                                      >  1:511643/1‑150 (MQ=255)
                                             aGCGTGATCCTGGCAGTGATGATGTATCGCTCGAAGGCCACGCAGAAGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTc                                                                                                     >  2:656847/1‑150 (MQ=255)
                                             aGCGTGATCCTGGCAGTGATGATGTATCGCTCGAAGGCCACGCAGAAGCTAAACAATATCGCTTACGCGCTGAGTATTGGAGGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTAAGta                                                                                                     >  2:48889/1‑149 (MQ=255)
                                              gCGTGATCCTGGCAGTGATGATGTATCGCTCGAAGGCCACGCAGAAGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCg                                                                                                    <  1:4419012/150‑1 (MQ=255)
                                                gTGATCCTGGCAGTGATGATGTATCGCTCGAAGGCCACGCAGAAGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGc                                                                                                  <  2:3873012/150‑1 (MQ=255)
                                                 tGATCCTGGCAGTGATGATGTATCGCTCGAAGGCCACGCAGAAGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCg                                                                                                 <  1:1640713/150‑1 (MQ=255)
                                                 tGATCCTGGCAGTGATGATGTATCGCTCGAAGGCCACGCAGAAGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCg                                                                                                 <  2:2139098/150‑1 (MQ=255)
                                                  gATCCTGGCAGTGCTGATGTATCGCTCGAAGGCCACGCAGAAGCTAAACAATCTCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCACCCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGa                                                                                                <  2:1354034/150‑1 (MQ=255)
                                                     ccTGGCAGTGATGATGTATCGCTCGAAGGCCACGCAGAAGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTg                                                                                             <  1:3528404/150‑1 (MQ=255)
                                                       tGGCAGTGATGATGTATCGCTCGAAGGCCACGCAGAAGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTGGc                                                                                           >  2:3282743/1‑150 (MQ=255)
                                                         gCAGTGATGATGTATCGCTCGAAGGCCACGCAGAAGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCGACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTGGCAc                                                                                         <  1:1246014/150‑1 (MQ=255)
                                                          cAGTGATGATGTATCGCTCGAAGGCCACGCAGAAGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTGGCACt                                                                                        >  2:3328942/1‑150 (MQ=255)
                                                          cAGTGATGATGTATCGCTCGAAGGCCACGCAGAAGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTGGCACt                                                                                        >  2:2993548/1‑150 (MQ=255)
                                                                 gatgTATCGCTCGAAGGCCACGCAGAAGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTGGCACTTCGCCAc                                                                                 <  1:3808137/150‑1 (MQ=255)
                                                                 gatgTATCGCTCGAAGGCCACGCAGAAGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTGGCACTTCGCCAc                                                                                 >  2:2022461/1‑150 (MQ=255)
                                                                    gTATCGCTCGAAGGCCACGCAGAAGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTGGCACTTCGCCACCtt                                                                              <  1:1030260/150‑1 (MQ=255)
                                                                    gTATCGCTCGAAGGCCACGCAGAAGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTGGCACTTCGCCACCtt                                                                              <  1:206331/150‑1 (MQ=255)
                                                                     tATCGCTCGAAGGCCACGCAGAAGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTGGCACTTCGCCACCTTc                                                                             >  1:1851324/1‑150 (MQ=255)
                                                                          cTCGAAGGCCACGCAGAAGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTGGCACTTCGCCACCTTCAACCt                                                                        >  2:4423589/1‑150 (MQ=255)
                                                                          cTCGAAGGCCACGCAGAAGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTGGCACTTCGCCACCTTCAACCt                                                                        >  2:1296927/1‑150 (MQ=255)
                                                                            cGAAGGCCACGCAGAAGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTGGCACTTCGCCACCTTCAACCTTg                                                                      <  1:1098375/150‑1 (MQ=255)
                                                                            cGAAGGCCACGCAGAAGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTGGCACTTCGCCACCTTCAACCTTg                                                                      <  1:1521228/150‑1 (MQ=255)
                                                                             gAAGGCCACGCAGAAGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTGGCACTTCGCCACCTTCAACCTTGc                                                                     >  2:3944710/1‑150 (MQ=255)
                                                                               aGGCCACGCAGAAGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTGGCACTTCGCCACCTTCAACCTTGCCg                                                                   >  1:3580983/1‑150 (MQ=255)
                                                                                 gCCACGCAGAAGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTGGCACTTCGCCACCTTCAACCTTGCCGAt                                                                 <  2:1878743/150‑1 (MQ=255)
                                                                                    aCGCAGAAGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTGGCACTTCGCCACCTTCAACCTTGCCGATACt                                                              <  2:3164804/150‑1 (MQ=255)
                                                                                     cGCAGAAGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTGGCACTTCGCCACCTTCAACCTTGCCGATACTg                                                             <  2:4620261/150‑1 (MQ=255)
                                                                                       cAGAAGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCGACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTGGCACTTCGCCACCTTCAACCTTGCCGATACTGcc                                                           <  1:3468224/150‑1 (MQ=255)
                                                                                        aGAAGCTAAACAATATCGCTTACGCGCTGATGATTGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCa                                                          <  2:1592200/150‑1 (MQ=255)
                                                                                          aaGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATc                                                        >  2:3471251/1‑150 (MQ=255)
                                                                                             cTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCtgt                                                     <  1:959148/150‑1 (MQ=255)
                                                                                             cTAAACAATATCGATTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGACGGCGACTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCtgt                                                     >  1:2439610/1‑150 (MQ=255)
                                                                                                aaCAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTc                                                  >  2:2303841/1‑150 (MQ=255)
                                                                                                   aaTATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGt                                               >  1:461905/1‑150 (MQ=255)
                                                                                                      atCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCg                                            >  2:2089151/1‑150 (MQ=255)
                                                                                                         gCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCa                                         <  2:4108359/150‑1 (MQ=255)
                                                                                                          cTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCAc                                        >  1:2914199/1‑150 (MQ=255)
                                                                                                          cTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCAc                                        >  1:484295/1‑150 (MQ=255)
                                                                                                              cgcgCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGAt                                    <  1:2649181/150‑1 (MQ=255)
                                                                                                              cgcgCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGAt                                    <  1:4367429/150‑1 (MQ=255)
                                                                                                               gcgcTGATTATTGGCGGCGCGCTGGGCAACCTGTTCGACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGAtt                                   <  2:638841/150‑1 (MQ=255)
                                                                                                                 gcTGATTATTGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATtgt                                 >  1:1463064/1‑150 (MQ=255)
                                                                                                                   tGATTATTGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGc                               >  1:757792/1‑150 (MQ=255)
                                                                                                                   tGATTATTGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGc                               <  1:906357/150‑1 (MQ=255)
                                                                                                                     attattGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTg                             >  1:1922365/1‑150 (MQ=255)
                                                                                                                      ttattGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTgg                            >  2:1944764/1‑150 (MQ=255)
                                                                                                                      ttattGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTgg                            >  2:3783297/1‑150 (MQ=255)
                                                                                                                      ttattGGCGGCGCGCTGGGCAACCTGTTCGACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTgg                            >  1:2314587/1‑150 (MQ=255)
                                                                                                                       tattGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGa                           >  2:3542060/1‑150 (MQ=255)
                                                                                                                        attGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGaa                          <  1:566400/150‑1 (MQ=255)
                                                                                                                          tGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAgg                        <  1:2986309/150‑1 (MQ=255)
                                                                                                                             cggcgCGCTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAGGttt                     <  1:1969943/150‑1 (MQ=255)
                                                                                                                              ggcgCGCTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAGGtttt                    <  2:2786816/150‑1 (MQ=255)
                                                                                                                               gcgcgcTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAGGttttt                   <  2:3461883/150‑1 (MQ=255)
                                                                                                                                cgcgcTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAGGtttttt                  <  2:397459/150‑1 (MQ=255)
                                                                                                                                cgcgcTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAGGtttttt                  <  1:1329947/150‑1 (MQ=255)
                                                                                                                                 gcgcTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAGGTTTTTTg                 <  2:1324535/150‑1 (MQ=255)
                                                                                                                                     tGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAGGTTTTTTGCCtt             >  2:1132933/1‑150 (MQ=255)
                                                                                                                                       ggCAACCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAGGTTTTTTGCCTTCt           <  2:1081414/150‑1 (MQ=255)
                                                                                                                                       ggCAACCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAGGTTTTTTGCCTTCt           <  1:4457839/150‑1 (MQ=255)
                                                                                                                                         cAACCTGTTCGACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAGGTTTTTTGCCTTCTag         >  1:1069119/1‑150 (MQ=255)
                                                                                                                                          aaCCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAGGTTTTTTGCCTTCTaga        >  2:2410571/1‑150 (MQ=255)
                                                                                                                                          aaCCTGTTCGACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTCTGCTGGAAGGTTTTTTGCCTTCTaga        >  2:3375894/1‑150 (MQ=255)
                                                                                                                                           aCCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAGGTTTTTTGCCTTCTagag       <  1:2530135/150‑1 (MQ=255)
                                                                                                                                              tGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAGGTTTTTTGCCTTCTAGAGCGa    >  2:3214350/1‑150 (MQ=255)
                                                                                                                                              tGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAGGTTTTTTGCCTTCTAGAGCGa    <  1:161056/150‑1 (MQ=255)
                                                                                                                                                ttCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAGGTTTTTTGCCTTCTAGAGCGaaa  <  2:1769313/150‑1 (MQ=255)
                                                                                                                                                ttCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAGGTTTTTTGCCTTCTAGAGCGaaa  <  1:4055788/150‑1 (MQ=255)
                                                                                                                                                ttCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAGGTTTTTTGCCTTCTAGAGCGaaa  <  1:2958545/150‑1 (MQ=255)
                                                                                                                                                ttCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAGGTTTTTTGCCTTCTAGAGCGaaa  <  2:613375/150‑1 (MQ=255)
                                                                                                                                                ttCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAGGTTTTTTGCCTTCTAGAGCGaaa  <  1:1841671/150‑1 (MQ=255)
                                                                                                                                                ttCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTAAGTCGGCGACTGGCACTTCGCCACCATCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAGGTTTTTTGCCTTCTAGAGCGaaa  >  2:4222952/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
GGCGGCTGGCAGCGTTGGTTCTTTGCCGGTATTGCGATTGGTATTAGCGTGATCCTGGCAGTGATGATGTATCGCTCGAAGGCCACGCAGAAGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCGACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAGGTTTTTTGCCTTCTAGAGCGAAA  >  CP009273/25399‑25692

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: