Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,289,374 C→A 92.0% V142V (GTG→GTT ccmC ← heme exporter subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,289,3740CA92.0% 1067.6 / 78.7 351V142V (GTG→GTTccmCheme exporter subunit
Reads supporting (aligned to +/- strand):  ref base C (12/16);  new base A (155/168);  total (167/184)
Fisher's exact test for biased strand distribution p-value = 6.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.21e-01

Only 100 of 357 total aligned reads displayed.

CTGCATCCGCGTTGATCCCTGATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATCACACCCACATACAAAAACAGCAGCACCAGTTCAGAAGTCAGACGTGCATCCCATACCCACCAGGTGCCCCACATCGGTTTTCCCCATGCAGAGCCGGTAACCAGGGCAATAAAGGTAAACACGGCACCAATGGGGGCCATCGCCG  >  CP009273/2289227‑2289519
                                                                                                                                                   |                                                                                                                                                 
 tgcATCCGCGTTGATCCCTGATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATaaca                                                                                                                                                <  2:2222856/150‑1 (MQ=255)
     tCCGCGTTGATCCCTGATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATAACACCca                                                                                                                                            <  2:3200004/150‑1 (MQ=255)
       cgcgTTGATCCCTGATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATAACACCcaca                                                                                                                                          >  1:3304652/1‑150 (MQ=255)
        gcgTTGATCCCTGATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATAACACCCACAt                                                                                                                                         <  1:3808202/150‑1 (MQ=255)
        gcgTTGATCCCTGATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATAACACCCACAt                                                                                                                                         <  1:355096/150‑1 (MQ=255)
        gcgTTGATCCCTGATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATAACACCCACAt                                                                                                                                         <  1:2746369/150‑1 (MQ=255)
         cgTTGATCCCTGATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATAACACCCACATa                                                                                                                                        >  1:2639099/1‑150 (MQ=255)
          gTTGATCCCTGATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATAACACCCACATAc                                                                                                                                       >  1:3028646/1‑150 (MQ=255)
          gTTGATCCCTGATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATAACACCCACATAc                                                                                                                                       >  2:3381955/1‑150 (MQ=255)
           ttGATCCCTGATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATCACACCCACATACa                                                                                                                                      >  1:2223253/1‑150 (MQ=255)
             gATCCCTGATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCGGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATAACACCCACATACaaa                                                                                                                                    <  2:433257/150‑1 (MQ=255)
                cccTGATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCGGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATAACACCCACATACAAAAAc                                                                                                                                 >  2:2910704/1‑150 (MQ=255)
                 ccTGATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATAACACCCACATACAAAAAca                                                                                                                                <  2:3566250/150‑1 (MQ=255)
                 ccTGATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATAACACCCACATACAAAAAca                                                                                                                                <  1:485282/150‑1 (MQ=255)
                   tGATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATCCCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATAACACCCACATACACAAAcagc                                                                                                                              >  2:3779006/1‑150 (MQ=255)
                   tGATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATAACACCCACATACAAAAAcagc                                                                                                                              >  2:3753747/1‑150 (MQ=255)
                        cAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATAACACCCACATACAAAAACAGCAGCAc                                                                                                                         <  1:4753319/150‑1 (MQ=255)
                           ggTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATAACACCCACATACAAAAACAGCAGCACCAg                                                                                                                      <  2:4448737/150‑1 (MQ=255)
                           ggTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATAACACCCACATACAAAAACAGCAGCACCAg                                                                                                                      <  2:1168223/150‑1 (MQ=255)
                            gtgtTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATAACACCCACATACAAAAACAGCAGCACCAGt                                                                                                                     >  2:3495512/1‑150 (MQ=255)
                             tgtTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATAACACCCACATACAAAAACAGCAGCACCAGtt                                                                                                                    >  1:939724/1‑150 (MQ=255)
                              gtTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATAACACCCACATACAAAAACAGCAGCACCAGTTc                                                                                                                   <  1:2985469/150‑1 (MQ=255)
                                tCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATAACACCCACATACAAAAACAGCAGCACCAGTTCAg                                                                                                                 >  1:4233875/1‑150 (MQ=255)
                                 ccaccaCTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATAACACCCACATACAAAAACAGCAGCACCAGTTCAGa                                                                                                                <  2:4611787/150‑1 (MQ=255)
                                   accacTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATAACACCCACATACAAAAACAGCAGCACCAGTTCAGAAg                                                                                                              >  1:3347264/1‑150 (MQ=255)
                                   accacTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATAACACCCACATACAAAAACAGCAGCACCAGTTCAGAAg                                                                                                              <  1:3450759/150‑1 (MQ=255)
                                    ccacTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATCACACCCACATACAAAAACAGCAGCACCAGTTCAGAAGt                                                                                                             <  1:1599951/150‑1 (MQ=255)
                                    ccacTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATAACACCCACATACAAAAACAGCAGCACCAGTTCAGAAGt                                                                                                             >  1:1000729/1‑150 (MQ=255)
                                    ccacTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATAACACCCACATACAAAAACAGCAGCACCAGTTCAGAAGt                                                                                                             <  2:3738070/150‑1 (MQ=255)
                                     cacTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATAACACCCACATACAAAAACAGCAGCACCAGTTCAGAAGTc                                                                                                            >  2:3880607/1‑150 (MQ=255)
                                       cTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATAACACCCACATACAAAAACAGCAGCACCAGTTCAGAAGTCAg                                                                                                          >  2:680068/1‑150 (MQ=255)
                                         ccACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATAACACCCACATACAAAAACAGCAGCACCAGTTCAGAAGTCAGAc                                                                                                        <  1:409234/150‑1 (MQ=255)
                                         ccACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATAACACCCACATACAAAAACAGCAGCACCAGTTCAGAAGTCAGAc                                                                                                        >  1:2023751/1‑150 (MQ=255)
                                          cACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATAACACCCACATACAAAAACAGCAGCACCAGTTCAGAAGTCAGACg                                                                                                       >  2:1876610/1‑150 (MQ=255)
                                            cGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATAACACCCACATACAAAAACAGCAGCACCAGTTCAGAAGTCAGACGTg                                                                                                     <  2:4713582/150‑1 (MQ=255)
                                                 tAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATAACACCCACATACAAAAACAGCAGCACCAGTTCAGAAGTCAGACGTGCATcc                                                                                                >  1:3413037/1‑150 (MQ=255)
                                                   aTGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATAACACCCACATACAAAAACAGCAGCACCAGTTCAGAAGTCAGACGTGCATCCCa                                                                                              >  1:3939532/1‑150 (MQ=255)
                                                        taatCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATAACACCCACATACAAAAACAGCAGCACCAGTTCAGAAGTCAGACGTGCATCCCATAccc                                                                                         <  2:3617515/150‑1 (MQ=255)
                                                        taatCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATAACACCCACATACAAAAACAGCAGCACCAGTTCAGAAGTCAGACGTGCATCCCATAccc                                                                                         >  1:3637976/1‑150 (MQ=255)
                                                        taatCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATAACACCCACATACAAAAACAGCAGCACCAGTTCAGAAGTCAGACGTGCATCCCATAccc                                                                                         >  1:462836/1‑150 (MQ=255)
                                                        taatCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATAACACCCACATACAAAAACAGCAGCACCAGTTCAGAAGTCAGACGTGCATCCCATAccc                                                                                         >  1:2855207/1‑150 (MQ=255)
                                                         aatCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATAACACCCACATACAAAAACAGCAGCACCAGTTCAGAAGTCAGACGTGCATCCCATACcca                                                                                        >  1:3977981/1‑150 (MQ=255)
                                                         aatCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATAACACCCACATACAAAAACAGCAGCACCAGTTCAGAAGTCAGACGTGCATCCCATACcca                                                                                        >  2:1835688/1‑150 (MQ=255)
                                                          atCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATAACACCCACATACAAAAACAGCAGCACCAGTTCAGAAGTCAGACGTGCATCCCATACccac                                                                                       >  2:3111926/1‑150 (MQ=255)
                                                            cGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATAACACCCACATACAAAAACAGCAGCACCAGTTCAGAAGTCAGACGTGCATCCCATACccacca                                                                                     <  2:1254758/150‑1 (MQ=255)
                                                            cGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATAACACCCACATACAAAAACAGCAGCACCAGTTCAGAAGTCAGACGTGCATCCCATACccacca                                                                                     <  2:603515/150‑1 (MQ=255)
                                                            cGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATAACACCCACATACAAAAACAGCAGCACCAGTTCAGAAGTCAGACGTGCATCCCATACccacca                                                                                     <  2:2993104/150‑1 (MQ=255)
                                                              gCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATAACACCCACATACAAAAACAGCAGCACCAGTTCAGAAGTCAGACGTGCATCCCATACCCACCAgg                                                                                   >  1:4503083/1‑150 (MQ=255)
                                                               cAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATAACACCCACATACAAAAACAGCAGCACCAGTTCAGAAGTCAGACGTGCATCCCATACCCACCAGGt                                                                                  >  1:2118518/1‑150 (MQ=255)
                                                                aGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATAACACCCACATACAAAAACAGCAGCACCAGTTCAGAAGTCAGACGTGCATCCCATACCCACCAGGTg                                                                                 <  2:2067777/150‑1 (MQ=255)
                                                                    tCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATAACACCCACATACAAAAACAGCAGCACCAGTTCAGAAGTCAGACGTGCATCCCAGACCCACCAGGTGcccc                                                                             >  2:4177468/1‑150 (MQ=255)
                                                                     caccacGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATAACACCCACATACAAAAACAGCAGCACCAGTTCAGAAGTCAGACGTGCATCCCATACCCACCAGGTGCCCca                                                                            <  2:231427/150‑1 (MQ=255)
                                                                        cacGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATAACACCCACATACAAAAACAGCAGCACCAGTTCAGAAGTCAGACGTGCATCCCATACCCACCAGGTGCCCCACAt                                                                         >  2:1298705/1‑150 (MQ=255)
                                                                                  aGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATAACACCCACATACAAAAACAGCAGCACCAGTTCAGAAGTCAGACGTGCATCCCATACCCACCAGGTGCCCCACATCGGTTTTccc                                                               <  2:495237/150‑1 (MQ=255)
                                                                                      ccAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATAACACCCACATACAAAAACAGCAGCACCAGTTCAGAAGTCAGACGTGCATCCCATACCCACCAGGTGCCCCACATCGGTTTTCCCCATg                                                           <  1:2198349/150‑1 (MQ=255)
                                                                                        aGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATAACACCCACATACAAAAACAGCAGCACCAGTTCAGAAGTCAGACGTGCATCCCATACCCACCAGGTGCCCCACATCGGTTTTCCCCATGCa                                                         >  2:1972127/1‑150 (MQ=255)
                                                                                         ggATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATAACACCCACATACAAAAACAGCAGCACCAGTTCAGAAGTCAGACGTGCATCCCATACCCACCAGGTGCCCCACATCGGTTTTCCCCATGCag                                                        <  1:4140341/150‑1 (MQ=255)
                                                                                         ggATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATAACACCCACATACAAAAACAGCAGCACCAGTTCAGAAGTCAGACGTGCATCCCATACCCACCAGGTGCCCCACATCGGTTTTCCCCATGCag                                                        <  2:4505961/150‑1 (MQ=255)
                                                                                         ggATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATAACACCCACATACAAAAACAGCAGCACCAGTTCAGAAGTCAGACGTGCATCCCATACCCACCAGGTGCCCCACATCGGTTTTCCCCATGCag                                                        >  2:334476/1‑150 (MQ=255)
                                                                                         ggATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATAACACCCACATACAAAAACAGCAGCACCAGTTCAGAAGTCAGACGTGCATCCCATACCCACCAGGTGCCCCACATCGGTTTTCCCCATGCag                                                        >  1:4522314/1‑150 (MQ=255)
                                                                                         ggATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATAACACCCACATACAAAAACAGCAGCACCAGTTCAGAAGTCAGACGTGCATCCCATACCCACCAGGTGCCCCACATCGGTTTTCCCCATGCag                                                        <  1:1736339/150‑1 (MQ=255)
                                                                                          gATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATAACACCCACATACAAAAACAGCAGCACCAGTTCAGAAGTCAGACGTGCATCCCATACCCACCAGGTGCCCCACATCGGTTTTCCCCATGCaga                                                       >  1:1488305/1‑150 (MQ=255)
                                                                                           aTACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATCACACCCACATACAAAAACAGCAGCACCAGTTCCGAAGTCAGACGTGCATCCCATACCCACCAGGTGCCCCACATCGGTTTTCCCCATGCagag                                                      >  2:4015688/1‑150 (MQ=255)
                                                                                           aTACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATAACACCCACATACAAAAACAGCAGCACCAGTTCAGAAGTCAGACGTGCATCCCATACCCACCAGGTGCCCCACATCGGTTTTCCCCATGCagag                                                      <  1:1954387/150‑1 (MQ=255)
                                                                                              ccTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATCACACCCACATACAAAAACAGCAGCACCAGTTCAGAAGTCAGACGTGCATCCCATACCCACCAGGTGCCCCACATCGGTTTTCCCCATGCAGAGCCg                                                   <  1:4480542/150‑1 (MQ=255)
                                                                                              ccTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATAACACCCACATACAAAAACAGCAGCACCAGTTCAGAAGTCAGACGTGCATCCCATACCCACCAGGTGCCCCACATCGGTTTTCCCCATGCAGAGCCg                                                   <  2:2786864/150‑1 (MQ=255)
                                                                                              ccTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATAACACCCACATACAAAAACAGCAGCACCAGTTCAGAAGTCAGACGTGCATCCCATACCCACCAGGTGCCCCACATCGGTTTTCCCCATGCAGAGCCg                                                   >  2:3144165/1‑150 (MQ=255)
                                                                                               cTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATAACACCCACATACAAAAACAGCAGCACCAGTTCAGAAGTCAGACGTGCATCCCATACCCACCAGGTGCCCCACATCGGTTTTCCCCATGCAGAGCCgg                                                  >  1:1222339/1‑150 (MQ=255)
                                                                                               cTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATAACACCCACATACAAAAACAGCAGCACCAGTTCAGAAGTCAGACGTGCATCCCATACCCACCAGGTGCCCCACATCGGTTTTCCCCATGCAGAGCCgg                                                  >  2:322804/1‑150 (MQ=255)
                                                                                                 gCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATAACACCCACATACAAAAACAGCAGCACCAGTTCAGAAGTCAGACGTGCATCCCATACCCACCAGGTGCCCCACATCGGTTTTCCCCATGCAGAGCCGGTa                                                <  1:613207/150‑1 (MQ=255)
                                                                                                    gCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATAACACCCACATACAAAAACAGCAGCACCAGTTCAGAAGTCAGACGTGCATCCCATACCCACCAGGTGCCCCACATCGGTTTTCCCCATGCAGAGCCGGTAAcc                                             <  2:1375148/150‑1 (MQ=255)
                                                                                                       cGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATAACACCCACATACAAAAACAGCAGCACCAGTTCAGAAGTCAGACGTGCATCCCATACCCACCAGGTGCCCCACATCGGTTTTCCCCATGCAGAGCCGGTAACCAgg                                          <  1:2445354/150‑1 (MQ=255)
                                                                                                       cGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATAACACCCACATACAAAAACAGCAGCACCAGTTCAGAAGTCAGACGTGCATCCCATACCCACCAGGTGCCCCACATCGGTTTTCCCCATGCAGAGCCGGTAACCAgg                                          <  1:3125563/150‑1 (MQ=255)
                                                                                                        ggCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATCACACCCACATACAAAAACAGCAGCACCAGTTCAGAAGTCAGACGTGCATCCCATACCCACCAGGTGCCCCACATCGGTTTTCCCCATGCAGAGCCGGTAACCAggg                                         >  2:2824996/1‑150 (MQ=255)
                                                                                                          cccGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATAACACCCACATACAAAAACAGCAGCACCAGTTCAGAAGTCAGACGTGCATCCCATTCCCACCAGGTGCCCCACATCGGTTTTCCCCATGCAGAGCCGGTAACCAGGGCa                                       <  1:1085819/150‑1 (MQ=255)
                                                                                                            cGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATAACACCCACATACAAAAACAGCAGCACCAGTTCAGAAGTCAGACGTGCATCCCATACCCACCAGGTGCCCCACATCGGTTTTCCCCATGCAGAGCCGGTAACCAGGGCAAt                                     <  2:2607617/150‑1 (MQ=255)
                                                                                                            cGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATAACACCCACATACAAAAACAGCAGCACCAGTTCAGAAGTCAGACGTGCATCCCATACCCACCAGGTGCCCCACATCGGTTTTCCCCATGCAGAGCCGGTAACCAGGGCAAt                                     <  1:2303071/150‑1 (MQ=255)
                                                                                                            cGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATAACACCCACATACAAAAACAGCAGCACCAGTTCAGAAGTCAGACGTGCATCCCATACCCACCAGGTGCCCCACATCGGTTTTCCCCATGCAGAGCCGGTAACCAGGGCAAt                                     <  1:2285752/150‑1 (MQ=255)
                                                                                                               cAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATAACACCCACATACAAAAACAGCAGCACCAGTTCAGAAGTCAGACGTGCATCCCATACCCACCAGGTGCCCCACATCGGTTTTCCCCATGCAGAGCCGGTAACCAGGGCAATaaa                                  >  2:2407732/1‑150 (MQ=255)
                                                                                                                aGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATCACACCCACATACAAAAACAGCAGCACCAGTTCAGAAGTCAGACGTGCATCCCATACCCACCAGGTGCCCCACATCGGTTTTCCCCATGCAGAGCCGGTAACCAGGGCAATAAAg                                 <  2:1419739/150‑1 (MQ=255)
                                                                                                                 gACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATAACACCCACATACAAAAACAGCAGCACCAGTTCAGAAGTCAGACGTGCATCCCATACCCACCAGGTGCCCCACATCGGTTTTCCCCATGCAGAGCCGGTAACCAGGGCAATAAAgg                                >  2:4272892/1‑150 (MQ=255)
                                                                                                                 gACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATAACACCCACATACAAAAACAGCAGCACCAGTTCAGAAGTCAGACGTGCATCCCATACCCACCAGGTGCCCCACATCGGTTTTCCCCATGCAGAGCCGGTAACCAGGGCAATAAAgg                                >  2:4273893/1‑150 (MQ=255)
                                                                                                                  aCGGCGGTCGTCGAAGGCGTGCCACAGGGCAATAACACCCACATACAAAAACAGCAGCACCAGTTCAGAAGTCAGACGTGCATCCCATACCCACCAGGTGCCCCACATCGGTTTTCCCCATGCAGAGCCGGTAACCAGGGCAATAAAGGt                               >  2:106523/1‑150 (MQ=255)
                                                                                                                    ggcggTCGTCGAAGGCGTGCCACAGGGCAATAACACCCACATACAAAAACAGCAGCACCAGTTCAGAAGTCAGACGTGCATCCCATACCCACCAGGTGCCCCACATCGGTTTTCCCCATGCAGAGCCGGTAACCAGGGCAATAAAGGTaa                             >  1:1033392/1‑150 (MQ=255)
                                                                                                                     gcggTCGTCGAAGGCGTGCCACAGGGCAATAACACCCACATACAAAAACAGCAGCACCAGTTCAGAAGTCAGACGTGCATCCCATACCCACCAGGTGCCCCACATCGGTTTTCCCCATGCAGAGCCGGTAACCAGGGCAATAAAGGTaaa                            <  1:1504374/150‑1 (MQ=255)
                                                                                                                      cggTCGTCGAAGGCGTGCCACAGGGCAATAACACCCACATACAAAAACAGCAGCACCAGTTCAGAAGTCAGACGTGCATCCCATACCCACCAGGTGCCCCACATCGGTTTTCCCCATGCAGAGCCGGTAACCAGGGCAATAAAGGTAAac                           <  2:4138362/150‑1 (MQ=255)
                                                                                                                       ggTCGTCGAAGGCGTGCCACAGGGCAATAACACCCACATACAAAAACAGCAGCACCAGTTCAGAAGTCAGACGTGCATCCCATACCCACCAGGTGCCCCACATCGGTTTTCCCCATGCAGAGCCGGTAACCAGGGCAATAAAGGTAAaca                          >  2:1504005/1‑150 (MQ=255)
                                                                                                                       ggTCGTCGAAGGCGTGCCACAGGGCAATAACACCCACATACAAAAACAGCAGCACCAGTTCAGAAGTCAGACGTGCATCCCATACCCACCAGGTGCCCCACATCGGTTTTCCCCATGCAGAGCCGGTAACCAGGGCAATAAAGGTAAaca                          >  1:4063997/1‑150 (MQ=255)
                                                                                                                           gtcgAAGGCGTGCCACAGGGCAATAACACCCACATACAAAAACAGCAGCACCAGTTCAGAAGTCAGACGTGCATCCCATACCCACCAGGTGCCCCACATCGGTTTTCCCCATGCAGAGCCGGTAACCAGGGCAATAAAGGTAAACACGGc                      >  1:784505/1‑150 (MQ=255)
                                                                                                                           gtcgAAGGCGTGCCACAGGGCAATAACACCCACATACAAAAACAGCAGCACCAGTTCAGAAGTCAGACGTGCATCCCATACCCACCAGGTGCCCCACATCGGTTTTCCCCATGCAGAGCCGGTAACCAGGGCAATAAAGGTAAACACGGc                      >  1:2672030/1‑150 (MQ=255)
                                                                                                                            tcgAAGGCGTGCCACAGGGCAATAACACCCACATACAAAAACAGCAGCACCAGTTCAGAAGTCAGACGTGCATCCCATACCCACCAGGTGCCCCACATCGGTTTTCCCCATGCAGAGCCGGTAACCAGGGCAATAAAGGTAAACACGGCa                     <  2:3028646/150‑1 (MQ=255)
                                                                                                                            tcgAAGGCGTGCCACAGGGCAATAACACCCACATACAAAAACAGCAGCACCAGTTCAGAAGTCAGACGTGCATCCCATACCCACCAGGTGCCCCACATCGGTTTTCCCCATGCAGAGCCGGTAACCAGGGCAATAAAGGTAAACACGGCa                     >  2:1288080/1‑150 (MQ=255)
                                                                                                                              gAAGGCGTGCCACAGGGCAATAACACCCACATACAAAAACAGCAGCACCAGTTCAGAAGTCAGACGTGCATCCCATACCCACCAGGTGCCCCACATCGGTTTTCCCCATGCAGAGCCGGTAACCAGGGCAATAAAGGTAAACACGGCAcc                   <  1:3787826/150‑1 (MQ=255)
                                                                                                                                aGGCGTGCCACAGGGCAATAACACCCACATACAAAAACAGCAGCACCAGTTCAGAAGTCAGACGTGCATCCCATACCCACCAGGTGCCCCACATCGGTTTTCCCCATGCAGAGCCGGTAACCAGGGCAATAAAGGTAAACACGGCACCaa                 <  1:657821/150‑1 (MQ=255)
                                                                                                                                   cGTGCCACAGGGCAATAACACCCACATACAAAAACAGCAGCACCAGTTCAGAAGTCAGACGTGCATCCCATACCCACCAGGTGCCCCACATCGGTTTTCCCCATGCAGAGCCGGTAACCAGGGCAATAAAGGTAAACACGGCACCAATgg              >  1:715996/1‑150 (MQ=255)
                                                                                                                                    gTGCCACAGGGCAATCACACCCACATACAAAAACAGCAGCACCAGTTCAGAAGTCAGACGTGCATCCCATACCCACCAGGTGCCCCACATCGGTTTTCCCCATGCAGAGCCGGTAACCAGGGCAATAAAGGTAAACACGGCACCAATggg             >  2:710421/1‑150 (MQ=255)
                                                                                                                                     tGCCACAGGGCACTAACACCCACATACAAAAACAGCAGCACCAGTTCAGAAGTCAGACGTGCATCCCATACCCACCAGGTGCCCCACATCGGTTTTCCCCATGCAGAGCCGGTAACCAGGGCAATAAAGGTAAACACGGCACCAATgggg            <  1:4368898/150‑1 (MQ=255)
                                                                                                                                     tGCCACAGGGCAATAACACCCACATACAAAAACAGCAGCACCAGTTCAGAAGTCAGACGTGCATCCCATACCCACCAGGTGCCCCACATCGGTTTTCCCCATGCAGAGCCGGTAACCAGGGCAATAAAGGTAAACACGGCACCAATgggg            <  1:3581014/150‑1 (MQ=255)
                                                                                                                                              gCAATAACACCCACATACAAAAACAGCAGCACCAGTTCAGAAGTCAGACGTGCATCCCATACCCACCAGGTGCCCCACATCGGTTTTCCCCATGCAGAGCCGGTAACCAGGGCAATAAAGGTAAACACGGCACCAATGGGGGCCATcgcc   >  2:837876/1‑150 (MQ=255)
                                                                                                                                               cAATCACACCCACATACAAAAACAGCAGCACCAGTTCAGAAGTCAGACGTGCATCCCATACCCACCAGGTGCCCCACATCGGTTTTCCCCATGCAGAGCCGGTAACCAGGGCAATAAAGGTAAACACGGCACCAATGGGGGCCATcgccg  <  2:1699486/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
CTGCATCCGCGTTGATCCCTGATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATCACACCCACATACAAAAACAGCAGCACCAGTTCAGAAGTCAGACGTGCATCCCATACCCACCAGGTGCCCCACATCGGTTTTCCCCATGCAGAGCCGGTAACCAGGGCAATAAAGGTAAACACGGCACCAATGGGGGCCATCGCCG  >  CP009273/2289227‑2289519

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: