Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,992,861 C→A 7.0% T130T (ACG→ACT ygeR ← LysM domain‑containing M23 family putative peptidase; septation lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,992,8610CA 7.0% 739.0 / 19.6 199T130T (ACG→ACTygeRLysM domain‑containing M23 family putative peptidase; septation lipoprotein
Reads supporting (aligned to +/- strand):  ref base C (87/98);  new base A (6/8);  total (93/106)
Fisher's exact test for biased strand distribution p-value = 7.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

Only 100 of 200 total aligned reads displayed.

AGATTACCGTAGCCACGCAGCTGGTTGCCCACATACACCACCTTTCCTGCACCCGCGGCGTAAATAGGTGTACCCCGTGGAGCTGAGATATCAATCCCTTTATTGCCGCCATCTGCTGTCGAATACGGCATGATAACTTTCCCTGTCGTTGGCCATAACCAACAACGTTGCCCTACTGGCGGCCAGGATGATTTCGGTACCGCTGATGACGGTGTAACCGATGCGGTTTTGGTCGTTGATTTGGCGGTTGATTTACGTGTAATACTGCTACTTTTCGCCCCACCCAGTTTTAGT  >  CP009273/2992715‑2993008
                                                                                                                                                  |                                                                                                                                                   
agaTTACCGTAGCCACGCAGCTGGTTGCCCACATACACCACCTTTCCTGCACCCGCGGCGTAAATAGGTGTACCCCGTGGAGCTGAGATATCAATCCCTTTATTGCCGCCATCTGCTGTCGAATACGGCATGATAACTTTCCCTGTCg                                                                                                                                                    >  1:2280638/1‑148 (MQ=255)
agaTTACCGTAGCCACGCAGCTGGTTGCCCACATACACCACCTTTCCTGCACCCGCGGCGTAAATAGGTGTACCCCGTGGAGCTGAGATATCAATCCCTTTATTGCCGCCATCTGCTGTCGAATACGGCATGATAACTTTCCCTGTCg                                                                                                                                                    >  1:619192/1‑148 (MQ=255)
agaTTACCGTAGCCACGCAGCTGGTTGCCCACATACACCACCTTTCCTGCACCCGCGGCGTAAATAGGTGTACCCCGTGGAGCTGAGATATCAATCCCTTTATTGCCGCCATCTGCTGTCGAATACGGCATGATAACTTTCCCTGTCg                                                                                                                                                    >  1:1125394/1‑148 (MQ=255)
  aTTACCGTAGCCACGCAGCTGGTTGCCCACATACACCACCTTTCCTGCACCCGCGGCGTAAATAGGTGTACCCCGTGGAGCTGAGATATCAATCCCTTTATTGCCGCCATCTGCTGTCGAATACGGCATGATAACTTTCCCTGTCGtt                                                                                                                                                  <  1:1214002/148‑1 (MQ=255)
  aTTACCGTAGCCACGCAGCTGGTTGCCCACATACACCACCTTTCCTGCACCCGCGGCGTAAATAGGTGTACCCCGTGGAGCTGAGATATCAATCCCTTTATTGCCGCCATCTGCTGTCGAATACGGCATGATAACTTTCCCTGTCGtt                                                                                                                                                  <  2:3065312/148‑1 (MQ=255)
     aCCGTAGCCACGCAGCTGGTTGCCCACATACACCACCTTTCCTGCACCCGCGGCGTAAATAGGTGTACCCCGTGGAGCTGAGATATCAATCCCTTTATTGCCGCCATCTGCTGTCGAATACGGCATGATAACTTTCCCTGTCGTTGGc                                                                                                                                               >  1:821342/1‑148 (MQ=255)
     aCCGTAGCCACGCAGCTGGTTGCCCACATACACCACCTTTCCTGCACCCGCGGCGTAAATAGGTGTACCCCGTGGAGCTGAGATATCAATCCCTTTATTGCCGCCATCTGCTGTCGAATACGGCATGATAACTTTCCCTGTCGTTGGc                                                                                                                                               <  2:1014526/148‑1 (MQ=255)
       cGTAGCCACGCAGCTGGTTGCCCACATACACCACCTTTCCTGCACCCGCGGCGTAAATAGGTGTACCCCGTGGAGCTGAGATATCAATCCCTTTATTGCCGCCATCTGCTGTCGAATACGGCATGATAACTTTCCCTGTCGTTGGCCa                                                                                                                                             <  1:1970767/148‑1 (MQ=255)
               cGCAGCTGGTTGCCCACATACACCACCTTTCCTGCACCCGCGGCGTAAATAGGTGTACCCCGTGGAGCTGAGATATCAATCCCTTTATTGCCGCCATCTGCTGTCGAATACGGCATGATAACTTTCCCTGTCGTTGGCCATAACcaac                                                                                                                                     >  1:431177/1‑148 (MQ=255)
                gCAGCTGGTTGCCCACATACACCACCTTTCCTGCACCCGCGGCGTAAATAGGTGTACCCCGTGGAGCTGAGATATCAATCCCTTTATTGCCGCCATCTGCTGTCGAATACGGCATGATAACTTTCCCTGTCGTTGGCCATAACcaaca                                                                                                                                    <  1:371946/148‑1 (MQ=255)
                gCAGCTGGTTGCCCACATACACCACCTTTCCTGCACCCGCGGCGTAAATAGGTGTACCCCGTGGAGCTGAGATATCAATCCCTTTATTGCCGCCATCTGCTGTCGAATACGGCATGATAACTTTCCCTGTCGTTGGCCATAACcaaca                                                                                                                                    >  1:288376/1‑148 (MQ=255)
                      ggTTGCCCACATACACCACCTTTCCTGCACCCGCGGCGTAAATAGGTGTACCCCGTGGAGCTGAGATATCAATCCCTTTATTGCCGCCATCTGCTGTCGAATACGGCATGATAACTTTCCCTGTCGTTGGCCATAACCAACAACGTTg                                                                                                                              >  2:2249060/1‑148 (MQ=255)
                        ttGCCCACATACACCACCTTTCCTGCACCCGCGGCGTAAATAGGTGTACCCCGTGGAGCTGAGATATCAATCCCTTTATTGCCGCCATCTGCTGTCGAATACGGCATGATAACTTTCCCTGTCGTTGGCCATAACCAACAACGTTGcc                                                                                                                            <  2:1464481/148‑1 (MQ=255)
                           cccACATACACCACCTTTCCTGCACCCGCGGCGTAAATAGGTGTACCCCGTGGAGCTGAGATATCAATCCCTTTATTGCCGCCATCTGCTGTCGAATACGGCATGATAACTTTCCCTGTCGTTGGCCATAACCAACAACGTTGCCCTa                                                                                                                         <  1:2646928/148‑1 (MQ=255)
                             cacaTACACCACCTTTCCTGCACCCGCGGCGTAAATAGGTGTACCCCGTGGAGCTGAGATATCAATCCCTTTATTGCCGCCATCTGCTGTCGAATACGGCATGATAACTTTCCCTGTCGTTGGCCATAACCAACAACGTTGCCCTACt                                                                                                                       <  2:2331203/148‑1 (MQ=255)
                                  acacCACCTTTCCTGCACCCGCGTCGTATATAGGTGTACCCCGTGGAGCTGAGATATCAATCCCTTTATTGCCGCCATCTGCTGTCGAATACGGCATGATAACTTTCCCTGTCGTTGTCCATAACCAACAACGTTGCCCTACTggcgg                                                                                                                  <  2:159026/148‑1 (MQ=255)
                                  acacCACCTTTCCTGCACCCGCGGCGTAAATAGGTGTACCCCGTGGAGCTGAGATATCAATCCCTTTATTGCCGCCATCTGCTGTCGAATACGGCATGATAACTTTCCCTGTCGTTGGCCATAACCAACAACGTTGCCCTACTggcgg                                                                                                                  >  1:255392/1‑148 (MQ=255)
                                   caccacCTTTCCTGCACCCGCGGCGTAAATAGGTGTACCCCGTGGAGCTGAGATATCAATCCCTTTATTGCCGCCATCTGCTGTCGAATACGGCATGATAACTTTCCCTGTCGTTGGCCATAACCAACCACGTTGCCCTACAggcggc                                                                                                                 >  1:1589265/1‑148 (MQ=255)
                                   caccacCTGTCCTGCACCCGCGGCGTAAATAGGTGTACCCCGTGGAGCTGAGATATCAATCCCTTTATTGCCGCCATCTGCTGTCGAATACGGCATGATAACTTTCCCTGTCGTTGGCCATAACCAACAACGTTGCCCTACTggcggc                                                                                                                 <  1:88928/148‑1 (MQ=255)
                                    accaccTTTCCTGCACCCGCGGCGTAAATAGGTGTACCCCGTGGAGCTGAGATATCAATCCCTTTATTGCCGCCATCTGCTGTCGAATACGGCATGATAACTTTCCCTGTCGTTGGCCATAACCAACAACGTTGCCCTACTGGCGGcc                                                                                                                <  2:1153065/148‑1 (MQ=255)
                                    accaccTGTCCTGCACCCGCGGCGTAAATAGGTGTACCCCGTGGAGCTGAGATATCAATCCCTTTATTGCCGCCATCTGCTGTCGAATACGGCATGATAACTTTCCCTGTCGTTGGCCATAACCAACAACGTTGCCCTACTGGCGGcc                                                                                                                <  2:2809387/148‑1 (MQ=255)
                                       accTTTCCTGCACCCGCGGCGTAAATAGGTGTACCCCGTGGAGCTGAGATATCAATCCCTTTATTGCCGCCATCTGCTGTCGAATACGGCATGATAACTTTCCCTGTCGTTGGCCATAACCAACAACGTTGCCCTACTGGCGGCCAgg                                                                                                             >  1:737532/1‑148 (MQ=255)
                                       accTTTCCTGCACCCGCGGCGTAAATAGGTGTACCCCGTGGAGCTGAGATATCAATCCCTTTATTGCCGCCATCTGCTGTCGAATACGGCATGATAACTTTCCCTGTAGTTGGCCATAACCAACAACGTTGCCCTACTGGCGGCCAgg                                                                                                             <  1:3051474/148‑1 (MQ=255)
                                        ccTTTCCTGCACCCGCGGCGTAAATAGGTGTACCCCGTGGAGCTGAGATATCAATCCCTTTATTGCCGCCATCTGCTGTCGAATACGGCATGATAACTTTCCCTGTCGTTGGCCATAACCAACAACGTTGCCCTACTGGCGGCCAGga                                                                                                            <  1:789508/148‑1 (MQ=255)
                                         cTTTCCTGCACCCGCGGCGTAAATAGGTGTACCCCGTGGAGCTGAGATATCAATCCCTTTATTGCCGCCATCTGCTGTCGAATACGGCATGATAACTTTCCCTGTCGTTGGCCATAACCAACAACGTTGCCCTACTGGCGGCCAGgat                                                                                                           >  2:2376492/1‑148 (MQ=255)
                                          tttCCTGCACCCGCGGCGTAAATAGGTGTACCCCGTGGAGCTGAGATATCAATCCCTTTATTGCCGCCATCTGCTGTCGAATACGGCATGATAACTTTCCCTGTCGTTGGCCATAACCAACAACGTTGCCCTACTGGCGGCCAGgatg                                                                                                          >  2:1951776/1‑148 (MQ=255)
                                          tttCCTGCACCCGCGGCGTAAATAGGTGTACCCCGTGGAGCTGAGATATCAATCCCTTTATTGCCGCCATCTGCTGTCGAATACGGCATGATAACTTTCCCGGTAGTTGGCCATAACCAACAACGTTGCCCTACTGGCGGCCAGgatg                                                                                                          <  2:3202791/148‑1 (MQ=255)
                                           ttCCTGCACCCGCGGCGTAAATAGGTGTACCCCGTGGAGCTGAGATATCAATCCCTTTATTGCCGCCATCTGCTGTCGAATACGGCATGATAACTTTCCCTGTCGTTGGCCATAACCAACAACGTTGCCCTACTGGCGGCCAGgatga                                                                                                         <  1:2095365/148‑1 (MQ=255)
                                           ttCCTGCACCCGCGGCGTAAATAGGTGTACCCCGTGGAGCTGAGATATCAATCCCTTTATTGCCGCCATCTGCTGTCGAATACGGCATGATAACTTTCCCTGTCGTTGGCCATAACCAACAACGTTGCCCTACTGGCGGCCAGgatga                                                                                                         <  2:1989337/148‑1 (MQ=255)
                                                 cACCCGCGGCGTAAATAGGTGTACCCCGTGGAGCTGAGATATCAATCCCTTTATTGCCGCCATCTGCTGTCGAATACGGCATGATAACTTTCCCTGTCGTTGGCCATAACCAACAACGTTGCCCTACTGGCGGCCAGGATGATTTCgg                                                                                                   >  2:209167/1‑148 (MQ=255)
                                                   cccGCGGCGTAAATAGGTGTACCCCGTGGAGCTGAGATATCAATCCCTTTATTGCCGCCATCTGCTGTCGAATACGGCATGATAACTTTCCCTGTCGTTGGCCATAACCAACAACGTTGCCCTACTGGCGGCCAGGATGATTTCGGTa                                                                                                 <  2:1410052/148‑1 (MQ=255)
                                                    ccGCGGCGTAAATAGGTGTACCCCGTGGAGCTGAGATATCAATCCCTTTATTGCCGCCATCTGCTGTCGAATACGGCATGATAACTTTCCCTGTAGTTGGCCATAACCAACAACGTTGCCCTACTGGCGGCCAGGATGATTTCGGTAc                                                                                                >  1:829290/1‑148 (MQ=255)
                                                    ccGCGGCGTAAATAGGTGTACCCCGTGGAGCTGAGATATCAATCCCTTTATTGCCGCCATCTGCTGTCGAATACGGCATGATAACTTTCCCTGTAGTTGGCCATAACCAACAACGTTGCCCTACTGGCGGCCAGGATGATTTCGGTAc                                                                                                <  1:1996038/148‑1 (MQ=255)
                                                         gcgTAAATAGGTGTACCCCGTGGAGCTGAGATATCAATCCCTTTATTGCCGCCATCTGCTGTCGAATACGGCATGATAACTTTCCCTGTCGTTGGCCATAACCAACAACGTTGCCCTACTGGCGGCCAGGATGATTTCGGTACCGCtg                                                                                           <  1:2729108/148‑1 (MQ=255)
                                                            tAAATAGGTGTACCCCGTGGAGCTGAGATATCAATCCCTTTATTGCCGCCATCTGCTGTCGAATACGGCATGATAACTTTCCCTGTCGTTGGCCATAACCAACAACGTTGCCCTACTGGCGGCCAGGATGATTTCGGTACCGCtgatg                                                                                        <  2:2475998/148‑1 (MQ=255)
                                                             aaaTAGGTGTACCCCGTGGAGCTGAGATATCAATCCCTTTATTGCCGCCATCTGCTGTCGAATACGGCATGATAACTTTCCCTGTCGTTGGCCATAACCAACAACGTTGCCCTACTGGCGGCCAGGATGATTTCGGTACCGCtgatga                                                                                       >  1:1619325/1‑148 (MQ=255)
                                                                 aGGTGTACCCCGTGGAGCTGAGATATCAATCCATTTATTGCCGCCATCTGCTGTCGAATACGGCATGATAACTTTCCCTGTCGTTGGCCATAACCAACAACGTTGCCCTACTGGCGGCCAGGATGATTTCGGTACCGCTGATGACGgt                                                                                   >  2:780599/1‑148 (MQ=255)
                                                                   gtgtACCCCGTGGAGCTGAGATATCAATCCCTTTATTGCCGCCATCTGCTGTCGAATACGGCATGATAACTTTCCCTGTCGTTGGCCATAACCAACAACGTTGCCCTACTGGCGGCCAGGCTGATTTCGGTACCGCTGATGACGgtgt                                                                                 >  1:2139216/1‑148 (MQ=255)
                                                                   gtgtACCCCGTGGAGCTGAGATATCAATCCCTTTATTGCCGCCATCTGCTGTCGAATACGGCATGATAACTTTCCCTGTCGTTGGCCATAACCAACAACGTTGCCCTACTGGCGGCCAGGATGATTTCGGTACCGCTGATGACGgtgt                                                                                 >  1:364236/1‑148 (MQ=255)
                                                                     gtACCCCGTGGAGCTGAGATATCAATCCCTTTATTGCCGCCATCTGCTGTCGAATACGGCATGATAACTTTCCCTGTCGTTGGCCATAACCAACAACGTTGCCCTACTGGCGGCCAGGATGATTTCGGTACCGCTGATGACGGTGTaa                                                                               >  2:2545129/1‑148 (MQ=255)
                                                                        ccccGTGGAGCTGAGATATCAATCCCTTTATTGCCGCCATCTGCTGTCGAATACGGCATGATAACTTTCCCTGTCGTTGGCCATAACCAACAACGTTGCCCTACTGGCGGCCAGGATGATTTCGGTACCGCTGATGACGGTGTAACCg                                                                            <  2:2121419/148‑1 (MQ=255)
                                                                        ccccGTGGAGCTGAGATATCAATCCCTTTATTGCCGCCATCTGCTGTCGAATACGGCATGATAACTTTCCCTGTCGTTGGCCATAACCAACAACGTTGCCCTACTGGCGGCCAGGATGATTTCGGTACCGCTGATGACGGTGTAACCg                                                                            >  2:3176389/1‑148 (MQ=255)
                                                                        ccccGTGGAGCTGAGATATCAATCCCTTTATTGCCGCCATCTGCTGTCGAATACGGCATGATAACTTTCCCTGTCGTTGGCCATAACCAACAACGTTGCCCTACTGGCGGCCAGGATGATTTCGGTACCGCTGATGACGGTGTAACCg                                                                            <  1:2366941/148‑1 (MQ=255)
                                                                         cccGTGGAGCTGAGATATCAATCCCTTTATTGCCGCCATCTGCTGTCGAATACGGCATGATAACTTTCCCTGTCGTTGGCCATAACCAACAACGTTGCCCTACTGGCGGCCAGGATGATTTCGGTACCGCTGATGACGGTGTAACCGa                                                                           <  1:1903008/148‑1 (MQ=255)
                                                                          ccGTGGAGCTGAGATATCAATCCCTTTATTGCCGCCATCTGCTGTCGAATACGGCATGATAACTTTCCCTGTCGTTGGCCATAACCAACAACGTTGCCCTACTGGCGGCCAGTATGATTTCGGTACCGCTGATGACGGTGTAACCGAt                                                                          <  2:2763471/148‑1 (MQ=255)
                                                                          ccGTGGAGCTGAGATATCAATCCCTTTATTGCCGCCATCTGCTGTCGAATACGGCATGATAACTTTCCCTGTCGTTGGCCATAACCAACAACGTTGCCCTACTGGCGGCCAGGATGATTTCGGTACCGCTGATGACGGTGTAACCGAt                                                                          <  2:1111402/148‑1 (MQ=255)
                                                                           cGTGGAGCTGAGATATCAATCCCTTTATTGCCGCCATCTGCTGTCGAATACGGCATGATAACTTTCCCTGTCGTTGGCCATAACCAACAACGTTGCCCTACTGGCGGCCAGGATGATTTCGGTACCGCTGATGACGGTGTAACCGATg                                                                         <  1:2420311/148‑1 (MQ=255)
                                                                              ggAGCTGAGATATCAATCCCTTTATTGCCGCCATCTGCTGTCGAATACGGCATGATAACTTTCCCTGTCGTTGGCCATAACCAACAACGTTGCCCTACTGGCGGCCAGGATGATTTCGGTACCGCTGATGACGGTGTAACCGATGCgg                                                                      <  1:1959511/148‑1 (MQ=255)
                                                                               gAGCTGAGATATCAATCCCTTTATTGCCGCCATCTGCTGTCGAATACGGCATGATAACTTTCCCTGTCGTTGGCCATAACCAACAACGTTGCCCTACTGGCGGCCAGGATGATTTCGGTACCGCTGATGACGGTGTAACCGATGCGGt                                                                     <  2:1984689/148‑1 (MQ=255)
                                                                               gAGCTGAGATATCAATCCCTTTATTGCCGCCATCTGCTGTCGAATACGGCATGATAACTTTCCCTGTCGTTGGCCATAACCAACAACGTTGCCCTACTGGCGGCCAGGATGATTTCGGTACCGCTGATGACGGTGTAACCGATGCGGt                                                                     <  2:319033/148‑1 (MQ=255)
                                                                                  cTGATATATCAATCCCTTTATTGCCGCCATCTGCTGTCGAATACGGCATGATAACTTTCCCTGTCGTTGGCCATAACCAACAACGTTGCCCTACTGGCGGCCAGGATGATTTCGGTACCGCTGATGACGGTGTAACCGATGCGGtttt                                                                  <  2:645482/148‑1 (MQ=255)
                                                                                    gagaTATCAATCCCTTTATTGCCGCCATCTGCTGTCGAATACGGCATGATAACTTTCCCTGTCGTTGGCCATAACCAACAACGTTGCCCTACTGGCGGCCAGGATGATTTCGGTACCGCTGATGACGGTGTAACCGATGCGGTTTTgg                                                                >  2:957482/1‑148 (MQ=255)
                                                                                     agaTATCAATCCCTTTATTGCCGCCATCTGCTGTCGAATACGGCATGATAACTTTCCCTGTCGTTGGCCATAACCAACAACGTTGCCCTACTGGCGGCCAGGCTGATTTCGGTACCGCTGATGACGGTGTAACCGATGCGGTTTTGGt                                                               >  1:1966118/1‑148 (MQ=255)
                                                                                     agaTATCAATCCCTTTATTGCCGCCATCTGCTGTCGAATACGGCATGATAACTTTCCCTGTCGTTGGCCATAACCAACAACGTTGCCCTACTGGCGGCCAGGATGATTTCGGTACCGCTGATGACGGTGTAACCGATGCGGTTTTGGt                                                               <  1:622659/148‑1 (MQ=255)
                                                                                     agaTATCAATCCCTTTATTGCCGCCATCTGCTGTCGAATACGGCATGATAACTTTCCCTGTCGTTGGCCATAACCAACAACGTTGCCCTACTGGCGGCCAGGATGATTTCGGTACCGCTGATGACGGTGTAACCGATGCGGTTTTGGt                                                               >  1:1237113/1‑148 (MQ=255)
                                                                                     agaTATCAATCCCTTTATTGCCGCCATCTGCTGTCGAATACGGCATGATAACTTTCCCTGTCGTTGGCCATAACCAACAACGTTGCCCTACTGGCGGCCAGGATGATTTCGGTACCGCTGATGACGGTGTAACCGATGCGGTTTTGGt                                                               <  1:2285085/148‑1 (MQ=255)
                                                                                     agaTATCAATCCCTTTATTGCCGCCATCTGCTGTCGAATACGGCATGATAACTTTCCCTGTCGTTGGCCATAACCAACAACGTTGCCCTACTGGCGGCCAGGATGATTTCGGTACCGCTGATGACGGTGTAACCGATGCGGTTTTGGt                                                               >  2:3313857/1‑148 (MQ=255)
                                                                                     agaTATCAATCCCTTTATTGCCGCCATCTGCTGTCGAATACGGCATGATAACTTTCCCTGTCGTTGGCCATAAACAACAACGTTGACCTAATGGCGGCCAGGATGATTTCGGTACCGCTGATGACGGTGTAACCGATGCGGTTTTGGt                                                               >  1:1238627/1‑148 (MQ=255)
                                                                                      gaTATCAATCCCTTTATTGCCGCCATCTGCTGTCGAATACGGCATGATAACTTTCCCTGTCGTTGGCCATAACCAACAACGTTGCCCTACTGGCGGCCAGGATGATTTCGGTACCGCTGATGACGGTGTAACCGATGCGGTTTTGGTc                                                              >  2:1757714/1‑148 (MQ=255)
                                                                                       atatCAATCCCTTTATTGCCGCCATCTGCTGTCGAATACGGCATGATAACTTTCCCTGTCGTTGGCCATAACCAACAACGTTGCCCTACTGGCGGCCAGGATGATTTCGGTACCGCTGATGACGGTGTAACCGATGCGGTTTTGGTCg                                                             <  2:1499893/148‑1 (MQ=255)
                                                                                           cAATCCCTTTATTGCCGCCATCTGCTGTCGAATACGGCATGATAACTTTCCCTGTCGTTGGCCATAACCAACAACGTTGCCCTACTGGCGGCCAGGATGATTTCGGTACCGCTGATGACGGTGTAACCGATGCGGTTTTGGTCGTTGa                                                         >  1:267648/1‑148 (MQ=255)
                                                                                           cAATCCCTTTATTGCCGCCATCTGCTGTCGAATACGGCATGATAACTTTCCCTGTCGTTGGCCATAACCAACAACGTTGCCCTACTGGCGGCCAGGATGATTTCGGTACCGCTGATGACGGTGTAACCGATGCGGTTTTGGTCGTTGa                                                         <  1:1870248/148‑1 (MQ=255)
                                                                                             aTCCCTTTATTGCCGCCATCTGCTGTCGAATACGGCATGATAACTTTCCCTGTCGTTGGCCATAACCAACAACGTTGCCCTACTGGCGGCCAGGATGATTTCGGTACCGCTGATGACGGTGTAACCGATGCGGTTTTGGTCGTTGAtt                                                       <  1:2415037/148‑1 (MQ=255)
                                                                                             aTCCCTTTATTGCCGCCATCTGCTGTCGAATACGGCATGATAACTTTCCCTGTCGTTGGCCATAACCAACAACGTTGCCCTACTGGCGGCCAGGATGATTTCGGTACCGCTGATGACGGTGTAACCGATGCGGTTTTGGTCGTTGAtt                                                       <  1:1344921/148‑1 (MQ=255)
                                                                                              tCCCTTTATTGCCGCCATCTGCTGTCGAATACGGCATGATAACTTTCCCTGTCGTTGGCCATAACCAACAACGTTGCCCTACTGGCGGCCAGGATGATTTCGGTACCGCTGATGACGGTGTAACCGATGCGGTTTTGGTCGTTGAttt                                                      >  1:699873/1‑148 (MQ=255)
                                                                                               cccTTTATTGCCGCCATCTGCTGTCGAATACGGCATGATAACTTTCCCTGTCGTTGGCCATAACCAACAACGTTGCCCTACTGGCGGCCAGGATGATTTCGGTACCGCTGATGACGGTGTAACCGATGCGGTTTTGGTCGTTGATTTg                                                     >  2:1054224/1‑148 (MQ=255)
                                                                                               cccTTTATTGCCGCCATCTGCTGTCGAATACGGCATGATAACTTTCCCTGTCGTTGGCCATAACCAACAACGTTGCCCTACTGGCGGCCAGGATGATTTCGGTACCGCTGATGACGGTGTAACCGATGCGGTTTTGGTCGTTGATTTg                                                     >  1:500064/1‑148 (MQ=255)
                                                                                                ccTTTATTGCCGCCATCTGCTGTCGAATACGGCATGATAACTTTCCCTGTCGTTGGCCATAACCAACAACGTTGCCCTACTGGCGGCCAGGATGATTTCGGTACCGCTGATGACGGTGTAACCGATGCGGTTTTGGTCGTTGATTTgg                                                    >  1:1689543/1‑148 (MQ=255)
                                                                                                  tttATTGCCGCCATCTGCTGTCGAATACGGCATGATAACTTTCCCTGTCGTTGGCCATAACCAACAACGTTGCCCTACTGGCGGCCAGGCTGATTTCGGTACCGCTGATGACGGTGTAACCGATGCGGTTTTGGTCGTTGATTTGGCg                                                  >  2:768361/1‑148 (MQ=255)
                                                                                                  tttATTGCCGCCATCTGCTGTCGAATACGGCATGATAACTTTCCCTGTCGTTGGCCATAACCAACAACGTTGCCCTACTGGCGGCCAGGATGATTTCGGTACCGCTGATGACGGTGTAACCGATGCGGTTTTGGTCGTTGATTTGGCg                                                  >  2:1845096/1‑148 (MQ=255)
                                                                                                   ttATTGCCGCCATCTGCTGTCGAATACGGCATGATAACTTTCCCTGTCGTTGGCCATAACCAACAACGTTGCCCTACTGGCGGCCAGGATGATTTCGGTACCGCTGATGACGGTGTAACCGATGCGGTTTTGGTCGTTGATTTGGCgg                                                 <  2:288376/148‑1 (MQ=255)
                                                                                                   ttATTGCCGCCATCTGCTGTCGAATACGGCATGATAACTTTCCCTGTCGTTGGCCATAACCAACAACGTTGCCCTACTGGCGGCCAGGATGATTTCGGTACCGCTGATGACGGTGTAACCGATGCGGTTTTGGTCGTTGATTTGGCgg                                                 >  2:544289/1‑148 (MQ=255)
                                                                                                      ttGCCGCCATCTGCTGTCGAATACGGCATGATAACTTTCCCTGTCGTTGGCCATAACCAACAACGTTGCCCTACTGGCGGCCAGGATGATTTCGGTACCGCTGATGACGGTGTAACCGATGCGGTTTTGGTCGTTGATTTGGCGGTTg                                              <  2:2678134/148‑1 (MQ=255)
                                                                                                       tGCCGCCATCTGCTGTCGAATACGGCATGATAACTTTCCCTGTCGTTGGCCATAACCAACAACGTTGCCCTACTGGCGGCCAGGATGATTTCGGTACCGCTGATGACGGTGTAACCGATGCGGTTTTGGTCGTTGATTTGGCGGTTGa                                             <  1:1935715/148‑1 (MQ=255)
                                                                                                         ccgccATCTGCTGTCGAATACGGCATGATAACTTTCCCTGTCGTTGGCCATAACCAACAACGTTGCCCTACTGGCGGCCAGGATGATTTCGGTACCGCTGATGACGGTGTAACCGATGCGGTTTTGGTCGTTGATTTGGCGGTTGAtt                                           <  1:1478293/148‑1 (MQ=255)
                                                                                                            ccATCTGCTGTCGAATACGGCATGATAACTTTCCCTGTCGTTGGCCATAACCAACAACGTTGCCCTACTGGCGGCCAGGATGATTTCGGTACCGCTGATGACGGTGTAACCGATGCGGTTTTGGTCGTTGATTTGGCGGTTGATTTAc                                        >  1:212410/1‑148 (MQ=255)
                                                                                                              aTCTGCTGTCGAATACGGCATGATAACTTTCCCTGTAGTTGGCCATAACCAACAACGTTGCCCTACTGGCGGCCAGGATGATTTCGGTACCGCTGATGACGGTGTAACCGATGCGGTTTTGGTCGTTGATTTGGCGGTTGATTTACgt                                      <  2:3163/148‑1 (MQ=255)
                                                                                                                  gctgTCGAATACGGCATGATAACTTTCCCTGTCGTTGGCCATAACCAACAACGTTGCCCTACTGGCGGCCAGGATGATTTCGGTACCGCTGATGACGGTGTAACCGATGCGGTTTTGGTCGTTGATTTGGCGGTTGATTTACGTGTaa                                  >  1:2312255/1‑148 (MQ=255)
                                                                                                                   ctgTCGAATACGGCATGATAACTTTCCCTGTCGTTGGCCATAACCAACAACGTTGCCCTACTGGCGGCCAGGATGATTTCGGTACCGCTGATGACGGTGTAACCGATGCGGTTTTGGTCGTTGATTTGGCGGTTGATTTACGTGTAAt                                 <  1:1289504/148‑1 (MQ=255)
                                                                                                                          aTACGGCATGATAACTTTCCCTGTAGTTGGCCATAACCAACAACGTTGCCCTACTGGCGGCCAGGATGATTTCGGTACCGCTGATGACGGTGTAACCGATGCGGTTTTGGTCGTTGATTTGGCGGTTGATTTACGTGTAATACTGCTa                          <  1:149893/148‑1 (MQ=255)
                                                                                                                           tACGGCATGATAACTGTCCCTGTCGTTGGCCATAACCAACAACGTTGCCCTACTGGCGGCCAGGATGATTTCGGTACCGCTGATGACGGTGTAACCGATGCGGTTTTGGTCGTTGATTTGGCGGTTGATTTACGTGTAATACTGCTAc                         <  1:3127794/148‑1 (MQ=255)
                                                                                                                             cGGCATGATAACTTTCCCTGTCGTTGGCCATAACCAACAACGTTGCCCTACTGGCGGCCAGGATGATTTCGGTACCGCTGATGACGGTGTAACCGATGCGGTTTTGGTCGTTGATTTGGCGGTTGATTTACGTGTAATACTGCTACtt                       <  1:966977/148‑1 (MQ=255)
                                                                                                                               gCATGATAACTTTCCCTGTCGTTGGCCATAACCAACAACGTTGCCCTACTGGCGGCCAGGATGATTTCGGTACCGCTGATGACGGTGTAACCGATGCGGTTTTGGTCGTTGATTTGGCGGTTGATTTACGTGTAATACTGCTACtttt                     <  2:821342/148‑1 (MQ=255)
                                                                                                                                cATGATAACTTTCCCTGTCGTTGGCCATAACCAACAACGTTGCCCTACTGGCGGCCAGGATGATTTCGGTACCGCTGATGACGGTGTAACCGATGCGGTTTTGGTCGTTGATTTGGCGGTTGATTTACGTGTAATACTGCTACTTTTc                    >  2:894869/1‑148 (MQ=255)
                                                                                                                                cATGATAACTTTCCCTGTCGTTGGCCATAACCAACAACGTTGCCCTACTGGCGGCCAGGATGATTTCGGTACCGCTGATGACGGTGTAACCGATGCGGTTTTGGTCGTTGATTTGGCGGTTGATTTACGTGTAATACTGCTACTTTTc                    >  1:1278131/1‑148 (MQ=255)
                                                                                                                                cATGATAACTTTCCCTGTCGTTGGCCATAACCAACAACGTTGCCCTACTGGCGGCCAGGATGATTTCGGTACCGCTGATGACGGTGTAACCGATGCGGTTTTGGTCGTTGATTTGGCGGTTGATTTACGTGTAATACTGCTACTTTTc                    >  2:2148025/1‑148 (MQ=255)
                                                                                                                                  tGATAACTTTCCCTGTCGTTGGCCATAACCAACAACGTTGCCCTACTGGCGGCCAGGATGATTTCGGTACCGCTGATGACGGTGTAACCGATGCGGTTTTGGTCGTTGATTTGGCGGTTGATTTACGTGTAATACTGCTACTTTTCGc                  >  2:1327593/1‑148 (MQ=255)
                                                                                                                                  tGATAACTTTCCCTGTCGTTGGCCATAACCAACAACGTTGCCCTACTGGCGGCCAGGATGATTTCGGTACCGCTGATGACGGTGTAACCGATGCGGTTTTGGTCGTTGATTTGGCGGTTGATTTACGTGTAATACTGCTACTTTTCGc                  >  2:458549/1‑148 (MQ=255)
                                                                                                                                  tGATAACTTTCCCTGTAGTTGGCCATAACCAACAACGTTGCCCTACTGGCGGCCAGGATGATTTCGGTACCGCTGATGACGGTGTAACCGATGCGGTTTTGGTCGTTGATTTGGCGGTTGATTTACGTGTAATACTGCTACTTTTCGc                  >  2:1111225/1‑148 (MQ=255)
                                                                                                                                    aTAACTTTCCCTGTCGTTGGCCATAACCAACAACGTTGCCCTACTGGCGGCCAGGATGATTTCGGTACCGCTGATGACGGTGTAACCGATGCGGTTTTGGTCGTTGATTTGGCGGTTGATTTACGTGTAATACTGCTACTTTTCGccc                <  1:2222540/148‑1 (MQ=255)
                                                                                                                                    aTAACTTTCCCTGTCGTTGGCCATAACCAACAACGTTGCCCTACTGGCGGCCAGGATGATTTCGGTACCGCTGATGACGGTGTAACCGATGCGGTTTTGGTCGTTGATTTGGCGGTTGATTTACGTGTAATACTGCTACTTTTCGccc                >  1:1262946/1‑148 (MQ=255)
                                                                                                                                     tAACTTTCCCTGTCGTTGGCCATAACCAACAACGTTGCCCTACTGGCGGCCAGGATGATTTCGGTACCGCTGATGACGGTGTAACCGATGCGGTTTTGGTCGTTGATTTGGCGGTTGATTTACGTGTAATACTGCTACTTTTCGcccc               >  2:2137551/1‑148 (MQ=255)
                                                                                                                                        cTTTCCCTGTCGTTGGCCATAACCAACAACGTTGCCCTACTGGCGGCCAGGATGATTTCGGTACCGCTGATGACGGTGTAACCGATGCGGTTTTGGTCGTTGATTTGGCGGTTGATTTACGTGTAATACTGCTACTTTTCGCcccacc            >  2:2979607/1‑148 (MQ=255)
                                                                                                                                         tttCCCTGTCGTTGGCCATAACCAACAACGTTGCCCTACTGGCGGCCAGGATGATTTCGGTACCGCTGATGACGGTGTAACCGATGCGGTTTTGGTCGTTGATTTGGCGGTTGATTTACGTGTAATACTGCTACTTTTCGCcccaccc           >  1:93659/1‑148 (MQ=255)
                                                                                                                                           tCCCTGTCGTTGGCCATAACCAACAACGTTGCCCTCCTGGCGGCCAGGATGATTTCGGTACCGCTGATGACGGTGTAACCGATGCGGTTTTGGTCGTTGATTTGGCGGTTGATTTACGTGTAATACTGCTACTTTTCGCCCCACCCAg         >  2:2265174/1‑148 (MQ=255)
                                                                                                                                            cccTGTCGTTGGCCATAACCAACAACGTTGCCCTACTGGCGGCCAGGATGATTTCGGTACCGCTGATGACGGTGTAACCGATGCGGTTTTGGTCGTTGATTTGGCGGTTGATTTACGTGTAATACTGCTACTTTTCGCCCCACCCAGt        <  1:918998/148‑1 (MQ=255)
                                                                                                                                            cccTGTCGTTGGCCATAACCAACAACGTTGCCCTACTGGCGGCCAGGATGATTTCGGTACCGCTGATGACGGTGTAACCGATGCGGTTTTGGTCGTTGATTTGGCGGTTGATTTACGTGTAATACTGCTACTTTTCGCCCCACCCAGt        <  1:2988545/148‑1 (MQ=255)
                                                                                                                                            cccTGTCGTTGGCCATAACCAACAACGTTGCCCTACTGGCGGCCAGGATGATTTCGGTACCGCTGATGACGGTGTAACCGATGCGGTTTTGGTCGTTGATTTGGCGGTTGATTTACGTGTAATACTGCTACTTTTCGCCCCACCCAGt        <  2:933972/148‑1 (MQ=255)
                                                                                                                                                gTCGTTGGCCATAACCAACAACGTTGCCCTACTGGCGGCCAGGATGATTTCGGTACCGCTGATGACGGTGTAACCGATGCGGTTTTGGTCGTTGATTTGGCGGTTGATTTACGTGTAATACTGCTACTTTTCGCCCCACCCAGTTTTa    <  2:2025500/148‑1 (MQ=255)
                                                                                                                                                  cGTTGGCCATAACCAACAACGTTGCCCTACTGGCGGCCAGGATGATTTCGGTACCGCTGATGACGGTGTAACCGATGCGGTTTTGGTCGTTGATTTGGCGGTTGATTTACGTGTAATACTGCTACTTTTCGCCCCACCCAGTTTTAGt  <  1:316317/148‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
AGATTACCGTAGCCACGCAGCTGGTTGCCCACATACACCACCTTTCCTGCACCCGCGGCGTAAATAGGTGTACCCCGTGGAGCTGAGATATCAATCCCTTTATTGCCGCCATCTGCTGTCGAATACGGCATGATAACTTTCCCTGTCGTTGGCCATAACCAACAACGTTGCCCTACTGGCGGCCAGGATGATTTCGGTACCGCTGATGACGGTGTAACCGATGCGGTTTTGGTCGTTGATTTGGCGGTTGATTTACGTGTAATACTGCTACTTTTCGCCCCACCCAGTTTTAGT  >  CP009273/2992715‑2993008

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: