Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,114,329 C→A 100% L58F (TTG→TTT)  yceJ ← putative cytochrome b561

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,114,3290CA87.9% 671.5 / 92.8 264L58F (TTG→TTT) yceJputative cytochrome b561
Reads supporting (aligned to +/- strand):  ref base C (16/16);  new base A (117/115);  total (133/131)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01
Rejected as polymorphism: Polymorphic base substitution creates a homopolymer stretch.

Only 100 of 266 total aligned reads displayed.

GCTGATGCCGATGGCGAAAAGCAATAGGTACAATGCGAGATGTCCAGCTCTTGCGGCAAGACGAGTCATGGGCGAATAACTCGGCAGCGGACCCGGCGGTGGAGATATGACACGCCACAGAACGCGAATAACCAGCCCCATCATTAACAAAATACCGATGCTTTTATGCAGTTCGGGTGCTTTGTGATACCAGCCATCGTAATAACTGAGCGTGACCATCCACAGACCTAAAGCAAACATGCCGTAGACAATAATTGCGCTTAGCCAATGGAAGGCCGCAGATATAACACCATAG  >  CP009273/1114182‑1114476
                                                                                                                                                   |                                                                                                                                                   
    aTGCCGATGGCGAAAAGCAATAGGTACACTGCGAGATGTCCAGCTCTTGCGGCAAGACGAGTCATGGGCGAATAACTCGGCAGCGGACCCGGCGGTGGAGATATGACACGCCACAGAACGCGAATAACCAGCCCCATCATTAAAAAAAt                                                                                                                                                <  2:2153061/149‑1 (MQ=255)
     tGCCGAGGGCGAAAAGCACTAGGTACAATGCGAGATGTCCAGCTCTTGCGGCAAGACGAGTCATGGGCGAATAACTCGGCAGCGGACCCGGCGGTGGAGATATGACACGCCACAGAACGCGAATAACCAG‑CCCATCATTAAAAAAATAc                                                                                                                                              <  2:2544045/149‑1 (MQ=255)
      gCCGATGGCGAAAAGCAATAGGTACAATGCGAGATGTCCAGCTCTTGCGGCAAGACGAGTCATGGGCGAATAACTCGGCAGCGGACCCGGCGGTGGAGATATGACACGCCACAGAACGCGAATAACCAGCCCCATCATTAAAAAAATAc                                                                                                                                              <  1:1279044/149‑1 (MQ=255)
        cGATGGCGAAAAGCAATAGGTACAATGCGAGATGTCCAGCTCTTGCGGCAAGACGAGTCATGGGCGAATAACTCGGCAGCGGACCCGGCGGTGGAGATATGACACGCCACAGAACGCGAATAACCAGCCCCATCATTAAAAAAATACCg                                                                                                                                            <  1:104863/149‑1 (MQ=255)
         gATGGCGAAAAGCAATAGGTACAATGCGAGATGTCCAGCTCTTGCGGCAAGACGAGTCATGGGCGAATAACTCGGCAGCGGACCCGGCGGTGGAGATATGACACGCCACAGAACGCGAATAACCAGCCCCATCATTAAAAAAATACCGa                                                                                                                                           >  2:2205266/1‑149 (MQ=255)
          aTGGCGAAAAGCAATAGGTACAATGCGAGATGTCCAGCTCTTGCGGCAAGACGAGTCATGGGCGAATAACTCGGCAGCGGACCCGGCGGTGGAGATATGACACGCCACAGAACGCGAATAACCAGCCCCATCATTAAAAAAATACCGAt                                                                                                                                          >  2:931519/1‑149 (MQ=255)
          aTGGCGAAAAGCAATAGGTACAATGCGAGATGTCCAGCTCTTGCGGCAAGACGAGTCATGGGCGAATAACTCGGCAGCGGACCCGGCGGTGGAGATATGACACGCCACAGAACGCGAATAACCAGCCCCATCATTAAAAAAATACCGAt                                                                                                                                          <  2:4751352/149‑1 (MQ=255)
            ggCGAAAAGCAATAGGTACAATGCGAGATGTCCAGCTCTTGCGGCAAGACGAGTCATGGGCGAATAACTCGGCAGCGGACCCGGCGGTGGAGATATGACACGCCACAGAACGCGAATAACCAGCCCCATCATTAACAAAATACCGATGc                                                                                                                                        <  1:4364243/149‑1 (MQ=255)
            ggCGAAAAGCAATAGGTACAATGCGAGATGTCCAGCTCTTGCGGCAAGACGAGTCATGGGCGAATAACTCGGCAGCGGACCCGGCGGTGGAGATATGACACGCCACAGAACGCGAATAACCAGCCCCATCATTAACAAAATACCGATGc                                                                                                                                        <  1:1625887/149‑1 (MQ=255)
              cGAAAAGCAATAGGTACAATGCGAGATGTCCAGCTCTTGCGGCAAGACGAGTCATGGGCGAATAACTCGGCAGCGGACCCGGCGGTGGAGATATGACACGCCACAGAACGCGAATAACCAGCCCCATCATTAACAAAATACCGATGCtt                                                                                                                                      <  1:2869846/149‑1 (MQ=255)
              cGAAAAGCAATAGGTACAATGCGAGATGTCCAGCTCTTGCGGCAAGACGAGTCATGGGCGAATAACTCGGCAGCGGACCCGGCGGTGGAGATATGACACGCCACAGAACGCGAATAACCAGCCCCATCATTAAAAAAATACCGATGCtt                                                                                                                                      <  1:1089404/149‑1 (MQ=255)
              cGAAAAGCAATAGGTACAATGCGAGATGTCCAGCTCTTGCGGCAAGACGAGTCATGGGCGAATAACTCGGCAGCGGACCCGGCGGTGGAGATATGACACGCCACAGAACGCGAATAACCAGCCCCATCATTAAAAAAATACCGATGCtt                                                                                                                                      <  1:4199010/149‑1 (MQ=255)
               gAAAAGCAATAGGTACAATGCGAGATGTCCAGCTCTTGCGGCAAGACGAGTCATGGGCGAATAACTCGGCAGCGGACCCGGCGGTGGAGATATGACACGCCACAGAACGCGAATAACCAGCCCCATCATTAAAAAAATACCGATGCttt                                                                                                                                     <  2:2839908/149‑1 (MQ=255)
                aaaaGCAATAGGTACAATGCGAGATGTCCAGCTCTTGCGGCAAGACGAGTCATGGGCGAATAACTCGGCAGCGGACCCGGCGGTGGAGATATGACACGCCACAGAACGCGAATAACCAGCCCCATCATTAAAAAAATACCGATGCtttt                                                                                                                                    >  1:3510740/1‑149 (MQ=255)
                aaaaGCAATAGGTACAATGCGAGATGTCCAGCTCTTGCGGCAAGACGAGTCATGGGCGAATAACTCGGCAGCGGACCCGGCGGTGGAGATATGACACGCCACAGAACGCGAATAACCAGCCCCATCATTAAAAAAATACCGATGCtttt                                                                                                                                    <  2:3619808/149‑1 (MQ=255)
                 aaaGCAATAGGTACAATGCGAGATGTCCAGCTCTTGCGGCAAGACGAGTCATGGGCGAATAACTCGGCAGCGGACCCGGCGGTGGAGATATGACACGCCACAGAACGCGAATAACCAGCCCCATCATTAAAAAAATACCGATGCTTTTa                                                                                                                                   <  1:536812/149‑1 (MQ=255)
                 aaaGCAATAGGTACAATGCGAGATGTCCAGCTCTTGCGGCAAGACGAGTCATGGGCGAATAACTCGGCAGCGGACCCGGCGGTGGAGATATGACACGCCACAGAACGCGAATAACCAGCCCCATCATTAAAAAAATACCGATGCTTTTa                                                                                                                                   <  2:1773963/149‑1 (MQ=255)
                  aaGCAATAGGTACAATGCGAGATGTCCAGCTCTTGCGGCAAGACGAGTCATGGGCGAATAACTCGGCAGCGGACCCGGCGGTGGAGATATGACACGCCACAGAACGCGAATAACCAGCCCCATCATTAAAAAAATACCGATGCTTTTAt                                                                                                                                  >  2:4612320/1‑149 (MQ=255)
                    gCAATAGGTACAATGCGAGATGTCCAGCTCTTGCGGCAAGACGAGTCATGGGCGAATAACTCGGCAGCGGACCCGGCGGTGGAGATATGACACGCCACAGAACGCGAATAACCAGCCCCATCATTAAAAAAATACCGATGCTTTTATGc                                                                                                                                <  1:3568901/149‑1 (MQ=255)
                      aaTAGGTACAATGCGAGATGTCCAGCTCTTGCGGCAAGACGAGTCATGGGCGAATAACTCGGCAGCGGACCCGGCGGTGGAGATATGACACGCCACAGAACGCGAATAACCAGCCCCATCATTAAAAAAATACCGATGCTTTTATGCAg                                                                                                                              >  1:3717223/1‑149 (MQ=255)
                      aaTAGGTACAATGCGAGATGTCCAGCTCTTGCGGCAAGACGAGTCATGGGCGAATAACTCGGCAGCGGACCCGGCGGTGGAGATATGACACGCCACAGAACGCGAATAACCAGCCCCATCATTAAAAAAATACCGATGCTTTTATGCAg                                                                                                                              <  2:350214/149‑1 (MQ=255)
                         aGGTACAATGCGAGATGTCCAGCTCTTGCGGCAAGACGAGTCATGGGCGAATAACTCGGCAGCGGACCCGGCGGTGGAGATATGACACGCCACAGAACGCGAATAACCAGCCCCATCATTAACAAAATACCGATGCTTTTATGCAGTTc                                                                                                                           <  2:4665298/149‑1 (MQ=255)
                         aGGTACAATGCGAGATGTCCAGCTCTTGCGGCAAGACGAGTCATGGGCGAATAACTCGGCAGCGGACCCGGCGGTGGAGATATGACACGCCACAGAACGCGAATAACCAGCCCCATCATTAAAAAAATACCGATGCTTTTATGCAGTTc                                                                                                                           >  1:1444924/1‑149 (MQ=255)
                             aCAATGCGAGATGTCCAGCTCTTGCGGCAAGACGAGTCATGGGCGAATAACTCGGCAGCGGACCCGGCGGTGGAGATATGACACGCCACAGAACGCGAATAACCAGCCCCATCATTAAAAAAATACCGATGCTTTTATGCAGTTCGGGt                                                                                                                       >  2:1337217/1‑149 (MQ=255)
                             aCAATGCGAGATGTCCAGCTCTTGCGGCAAGACGAGTCATGGGCGAATAACTCGGCAGCGGACCCGGCGGTGGAGATATGACACGCCACAGAACGCGAATAACCAGCCCCATCATTAAAAAAATACCGATGCTTTTATGCAGTTCGGGt                                                                                                                       <  2:169675/149‑1 (MQ=255)
                               aaTGCGAGATGTCCAGCTCTTGCGGCAAGACGAGTCATGGGCGAATAACTCGGCAGCGGACCCGGCGGTGGAGATATGACACGCCACAGAACGCGAATAACCAGCCCCATCATTAAAAAAATACCGATGCTTTTATGCAGTTCGGGTGc                                                                                                                     >  1:4724641/1‑149 (MQ=255)
                               aaTGCGAGATGTCCAGCTCTTGCGGCAAGACGAGTCATGGGCGAATAACTCGGCAGCGGACCCGGCGGTGGAGATATGACACGCCACAGAACGCGAATAACCAGCCCCATCATTAAAAAAATACCGATGCTTTTATGCAGTTCGGGTGc                                                                                                                     >  1:752305/1‑149 (MQ=255)
                                 tGCGAGATGTCCAGCTCTTGCGGCAAGACGAGTCATGGGCGAATAACTCGGCAGCGGACCCGGCGGTGGAGATATGACACGCCACAGAACGCGAATAACCAGCCCCATCATTAAAAAAATACCGATGCTTTTATGCAGTTCGGGTGCtt                                                                                                                   >  2:4587439/1‑149 (MQ=255)
                                 tGCGAGATGTCCAGCTCTTGCGGCAAGACGAGTCATGGGCGAATAACTCGGCAGCGGACCCGGCGGTGGAGATATGACACGCCACAGAACGCGAATAACCAGCCCCATCATTAAAAAAATACCGATGCTTTTATGCAGTTCGGGTGCtt                                                                                                                   <  2:814503/149‑1 (MQ=255)
                                  gCGAGATGTCCAGCTCTTGCGGCAAGACGAGTCATGGGCGAATAACTCGGCAGCGGACCCGGCGGTGGAGATATGACACGCCACAGAACGCGAATAACCAGCCCCATCATTAAAAAAATACCGATGCTTTTATGCAGTTCGGGTGCttt                                                                                                                  >  2:1081236/1‑149 (MQ=255)
                                    gagaTGTCCAGCTCTTGCGGCAAGACGAGTCATGGGCGAATAACTCGGCAGCGGACCCGGCGGTGGAGATATGACACGCCACAGAACGCGAATAACCAGCCCCATCATTAAAAAAATACCGATGCTTTTATGCAGTTCGGGTGCTTtgt                                                                                                                >  1:2966872/1‑149 (MQ=255)
                                       aTGTCCAGCTCTTGCGGCAAGACGAGTCATGGGCGAATAACTCGGCAGCGGACCCGGCGGTGGAGATATGACACGCCACAGAACGCGAATAACCAGCCCCATCATTAACAAAATACCGATGCTTTTATGCAGTTCGGGTGCTTTGTGAt                                                                                                             >  1:3972746/1‑149 (MQ=255)
                                       aTGTCCAGCTCTTGCGGCAAGACGAGTCATGGGCGAATAACTCGGCAGCGGACCCGGCGGTGGAGATATGACACGCCACAGAACGCGAATAACCAGCCCCATCATTAACAAAATACCGATGCTTTTATGCAGTTCGGGTGCTTTGTGAt                                                                                                             >  1:1583577/1‑149 (MQ=255)
                                          tCCAGCTCTTGCGGCAAGACGAGTCATGGGCGAATAACTCGGCAGCGGACCCGGCGGTGGAGATATGACACGCCACAGAACGCGAATAACCAGCCCCATCATTAAAAAAATACCGATGCTTTTATGCAGTTCGGGTGCTTTGTGATAcc                                                                                                          <  2:150987/149‑1 (MQ=255)
                                           ccAGCTCTTGCGGCAAGACGAGTCATGGGCGAATAACTCGGCAGCGGACCCGGCGGTGGAGATATGACACGCCACAGAACGCGAATAACCAGCCCCATCATTAAAAAAATACCGATGCTTTTATGCAGTTCGGGTGCTTTGTGATACCa                                                                                                         <  2:39133/149‑1 (MQ=255)
                                           ccAGCTCTTGCGGCAAGACGAGTCATGGGCGAATAACTCGGCAGCGGACCCGGCGGTGGAGATATGACACGCCACAGAACGCGAATAACCAGCCCCATCATTAAAAAAATACCGATGCTTTTATGCAGTTCGGGTGCTTTGTGATACCa                                                                                                         >  1:1786646/1‑149 (MQ=255)
                                            cAGCTCTTGCGGCAAGACGAGTCATGGGCGAATAACTCGGCAGCGGACCCGGCGGTGGAGATATGACACGCCACAGAACGCGAATAACCAGCCCCATCATTAAAAAAATACCGATGCTTTTATGCAGTTCGGGTGCTTTGTGATACCAg                                                                                                        <  2:4204700/149‑1 (MQ=255)
                                              gCTCTTGCGGCAAGACGAGTCATGGGCGAATAACTCGGCAGCGGACCCGGCGGTGGAGATATGACACGCCACAGAACGCGAATAACCAGCCCCATCATTAAAAAAATACCGATGCTTTTATGCAGTTCGGGTGCTTTGTGGTACCAGcc                                                                                                      <  2:295120/149‑1 (MQ=255)
                                              gCTCTTGCGGCAAGACGAGTCATGGGCGAATAACTCGGCAGCGGACCCGGCGGTGGAGATATGACACGCCACAGAACGCGAATAACCAGCCCCATCATTAAAAAAATACCGATGCTTTTATGCAGTTCGGGTGCTTTGTGATACCAGcc                                                                                                      >  2:3798390/1‑149 (MQ=255)
                                              gCTCTTGCGGCAAGACGAGTCATGGGCGAATAACTCGGCAGCGGACCCGGCGGTGGAGATATGACACGCCACAGAACGCGAATAACCAGCCCCATCATTAAAAAAATACCGATGCTTTTATGCAGTTCGGGTGCTTTGTGATACCAGcc                                                                                                      >  2:40288/1‑149 (MQ=255)
                                                     cGGCAAGACGAGTCATGGGCGAATAACTCGGCAGCGGACCCGGCGGTGGAGATATGACACGCCACAGAACGCGAATAACCAGCCCCATCATTAAAAAAATACCGATGCTTTTATGCAGTTCGGGTGCTTTGTGATACCAGCCATCGtaa                                                                                               <  2:594923/149‑1 (MQ=255)
                                                      ggCAAGACGAGTCATGGGCGAATAACTCGGCAGCGGACCCGGCGGTGGAGATATGACACGCCACAGAACGCGAATAACCAGCCCCATCATTAAAAAAATACCGATGCTTTTATGCAGTTCGGGTGCTTTGTGATACCAGCCATCGtaat                                                                                              >  2:3222863/1‑149 (MQ=255)
                                                        cAAGACGAGTCATGGGCGAATAACTCGGCAGCGGACCCGGCGGTGGAGATATGACACGCCACAGAACGCGAATAACCAGCCCCATCATTAAAAAAATACCGATGCTTTTATGCAGTTCGGGTGCTTTGTGATACCAGCCATCGtaataa                                                                                            <  1:3898196/149‑1 (MQ=255)
                                                         aaGACGAGTCATGGGCGAATAACTCGGCAGCGGACCCGGCGGTGGAGATATGACACGCCACAGAACGCGAATAACCAGCCCCATCATTAAAAAAATACCGATGCTTTTATGCAGTTCGGGTGCTTTGTGATACCAGCCATCGTAATAAc                                                                                           >  2:3991390/1‑149 (MQ=255)
                                                             cGAGTCATGGGCGAATAACTCGGCAGCGGACCCGGCGGTGGAGATATGACACGCCACAGAACGCGAATAACCAGCCCCATCATTAAAAAAATACCGATGCTTTTATGCAGTTCGGGTGCTTTGTGATACCAGCCATCGTAATAACTGAg                                                                                       <  1:3343314/149‑1 (MQ=255)
                                                             cGAGTCATGGGCGAATAACTCGGCAGCGGACCCGGCGGTGGAGATATGACACGCCACAGAACGCGAATAACCAGCCCCATCATTAAAAAAATACCGATGCTTTTATGCAGTTCGGGTGCTTTGTGATACCAGCCATCGTAATAACTGAg                                                                                       >  2:2301614/1‑149 (MQ=255)
                                                               aGTCATGGGCGAATAACTCGGCAGCGGACCCGGCGGTGGAGATATGACACGCCACAGAACGCGAATAACCAGCCCCATCATTAAAAAAATACCGATGCTTTTATGCAGTTCGGGTGCTTTGTGATACCAGCCATCGTAATAACTGAGCg                                                                                     <  1:1140581/149‑1 (MQ=255)
                                                                gTCATGGGCGAATAACTCGGCAGCGGACCCGGCGGTGGAGATATGACACGCCACAGAACGCGAATCACCAGCCCCATCATTAAAAAAATACCGATGCTTTTATGCAGTTCGGGTGCTTTGTGATACCAGCCATCGTAATAACTGAGCGt                                                                                    <  2:1200488/149‑1 (MQ=255)
                                                                gTCATGGGCGAATAACTCGGCAGCGGACCCGGCGGTGGAGATATGACACGCCACAGAACGCGAATAACCAGCCCCATCATTAACAAAATACCGATGCTTTTATGCAGTTCGGGTGCTTTGTGATACCAGCCATCGTAATAACTGAGCGt                                                                                    >  2:1206822/1‑149 (MQ=255)
                                                                    tGGGCGAATAACTCGGCAGCGGACCCGGCGGTGGAGATATGACACGCCACAGAACGCGAATAACCAGCCCCATCATTAAAAAAATACCGATGCTTTTATGCAGTTCGGGTGCTTTGTGATACCAGCCATCGTAATAACTGAGCGTGAcc                                                                                >  2:1420490/1‑149 (MQ=255)
                                                                     gggCGAATAACTCGGCAGCGGACCCGGCGGTGGAGATAAGACACGCCACAGAACGCGAATAACCAGCCCCATCATTAAAAAAATACCGATGCTTTTATGCAGTTCGGGTGCTTTGTGATACCAGCCATCGTAATAACTGAGCGTGAccc                                                                               >  1:4428312/1‑148 (MQ=255)
                                                                         gAATAACTCGGCAGCGGACCCGGCGGTGGAGATATGACACGCCACAGAACGCGAATAACCAGCCCCATCATTAAAAAAATACCGATGCTTTTATGCAGTTCGGGTGCTTTGTGATGCCAGCCATCGTAATAACTGAGCGTGACCATCca                                                                           <  1:2486475/149‑1 (MQ=255)
                                                                         gAATAACTCGGCAGCGGACCCGGCGGTGGAGATATGACACGCCACAGAACGCGAATAACCAGCCCCATCATTAAAAAAATACCGATGCTTTTATGCAGTTCGGGTGCTTTGTGATACCAGCCATCGTAATAACTGAGCGTGACCATCca                                                                           <  2:66044/149‑1 (MQ=255)
                                                                          aaTAACTCGGCAGCGGACCCGGCGGTGGAGATATGACACGCCACAGAACGCGAATAACCAGCCCCATCATTAACAAAATACCGATGCTTTTATGCAGTTCGGGTGCTTTGTGATACCAGCCATCGTAATAACTGAGCGTGACCATCcac                                                                          <  1:3149124/149‑1 (MQ=255)
                                                                           aTAACTCGGCAGCGGACCCGGCGGTGGAGATATGACACGCCACAGAACGCGAATAACCAGCCCCATCATTAAAAAAATACCGATGCTTTTATGCAGTTCGGGTGCTTTGTGATACCAGCCATCGTAATAACTGAGCGTGACCATCcaca                                                                         <  2:547103/149‑1 (MQ=255)
                                                                             aaCTCGGCAGCGGACCCGGCGGTGGAGATATGACACGCCACAGAACGCGAATAACCAGCCCCATCATTAAAAAAATACCGATGCTTTTATGCAGTTCGGGTGCTTTGTGATACCAGCCATCGTAATAACTGAGCGTGACCATCCACTGa                                                                       <  2:4282226/149‑1 (MQ=255)
                                                                              aCTCGGCAGCGGACCCGGCGGTGGAGATATGACACGCCACAGAACGCGAATAACCAGCCCCATCATTAAAAAAATACCGATGCTTTTATGCAGTTCGGGTGCTTTGTGATACCAGCCATCGTAATAACTGAGCGTGACCATCCACAGAc                                                                      >  1:825148/1‑149 (MQ=255)
                                                                               cTCGGCAGCGGACCCGGCGGTGGAGATATGACACGCCACAGAACGCGAATAACCAGCCCCATCATTAAAAAAATACCGATGCTTTTATGCAGTTCGGGTGCTTTGTGATACCAGCCATCGTAATAACTGAGCGTGACCATCCACAGAcc                                                                     <  2:3393127/149‑1 (MQ=255)
                                                                                tCGGCAGCGGACCCGGCGGTGGAGATATGACACGCCACAGAACGCGAATAACCAGCCCCATCATTAAAAAAATACCGATGCTTTTATGCAGTTCGGGTGCTTTGTGATACCAGCCATCGTAATAACTGAGCGTGACCATCCACAGACCt                                                                    >  1:2580831/1‑149 (MQ=255)
                                                                                tCGGCAGCGGACCCGGCGGTGGAGATATGACACGCCACAGAACGCGAATAACCAGCCCCATCATTAAAAAAATACCGATGCTTTTATGCAGTTCGGGTGCTTTGTGATACCAGCCATCGTAATAACTGAGCGTGACCATCCACAGACCt                                                                    <  1:1389340/149‑1 (MQ=255)
                                                                                  ggCAGCGGACCCGGCGGTGGAGATATGACACGCCACAGAACGCGAATAACCAGCCCCATCATTAACAAAATACCGATGCTTTTATGCAGTTCGGGTGCTTTGTGATACCAGCCATCGTAATAACTGAGCGTGACCATCCACAGACCTaa                                                                  <  1:1941648/149‑1 (MQ=255)
                                                                                   gCAGCGGACCCGGCGGTGGAGATATGACACGCCACAGAACGCGAATAACCAGCCCCATCATTAAAAAAATACCGATGCTTTTATGCAGTTCGGGTGCTTTGTGATACCAGCCATCGTAATAACTGAGCGTGACCATCCACAGACCTaaa                                                                 >  1:2730469/1‑149 (MQ=255)
                                                                                    cAGTGGACCCGGCGGTGGAGATATGACACGCCACAGAACGCGAATAACCAGCCCCATCATTAAAAAAATACCGATGCTTTTATGCAGTTCGGGTGCTTTGTGATACCAGCCATCGTAATAACTGAGCGTGACCATCCACAGACCTAAAg                                                                <  2:1603454/149‑1 (MQ=255)
                                                                                    cAGCGGACCCGGCGGTGGAGATATGACACGCCACAGAACGCGAATAACCAGCCCCATCATTAAAAAAATACCGATGCTTTTATGCAGTTCGGGTGCTTTGTGATACCAGCCATCGTAATAACTGAGCGTGACCATCCACAGACCTAAAg                                                                <  1:1400214/149‑1 (MQ=255)
                                                                                        ggACCCGGCGGTGGAGATATGACACGCCACAGAACGCGAATAACCAGCCCCATCATTAAAAAAATACCGATGCTTTTATGCAGTTCGGGTGCTTTGTGATACCAGCCATCGTAATAACTGAGCGTGACCATCCACAGACCTAAAGCaaa                                                            <  1:4294335/149‑1 (MQ=255)
                                                                                         gACCCGGCGGTGGAGATATGACACGCCACAGAACGCGAATAACCAGCCCCATCATTAAAAAAATACCGATGCTTTTATGCAGTTCGGGTGCTTTGTGATACCAGCCATCGTAATAACTGAGCGTGACCATCCACAGACCTAAAGCAAAc                                                           <  2:1908588/149‑1 (MQ=255)
                                                                                           cccGGCGGTGGAGATATGACACGCCACAGAACGCGAATAACCAGCCCCATCATTAAAAAAATACCGATGCTTTTATGCAGTTCGGGTGCTTTGTGATACCAGCCATCGTAATAACTGAGCGTGACCATCCACAGACCTAAAGCAAACAt                                                         >  1:3259391/1‑149 (MQ=255)
                                                                                           cccGGCGGTGGAGATAAGACACGCCACAGAACGCGAATAACCAGCCCCATCATTAAAAAAATACCGATGCTTTTATGCAGTTAGGGTGCTTTTTGATACCAGCCATCGTAATAACAGAGCGTGACCATCAACAGAACTAAAGCAAACaa                                                         >  1:870515/1‑148 (MQ=255)
                                                                                                        atatGACACGCCACAGAACGCGAATAACCAGCCCCATCATTAAAAAAATACCGATGCTTTTATGCAGTTCGGGTGCTTTGTGATACCAGCCATCGTAATAACTGAGCGTGACCATCCACAGACCTAAAGCAAACATGCCGTAGACaata                                            <  2:4453468/149‑1 (MQ=255)
                                                                                                          atGACACGCCACAGAACGCGAATAACCAGCCCCATCATTAAAAAAATACCGATGCTTTTATGCAGTTCGGGTGCTTTGTGATACCAGCCATCGTAATAACTGAGCGTGACCATCCACAGACCTAAAGCAAACATGCCGTAGACaataat                                          <  2:3075950/149‑1 (MQ=255)
                                                                                                            gACACGCCACAGAACGCGAATAACCAGCCCCATCATTAACAAAATACCGATGCTTTTATGCAGTTCGGGTGCTTTGTGATACCAGCCATCGTAATAACTGAGCGTGACCATCCACAGACCTAAAGCAAACATGCCGTAGACAATAATTg                                        >  2:645254/1‑149 (MQ=255)
                                                                                                              cacGCCACAGAACGCGAATAACCAGCCCCATCATTAAAAAAATACCGATGCTTTTATGCAGTTCGGGTGCTTTGTGATACCAGCCATCGTAATAACTGAGCGTGACCATCCACAGACCTAAAGCAAACATGCCGTAGACAATAATTgcg                                      <  1:3439083/149‑1 (MQ=255)
                                                                                                               acGCCACAGAACGCGAATAACCAGCCCCATCATTAAAAAAATACCGATGCTTTTATGCAGTTCGGGTGCTTTGTGATACCAGCCATCGTAATAACTGAGCGTGACCATCCACAGACCTAAAGCAAACATGCCGTAGACAATAATTgcgc                                     <  1:3259630/149‑1 (MQ=255)
                                                                                                                 gCCACAGAACGCGAATAACCAGCCCCATCATTAAAAAAATACCGATGCTTTTATGCAGTTCGGGTGCTTTGTGATACCAGCCATCGTAATAACTGAGCGTGACCATCCACAGACCTAAAGCAAACATGCCGTAGACAATAATTGCGCtt                                   >  2:2084464/1‑149 (MQ=255)
                                                                                                                 gCCACAGAACGCGAATAACCAGCCCCATCATTAAAAAAATACCGATGCTTTTATGCAGTTCGGGTGCTTTGTGATACCAGCCATCGTAATAACTGAGCGTGACCATCCACAGACCTAAAGCAAACATGCCGTAGACAATAATTGCGCtt                                   <  2:2054204/149‑1 (MQ=255)
                                                                                                                    acaGAACGCGAATAACCAGCCCCATCATTAAAAAAATACCGATGCTTTTATGCAGTTCGGGTGCTTTGTGATACCAGCCATCGTAATAACTGAGCGTGACCATCCACAGACCTAAAGCAAACATGCCGTAGACAATAATTGCGCTTAGc                                <  1:1337126/149‑1 (MQ=255)
                                                                                                                    acaGAACGCGAATAACCAGCCCCATCATTAAAAAAATACCGATGCTTTTATGCAGTTCGGGTGCTTTGTGATACCAGCCATCGTAATAACTGAGCGTGACCATCCACAGACCTAAAGCAAACATGCCGTAGACAATAATTGCGCTTAGc                                >  2:3532563/1‑149 (MQ=255)
                                                                                                                    acaGAACGCGAATAACCAGCCCCATCATTAAAAAAATACCGATGCTTTTATGCAGTTCGGGTGCTTTGTGATACAAGCCATCGTAATAACTGAGCGTGACCATCCACAGACCTAAAGCAAACATGCCGTAGACAATAATTGCGCTTAGc                                >  1:27939/1‑149 (MQ=255)
                                                                                                                     caGAACGCGAATAACCAGCCCCATCATTAAAAAAATACCGATGCTTTTATGCAGTTCGGGTGCTTTGTGATACCAGCCATCGTAATAACTGAGCGTGACCATCCACAGACCTAAAGCAAACATGCCGTAGACAATAATTGCGCTTAGcc                               >  2:4024004/1‑149 (MQ=255)
                                                                                                                       gAACGCGAATAACCAGCCCCATCATTAAAAAAATACCGATGCTTTTATGCAGTTCGGGTGCTTTGTGATACCAGCCATCGTAATAACTGAGCGTGACCATCCACAGACCTAAAGCAAACATGCCGTAGACAATAATTGCGCTTAGCCaa                             >  2:524592/1‑149 (MQ=255)
                                                                                                                       gAACGCGAATAACCAGCCCCATCATTAAAAAAATACCGATGCTTTTATGCAGTTCGGGTGCTTTGTGATACCAGCCATCGTAATAACTGAGCGTGACCATCCACAGACCTAAAGCAAACATGCCGTAGACAATAATTGCGCTTAGCCaa                             <  2:332383/149‑1 (MQ=255)
                                                                                                                         aCGCGAATAACCAGCCCCATCATTAAAAAAATACCGATGCTTTTATGCAGTTCGGGTGCTTTGTGATACCAGCCATCGTAATAACTGAGCGTGACCATCCACAGACCTAAAGCAAACATGCCGTAGACAATAATTGCGCTTAGCCAATg                           >  1:2673415/1‑149 (MQ=255)
                                                                                                                           gcgAATAACCAGCCCCATCATTAACAAAATACCGATGCTTTTATGCAGTTCGGGTGCTTTGTGATACCAGCCATCGTAATAACTGAGCGTGACCATCCACAGACCTAAAGCAAACATGCCGTAGACAATAATTGCGCTTAGCCAATGGa                         <  2:1805982/149‑1 (MQ=255)
                                                                                                                           gcgAATAACCAGCCCCATCATTAAAAAAATACCGATGCTTTTATGCAGTTCGGGTGCTTTGTGATACCAGCCATCGTAATAACTGAGCGTGACCATCCACAGACCTAAAGCAAACATGCCGTAGACAATAATTGCGCTTAGCCAATGGa                         <  2:4095862/149‑1 (MQ=255)
                                                                                                                             gAATAACCAGCCCCATCATTAAAAAAATACCGATGCTTTTATGCAGTTCGGGTGCTTTGTGATACCAGCCATCGTAATAACTGAGCGTGACCATCCACAGACCTAAAGCAAACATGCCGTAGACAATAATTGCGCTTAGCCAATGGAAg                       >  1:3052427/1‑149 (MQ=255)
                                                                                                                              aaTAACCAGCCCCATCATTAACAAAATACCGATGCTTTTATGCAGTTCGGGTGCTTTGTGATACCAGCCATCGTAATAACTGAGCGTGACCATCCACAGACCTAAAGCAAACATGCCGTAGACAATAATTGCGCTTAGCCAATGGAAgg                      >  2:112551/1‑149 (MQ=255)
                                                                                                                               aTAACCAGCCCCATCATTAAAAAAATACCGATGCTTTTATGCAGTTCGGGTGCTTTTTGATACCAGCCATAGTAATAACTGAGAGAGACCATCCACAGACCTAAAGCAAACATGCCGTAGACAATAATTTCGCTTAGCCAATGGAAGGc                     >  1:3121880/1‑149 (MQ=255)
                                                                                                                                tAACCAGCCCCATCATTAAAAAAATACCGATGCTTTTATGCAGTTCGGGTGCTTTGTGATACCAGCCATCGTAATAACTGAGCGTGACCATCCACAGACCTAAAGCAAACATGCCGTAGACAATAATTGCGCTTAGCCAATGGAAGGcc                    <  2:1862967/149‑1 (MQ=255)
                                                                                                                                    cAGCCCCATCATTAAAAAAATACCGATGCTTTTATGCAGTTCGGGTGCTTTGTGATACCAGCCATCGTAATAACTGAGCGTGACCATCCACAGACCTAAAGCAAACATGCCGTAGACAATAATTGCGCTTAGCCAATGGAAGGCCGCAg                >  1:2079090/1‑149 (MQ=255)
                                                                                                                                      gCCCCATCATTAAAAAAATACCGATGCTTTTATGCAGTTCGGGTGCTTTGTGATACCAGCCATCTTAATAACTGAGCGTGACCATCCACAGACCTAAAGCAAACATGCCGTAGACAATAATTGCGCTTAGCCAATGGAAGGCCGCAGat              >  1:1168586/1‑149 (MQ=255)
                                                                                                                                      gCCCCATCATTAAAAAAATACCGATGCTTTTATGCAGTTCGGGTGCTTTGTGATACCAGCCATCGTAATAACTGAGCGTGACCATCCACAGACCTAAAGCAAACATGCCGTAGACAATAATTGCGCTTAGCCAATGGAAGGCCGCAGat              >  2:3386289/1‑149 (MQ=255)
                                                                                                                                       ccccATCATTAAAAAAATACCGATGCTTTTATGCAGTTCGGGTGCTTTGTGATACCAGCCATCGTAATAACTGAGCGTGACCATCCACAGACCTAAAGCAAACATGCCGTAGACAATAATTGCGCTTAGCCAATGGAAGGCCGCAGata             >  1:239181/1‑149 (MQ=255)
                                                                                                                                          catcatTAAAAAAATACCGATGCTTTTATGCAGTTCGGGTGCTTTGTGATACCAGCCATCGTAATAACTGAGCGTGACCATCCACAGACCTAAAGCAAACATGCCGTAGACAATAATTGCGCTTAGCCAATGGAAGGCCGCAGATATaa          >  1:2554298/1‑149 (MQ=255)
                                                                                                                                          catcatTAAAAAAATACCGATGCTTTTATGCAGTTCGGGTGCTTTGTGATACCAGCCATCGTAATAACTGAGCGTGACCATCCACAGACCTAAAGCAAACATGCCGTAGACAATAATTGCGCTTAGCCAATGGAAGGCCGCAGATATaa          >  2:1328518/1‑149 (MQ=255)
                                                                                                                                           atcatTAAAAAAATACCGATGCTTTTATGCAGTTCGGGTGCTTTGTGATACCAGCCATCGTAATAACTGAGCGTGACCATCCACAGACCTAAAGCAAACATGCCGTAGACAATAATTGCGCTTAGCCCATGGAAGGCCGCAGATATAac         >  1:623741/1‑149 (MQ=255)
                                                                                                                                            tcatTAAAAAAATACCGATGCTTTTATGCAGTTCGGGTGCTTTGTGATACCAGCCATCGTAATAACTGAGCGTGACCATCCACAGACCTAAAGCAAACATGCCGTAGACAATAATTGCGCTTAGCCAATGGAAGGCCGCAGATATAaca        >  2:3267063/1‑149 (MQ=255)
                                                                                                                                            tcatTAAAAAAATACCGATGCTTTTATGCAGTTCGGGTGCTTTGTGATACCAGCCATCGTAATAACTGAGCGTGACCATCCACAGACCTAAAGCAAACATGCCGTAGACAATAATTGCGCTTAGCCAATGGAAGGCCGCAGATATAaca        >  2:776221/1‑149 (MQ=255)
                                                                                                                                             catTAAAAAAATACCGCTGCTTTTATGCAGTTCGGGTGCTTTGTGATACCAGCCATCGTAATAACTGAGCGTGACCATCCACAGACCTAAAGCAAACATGCCGTAGACAATAATTGCGCTTAGCCAATGGAAGGCCGCAGATATAacac       <  2:2764358/149‑1 (MQ=255)
                                                                                                                                             catTAAAAAAATACCGATGCTTTTATGCAGTTCGGGTGCTTTGTGATACCAGCCATCGTAATAACTGAGCGTGACCATCCACAGACCTAAAGCAAACATGCCGTAGACAATAATTGCGCTTAGCCAATGGAAGGCCGCAGATATAacac       >  1:4606394/1‑149 (MQ=255)
                                                                                                                                               ttAACAAAATACCGATGCTTTTATGCAGTTCGGGTGCTTTGTGATACCAGCCATCGTAATAACTGAGCGTGACCATCCACAGACCTAAAGCAAACATGCCGTAGACAATAATTGCGCTTAGCCAATGGAAGGCCGCAGATATAACACCa     >  1:952712/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
GCTGATGCCGATGGCGAAAAGCAATAGGTACAATGCGAGATGTCCAGCTCTTGCGGCAAGACGAGTCATGGGCGAATAACTCGGCAGCGGACCCGGCGGTGGAGATATGACACGCCACAGAACGCGAATAACCAGCCCCATCATTAACAAAATACCGATGCTTTTATGCAGTTCGGGTGCTTTGTGATACCAGCCATCGTAATAACTGAGCGTGACCATCCACAGACCTAAAGCAAACATGCCGTAGACAATAATTGCGCTTAGCCAATGGAAGGCCGCAGATATAACACCATAG  >  CP009273/1114182‑1114476

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: ‑