Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,287,028 C→A 5.8% T452T (ACG→ACT)  ccmF ← heme lyase, CcmF subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,287,0280CA 5.8% 976.0 / 21.0 267T452T (ACG→ACT) ccmFheme lyase, CcmF subunit
Reads supporting (aligned to +/- strand):  ref base C (108/142);  new base A (6/10);  total (114/153)
Fisher's exact test for biased strand distribution p-value = 7.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 267 total aligned reads displayed.

AATGCCAACAATTGTCACTGCCAGCCCAAGGTGAGCCGCCACCATCCCCCAATAACTGAAGGTGGTTTTCGTGCCGCGTGAAATACGTAGCGCAGCTTCCGCAATTGCCAGCACCGCAATCCAGCAGGCCATTGCCAGGCCGAGCACCGTCATCGCCACAACTTTGCTTTCGAACAGCCACGGCAACAGCAGCGACAGCACCAGCGTAGAGATGAAGGCGATAATCAATAAATTGCGGATCTTACGCGGGCGATCCCGCCCCCAGCGCACCAGAGGACCGACACCAAGCAGTAG  >  CP009273/2286881‑2287174
                                                                                                                                                   |                                                                                                                                                  
    ccAACAATTGTCACTGCCAGCCCAAGGTGAGCCGCCACCATCCCCCAATAACTGCAGGTGGTTTTCGTGCCGCGTGAAATACGTAGCGCAGCTTCCGCAATTGCCAGCACCGCAATCCAGCAGGCCATTGCCAGGCCGAGCACCGTCAt                                                                                                                                               <  2:238993/149‑1 (MQ=255)
     caacaaTTGTCACTGCCAGCCCAAGGTGAGCCGCCACCATCCCCCAATAACTGAAGGTGGTTTTCGTGCCGCGTGAAATACGTAGCGCAGCTTCCGCAATTGCCAGCACCGCAATCCAGCAGGCCATTGCCAGGCCGAGCACCGTCATc                                                                                                                                              <  1:562406/149‑1 (MQ=255)
       acaaTTGTCACTGCCAGCCCAAGGTGAGCCGCCACCATCCCCCAATAACTGAAGGTGGTTTTCGTGCCGCGTGAAATACGTAGCGCAGCTTCCGCAATTGCCAGCACCGCAATCCAGCAGGCCATTGCCAGGCCGAGCACAGTCATCGc                                                                                                                                            >  2:3793912/1‑149 (MQ=255)
        caaTTGTCACTGCCAGCCCAAGGTGAGCCGCCACCATCCCCCAATAACTGAAGGTGGTTTTCGTGCCGCGTGAAATACGTAGCGCAGCTTCCGCAATTGCCAGCACCGCAATCCAGCAGGCCATTGCCAGGCCGAGCACCGTCATCGcc                                                                                                                                           <  1:1993207/149‑1 (MQ=255)
         aaTTGTCACTGCCAGCCCAAGGTGAGCCGCCACCATCCCCCAATAACTGAAGGTGGTTTTCGTGCCGCGTGAAATACGTAGCGCAGCTTCCGCAATTGCCAGCACCGCAATCCAGCAGGCCATTGCCAGGCCGAGCACCGTCATCGCca                                                                                                                                          >  1:3024882/1‑149 (MQ=255)
         aaTTGTCACTGCCAGCCCAAGGTGAGCCGCCACCATCCCCCAATAACTGAAGGTGGTTTTCGTGCCGCGTGAAATACGTAGCGCAGCTTCCGCAATTGCCAGCACCGCAATCCAGCAGGCCATTGCCAGGCCGAGCACCGTCATCGCca                                                                                                                                          >  1:2365062/1‑149 (MQ=255)
          aTTGTCACTGCCAGCCCAAGGTGAGCCGCCACCATCCCCCAATAACTGAAGGTGGTTTTCGTGCCGCGTGAAATACGTAGCGCAGCTTCCGCAATTGCCAGCACCGCAATCCAGCAGGCCATTGCCAGGCCGAGCACCGTCATCGCcac                                                                                                                                         <  1:4197147/149‑1 (MQ=255)
           ttGTCACTGCCAGCCCAAGGTGAGCCGCCACCATCCCCCAATAACTGAAGGTGGTTTTCGTGCCGCGTGAAATACGTAGCGCAGCTTCCGCAATTGCCAGCACCGCAATCCAGCAGGCCATTGCCAGGCCGAGCACCGTCATCGCcaca                                                                                                                                        >  2:401511/1‑149 (MQ=255)
            tGTCACTGCCAGCCCAAGGTGAGCCGCCACCATCCCCCAATAACTGAAGGTGGTTTTCGTGCCGCGTGAAATACGTAGCGCAGCTTCCGCAATTGCCAGCACCGCAATCCAGCAGGCCATTGCCAGGCCGAGCACCGTCATCGCCACaa                                                                                                                                       <  1:2505360/149‑1 (MQ=255)
             gTCACTGCCAGCCCAAGGTGAGCCGCCACCATCCCCCAATAACTGAAGGTGGTTTTCGTGCCGCGTGAAATACGTAGCGCAGCTTCCGCAATTGCCAGCACCGCAATCCAGCAGGCCATTGCCAGGCCGAGCACCGTCATCGCCACAAc                                                                                                                                      >  1:3501851/1‑149 (MQ=255)
             gTCACTGCCAGCCCAAGGTGAGCCGCCACCATCCCCCAATAACTGAAGGTGGTTTTCGTGCCGCGTGAAATACGTAGCGCAGCTTCCGCAATTGCCAGCACCGCAATCCAGCAGGCCATTGCCAGGCCGAGCACCGTCATCGCCACAAc                                                                                                                                      >  2:1224860/1‑149 (MQ=255)
               cACTGCCAGCCCAAGGTGAGCCGCCACCATCCCCCAATAACTGAAGGTGGTTTTCGTGCCGCGTGAAATACGTAGCGCAGCTTCCGCAATTGCCAGCACCGCAATCCAGCAGGCCATTGCCAGGCCGAGCACCGTCATCGCCACAACtt                                                                                                                                    >  2:2175724/1‑149 (MQ=255)
                aCTGCCAGCCCAAGGTGAGCCGCCACCATCCCCCAATAACTGAAGGTGGTTTTCGTGCCGCGTGAAATACGTAGCGCAGCTTCCGCAATTGCCAGCACCGCAATCCAGCAGGCCATTGCCAGGCCGAGCACCGTCATCGCCACAACttt                                                                                                                                   >  2:3698290/1‑149 (MQ=255)
                aCTGCCAGCCCAAGGTGAGCCGCCACCATCCCCCAATAACTGAAGGTGGTTTTCGTGCCGCGTGAAATACGTAGCGCAGCTTCCGCAATTGCCAGCACCGCAATCCAGCAGGCCATTGCCAGGCCGAGCACAGTCATCGCCACAACttt                                                                                                                                   >  1:2654602/1‑149 (MQ=255)
                  tGCCAGCCCAAGGTGAGCCGCCACCATCCCCCAATAACTGAAGGTGGTTTTCGTGCCGCGTGAAATACGTAGCGCAGCTTCCGCAATTGCCAGCACCGCAATCCAGCAGGCCATTGCCAGGCCGAGCACCGTCATCGCCACAACTTTGc                                                                                                                                 <  2:840759/149‑1 (MQ=255)
                  tGCCAGCCCAAGGTGAGCCGCCACCATCCCCCAATAACTGAAGGTGGTTTTCGTGCCGCGTGAAATACGTAGCGCAGCTTCCGCAATTGCCAGCACCGCAATCCAGCAGGCCATTGCCAGGCCGAGCACCGTCATCGCCACAACTTTGc                                                                                                                                 <  2:555763/149‑1 (MQ=255)
                   gCCAGCCCAAGGTGAGCCGCCACCATCCCCCAATAACTGAAGGTGGTTTTCGTGCCGCGTGAAATACGTAGCGCAGCTTCCGCAATTGCCAGCACCGCAATCCAGCAGGCCATTGCCAGGCCGAGCACCGTCATCGCCACATCTTTGCt                                                                                                                                <  1:424973/149‑1 (MQ=255)
                    ccAGCCCAAGGTGAGCCGCCACCATCCCCCAATAACTGAAGGTGGTTTTCGTGCCGCGTGAAATACGTAGCGCAGCTTCCGCAATTGCCAGCACCGCAATCCAGCAGGCCATTGCCAGGCCGAGCACCGTCATCGCCACAACTTTGCtt                                                                                                                               <  1:27533/149‑1 (MQ=255)
                      aGCCCAAGGTGAGCCGCCACCATCCCCCAATAACTGAAGGTGGTTTTCGTGCCGCGTGAAATACGTAGCGCAGCTTCCGCAATTGCCAGCACCGCAATCCAGCAGGCCATTGCCAGGCCGAGCACCGTCATCGCCACAACTTTGCTTTc                                                                                                                             >  2:2658949/1‑149 (MQ=255)
                      aGCCCAAGGTGAGCCGCCACCATCCCCCAATAACTGAAGGTGGTTTTCGTGCCGCGTGAAATACGTAGCGCAGCTTCCGCAATTGCCAGCACCGCAATCCAGCAGGCCATTGCCAGGCCGAGCACCGTCATCGCCACAACTTTGCTTTc                                                                                                                             <  1:3452187/149‑1 (MQ=255)
                         ccAAGGTGAGCCGCCACCATCCCCCAATAACTGAAGGTGGTTTTCGTGCCGCGTGAAATACGTAGCGCAGCTTCCGCAATTGCCAGCACCGCAATCCAGCAGGCCATTGCCAGGCCGAGCACCGTCATCGCCACAACTTTGCTTTCGaa                                                                                                                          >  2:4429858/1‑149 (MQ=255)
                           aaGGTGAGCCGCCACCATCCCCCAATAACTGAAGGTGGTTTTCGTGCCGCGTGAAATACGTAGCGCAGCTTCCGCAATTGCCAGCACCGCAATCCAGCAGGCCATTGCCAGGCCGAGCACAGTCATCGCCACAACTTTGCTTTCGAACa                                                                                                                        <  2:3795042/149‑1 (MQ=255)
                              gTGAGCCGCCACCATCCCCCAATAACTGAAGGTGGTTTTCGTGCCGCGTGAAATACGTAGCGCAGCTTCCGCAATTGCCAGCACCGCAATCCAGCAGGCCATTGCCAGGCCGAGCACCGTCATCGCCACAACTTTGCTTTCGAACAGcc                                                                                                                     >  1:3645650/1‑149 (MQ=255)
                              gTGAGCCGCCACCATCCCCCAATAACTGAAGGTGGTTTTCGTGCCGCGTGAAATACGTAGCGCAGCTTCCGCAATTGCCAGCACCGCAATCCAGCAGGCCATTGCCAGGCCGAGCACCGTCATCGCCACAACTTTGCTTTCGAACAGcc                                                                                                                     <  1:1097995/149‑1 (MQ=255)
                                gAGCCGCCACCATCCCCCAATAACTGAAGGTGGTTTTCGTGCCGCGTGAAATACGTAGCGCAGCTTCCGCAATTGCCAGCACCGCAATCCAGCAGGCCATTGCCAGGCCGAGCACCGTCATCGCCACAACTTTGCTTTCGAACAGCCAc                                                                                                                   <  2:1238206/149‑1 (MQ=255)
                                      ccaccaTCCCCCAATAACTGAAGGTGGTTTTCGTGCCGCGTGAAATACGTAGCGCAGCTTCCGCAATTGCCAGCACCGCAATCCAGCAGGCCATTGCCAGGCCGAGCACCGTCATCGCCACAACTTTGCTTTCGAACAGCCACGGCAAc                                                                                                             >  2:1186690/1‑149 (MQ=255)
                                      ccaccaTCCCCCAATAACTGAAGGTGGTTTTCGTGCCGCGTGAAATACGTAGCGCAGCTTCCGCAATTGCCAGCACCGCAATCCAGCAGGCCATTGCCAGGCCGAGCACCGTCATCGCCACAACTTTGCTTTCGAACAGCCACGGCAAc                                                                                                             <  1:4563502/149‑1 (MQ=255)
                                       caccaTCCCCCAATAACTGAAGGTGGTTTTCGTGCCGCGTGAAATACGTAGCGCAGCTTCCGCAATTGCCAGCACCGCAATCCAGCAGGCCATTGCCAGGCCGAGCACCGTCATCGCCACAACTTTGCTTTCGAACAGCCACGGCAAca                                                                                                            >  1:866569/1‑149 (MQ=255)
                                        accaTCCCCCAATAACTGAAGGTGGTTTTCGTGCCGCGTGAAATACGTAGCGCAGCTTCCGCAATTGCCAGCACCGCAATCCAGCAGGCCATTGCCAGGCCGAGCACCGTCATCGCCACAACTTTGCTTTCGAACAGCCACGGCAAcag                                                                                                           >  1:1547160/1‑149 (MQ=255)
                                        accaTCCCCCAATAACTGAAGGTGGTTTTCGTGCCGCGTGAAATACGTAGCGCAGCTTCCGCAATTGCCAGCACCGCAATCCAGCAGGCCATTGCCAGGCCGAGCACCGTCATCGCCACAACTTTGCTTTCGAACAGCCACGGCAAcag                                                                                                           >  2:1963950/1‑149 (MQ=255)
                                               cccAATAACTGAAGGTGGTTTTCGTGCCGCGTGAAATACGTAGCGCAGCTTCCGCAATTGCCAGCACCGCAATCCAGCAGGCCATTGCCAGGCCGAGCACCGTCATCGCCACAACTTTGCTTTCGAACAGCCACGGCAACAGCAGCGAc                                                                                                    <  2:625509/149‑1 (MQ=255)
                                               cccAATAACTGAAGGTGGTTTTCGTGCCGCGTGAAATACGTAGCGCAGCTTCCGCAATTGCCAGCACCGCAATCCAGCAGGCCATTGCCAGGCCGAGCACCGTCATCGCCACAACTTTGCTTTCGAACAGCCACGGCAACAGCAGCGAc                                                                                                    <  2:498813/149‑1 (MQ=255)
                                                 cAATAACTGAAGGTGGTTTTCGTGCCGCGTGAAATACGTAGCGCAGCTTCCGCAATTGCCAGCACCGCAATCCAGCAGGCCATTGCCAGGCCGAGCACCGTCATCGCCACAACTTTGCTTTCGAACAGCCACGGCAACAGCAGCGACAg                                                                                                  <  2:4131532/149‑1 (MQ=255)
                                                   aTAACTGAAGGTGGTTTTCGTGCCGCGTGAAATACGTAGCGCAGCTTCCGCAATTGCCAGCACCGCAATCCAGCAGGCCATTGCCAGGCCGAGCACCGTCATCGCCACAACTTTGCTTTCGAACAGCCACGGCAACAGCAGCGACAGCa                                                                                                >  1:1852624/1‑149 (MQ=255)
                                                    tAACTGAAGGTGGTTTTCGTGCCGCGTGAAATACGTAGCGCAGCTTCCGCAATTGCCAGCACCGCAATCCAGCAGGCCATTGCCAGGCCGAGCACCGTCATCGCCACAACTTTGCTTTCGAACAGCCACGGCAACAGCAGCGACAGCAc                                                                                               >  2:2082446/1‑149 (MQ=255)
                                                     aaCTGAAGGTGGTTTTCGTGCCGCGTGAAATACGTAGCGCAGCTTCCGCAATTGCCAGCACCGCAATCCAGCAGGCCATTGCCAGGCCGAGCACCGTCATCGCCACAACTTTGCTTTCGAACAGCCACGGCAACAGCAGCGACAGCAcc                                                                                              >  1:1540785/1‑149 (MQ=255)
                                                      aCTGAAGGTGGTTTTCGTGCCGCGTGAAATACGTAGCGCAGCTTCCGCAATTGCCAGCACCGCAATCCAGCAGGCCATTGCCAGGCCGAGCACCGTCATCGCCACAACTTTGCTTTCGAACAGCCACGGCAACAGCAGCGACAGCACCa                                                                                             <  2:1167921/149‑1 (MQ=255)
                                                      aCTGAAGGTGGTTTTCGTGCCGCGTGAAATACGTAGCGCAGCTTCCGCAATTGCCAGCACCGCAATCCAGCAGGCCATTGCCAGGCCGAGCACCGTCATCGCCACAACTTTGCTTTCGAACAGCCACGGCAACAGCAGCGACAGCACCa                                                                                             <  2:1634710/149‑1 (MQ=255)
                                                      aCTGAAGGTGGTTTTCGTGCCGCGTGAAATACGTAGCGCAGCTTCCGCAATTGCCAGCACCGCAATCCAGCAGGCCATTGCCAGGCCGAGCACCGTCATCGCCACAACTTTGCTTTCGAACAGCCACGGCAACAGCAGCGACAGCACCa                                                                                             >  2:1319512/1‑149 (MQ=255)
                                                      aCTGAAGGTGGTTTTCGTGCCGCGTGAAATACGTAGCGCAGCTTCCGCAATTGCCAGCACCGCAATCCAGCAGGCCATTGCCAGGCCGAGCACCGTCATCGCCACAACTTTGCTTTCGAACAGCCACGGCAACAGCAGCGACAGCACCa                                                                                             >  2:4622431/1‑149 (MQ=255)
                                                       cTGAAGGTGGTTTTCGTGCCGCGTGAAATACGTAGCGCAGCTTCCGCAATTGCCAGCACCGCAATCCAGCAGGCCATTGCCAGGCCGAGCACCGTCATCGCCACAACTTTGCTTTCGAACAGCCACGGCAACAGCAGCGACAGCACCAg                                                                                            <  1:536119/149‑1 (MQ=255)
                                                        tGAAGGTGGTTTTCGTGCCGCGTGAAATACGTAGCGCAGCTTCCGCAAGTGCCAGCACCGCAATCCAGCAGGCCATTGCCAGGCCGAGCACCGTCATCGCCACAACTTTGCTTTCGAACAGCCACGGCAACAGCAGCGACAGCACCAGc                                                                                           <  2:4088142/149‑1 (MQ=255)
                                                           aGGTGGTTTTCGTGCCGCGTGAAATACGTAGCGCAGCTTCCGCAATTGCCAGCACCGCAATCCAGCAGGCCATTGCCAGGCCGAGCACCGTCATCGCCACAACTTTGCTTTCGAACAGCCACGGCAACAGCAGCGACAGCACCAGCGTa                                                                                        <  1:1337150/149‑1 (MQ=255)
                                                              tggtTTTCGTGCCGCGTGAAATACGTAGCGCAGCTTCCGCAATTGCCAGCACCGCAATCCAGCAGGCCATTGCCAGGCCGAGCACCGTCATCGCCACAACTTTGCTTTCGAACAGCCACGGCAACAGCAGCGACAGCACCAGCGTagag                                                                                     >  2:2200604/1‑149 (MQ=255)
                                                               ggtTTTCGTGCCGCGTGAAATACGTAGCGCAGCTTCCGCAATTGCCAGCACCGCAATCCAGCAGGCCATTGCCAGGCCGAGCACCGTCATCGCCACAACTTTGCTTTCGAACAGCCACGGCAACAGCAGCGACAGCACCAGCGTAgaga                                                                                    <  2:1351213/149‑1 (MQ=255)
                                                                gtTTTCGTGCCGCGTGAAATACGTAGCGCAGCTTCCGCAATTGCCAGCACCGCAATCCAGCAGGCCATTGCCAGGCCGAGCACCGTCATCGCCACAACTTTGCTTTCGAACAGCCACGGCAACAGCAGCGACAGCACCAGCGTAGAGAt                                                                                   <  2:13178/149‑1 (MQ=255)
                                                                   ttCGTGCCGCGTGAAATACGTAGCGCAGCTTCCGCAATTGCCAGCACCGCAATCCAGCAGGCCATTGCCAGGCCGAGCACCGTCATCGCCACAACTTTGCTTTCGAACAGCCACGGCAACAGCAGCGACAGCACCAGCGTAGAGATGaa                                                                                <  1:817849/149‑1 (MQ=255)
                                                                    tCGTGCCGCGTGAAATACGTAGCGCAGCTTCCGCAATTGCCAGCACCGCAATCCAGCAGACCATTGCCAGGCCGAGCACCGTCATCGCCACAACTTTGCTTTCGAACAGCCACGGCAACAGCAGCGACAGCACCAGCGTAGAGATGAAg                                                                               <  2:2446884/149‑1 (MQ=255)
                                                                      gTGCCGCGTGAAATACGTAGCGCAGCTTCCGCAATTGCCAGCACCGCAATCCAGCAGGCCATTGCCAGGCCGAGCACCGTCATCGCCACAACTTTGCTTTCGAACAGCCACGGCAACAGCAGCGACAGCACCAGCGTAGAGATGAAGGc                                                                             >  2:4498685/1‑149 (MQ=255)
                                                                       tGCCGCGTGAAATACGTAGCGCAGCTTCCGCAATTGCCAGCACCGCAATCCAGCAGGCCATTGCCAGGCCGAGCACCGTCATCGCCACAACTTTGCTTTCGAACAGCCACGGCAACAGCAGCGACAGCACCAGCGTAGAGATGAAGGCg                                                                            <  1:2372875/149‑1 (MQ=255)
                                                                         ccGCGTGAAATACGTAGCGCAGCTTCCGCAATTGCCAGCACCGCAATCCAGCAGGCCATTGCCAGGCCGAGCACGGTCATCGCCACATCTTTGCTTTCGAACCGCCACGGAAACTGCAGCGACAGCACCAGCGTAGAGATGAAGGCGAt                                                                          <  1:4377698/149‑1 (MQ=255)
                                                                           gcgTGAAATACGTAGCGCAGCTTCCGCAATTGCCAGCACCGCAATCCAGCAGGCCATTGCCAGGCCGAGCACCGTCATCGCCACAACTTTGCTTTCGAACAGCCACGGCAACAGCAGCGACAGCACCAGCGTAGAGATGAAGGCGATaa                                                                        <  1:1308453/149‑1 (MQ=255)
                                                                           gcgTGAAATACGTAGCGCAGCTTCCGCAATTGCCAGCACCGCAATCCAGCAGGCCATTGCCAGGCCGAGCACCGTCATCGCCACAACTTTGCTTTCGAACAGCCACGGCAACAGCAGCGACAGCACCAGCGTAGAGATGAAGGCGATaa                                                                        <  2:2488813/149‑1 (MQ=255)
                                                                           gcgTGAAATACGTAGCGCAGCTTCCGCAATTGCCAGCACCGAAATCCAGCAGGCCATTGCCAGGCCGAGCACAGTCATCGCCACAACTTTGCTTTCGAACAGCCACGGCAACAGCAGCGACAGCACCAGCGTAGAGATGAAGGCGATaa                                                                        <  1:2057062/149‑1 (MQ=255)
                                                                               gAAATACGTAGCGCAGCTTCCGCAATTGCCAGCACCGCAATCCAGCAGGCCATTGCCAGGCCGAGCACCGTCATCGCCACAACTTTGCTTTCGAACAGCCACGGCAACAGCAGCGACAGCACCAGCGTAGAGATGAAGGCGATAATCaa                                                                    >  2:4192432/1‑149 (MQ=255)
                                                                               gAAATACGTAGCGCAGCTTCCGCAATTGCCAGCACCGCAATCCAGCAGGCCATTGCCAGGCCGAGCACCGTCATCGCCACAACTTTGCTTTCGAACAGCCACGGCAACAGCAGCGACAGCACCAGCGTAGAGATGAAGGCGATAATCaa                                                                    >  1:2820526/1‑149 (MQ=255)
                                                                                aaaTACGTAGCGCAGCTTCCGCAATTGCCAGCACCGCAATCCAGCAGGCCATTGCCAGGCCGAGCACCGTCATCGCCACAACTTTGCTTTCGAACAGCCACGGCAACAGCAGCGACAGCACCAGCGTAGAGATGAAGGCGATAATCAAt                                                                   <  1:3365865/149‑1 (MQ=255)
                                                                                aaaTACGTAGCGCAGCTTCCGCAATTGCCAGCACCGCAATCCAGCAGGCCATTGCCAGGCCGAGCACCGTCATCGCCACAACTTTGCTTTCGAACAGCCACGGCAACAGCAGCGACAGCACCAGCGTAGAGATGAAGGCGATAATCAAt                                                                   >  2:4530671/1‑149 (MQ=255)
                                                                                 aaTACGTAGCGCAGCTTCCGCAATTGCCAGCACCGCAATCCAGCAGGCCATTGCCAGGCCGAGCACCGTCATCGCCACAACTTTGCTTTCGAACAGCCACGGCAACAGCAGCGACAGCACCAGCGTAGAGATGAAGGCGATAATCAATa                                                                  >  2:4354092/1‑149 (MQ=255)
                                                                                 aaTACGTAGCGCAGCTTCCGCAATTGCCAGCACCGCAATCCAGCAGGCCATTGCCAGGCCGAGCACCGTCATCGCCACAACTTTGCTTTCGAACAGCCACGGCAACAGCAGCGACAGCACCAGCGTAGAGATGAAGGCGATAATCAATa                                                                  <  1:2785326/149‑1 (MQ=255)
                                                                                  aTACGTAGCGCAGCTTCCGCAATTGCCAGCACCGCAATCCAGCAGGCCATTGCCAGGCCGAGCACCGTCATCGCCACAACTTTGCTTTCGAACAGCCACGGCAACAGCAGCGACAGCACCAGCGGAGAGATGAAGGCGAGAATCAATaa                                                                 >  1:3869577/1‑149 (MQ=255)
                                                                                    aCGTAGCGCAGCTTCCGCAATTGCCAGCACCGCAATCCCGCAGGCCATTGCCAGGCCGAGCACCGTCATCGCCACAACTTTGCTTTCGAACAGCCACGGCAACAGCAGCGACAGCACCAGCGTAGAGATGAAGGCGATAATCAATAAAt                                                               <  1:1244008/149‑1 (MQ=255)
                                                                                    aCGTAGCGCAGCTTCCGCAATTGCCAGCACCGCAATCCAGCAGGCCATTGCCAGGCCGAGCACCGTCATCGCCACAACTTTGCTTTCGAACAGCCACGGCAACAGCAGCGACAGCACCAGCGTAGAGATGAAGGCGATAATaaataaat                                                               >  1:4007722/1‑149 (MQ=255)
                                                                                      gTAGCGCAGCTTCCGCAATTGCCAGCACCGCAATCCAGCAGGCCATTGCCAGGCCGAGCACCGTCATCGCCACAACTTTGCTTTCGAACAGCCACGGCAACAGCCGCGACAGCACCAGCGTCGAGATGAAGGCGATAATCAATAAATTg                                                             >  2:4228575/1‑149 (MQ=255)
                                                                                       tAGCGCAGCTTCCGCAATTGCCAGCACCGCAATCCAGCAGGCCATTGCCAGGCCGAGCACCGTCATCGCCACAACTTTGCTTTCGAACAGCCACGGCAACAGCAGCGACAGCACCAGCGTAGAGATGAAGGCGATAATCAATAAATTGc                                                            <  2:680949/149‑1 (MQ=255)
                                                                                             aGCTTCCGCAATTGCCAGCACCGCAATCCAGCAGGCCATTGCCAGGCCGAGCACCGTCATCGCCACAACTTTGCTTTCGAACAGCCACGGCAACAGCAGCGACAGCACCAGCGTAGAGATGAAGGCGATAATCAATAAATTGCGGATCt                                                      <  1:3632291/149‑1 (MQ=255)
                                                                                                ttCCGCAATTGCCAGCACCGCAATCCAGCAGGCCATTGCCAGGCCGAGCACCGTCATCGCCACAACTTTGCTTTCGAACAGCCACGGCAACAGCAGCGACAGCACCAGCGTAGAGATGAAGGCGATAATCAATAAATTGCGGATCTTAc                                                   >  2:40539/1‑149 (MQ=255)
                                                                                                  ccGCAATGGCCAGCACCGCAATCCAGCAGGCCATTGCCAGGCCGAGCACCGTCATCGCCACAACTTTGCTTTCGAACAGCCACGGCAACAGCAGCGACAGCACCAGCGTAGAGATGAAGGCGATAATCAATAAATTGCGGATCTTAcgc                                                 <  1:4549254/149‑1 (MQ=255)
                                                                                                       aTTGCCAGCACCGCAATCCAGCAGGCCATTGCCAGGCCGAGCACCGTCATCGCCACAACTTTGCTTTCGAACAGCCACGGCACCAGCAGCGACAGCACCAGCGTAGAGATGAAGGCGATAATCAATAAATTGCGGATCTTACGCGGGCg                                            >  1:3135417/1‑149 (MQ=255)
                                                                                                        ttGCCAGCACCGCAATCCAGCAGGCCATTGCCAGGCCGAGCACCGTCATCGCCACAACTTTGCTTTCGAACAGCCACGGCAACAGCAGCGACAGCACCAGCGTAGAGATGAAGGCGATAATCAATAAATTGCGGATCTTACGCGGGCGa                                           <  1:3162220/149‑1 (MQ=255)
                                                                                                        ttGCCAGCACCGCAATCCAGCAGGCCAGTGCCAGGCCGAGCACCGTCATCGCCACAACTTTGCTTTCGAACCGCCACGGCAACAGCAGCGACAGCACCAGCGTAGAGATGAAGGCGATAATCAATAAATTGCGGATCTTACGCGGGCGa                                           <  1:102005/149‑1 (MQ=255)
                                                                                                          gCCAGCACCGCAATCCAGCAGGCCATTGCCAGGCCGAGCACCGTCATCGCCACAACTTTGCTTTCGAACAGCCACGGCAACAGCAGCGACAGCACCAGCGTAGAGATGAAGGCGATAATCAATAAATTGCGGATCTTACGCGGGCGATc                                         >  1:2968222/1‑149 (MQ=255)
                                                                                                           ccAGCACCGCAATCCAGCAGGCCATTGCCAGGCCGAGCACCGTCATCGCCACAACTTTGCTTTCGAACAGCCACGGCAACAGCAGCGACAGCACCAGCGTAGAGATGAAGGCGATAATCAATAAATTGCGGATCTTACGCGGGCGATcc                                        <  2:710782/149‑1 (MQ=255)
                                                                                                            cAGCACCGCAATCCCGCAGGCCATTGCCAGGCCGAGCACCGTCATCGCCACAACTTTGCTTTCGAACAGCCACGGCAACAGCAGCGACAGCACCAGCGTAGAGATGAAGGCGATAATCAATAAATTGCGGATCTTACGCGGGCGATccc                                       <  1:1676625/149‑1 (MQ=255)
                                                                                                            cAGCACCGCAATCCAGCAGGCCATTGCCAGGCCGAGCACCGTCATCGCCACAACTTAGCTTTCGAACAGCCACGGCAACAGCAGCGACAGCACCAGCGTAGAGATGAAGGAGATAATCAATAAATTGCGGATCTTACGCGGGCGATccc                                       >  2:3279696/1‑149 (MQ=255)
                                                                                                              gCACCGCAATCCAGCAGGCCATTGCCAGGCCGAGCACCGTCATCGCCACAACTTTGCTTTCGAACAGCCACGGCAACAGCAGCGACAGCACCAGCGTAGAGATGAAGGCGATAATCAATAAATTGCGGATCTTACGCGGGCGATCCCGc                                     >  2:4709815/1‑149 (MQ=255)
                                                                                                                cccGCAATCCAGCAGGCCATTGCCAGGCCGAGCACAGTCATCGCCACAACTTGGCTTTCGAACAGCCACGGCACCAGCAGCGACAGCACCAGCGTAGAGATGAAGGCGATAATCAATAAATTGCGGATCTTACGCGGGCGATCCCGccc                                   <  1:4129120/148‑1 (MQ=255)
                                                                                                                aCCGCAATCCAGCAGGCCATGGCCAGGCCGAGCACAGTCATCGCCACAACTTTGCTTTCGAACCGCCACGGCAACAGCAGCGACAGCACCAGCGTAGAGATGAAGGCGATAATCAATAAATTGCGGATCTTACGCGGGCGATCCCGccc                                   <  1:1739548/149‑1 (MQ=255)
                                                                                                                 ccGCAATCCAGCAGGCCATTGCCAGGCCGAGCACCGTCATCGCCACAACTTTGCTTTCGAACAGCCACGGCAACAGCAGCGACAGCACCAGCGTAGAGATGAAGGCGATAATCAATAAATTGCGGATCTTACGCGGGCGATCCCGcccc                                  >  1:4640746/1‑149 (MQ=255)
                                                                                                                 ccGCAATCCAGCAGGCCATTGCCAGGCCGAGCACCGTCATCGCCACAACTTTGCTTTCGAACAGCCACGGCAACAGCAGCGACAGCACCAGCGTAGAGATGAAGGCGATAATCAATAAATTGCGGATCTTACGCGGGCGATCCCGcccc                                  >  1:1270167/1‑149 (MQ=255)
                                                                                                                  cGCAATCCAGCAGGCCATTGCCAGGCCGAGCACCGTCATCGCCACAACTTTGCTTTCGAACAGCCACGGCAACAGCAGCGACAGCACCAGCGTAGAGATGAAGGCGATAATCAATAAATTGCGGATCTTACGCGGGCGATCCCGccccc                                 <  2:1109477/149‑1 (MQ=255)
                                                                                                                         cagcagGCCATTGCCAGGCCGAGCACCGTCATCGCCACAACTTTGCTTTCGAACAGCCACGGCAACAGCAGCGACAGCACCAGCGTAGAGATGAAGGCGATAATCAATAAATTGCGGATCTTACGCGGGCGATCCCGCCCCCAGCGCAc                          >  2:1667570/1‑149 (MQ=255)
                                                                                                                           gcagGGCATTGCCAGGCCGAGCACCGTCATCGCCACAACTTTGCTTTCGAACAGCCACGGCAACAGCAGCGACAGCACCAGCGTAGAGATGAAGGCGATAATCAATAAATTGCGGATCTTACGCGGGCGATCCCGCCCCCAGCGCACCa                        >  2:4541310/1‑149 (MQ=255)
                                                                                                                             agGCCATTGCCAGGCCGAGCACCGTCATCGCCACAACTTTGCTTTCGAACAGCCACGGCAACAGCAGCGACAGCACCAGCGTAGAGATGAAGGCGATAATCAATAAATTGCGGATCTTACGCGGGCGATCCCGCCCCCAGCGCACCaga                      <  2:4501816/149‑1 (MQ=255)
                                                                                                                              ggCCATTGCCAGGCCGAGCACCGTCATCGCCACAACTTTGCTTTCGAACAGCCACGGCAACAGCAGCGACAGCACCAGCGTAGAGATGAAGGCGATAATCAATAAATTGCGGATCTTACGCGGGCGATCCCGCCCCCAGCGCACCagag                     <  1:864924/149‑1 (MQ=255)
                                                                                                                                 cATTGCCAGGCCGAGCACCGTCATCGCCACAACTTTGCTTTCGAACAGCCACGGCAACAGCAGCGACAGCACCAGCGTAGAGATGAAGGCGATAATCAATAAATTGCGGATCTTACGCGGGCGATCCCGCCCCCAGCGCACCAGAGGAc                  <  2:1830948/149‑1 (MQ=255)
                                                                                                                                 cATTGCCAGGCCGAGCACCGTCATCGCCACAACTTTGCTTTCGAACAGCCACGGCAACAGCAGCGACAGCACCAGCGTAGAGATGAAGGCGATAATCAATAAATTGCGGATCTTACGCGGGCGATCCCGCCCCCAGCGCACCAGAGGAc                  <  2:335720/149‑1 (MQ=255)
                                                                                                                                   ttGCCAGGCCGAGCACCGTCATCGCCACAACTTTGCTTTCGAACAGCCACGGCAACAGCAGCGACAGCACCAGCGTAGAGATGAAGGCGATAATCAATAAATTGCGGATCTTACGCGGGCGATCCCGCCCCCAGCGCACCAGAGGACCg                <  2:1929185/149‑1 (MQ=255)
                                                                                                                                   ttGCCAGGCCGAGCACCGTCATCGCCACAACTTTGCTTTCGAACAGCCACGGCAACAGCAGCGACAGCACCAGCGTAGAGATGAAGGCGATAATCAATAAATTGCGGATCTTACGCGGGCGATCCCGCCCCCAGCGCACCAGAGGACCg                <  2:2850907/149‑1 (MQ=255)
                                                                                                                                   ttGCCAGGCCGAGCACCGTCATCGCCACAACTTTGCTTTCGAACAGCCACGGCAACAGCAGCGACAGCACCAGCGTAGAGATGAAGGCGATAATCAATAAATTGCGGATCTTACGCGGGCGATCCCGCCCCCAGCGCACCAGAGGACCg                >  2:3531555/1‑149 (MQ=255)
                                                                                                                                        aGGCCGAGCACCGTCATCGCCACAACTTTGCTTTCGAACAGCCACGGCAACAGCAGCGACAGCACCAGCGTAGAGATGAAGGCGATAATCAATAAATTGCGGATCTTACGCGGGCGATCCCGCCCCCAGCGCACCAGAGGACCGACAcc           <  1:1738860/149‑1 (MQ=255)
                                                                                                                                         ggCCGAGCACCGTCATCGCCACAACTTTGCTTTCGAACAGCCACGGCAACAGCAGCGACAGCACCAGCGTAGAGATGAAGGCGATAATCAATAAATTGCGGATCTTACGCGGGCGATCCCGCCCCCAGCGCACCAGAGGACCGACACCa          <  2:1066914/149‑1 (MQ=255)
                                                                                                                                         ggCCGAGCACCGTCATCGCCACAACTTTGCTTTCGAACAGCCACGGCAACAGCAGCGACAGCACCAGCGTAGAGATGAAGGCGATAATCAATAAATTGCGGATCTTACGCGGGCGATCCCGCCCCCAGCGCACCAGAGGACCGACACCa          <  2:3234874/149‑1 (MQ=255)
                                                                                                                                         ggCCGAGCACCGTCATCGCCACAACTTTGCTTTCGAACAGCCACGGCAACAGCAGCGACAGCACCAGCGTAGAGATGAAGGCGATAATCAATAAATTGCGGATCTTACGCGGGCGATCCCGCCCCCAGCGCACCAGAGGACCGACACCa          >  2:3113201/1‑149 (MQ=255)
                                                                                                                                          gCCGAGCACCGTCATCGCCACAACTTTGCTTTCGAACAGCCACGGCAACAGCAGCGACAGCACCAGCGTAGAGATGAAGGCGATAATCAATAAATTGCGGATCTTACGCGGGCGATCCCGCCCCCAGCGCACCAGAGGACCGACACCaa         <  2:1006900/149‑1 (MQ=255)
                                                                                                                                          gCCGAGCACCGTCATCGCCACAACTTTGCTTTCGAACAGCCACGGCAACAGCAGCGACAGCACCAGCGTAGAGATGAAGGCGATAATCAATAAATTGCGGATCTTACGCGGGCGATCCCGCCCCCAGCGCACCAGAGGACCGACACCaa         <  2:3396470/149‑1 (MQ=255)
                                                                                                                                             gAGCACCGTCATCGCCACAACTTTGCTTTCGAACAGCCACGGCAACAGCAGCGACAGCACCAGCGTAGAGATGAAGGCGATAATCAATAAATTGCGGATCTTACGCGGGCGATCCCGCCCCCAGCGCACCAGAGGACCGACACCAAGCa      >  1:1995200/1‑149 (MQ=255)
                                                                                                                                              aGCACCGTCATCGCCACAACTTTTCTTTCGAACAGCCACGGCAACAGCAGCGACAGCACCAGCGTAGAGATGAAGGCGATAATCAATAAATTGCGGATCTTACGAGGGCGATCCCGCCCCCAGCGCACCAGAGGACCGACACCAAGCAg     >  1:3816846/1‑149 (MQ=255)
                                                                                                                                               gCACCGTCATCGCCACAACTTTGCTTTCGAACAGCCACGGCAACAGCAGCGACAGCACCAGCGTAGAGATGAAGGCGATAATCAATAAATTGCGGATCTTACGCGGGCGATCCCGCCCCCAGCGCACCAGAGGACCGACACCAAGCAGt    >  2:2971048/1‑149 (MQ=255)
                                                                                                                                               gCACCGTCATCGCCACAACTTTGCTTTCGAACAGCCACGGCAACAGCAGCGACAGCACCAGCGTAGAGATGAAGGCGATAATCAATAAATTGCGGATCTTACGCGGGCGATCCCGCCCCCAGCGCACCAGAGGACCGACACCAAGCAGt    <  1:3622681/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
AATGCCAACAATTGTCACTGCCAGCCCAAGGTGAGCCGCCACCATCCCCCAATAACTGAAGGTGGTTTTCGTGCCGCGTGAAATACGTAGCGCAGCTTCCGCAATTGCCAGCACCGCAATCCAGCAGGCCATTGCCAGGCCGAGCACCGTCATCGCCACAACTTTGCTTTCGAACAGCCACGGCAACAGCAGCGACAGCACCAGCGTAGAGATGAAGGCGATAATCAATAAATTGCGGATCTTACGCGGGCGATCCCGCCCCCAGCGCACCAGAGGACCGACACCAAGCAGTAG  >  CP009273/2286881‑2287174

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑