Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 25,546 G→T 89.1% 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.1% 762.8 / 81.8 283D114Y (GAC→TAC) lspAprolipoprotein signal peptidase (signal peptidase II)
Reads supporting (aligned to +/- strand):  ref base G (18/14);  new base T (106/145);  total (124/159)
Fisher's exact test for biased strand distribution p-value = 1.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 283 total aligned reads displayed.

GCGGCTGGCAGCGTTGGTTCTTTGCCGGTATTGCGATTGGTATTAGCGTGATCCTGGCAGTGATGATGTATCGCTCGAAGGCCACGCAGAAGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCGACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAGGTTTTTTGCCTTCTAGAGCGAAA  >  CP009273/25400‑25692
                                                                                                                                                  |                                                                                                                                                  
gcggcTGGCAGCGTTGGTTCTTTGCCGGTATTGCGATTGGTATTAGCGTGATCCTGGCAGTGATGATGTATCGCTCGAAGGCCACGCAGAAGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCTAcc                                                                                                                                                 <  2:2949376/150‑1 (MQ=255)
gcggcTGGCAGCGTTGGTTCTTTGCCGGTATTGCGATTGGTATTAGCGTGATCCTGGCAGTGATGATGTATCGCTCGAAGGCCACGCAGAAGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCTAcc                                                                                                                                                 <  2:670191/150‑1 (MQ=255)
 cggcTGGCAGCGTTGGTTCTTTGCCGGTATTGCGATTGGTATTTGCGTGATACTGGCAGTGATGATGTATCGCTCGAAGGCCACGCAGAAGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCTACCg                                                                                                                                                <  1:4064907/150‑1 (MQ=255)
     tGGCAGCGTTGGTTCTTTGCCGGTATTGCGATTGGTATTAGCGTGATCCTGGCAGTGATGATGTATCGCTCGAAGGCCACGCAGAAGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCTACCGCCtg                                                                                                                                            >  2:44603/1‑150 (MQ=255)
       gCAGCGTTGGTTCTTTGCCGGTATTGCGATTGGTATTAGCGTGATCCTGGCAGTGATGATGTATCGCTCGAAGGCCACGCAGAAGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCTACCGCCtgtg                                                                                                                                          <  1:1111181/150‑1 (MQ=255)
        cAGCGTTGGTTCTTTGCCGGTATTGCGATTGGTATTAGCGTGATCCTGGCAGTGATGATGTATCGCTCGAAGGCCACGCAGAAGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTgg                                                                                                                                         <  2:2444554/150‑1 (MQ=255)
         aGCGTTGGTTCTTTGCCGGTATTGCGATTGGTATTAGCGTGATCCTGGCAGTGATGATGTATCGCTCGAAGGCCACGCAGAAGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGc                                                                                                                                        >  2:2311133/1‑150 (MQ=255)
         aGCGTTGGTTCTTTGCCGGTATTGCGATTGGTATTAGCGTGATCCTGGCAGTGATGATGTATCGCTCGAAGGCCACGCAGAAGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGc                                                                                                                                        >  2:1951412/1‑150 (MQ=255)
         aGCGTTGGTTCTTTGCCGGTATTGCGATTGGTATTAGCGTGATCCTGGCAGTGATGATGTATCGCTCGAAGGCCACGCAGAAGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGc                                                                                                                                        >  2:1493709/1‑150 (MQ=255)
           cGTTGGTTCTTTGCCGGTATTGCGATTGGTATTAGCGTGATCCTGGCAGTGATGATGTATCGCTCGAAGGCCACGCAGAAGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCGACCGCCTGTGGCAc                                                                                                                                      >  1:2564824/1‑150 (MQ=255)
            gTTGGTTCTTTGCCGGTATTGCGATTGGTATTAGCGTGATCCTGGCAGTGATGATGTATCGCTCGAAGGCCACGCAGAAGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGCACg                                                                                                                                     >  2:3359227/1‑150 (MQ=255)
            gTTGGTTCTTTGCCGGTATTGCGATTGGTATTAGCGTGATCCTGGCAGTGATGATGTATCGCTCGAAGGCCACGCAGAAGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGCACg                                                                                                                                     >  1:1889489/1‑150 (MQ=255)
            gTTGGTTCTTTGCCGGTATTGCGATTGGTATTAGCGTGATCCTGGCAGTGATGATGTATCGCTCGAAGGCCACGCAGAAGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGCACg                                                                                                                                     <  1:4144631/150‑1 (MQ=255)
            gTTGGTTCTTTGCCGGTATTGCGATTGGTATTAGCGTGATCCTGGCAGTGATGATGTATCGCTCGAAGGCCACGCAGAAGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGCACg                                                                                                                                     <  2:3350535/150‑1 (MQ=255)
                gTTCTTTGCCGGTATTGCGATTGGTATTAGTGTGATCCTGGCAGTGATGATGTATCGCTCGAAGGCCACGCAGAAGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGCACGGCtt                                                                                                                                 >  1:1504750/1‑150 (MQ=255)
                 ttCTTTGCCGGTATTGCGATTGGTATTAGCGTGATCCTGGCAGTGATGATGTATCGCTCGAAGGCCACGCAGAAGCTAAACAATATCGCTTACGCGCTGATTATTGTCGGCGCGCTGGGCAACCTGTTCGACCGCCTGTGGCACGGCTTc                                                                                                                                >  2:2388877/1‑150 (MQ=255)
                 ttCTTTGCCGGTATTGCGATTGGTATTAGCGTGATCCTGGCAGTGATGATGTATCGCTCGAAGGCCACGCAGAAGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTc                                                                                                                                <  1:4293173/150‑1 (MQ=255)
                  tCTTTGCCGGTATTGCGATTGGTATTAGCGTGATCCTGGCAGTGATGATGTATCGCTCGAAGGCCACGCAGAAGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCg                                                                                                                               <  1:2383455/150‑1 (MQ=255)
                   cTTTGCCGGTATTGCGATTGGTATTAGCGTGATCCTGGCAGTGATGATCTATCGCTCGAAGGCCACGCAGAAGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCGACCGCCTGTGGCACGGCTTCGt                                                                                                                              >  1:2940246/1‑150 (MQ=255)
                     ttGCCGGTTTTGCGATTGGTATTAGCGTGATCCTGGCAGTGATGATGTATCGCTCGAAGGCCACGCAGAAGCTATTCAATATCGCTTACGCGCTGATTATTGGAGGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCGTTg                                                                                                                            <  2:3249523/150‑1 (MQ=255)
                     ttGCCGGTATTGCGATTGGTATTAGCGTGATCCTGGCAGTGATGATGTATCGCTCGAAGGCCACGCAGAAGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCGTTg                                                                                                                            <  2:3250518/150‑1 (MQ=255)
                       gCCGGTATTGCGATTGGTATTAGCGTGATCCTGGCAGTGATGATGTATCGCTCGAAGGCCACGCAGAAGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTc                                                                                                                          <  2:4511933/150‑1 (MQ=255)
                         cGGTATTGCGATTGGTATTAGCGTGATCCTGGCAGTGATGATGTATCGCTCGAAGGCCACGCAGAAGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGa                                                                                                                        <  2:1574808/150‑1 (MQ=255)
                         cGGTATTGCGATTGGTATTAGCGTGATCCTGGCAGTGATGATGTATCGCTCGAAGGCCACGCAGAAGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGa                                                                                                                        >  1:3058745/1‑150 (MQ=255)
                          ggTATTGCGATTGGTATTAGCGTGATCCTGGCAGTGATGATGTATCGCTCGAAGGCCACGCAGAAGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGat                                                                                                                       >  2:2012143/1‑150 (MQ=255)
                              ttGCGATTGGTATTAGCGTGATCCTGGCAGTGATGATGTATCGCTCGAAGGCCACGCAGAAGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGa                                                                                                                   <  1:1179926/150‑1 (MQ=255)
                                     tGGTATTAGCGTGATCCTGGCAGTGATGATGTATCGCTCGAAGGCCACGCAGAAGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACtt                                                                                                            <  1:543673/150‑1 (MQ=255)
                                       gTATTAGCGTGATCCTGGCAGTGATGATGTATCGCTCGAAGGCCACGCAGAAGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCGACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCt                                                                                                          <  1:405698/150‑1 (MQ=255)
                                        tATTAGCGTGATCCTGGCAGTGATGATGTATCGCTCGAAGGCCACGCAGAAGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTa                                                                                                         <  1:1672745/150‑1 (MQ=255)
                                        tATTAGCGTGATCCTGGCAGTGATGATGTATCGCTCGAAGGCCACGCAGAAGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTa                                                                                                         <  1:1672037/150‑1 (MQ=255)
                                             gCGTGATCCTGGCAGTGATGATGTATCGCTCGAAGGCCACGCAGAAGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCg                                                                                                    <  2:613505/150‑1 (MQ=255)
                                             gCGTGATCCTGGCAGTGATGATGTATCGCTCGAAGGCCACGCAGAAGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCg                                                                                                    <  2:614162/150‑1 (MQ=255)
                                              cGTGATCCTGGCAGTGATGATGTATCGCTCGAAGGCCACGCAGAAGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCgg                                                                                                   >  1:4321332/1‑150 (MQ=255)
                                              cGTGATCCTGGCAGTGATGATGTATCGCTCGAAGGCCACGCAGAAGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCgg                                                                                                   <  2:1490366/150‑1 (MQ=255)
                                              cGTGATCCTGGCAGTGATGATGTATCGCTCGAAGGCCACGCAGAAGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCgg                                                                                                   >  2:2687278/1‑150 (MQ=255)
                                              cGTGATCCTGGCAGTGATGATGTATCGCTCGAAGGCCACGCAGAAGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCgg                                                                                                   >  2:401072/1‑150 (MQ=255)
                                                tGATCCTGGCAGTGATGATGTATCGCTCGAAGGCCACGCAGAAGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCg                                                                                                 <  2:2351975/150‑1 (MQ=255)
                                                   tCCTGGCAGTGATGATGTATCGCTCGAAGGCCACGCAGAAGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACt                                                                                              <  1:2788900/150‑1 (MQ=255)
                                                   tCCTGGCAGTGATGATGTATCGCTCGAAGGCCACGCAGAAGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACt                                                                                              >  1:1979610/1‑150 (MQ=255)
                                                     cTGGCAGTGATGATGTATCGCTCGAAGGCCACGCAGAAGCTGAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTgg                                                                                            <  2:4245043/150‑1 (MQ=255)
                                                      tGGCAGTGATGATGTATCGCTCGAAGGCCACGCAGAAGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTGGc                                                                                           >  1:1880031/1‑150 (MQ=255)
                                                       ggCAGTGATGATGTATCGCTCGAAGGCCACGCAGAAGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTGGCa                                                                                          >  1:4210188/1‑150 (MQ=255)
                                                          aGTGATGATGTATCGCTCGAAGGCCACGCAGAAGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTGGCACtt                                                                                       <  1:148637/150‑1 (MQ=255)
                                                           gTGATGATGTATCGCTCGAAGGCCACGCAGCAGCTAAACAATATCGCGTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTGGCACTTc                                                                                      <  2:2593584/150‑1 (MQ=255)
                                                           gTGATGATGTATCGCTCGAAGGCCACGCAGAAGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTGGCACTTc                                                                                      <  1:3358978/150‑1 (MQ=255)
                                                           gTGATGATGTATCGCTCGAAGGCCACGCAGAAGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTGGCACTTc                                                                                      <  1:1916293/150‑1 (MQ=255)
                                                           gTGATGATGTATCGCTCGAAGGCCACGCAGAAGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTGGCACTTc                                                                                      >  1:530127/1‑150 (MQ=255)
                                                              atgatgTATCGCTCGAAGGCCACGCAGAAGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTGGCACTTCGcc                                                                                   <  2:1277389/150‑1 (MQ=255)
                                                                gatgTATCGCTCGAAGGCCACGCAGAAGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCCACCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTGGCACTTCGCCAc                                                                                 >  1:2883695/1‑150 (MQ=255)
                                                                gatgTATCGCTCGAAGGCCACGCAGAAGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTGGCACTTCGCCAc                                                                                 <  2:3443022/150‑1 (MQ=255)
                                                                  tgTATCGCTCGAAGGCCACGCAGAAGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTGGCACTTCGCCACCt                                                                               >  2:4055024/1‑150 (MQ=255)
                                                                   gTATCGCTCGAAGGCCACGCAGAAGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTGGCACTTCGCCACCtt                                                                              <  1:2721815/150‑1 (MQ=255)
                                                                   gTATCGCTCGAAGGCCACGCAGAAGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTGGCACTTCGCCACCtt                                                                              <  1:679659/150‑1 (MQ=255)
                                                                    tATCGCTCGAAGGCCACGCAGAAGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTGGCACTTCGCCACCTTc                                                                             <  2:3705299/150‑1 (MQ=255)
                                                                        gCTCGAAGGCCACGCAGAAGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTGGCACTTCGCCACCTTCAAcc                                                                         >  2:3488449/1‑150 (MQ=255)
                                                                         cTCGAAGGCCACGCAGAAGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTGGCACTTCGCCACCTTCAACCt                                                                        >  1:2678861/1‑150 (MQ=255)
                                                                         cTCGAAGGCCACGCAGAAGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTGGCACTTCGCCACCTTCAACCt                                                                        <  1:1540651/150‑1 (MQ=255)
                                                                              aGGCCACGCAGAAGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTGGCACTTCGCCACCTTCAACCTTGCCg                                                                   >  2:1663387/1‑150 (MQ=255)
                                                                                 ccACGCAGAAGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTGGCACTTCGCCACCTTCAACCTTGCCGATa                                                                >  1:3627749/1‑150 (MQ=255)
                                                                                    cGCAGAAGCTAAACAATATCGCTTACGCGCTGATTATAGGCTGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTGGCACTTCGCCACCTTCAACCTTGCCGATACTg                                                             >  1:3689346/1‑150 (MQ=255)
                                                                                      cAGAAGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTGGCACTTCGCCACCTTCAACCTTGCCGATACTGcc                                                           >  2:77570/1‑150 (MQ=255)
                                                                                      cAGAAGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTGGCACTTCGCCACCTTCAACCTTGCCGATACTGcc                                                           >  2:2381720/1‑150 (MQ=255)
                                                                                      cAGAAGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTGGCACTTCGCCACCTTCAACCTTGCCGATACTGcc                                                           >  1:2837494/1‑150 (MQ=255)
                                                                                        gAAGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCAt                                                         <  2:448203/150‑1 (MQ=255)
                                                                                         aaGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATc                                                        >  2:3076712/1‑150 (MQ=255)
                                                                                         aaGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCGACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATc                                                        <  1:3280858/150‑1 (MQ=255)
                                                                                          aGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCt                                                       >  1:2271171/1‑150 (MQ=255)
                                                                                           gCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCtg                                                      <  1:556613/150‑1 (MQ=255)
                                                                                              aaaCAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTgtgt                                                   >  1:975653/1‑150 (MQ=255)
                                                                                              aaaCAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCGACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTgtgt                                                   >  1:1090607/1‑150 (MQ=255)
                                                                                                  aaTATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGt                                               <  1:2196002/150‑1 (MQ=255)
                                                                                                   atatCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCGACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTg                                              >  1:2091143/1‑150 (MQ=255)
                                                                                                         cTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCAc                                        <  2:4032445/150‑1 (MQ=255)
                                                                                                         cTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCAc                                        >  2:1843910/1‑150 (MQ=255)
                                                                                                          ttACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACt                                       <  1:3553567/150‑1 (MQ=255)
                                                                                                                gcTGATTATTGGCGGCGCGCTGGGCAACCTGTTCGACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATtgt                                 >  2:3423898/1‑150 (MQ=255)
                                                                                                                 cTGATTATTGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATtgtg                                <  2:491951/150‑1 (MQ=255)
                                                                                                                   gATTATTGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCt                              <  1:2552705/150‑1 (MQ=255)
                                                                                                                   gATTATTGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCt                              <  2:2966844/150‑1 (MQ=255)
                                                                                                                   gATTATTGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCt                              <  2:2201024/150‑1 (MQ=255)
                                                                                                                    attattGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTg                             <  2:2864778/150‑1 (MQ=255)
                                                                                                                    attattGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTg                             <  2:1829127/150‑1 (MQ=255)
                                                                                                                     ttattGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTgg                            <  1:1061398/150‑1 (MQ=255)
                                                                                                                       attGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGaa                          >  2:3893502/1‑150 (MQ=255)
                                                                                                                        ttGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAg                         <  1:4507391/150‑1 (MQ=255)
                                                                                                                          ggcggcGCGCTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAGGt                       >  2:2280337/1‑150 (MQ=255)
                                                                                                                            cggcgCGCTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAGGttt                     >  2:4218558/1‑150 (MQ=255)
                                                                                                                            cggcgCGCTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAGGttt                     >  2:980757/1‑150 (MQ=255)
                                                                                                                            cggcgCGCTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAGGttt                     <  1:4376224/150‑1 (MQ=255)
                                                                                                                            cggcgCGCTGGGCAACCTGTTCGACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAGGttt                     >  2:1161793/1‑150 (MQ=255)
                                                                                                                                gcgcTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAGGTTTTTTg                 >  1:4419118/1‑150 (MQ=255)
                                                                                                                                 cgcTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAGGTTTTTTGc                <  1:763979/150‑1 (MQ=255)
                                                                                                                                      ggCAACCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAGGTTTTTTGCCTTCt           >  2:1369929/1‑150 (MQ=255)
                                                                                                                                      ggCAACCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAGGTTTTTTGCCTTCt           >  2:692314/1‑150 (MQ=255)
                                                                                                                                       gCAACCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAGGTTTTTTGCCTTCTa          <  2:3175103/150‑1 (MQ=255)
                                                                                                                                         aaCCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAGGTTTTTTGCCTTCTaga        >  1:1144686/1‑150 (MQ=255)
                                                                                                                                          aCCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAGGTTTTTTGCCTTCTagag       <  1:401072/150‑1 (MQ=255)
                                                                                                                                              gTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAGGTTTTTTGCCTTCTAGAGCGaa   <  1:3113514/150‑1 (MQ=255)
                                                                                                                                              gTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAGGTTTTTTGCCTTCTAGAGCGaa   <  1:253799/150‑1 (MQ=255)
                                                                                                                                               ttCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAGGTTTTTTGCCTTCTAGAGCGaaa  >  2:313480/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
GCGGCTGGCAGCGTTGGTTCTTTGCCGGTATTGCGATTGGTATTAGCGTGATCCTGGCAGTGATGATGTATCGCTCGAAGGCCACGCAGAAGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCGACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAGGTTTTTTGCCTTCTAGAGCGAAA  >  CP009273/25400‑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: