Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,289,374 C→A 62.6% 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,3740CA62.6% 238.6 / inf 332V142V (GTG→GTTccmCheme exporter subunit
Reads supporting (aligned to +/- strand):  ref base C (69/55);  new base A (89/118);  total (158/174)
Fisher's exact test for biased strand distribution p-value = 3.08e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.30e-01

Only 100 of 337 total aligned reads displayed.

TGCTGCATCCGCGTTGATCCCTGATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATCACACCCACATACAAAAACAGCAGCACCAGTTCAGAAGTCAGACGTGCATCCCATACCCACCAGGTGCCCCACATCGGTTTTCCCCATGCAGAGCCGGTAACCAGGGCAATAAAGGTAAACACGGCACCAATGGGGGCCATCGCCGCC  >  CP009273/2289225‑2289521
                                                                                                                                                     |                                                                                                                                                   
   tgcATCCGCGTTGATCCCTGATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATcaca                                                                                                                                                  <  2:2829543/150‑1 (MQ=255)
   tgcATCCGCGTTGATCCCTGATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATaaca                                                                                                                                                  <  1:1373666/150‑1 (MQ=255)
     cATCCGCGTTGATCCCTGATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATAACAcc                                                                                                                                                <  2:1089820/150‑1 (MQ=255)
     cATCCGCGTTGATCCCTGATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATAACAcc                                                                                                                                                <  1:2497909/150‑1 (MQ=255)
      aTCCGCGTTGATCCCTGATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATAACAccc                                                                                                                                               >  1:1541638/1‑150 (MQ=255)
       tCCGCGTTGATCCCTGATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATCACACCca                                                                                                                                              <  2:4956450/150‑1 (MQ=255)
       tCCGCGTTGATCCCTGATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATCACACCca                                                                                                                                              >  1:690703/1‑150 (MQ=255)
         cgcgTTGATCCCTGCTGCAGGGTGTTCCACCACTCCACGGAGTAATGAATACTCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATAACACCcaca                                                                                                                                            <  2:1504369/150‑1 (MQ=255)
            gTTGATCCCTGATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATAACACCCACATAc                                                                                                                                         <  1:278359/150‑1 (MQ=255)
            gTTGATCCCTGATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATAACACCCACATAc                                                                                                                                         <  1:2323392/150‑1 (MQ=255)
              tGATCCCTGATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATAACACCCACATACaa                                                                                                                                       >  1:4739500/1‑150 (MQ=255)
               gATCCCTGATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATAGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATCACACCCACATACaaa                                                                                                                                      >  1:2523520/1‑150 (MQ=255)
                    cTGATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATAACACCCACATACAAAAAcag                                                                                                                                 >  1:3926958/1‑150 (MQ=255)
                      gATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATAACACCCACATACAAAAAcagca                                                                                                                               <  1:519412/150‑1 (MQ=255)
                         gCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATCACACCCACATACAAAAACAgcagca                                                                                                                            <  1:4466445/150‑1 (MQ=255)
                          cAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATCACACCCACATACAAAAACAGCAGCAc                                                                                                                           >  1:2870081/1‑150 (MQ=255)
                           aGGGTGTTCCACCACTCCACGGAGTAATGAATAATAGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATCACACCCACATACAAAAACAGCAGCAcc                                                                                                                          <  1:2085786/150‑1 (MQ=255)
                             ggTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATCACACCCACATACAACAACAGCAGCACCAg                                                                                                                        >  2:4538785/1‑150 (MQ=255)
                             ggTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATAACACCCACATACAAAAACAGCAGCACCAg                                                                                                                        <  2:1163923/150‑1 (MQ=255)
                                gtTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGTCGTGCCACATGTCAATAACACCCACATTCAAAAACAGCAGCACCAGTTc                                                                                                                     <  2:4435717/150‑1 (MQ=255)
                                gtTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATCACACCCACATACAAAAACAGCAGCACCAGTTc                                                                                                                     <  1:2409383/150‑1 (MQ=255)
                                gtTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATCACACCCACATACAAAAACAGCAGCACCAGTTc                                                                                                                     >  2:766079/1‑150 (MQ=255)
                                gtTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATAACACCCACATACAAAAACAGCAGCACCAGTTc                                                                                                                     <  1:4241294/150‑1 (MQ=255)
                                gtTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATAACACCCACATACAAAAACAGCAGCACCAGTTc                                                                                                                     >  1:4952702/1‑150 (MQ=255)
                                  tCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATAACACCCACATACAAAAACAGCAGCACCAGTTCAg                                                                                                                   >  1:3046664/1‑150 (MQ=255)
                                  tCCACCACTCCACGGAGTAATGAATAATAGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATCACACCCACATACAAAAACAGCAGCACCAGTTCAg                                                                                                                   >  2:2754408/1‑150 (MQ=255)
                                    caccacTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATAACACCCACATACAAAAACAGCAGCACCAGTTCAGaa                                                                                                                 >  2:4267664/1‑150 (MQ=255)
                                      ccacTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATAACACCCACATACAAAAACAGCAGCACCAGTTCAGAAGt                                                                                                               <  2:1417184/150‑1 (MQ=255)
                                       cacTCCACGGAGTAATGAATAATCGGCAGAGTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATAACACCCACATACAAAAACAGCAGCACCAGTTCAGAAGTc                                                                                                              <  2:5067094/150‑1 (MQ=255)
                                         cTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATAACACCCACATACAAAAACAGCAGCACCAGTTCAGAAGTCAg                                                                                                            <  2:31093/150‑1 (MQ=255)
                                          tCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATAACACCCACATACAAAAACAGCCGCACCAGTTCAGAAGTCAGa                                                                                                           >  2:2933082/1‑150 (MQ=255)
                                            cACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGTCGGCGGTCGTCGAAGGCGTGCCACAGGGCAATAACACCCACATACAAAAACAGCAGCACCAGTTCAGAAGTCAGACg                                                                                                         <  2:4784934/150‑1 (MQ=255)
                                             aCGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATCACACCCACATACAAAAACAGCAGCACCAGTTCAGAAGTCAGACgn                                                                                                        <  2:173012/150‑2 (MQ=255)
                                              cGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATAACACCCACATACAAAAACAGCAGCACCAGTTCAGAAGTCAGACGTg                                                                                                       <  1:1877393/150‑1 (MQ=255)
                                                gAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATAACACCCACATACAAAAACAGCAGCACCAGTTCAGAAGTCAGACGTGCa                                                                                                     <  2:3975119/150‑1 (MQ=255)
                                                gAGTAATGAATAATAGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATCACACCCACATACAAAAACAGCAGCACCAGTTCAGAAGTCAGACGTGCa                                                                                                     <  1:2840285/150‑1 (MQ=255)
                                                   tAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATCACCCCCACATACAAAAACAGCAGCACCAGTTCAGAAGTCAGACGTGCATcc                                                                                                  <  2:664150/150‑1 (MQ=255)
                                                   tAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATCACACCCACATACAAAAACAGCAGCACCAGTTCAGAAGTCAGACGTGCATcc                                                                                                  <  2:3631939/150‑1 (MQ=255)
                                                     aTGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATAACACCCACATACAAAAACAGCAGCACCAGTTCAGAAGTCAGACGTGCATCCCa                                                                                                >  1:3662359/1‑150 (MQ=255)
                                                      tGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATAACACCCACATACAAAAACAGCAGCACCAGTTCAGAAGTCAGACGTGCATCCCAt                                                                                               >  1:1787772/1‑150 (MQ=255)
                                                        aataatCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACCGCCAGCCAGACGGCGGTCGTCGAAGGCGTGACACAGGGCAATAACACCCACATACAAAAACAGCAGCACCAGTTCAGAAGTCAGACGTGCATCCCATaa                                                                                             >  1:3159717/1‑149 (MQ=255)
                                                         ataatCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATAACACCCACATACAAAAACAGCAGCACCAGTTCAGAAGTCAGACGTGCATCCCATAcc                                                                                            >  1:294668/1‑150 (MQ=255)
                                                         ataatCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATAACACCCACATACAAAAACAGCAGCACCAGTTCAGAAGTCAGACGTGCATCCCATAcc                                                                                            <  1:3303412/150‑1 (MQ=255)
                                                           aatCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATAACACCCACATACAAAAACAGCAGCACCAGTTCAGAAGTCAGACGTGCATCCCATACcca                                                                                          <  2:410112/150‑1 (MQ=255)
                                                           aatCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATAACACCCACATACAAAAACAGCAGCACCAGTTCAGAAGTCAGACGTGCATCCCATACcca                                                                                          >  1:1099574/1‑150 (MQ=255)
                                                           aatCGGCAGAGTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATTACACCCACATACAAATACAGCAGCACCAGTTCAGTAGACAGACGTGCATCCCATACcca                                                                                          <  2:3197523/150‑1 (MQ=255)
                                                              cGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCACTCACACCCACATACAAAAACAGCAGCACCAGTTCAGAAGTCAGACGTGCATCCCATACccacca                                                                                       >  2:4597672/1‑150 (MQ=255)
                                                              cGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATAACACCCACATACAAAAACAGCAGCACCAGTTCAGAAGTCAGACGTGCATCCCATACccacca                                                                                       <  2:3731179/150‑1 (MQ=255)
                                                               ggCAGATTCACCACGCCATTCATCACCAGGATACGTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATCACACCCACTTTCAAAAACAGCAGCTCCAGTTCAGAAGTCAGACGTGCTTCCCATACCCACCAg                                                                                      <  2:4342953/150‑1 (MQ=255)
                                                               ggCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATCACACCCACATACAAAAACAGCAGCACCAGTTCAGAAGTCAGACGTGCATCCCATACCCACCAg                                                                                      <  2:1735521/150‑1 (MQ=255)
                                                                gCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATCACACCCACATACAAAAACAGCAGCACCAGTTCAGAAGTCAGACGTGCATCCCATACCCACCAgg                                                                                     >  1:3856944/1‑150 (MQ=255)
                                                                  aGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATCACACCCACATACAAAAACAGCAGCACCAGTTCAGAAGTCAGACGTGCATCCCATACCCACCAGGTg                                                                                   >  1:5058780/1‑150 (MQ=255)
                                                                  aGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATCACACCCACATACAAAAACAGCAGCACCAGTTCAGAAGTCAGACGTGCATCCCATACCCACCAGGTg                                                                                   >  1:4888463/1‑150 (MQ=255)
                                                                   gATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCACTCACACCCACATACAAAAACAGCAGCACCAGTTCAGAAGTCAGACGTGCATCCCATACCCACCAGGTGc                                                                                  >  2:1356082/1‑150 (MQ=255)
                                                                    aTTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATAACACCCACATACAAAAACAGCAGCACCAGTTCAGAAGTCAGACGTGCATCCCATACCCACCAGGTGcc                                                                                 <  1:927954/150‑1 (MQ=255)
                                                                      tCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATCACACCCACATACAAAAACAGCAGCACCAGTTCAGAAGTCAGACGTGCATCCCATACCCACCAGGTGcccc                                                                               >  1:825675/1‑150 (MQ=255)
                                                                        accacGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATCACACCCACATACAAAAACAGCAGCACCAGTTCAGAAGTCAGACGTGCATCCCATACCCACCAGGTGCCCcac                                                                             >  1:2256554/1‑150 (MQ=255)
                                                                         ccacGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATCACACCCACATACAAAAACAGCAGCACCAGTTCAGAAGTCAGACGTGCATCCCATACCCACCAGGTGCCCcaca                                                                            >  1:2982960/1‑150 (MQ=255)
                                                                          cacGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATAACACCCACATACAAAAACAGCAGCACCAGTTCAGAAGTCAGACGTGCATCCCATACCCACCAGGTGCCCCACAt                                                                           >  1:1643099/1‑150 (MQ=255)
                                                                           acGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATAACACCCACATACAAAAACAGCAGCACCAGTTCAGAAGTCAGACGTGCATCCCATACCCACCAGGTGCCCCACATc                                                                          <  2:4885999/150‑1 (MQ=255)
                                                                            cGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATAACACCCACATACAAAAACAGCAGCACCAGTTCAGAAGTCAGACGTGCATCCCATACCCACCAGGTGCCCCACATCg                                                                         >  1:336935/1‑150 (MQ=255)
                                                                             gCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATAACACCCACATACAAAAACAGCAGCACCAGTTCAGAAGTCAGACGTGCATCCCATCCCCACCAGGTGCCCCACATCgg                                                                        >  2:3837135/1‑150 (MQ=255)
                                                                                aaTCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATAACACCCACATACAAAAACAGCAGCACCAGTTCAGAAGTCAGACGTGCATCCCATACCCACCAGGTGCCCCACATCGGttt                                                                     <  1:4752840/150‑1 (MQ=255)
                                                                                 aTCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATAACACCCACATACAAACACAGCAGCACCAGTTCAGAAGTCAGACGTGCATCCCATACCCACCAGGTGCCCCACATCGGtttt                                                                    <  2:3895210/150‑1 (MQ=255)
                                                                                   cAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATCACACCCACATACAAAAACAGCAGCACCAGTTCAGAAGTCAGACGTGCATCCCATACCCACCAGGTGCCCCACATCGGTTTTcc                                                                  >  1:204323/1‑150 (MQ=255)
                                                                                      cACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATAACACCCACATACAAAAACAGCAGCACCAGTTCAGAAGTCAGACGTGCATCCCATACCCACCAGGTGCCCCACATCGGTTTTCCCCa                                                               <  1:1764200/150‑1 (MQ=255)
                                                                                         cAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATAACACCCACATACAAAAACAGCAGCACCAGTTCAGAAGTCAGACGTGCATCCCATACCCACCAGGTGCCCCACATCGGTTTTCCCCATga                                                            >  2:1469389/1‑149 (MQ=255)
                                                                                           ggATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATAACACCCACATACAAAAACAGCAGCACCAGTTCAGAAGTCAGACGTGCATCCCATACCCACCAGGTGCCCCACATCGGTTTTCCCCATGCag                                                          >  1:1035052/1‑150 (MQ=255)
                                                                                            gATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATAACACCCACATACAAAAACAGCAGCACCAGTTCAGAAGTCAGACGTGCATCCCATACCCACCAGGTGCCCCACATCGGTTTTCCCCATGCaga                                                         >  1:3516413/1‑150 (MQ=255)
                                                                                              tACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATAACACCCACATACAAAAACAGCAGCACCAGTTCAGAAGTCAGACGTGCATCCCATACCCACCAGGTGCCCCACATCGGTTTTCCCCATGCAGAGc                                                       >  2:4261632/1‑150 (MQ=255)
                                                                                               aCCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATCACACCCACATACAAAAACAGCAGCACCAGTTCAGAAGTCAGACGTGCATCCCATACCCACCAGGTGCCCCACATCGGTTTTCCCCATGCAGAGcc                                                      <  2:5004160/150‑1 (MQ=255)
                                                                                               aCCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATAACACCCACATACAAAAACAGCAGCACCAGTTCAGAAGTCAGACGTGCATCCCATACCCACCAGGTGCCCCACATCGGTTTTCCCCATGCAGAGcc                                                      <  1:3409965/150‑1 (MQ=255)
                                                                                               aCCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATAACACCCACATACAAAAACAGCAGCACCAGTTCAGAAGTCAGACGTGCATCCCATACCCACCAGGTGCCCCACATCGGTTTTCCCCATGCAGAGcc                                                      <  1:785254/150‑1 (MQ=255)
                                                                                                ccTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATAACACCCACATACAAAAACAGCAGCACCAGTTCAGAAGTCAGACGTGCATCCCATACCCACCAGGTGCCCCACATCGGTTTTCCCCATGCAGAGCCg                                                     <  2:1643099/150‑1 (MQ=255)
                                                                                                   gCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATAACACCCACATACAAACACAGCAGCACCAGTTCAGAAGTCAGACGTGCATCCCATACCCACCAGGTGCCCCACATCGGTTTTCCCCATGCAGAGCCGGTa                                                  <  2:102466/150‑1 (MQ=255)
                                                                                                   gCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATAACACCCACATACAAAAACAGCAGCACCAGTTCAGAAGTCAGACGTGCATCCCATACCCACCAGGTGCCCCACATCGGTTTTCCCCATGCAGAGCCGGTa                                                  <  2:1847962/150‑1 (MQ=255)
                                                                                                        aCGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATCACACCCACATACAAAAACAGCAGCACCAGTTCAGAAGTCAGACGTGCATCCCATACCCACCAGGTGCCCCACATCGGTTTTCCCCATGCAGAGCCGGTAACCAg                                             <  2:3458384/150‑1 (MQ=255)
                                                                                                           gCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATCACACCCACATACAAAAACAGCAGCACCAGTTCAGAAGTCAGACGTGCATCCCATACCCACCAGGTGCCCCACATCGGTTTTCCCCATGCAGAGCCGGTAACCAGGGc                                          <  2:4385944/150‑1 (MQ=255)
                                                                                                           gCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATAACACCCACATACAAAAACAGCAGCACCAGTTCAGAAGTCAGACGTGCATCCCATACCCACCAGGTGCCCCACATCGGTTTTCCCCATGCAGAGCCGGTAACCAGGGc                                          <  1:3254092/150‑1 (MQ=255)
                                                                                                                 cAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATCACACCCACATACAAAAACAGCAGCACCAGTTCAGAAGTCAGACGTGCATCCCATACCCACCAGGTGCCCCACATCGGTTTTCCCCATGCAGAGCCGGTAACCAGGGCAAGaaa                                    >  2:3081525/1‑150 (MQ=255)
                                                                                                                     cggcggTCGTCGAAGGCGTGCCACAGGGCAATCACACCCACATACAAAAACAGCAGCACCAGTTCAGAAGTCAGACGTGCATCCCATACCCACCAGGTGCCCCACATCGGTTTTCCCCATGCAGAGCCGGTAACCAGGGCAATAAAGGTa                                <  1:4024825/150‑1 (MQ=255)
                                                                                                                           tcgtcgAAGGCGTGCCACAGGGCAATCACACCCACATACAAAAACAGCAGCACCAGTTCAGAAGTCAGACGTGCATCCCATACCCACCAGGTGCCCCACATCGGTTTTCCCCATGCAGAGCCGGTAACCAGGGCAATAAAGGTAAACACg                          <  1:4110381/150‑1 (MQ=255)
                                                                                                                             gtcgAAGGCGTGCCACAGGGCAATCACACCCACATACAAAAACAGCAGCACCAGTTCAGAAGTCAGACGTGCATCCCATACCCACCAGGTGCCCCACATCGGTTTTCCCCATGCAGAGCCGGTAACCAGGGCAATAAAGGTAAACACGGc                        >  2:2269137/1‑150 (MQ=255)
                                                                                                                             gtcgAAGGCGTGCCACAGGGCAATAACACCCACATACAAAAACAGCAGCACCAGTTCAGAAGTCAGACGTGCATCCCATACCCACCAGGTGCCCCACATCGGTTTTCCCCATGCAGAGCCGGTAACCAGGGCAATAAAGGTAAACACGGc                        >  2:2214480/1‑150 (MQ=255)
                                                                                                                              tcgAAGGCGTGCCACAGGGCAATCACACCCACATACAAAAACAGCAGCACCAGTTCAGAAGTCAGACGTGCATCCCATACCCACCAGGTGCCCCACATCGGTTTTCCCCATGCAGAGCCGGTAACCAGGGCAATAAAGGTAAACACGGCa                       >  1:1432214/1‑150 (MQ=255)
                                                                                                                              tcgAAGGCGTGCCACAGGGCAATAACACCCACATACAAAAACAGCAGCACCAGTTCAGAAGTCAGACGTGCATCCCATACCCACCAGGTGCCCCACATCGGTTTTCCCCATGCAGAGCCGGTAACCAGGGCAATAAAGGTAAACACGGCa                       <  2:3357636/150‑1 (MQ=255)
                                                                                                                               cgAAGGCGTGCCACAGGGCAATCACACCCACATACAAAAACAGCAGCACCAGTTCAGAAGTCAGACGTGCATCCCATACCCACCAGGTGCCCCACATCGGTTTTCCCCATGCAGAGCCGGTAACCAGGGCAATAAAGGTAAACACGGCAc                      >  1:1281391/1‑150 (MQ=255)
                                                                                                                                gAAGGCGTGCCACAGGGCAATAACACCCACATACAAAAACAGCAGCACCAGTTCAGAAGTCAGACGTGCATCCCATACCCACCAGGTGCCCCACATCGGTTTTCCCCATGCAGAGCCGGTAACCAGGGCAATAAAGGTAAACACGGCAcc                     <  1:1957544/150‑1 (MQ=255)
                                                                                                                                 aaGGCGTGCCACAGGGCAATCACACCCACATACAAAAACAGCAGCACCAGTTCAGAAGTCAGACGTGCATCCCATACCCACCAGGTGCCCCACATCGGTTTTCCCCATGCAGAGCCGGTAACCAGGGCAAGAAAGGTAAACACGGCACCa                    >  2:262196/1‑150 (MQ=255)
                                                                                                                                  aGGCGTGCCACAGGGCAATCACACCCACATACAAAAACAGCAGCACCAGTTCAGAAGTCAGACGTGCATCCCATACCCACCAGGTGCCCCACATCGGTTTTCCCCATGCAGAGCCGGTAACCAGGGCAATAAAGGTAAACACGGCACCaa                   <  2:3261346/150‑1 (MQ=255)
                                                                                                                                     cGTGCCACAGGGCAATAACACCCACATACAAAAACAGCAGCACCAGTTCAGAAGTCAGACGTGCATCCCATACCCACCAGGTGCCCCACATCGGTTTTCCCCATGCAGAGCCGGTAACCAGGGCAATAAAGGTAAACACGGCACCAATgg                >  2:4881247/1‑150 (MQ=255)
                                                                                                                                        gCCACAGGGCAATCACACCCACATACAAAAACAGCAGCACCAGTTCAGAAGTCAGACGTGCATCCCATACCCACCAGGTGCCCCACATCGGTTTTCCCCATGCAGAGCCGGTAACCAGGGCAATAAAGGTAAACACGGCACCAATggggg             >  2:1705828/1‑150 (MQ=255)
                                                                                                                                          cacaGGGCAATCACACCCACATACAAAAACAGCAGCACCAGTTCAGAAGTCAGACGTGCATCCCATACCCACCAGGTGCCCCACATCGGTTTTCCCCATGCAGAGCCGGTAACCAGGGCAATAAAGGTAAACACGGCACCAATGGGGGcc           <  2:2035412/150‑1 (MQ=255)
                                                                                                                                           acaGGGCAATCACACCCACATACAAAAACAGCAGCACCAGTTCAGAAGTCAGACGTGCATCCCATACCCACCAGGTGCCCCACATCGGTTTTCCCCATGCAGAGCCGGTAACCAGGGCAATAAAGGTAAACACGGCACCAATGGGGGCCa          >  1:4013568/1‑150 (MQ=255)
                                                                                                                                           acaGGGCAATCACACCCACATACAAAAACAGCAGCACCAGTTCAGAAGTCAGACGTGCATCCCATACCCACCAGGTGCCCCACATCGGTTTTCCCCATGCAGAGCCGGTAACCAGGGCAATAAAGGTAAACACGGCACCAATGGGGGCCa          >  1:1742032/1‑150 (MQ=255)
                                                                                                                                              gggCAATCACACCCACATACAAAAACAGCAGCACCAGTTCAGAAGTCAGACGTGCATCCCATACCCACCAGGTGCCCCACATCGGTTTTCCCCATGCAGAGCCGGTAACCAGGGCAATAAAGGTAAACACGGCACCAATGGGGGCCATcg       >  2:1569004/1‑150 (MQ=255)
                                                                                                                                              gggCAATAACACCCACATACAAAAACAGCAGCACCAGTTCAGAAGTCAGACGTGCATCCCATACCCACCAGGTGCCCCACATCGGTTTTCCCCATGCAGAGCCGGTAACCAGGGCAATAAAGGTAAACACGGCACCAATGGGGGCCATcg       <  1:2574075/150‑1 (MQ=255)
                                                                                                                                               ggCAATCACACCCACATACAAAAACAGCAGCACCAGTTCAGAAGTCAGACGTGCATCCCATACCCACCAGGTGCCCCACATCGGTTTTCCCCATGCAGAGCCGGTAACCAGGGCAATAAAGGTAAACACGGCACCAATGGGGGCCATcgc      >  1:4608220/1‑150 (MQ=255)
                                                                                                                                               ggCAATAACACCCACATACAAAAACAGCAGCACCAGTTCAGAAGTCAGACGTGCATCCCATACCCACCAGGTGCCCCACATCGGTTTTCCCCATGCAGAGCCGGTAACCAGGGCAATAAAGGTAAACACGGCACCAATGGGGGCCATcgc      >  2:2394119/1‑150 (MQ=255)
                                                                                                                                                   aTCACACCCACATACAAAAACAGCAGCACCAGTTCAGAAGTCAGACGTGCATCCCATACCCACCAGGAGCCCCACAACGGTTTTCCCCATGCAGAGCCGGTAACCAGGGCACTAAAGGTAAACACGGCACCAATGGGGGCCATCgccgcc  >  1:2638906/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
TGCTGCATCCGCGTTGATCCCTGATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATCACACCCACATACAAAAACAGCAGCACCAGTTCAGAAGTCAGACGTGCATCCCATACCCACCAGGTGCCCCACATCGGTTTTCCCCATGCAGAGCCGGTAACCAGGGCAATAAAGGTAAACACGGCACCAATGGGGGCCATCGCCGCC  >  CP009273/2289225‑2289521

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: