Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 25,546 G→T 100% 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,5460GT97.2% 975.4 / 15.5 288D114Y (GAC→TAC) lspAprolipoprotein signal peptidase (signal peptidase II)
Reads supporting (aligned to +/- strand):  ref base G (1/7);  new base T (129/150);  total (131/157)
Fisher's exact test for biased strand distribution p-value = 7.62e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.76e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 288 total aligned reads displayed.

GCGGCTGGCAGCGTTGGTTCTTTGCCGGTATTGCGATTGGTATTAGCGTGATCCTGGCAGTGATGATGTATCGCTCGAAGGCCACGCAGAAGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCGACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAGGTTTTTTGCCTTCTAGAGCGAAA  >  CP009273/25400‑25692
                                                                                                                                                  |                                                                                                                                                  
gcggcTGGCAGCGTTGGTTCTTTGCCGGTATTGCGATTGGTATTAGCGTGATCCTGGCAGTGATGATGTATCGCTCGAAGGCCACGCAGAAGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCTAcc                                                                                                                                                 <  1:3083658/150‑1 (MQ=255)
gcggcTGGCAGCGTTGGTTCTTTGCCGGTATTGCGATTGGTATTAGCGTGATCCTGGCAGTGATGATGTATCGCTCGAAGGCCACGCAGAAGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCTAcc                                                                                                                                                 <  2:274824/150‑1 (MQ=255)
    cTGGCAGCGTTGGTTCTTTGCCGGTATTGCGATTGGTATTAGCGTGATCCTGGCAGTGATGATGTATCGCTCGAAGGCCACGCAGAAGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCTACCGCCt                                                                                                                                             <  2:2129035/150‑1 (MQ=255)
      ggCAGCGTTGGTTCTTTGCCGGTATTGCGATTGGTATTAGCGTGATCCTGGCAGTGATGATGTATCGCTCGAAGGCCACGCAGAAGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCTACCGCCtgt                                                                                                                                           >  2:1102092/1‑150 (MQ=255)
      ggCAGCGTTGGTTCTTTGCCGGTATTGCGATTGGTATTAGCGTGATCCTGGCAGTGATGATGTATCGCTCGAAGGCCACGCAGAAGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCTACCGCCtgt                                                                                                                                           >  2:1873084/1‑150 (MQ=255)
       gCAGCGTTGGTTCTTTGCCGGTATTGCGATTGGTATTAGCGTGATCCTGGCAGTGATGATGTATCGCTCGAAGGCCACGCAGAAGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCTACCGCCtgtg                                                                                                                                          <  2:913241/150‑1 (MQ=255)
       gCAGCGTTGGTTCTTTGCCGGTATTGCGATTGGTATTAGCGTGATCCTGGCAGTGATGATGTATCGCTCGAAGGCCACGCAGAAGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCTACCGCCtgtg                                                                                                                                          <  2:651821/150‑1 (MQ=255)
         aGCGTTGGTTCTTTGCCGGTATTGCGATTGGTATTAGCGTGATCCTGGCAGTGATGATGTATCGCTCGAAGGCCACGCAGAAGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGc                                                                                                                                        <  1:2134554/150‑1 (MQ=255)
         aGCGTTGGTTCTTTGCCGGTATTGCGATAGGTATTAGCGTGATCCTGGCAGTGATGATGTATCGCTCGAAGGCCACGCAGAAGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGc                                                                                                                                        >  2:116947/1‑150 (MQ=255)
           cGTTGGTTCTTTGCCGGTATTGCGATTGGTATTAGCGTGATCCTGGCAGTGATGATGTATCGCTCGAAGGCCACGCAGAAGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGCAc                                                                                                                                      >  2:79352/1‑150 (MQ=255)
           cGTTGGTTCTTTGCCGGTATTGCGATTGGTATTAGCGTGATCCTGGCAGTGATGATGTATCGCTCGAAGGCCACGCAGAAGCTAAACAATATCGCTTACGCGCTGATTATCGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGCAc                                                                                                                                      >  1:1840569/1‑150 (MQ=255)
                   cTTTGCCGGTATTGCGATTGGTATTAGCGTGATCCTGGCAGTGATGATGTATCGCTCGAAGGCCACGCAGAAGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCGt                                                                                                                              >  1:3115707/1‑150 (MQ=255)
                    tttGCCGGTATTGCGATTGGTATTAGCGTGATCCTGGCAGTGATGATGTATCGCTCGAAGGCCACGCAGAAGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCGtt                                                                                                                             <  2:2003531/150‑1 (MQ=255)
                    tttGCCGGTATTGCGATTGGTATTAGCGTGATCCTGGCAGTGATGATGTATCGCTCGAAGGCCACGCAGAAGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCGtt                                                                                                                             <  2:3938612/150‑1 (MQ=255)
                     ttGCCGGTATTGCGATTGGTATTAGCGTGATCCTGGAAGTGATGATGTATCGCTAGAAGGCCACGCAGAAGCTAAACAATATCTCTTACGCGCTGATTATTGGAGGCGAGCTGGGCAACCTGTACTACCGCCTGTGGCACGTCTTCGTTg                                                                                                                            >  1:3699177/1‑150 (MQ=255)
                       gCCGGTATTGCGATTGGTATTAGCGTGATCCTGGCAGTGATGATGTATCGCTCGAAGGCCACGCAGAAGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTc                                                                                                                          >  2:319575/1‑150 (MQ=255)
                        ccGGTATTGCGATTGGTATTAGCGTGATCCTGGCAGTGATGATGTATCGCTCGAAGGCCACGCAGAAGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCAGTTCTACCGCCTGTGGCACGGCTTCGTTGTCg                                                                                                                         >  2:42375/1‑150 (MQ=255)
                         cGGTATTGCGATTGGTATTAGCGTGATCCTGGCAGTGATGATGTATCGCTCGAAGGCCACGCAGAAGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGa                                                                                                                        >  1:1205715/1‑150 (MQ=255)
                            tATTGCGATTGGTATTAGCGTGATCCTGGCAGTGATGATGTATCGCTCGAAGGCCACGCAGAAGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGatat                                                                                                                     >  1:2324600/1‑150 (MQ=255)
                             aTTGCGATTGGTATTAGCGTGATCCTGGCAGTGATGATGTATCGCTCGAAGGCCACGCAGAAGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATg                                                                                                                    >  1:706013/1‑150 (MQ=255)
                             aTTGCGATTGGTATTAGCGTGATCCTGGCAGTGATGATGTATCGCTCGAAGGCCACGCAGAAGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATg                                                                                                                    <  1:448017/150‑1 (MQ=255)
                             aTTGCGATTGGTATTAGCGTGATCCTGGCAGTGATGATGTATCGCTCGAAGGCCACGCAGAAGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATg                                                                                                                    >  1:28301/1‑150 (MQ=255)
                                gCGATTGGTATTAGCGTGATCCTGGCAGTGATGATGTATCGCTCGAAGGCCACGCAGAAGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATc                                                                                                                 >  2:2960449/1‑150 (MQ=255)
                                 cGATTGGTATTAGCGTGATCCTGGCAGTGATGATGTATCGCTCGAAGGCCACGCAGAAGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCg                                                                                                                >  2:124501/1‑150 (MQ=255)
                                 cGATTGGTATTAGCGTGATCCTGGCAGTGATGATGTATCGCTCGAAGGCCACGCAGAAGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCg                                                                                                                <  1:3072825/150‑1 (MQ=255)
                                 cGATTGGTATTAGCGTGATCCTGGCAGTGATGATGTATCGCTCGAAGGCCACGCAGAAGATAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCg                                                                                                                <  1:112563/150‑1 (MQ=255)
                                   aTTGGTATTAGCGTGATCCTGGCAGTGATGATGTATCGCTCGAAGGCCACGCAGAAGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGAc                                                                                                              <  1:859397/150‑1 (MQ=255)
                                    ttGGTATTAGCGTGATCCTGGCAGTGATGATGTATCGCTCGAAGGCCACGCAGAAGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACt                                                                                                             >  2:4154000/1‑150 (MQ=255)
                                      ggTATTAGCGTGATCCTGGCAGTGATGATGTATCGCTCGAAGGCCACGCAGAAGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTc                                                                                                           >  2:2474999/1‑150 (MQ=255)
                                           tAGCGTGATCCTGGCAGTGATGATGTATCGCTCGAAGGCCACGCAGAAGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGt                                                                                                      <  1:2758524/150‑1 (MQ=255)
                                           tAGCGTGATCCTGGCAGTGATGATGTATCGCTCGAAGGCCACGCAGAAGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGt                                                                                                      >  1:1357026/1‑150 (MQ=255)
                                           tAGCGTGATCCTGGCAGTGATGATGTATCGCTCGAAGGCCACGCAGAAGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGt                                                                                                      >  2:2461301/1‑150 (MQ=255)
                                               gTGATCCTGGCAGTGATGATGTATCGCTCGAAGGCCACGCAGAAGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGc                                                                                                  >  2:684101/1‑150 (MQ=255)
                                                 gATCCTGGCAGTGATGATGTATCGCTCGAAGGCCACGCAGAAGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGa                                                                                                >  1:358282/1‑150 (MQ=255)
                                                     cTGGCAGTGATGATGTATCGCTCGAAGGCCACGCAGAAGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTgg                                                                                            >  1:1375567/1‑150 (MQ=255)
                                                      tGGCAGTGATGATGTATCGCTCGAAGGCCACGCAGAAGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTGGc                                                                                           >  1:3918357/1‑150 (MQ=255)
                                                      tGGCAGTGATGATGTATCGCTCGAAGGCCACGCAGAAGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTGGc                                                                                           >  2:1880212/1‑150 (MQ=255)
                                                         cAGTGATGATGTATCGCTCGAAGGCCACGCAGAAGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTGGCACt                                                                                        <  1:107838/150‑1 (MQ=255)
                                                          aGTGATGATGTATCGCTCGAAGGCCACGCAGAAGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTGGCACtt                                                                                       <  1:1062138/150‑1 (MQ=255)
                                                          aGTGATGATGTATCGCTCGAAGGCCACGCAGAAGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTGGCACtt                                                                                       <  1:12937/150‑1 (MQ=255)
                                                          aGTGATGATGTATCGCTCGAAGGCCACGCAGAAGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTGGCACtt                                                                                       >  2:531614/1‑150 (MQ=255)
                                                           gTGATGATGTATCGCTCGAAGGCCACGCAGAAGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTGGCACTTc                                                                                      >  2:195084/1‑150 (MQ=255)
                                                            tgatgaTGTATCGCTCGAAGGCCACGCAGAAGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTGGCACTTCg                                                                                     <  2:470339/150‑1 (MQ=255)
                                                            tgatgaTGTATCGCTCGAAGGCCACGCAGAAGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTGGCACTTCg                                                                                     <  2:56592/150‑1 (MQ=255)
                                                               tgatgTATCGCTCGAAGGCCACGCAGAAGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTGGCACTTCGCCa                                                                                  <  1:2434958/150‑1 (MQ=255)
                                                                    tATCGCTCGAAGGCCACGCAGAAGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTGGCACTTCGCCACCTTc                                                                             >  1:274679/1‑150 (MQ=255)
                                                                     aTCGCTCGAAGGCCACGCAGAAGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTGGCACTTCGCCACCTTCa                                                                            <  1:1748829/150‑1 (MQ=255)
                                                                      tCGCTCGAAGGCCACGCAGAAGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTGGCACTTCGCCACCTTCaa                                                                           <  2:835004/150‑1 (MQ=255)
                                                                       cGCTCGAAGGCCACGCAGAAGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTGGCACTTCGCCACCTTCAAc                                                                          <  2:1738112/150‑1 (MQ=255)
                                                                        gCTCGAAGGCCACGCAGAAGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTGGCACTTCGCCACCTTCAAcc                                                                         <  2:604971/150‑1 (MQ=255)
                                                                          tCGAAGGCCACGCAGAAGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTGGCACTTCGCCACCTTCAACCtt                                                                       >  1:53029/1‑150 (MQ=255)
                                                                           cGAAGGCCACGCAGAAGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTGGCACTTCGCCACCTTCAACCTTg                                                                      >  1:1954397/1‑150 (MQ=255)
                                                                           cGAAGGCCACGCAGAAGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTGGCACTTCGCCACCTTCAACCTTg                                                                      >  1:100001/1‑150 (MQ=255)
                                                                            gAAGGCCACGCGGAAGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTGGCACTTCGCCACCTTCAACCTTGc                                                                     >  1:1015452/1‑150 (MQ=255)
                                                                              aGGCCACGCAGAAGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTGGCACTTCGCCACCTTCAACCTTGCCg                                                                   <  2:3347733/150‑1 (MQ=255)
                                                                                gCCACGCAGAAGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTGGCACTTCGCCACCTTCAACCTTGCCGAt                                                                 <  1:955021/150‑1 (MQ=255)
                                                                                  cACGCAGAAGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTGGCACTTCGCCACCTTCAACCTTGCCGATAc                                                               >  1:4103894/1‑150 (MQ=255)
                                                                                   aCGCAGAAGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTGGCACTTCGCCACCTTCAACCTTGCCGATACt                                                              <  1:1536873/150‑1 (MQ=255)
                                                                                       aGAAGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCa                                                          >  1:3759295/1‑150 (MQ=255)
                                                                                            cTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCtgt                                                     >  2:4119702/1‑150 (MQ=255)
                                                                                            cTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCtgt                                                     >  1:787602/1‑150 (MQ=255)
                                                                                             tAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCtgtg                                                    >  1:821395/1‑150 (MQ=255)
                                                                                                 cAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCCTCTGTGTCgg                                                >  2:3655177/1‑150 (MQ=255)
                                                                                                  aaTATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGt                                               <  2:3647649/150‑1 (MQ=255)
                                                                                                   atatCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTg                                              <  2:4071612/150‑1 (MQ=255)
                                                                                                   atatCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTg                                              >  1:1469865/1‑150 (MQ=255)
                                                                                                      tCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCgg                                           >  1:3285419/1‑150 (MQ=255)
                                                                                                         cTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCAc                                        >  2:2385849/1‑150 (MQ=255)
                                                                                                         cTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCAc                                        >  2:3579953/1‑150 (MQ=255)
                                                                                                          ttACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACt                                       <  1:3974936/150‑1 (MQ=255)
                                                                                                          ttACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACt                                       <  1:3948431/150‑1 (MQ=255)
                                                                                                           tACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTg                                      <  2:1414536/150‑1 (MQ=255)
                                                                                                              gcgcTGATTATTGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGAtt                                   <  1:2925329/150‑1 (MQ=255)
                                                                                                                gcTGATTATTGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATtgt                                 >  2:220526/1‑150 (MQ=255)
                                                                                                                 cTGATTATTGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATtgtg                                >  1:1344636/1‑150 (MQ=255)
                                                                                                                 cTGATTATTGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATtgtg                                <  2:1502547/150‑1 (MQ=255)
                                                                                                                   gATTATTGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCt                              >  2:1042650/1‑150 (MQ=255)
                                                                                                                   gATTATTGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCt                              >  2:3743929/1‑150 (MQ=255)
                                                                                                                     ttattGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTgg                            >  1:2883672/1‑150 (MQ=255)
                                                                                                                       attGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGaa                          >  2:721879/1‑150 (MQ=255)
                                                                                                                       attGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGaa                          >  1:701346/1‑150 (MQ=255)
                                                                                                                        ttGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAg                         >  1:2844339/1‑150 (MQ=255)
                                                                                                                         tGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAgg                        <  2:2898844/150‑1 (MQ=255)
                                                                                                                         tGGCGGCGCGCTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAgg                        <  2:417771/150‑1 (MQ=255)
                                                                                                                           gcggcgCGCTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAGGtt                      <  2:538415/150‑1 (MQ=255)
                                                                                                                           gcggcgCGCTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAGGtt                      <  1:2519923/150‑1 (MQ=255)
                                                                                                                            cggcgCGCTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAGGttt                     <  1:3170495/150‑1 (MQ=255)
                                                                                                                             ggcgCGCTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTGGCACTTCGCCACCTTCAACCTGGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAGGtttt                    <  2:335847/150‑1 (MQ=255)
                                                                                                                              gcgcgcTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTGGCACGTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAGGttttt                   >  2:1489904/1‑150 (MQ=255)
                                                                                                                                gcgcTGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAGGTTTTTTg                 <  1:1239503/150‑1 (MQ=255)
                                                                                                                                    tGGGCAACCTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTGACACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAGGTTTTTTGCCtt             >  1:1679790/1‑150 (MQ=255)
                                                                                                                                     gggCAACCTGTTCGACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAGGTTTTTTGCCTTc            <  2:3131974/150‑1 (MQ=255)
                                                                                                                                          aCCTGTTCAACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACGGATTGTGCTGGAAGGTTTTTTGCCTTCTagag       >  1:3450024/1‑150 (MQ=255)
                                                                                                                                           ccTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAGGTTTTTTGCCTTCTAGAGc      >  1:2944461/1‑150 (MQ=255)
                                                                                                                                           ccTGTTCGACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAGGTTTTTTGCCTTCTAGAGc      <  1:2588997/150‑1 (MQ=255)
                                                                                                                                            cTGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAGGTTTTTTGCCTTCTAGAGCg     >  1:2159161/1‑150 (MQ=255)
                                                                                                                                             tGTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAGGTTTTTTGCCTTCTAGAGCGa    >  1:3104869/1‑150 (MQ=255)
                                                                                                                                              gTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAGGTTTTTTGCCTTCTAGAGCGaa   >  1:2264792/1‑150 (MQ=255)
                                                                                                                                              gTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAGGTTTTTTGCCTTCTAGAGCGaa   <  2:3471372/150‑1 (MQ=255)
                                                                                                                                               ttCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAGGTTTTTTGCCTTCTAGAGCGaaa  <  2:4045942/150‑1 (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: