Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 768,068 G→T 5.3% P385P (CCG→CCT cydA → cytochrome d terminal oxidase, subunit I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273768,0680GT 5.3% 1104.6 / 23.2 291P385P (CCG→CCTcydAcytochrome d terminal oxidase, subunit I
Reads supporting (aligned to +/- strand):  ref base G (129/146);  new base T (10/6);  total (139/152)
Fisher's exact test for biased strand distribution p-value = 3.04e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 303 total aligned reads displayed.

ACCGACCAGGCGGTTCGTGACCAGTTCAATAGCATGAAGAAAGACCTCGGTTACGGTCTGCTGCTGAAACGCTATACGCCAAACGTGGCTGATGCGACTGAAGCGCAGATTCAACAGGCAACCAAAGACTCCATCCCGCGTGTAGCGCCGCTGTACTTTGCGTTCCGTATCATGGTGGCGTGTGGCTTCCTGCTTCTGGCAATCATCGCGCTCTCTTTCTGGAGTGTCATCCGCAACCGCATTGGCGAGAAAAAATGGCTTCTGCGCGCCGCGCTGTACGGTATTCCGCTGCCGTGGAT  >  CP009273/767919‑768217
                                                                                                                                                     |                                                                                                                                                     
aCCGACCAGGCGGTTCGTGACCAGTTCAATAGCATGAAGAAAGACCTCGGTTACGGTCTGCTGCTGAAACGCTATACGCCAAACGTGGCTGATGCGACTGAAGCGCAGATTCAACAGGCAACAAAAGACTCCATCCCGCGTGTAGcgccg                                                                                                                                                       >  2:941771/1‑150 (MQ=255)
      cAGGCGGTTCGTGACCAGTTCAATAGCATGAAGAAAGACCTCGGTTACGGTCTGCTGCTGAAACGCTATACGCCAAACGTGGCTGATGCGACTGAAGCGCAGATTCAACAGGCAACCAAAGACTCCATCCCGCGTGTAGCGCCGCTGTAc                                                                                                                                                 >  1:4399421/1‑150 (MQ=255)
       aGGCGGTTCGTGACCAGTTCAATAGCATGAAGAAAGACCTCGGTTACGGTCTGCTGCTGAAACGCTATACGCCAAACGTGGCTGATGCGACTGAAGCGCAGATTCAACAGGCAACCAAAGACTCCATCCCGCGTGTAGCGCCGCTGTACt                                                                                                                                                >  1:1114650/1‑150 (MQ=255)
       aGGCGGTTCGTGACCAGTTCAATAGCATGAAGAAAGACCTCGGTTACGGTCTGCTGCTGAAACGCTATACGCCAAACGTGGCTGATGCGACTGAAGCGCAGATTCAACAGGCAACCAAAGACTCCATCCCGCGTGTAGCGCCGCTGTACt                                                                                                                                                <  1:83535/150‑1 (MQ=255)
       aGGCGGTTCGTGACCAGTTCAATAGCATGAAGAAAGACCTCGGTTACGGTCTGCTGCTGAAACGCTATACGCCAAACGTGGCTGATGCGACTGAAGCGCAGATTCAACAGGCAACCAAAGACTCCATCCCGCGTGTAGCGCCGCTGTACt                                                                                                                                                <  1:2295751/150‑1 (MQ=255)
          cGGTTCGTGACCAGTTCAATAGCATGAAGAAAGACCTCGGTTACGGTCTGCTGCTGAAACGCTATACGCCAAACGTGGCTGATGCGACTGAAGCGCAGATTCAACAGGCAACCAAAGACTCCATCCCGCGTGTAGCGCCGCTGTACTTTg                                                                                                                                             <  2:3785645/150‑1 (MQ=255)
           ggTTCGTGACCAGTTCAATAGCATGAAGAAAGACCTCGGTTACGGTCTGCTGCTGAAACGCTATACGCCAAACGTGGCTGATGCGACTGAAGCGCAGATTCAACAGGCAACCAAAGACTCCATCCCGCGTGTAGCGCCGCTGTACTTTGc                                                                                                                                            >  1:102165/1‑150 (MQ=255)
            gTTCGTGACCAGTTCAATAGCATGAAGAAAGACCTCGGTTACGGTCTGCTGCTGAAACGCTATACGCCAAACGTGGCTGATGCGACTGAAGCGCAGATTCAACAGGCAACCAAAGACTCCATCCCGCGTGTAGCGCCGCTGTACTTTGCg                                                                                                                                           >  2:2604615/1‑150 (MQ=255)
             ttCGTGACCAGTTCAATAGCATGAAGAAAGACCTCGGTTACGGTCTGCTGCTGAAACGCTATACGCCAAACGTGGCTGATGCGACTGAAGCGCAGATTCAACAGGCAACCAAAGACTCCATCCCGCGTGTAGCGCCGCTGTACTTTGCGt                                                                                                                                          <  1:1177165/150‑1 (MQ=255)
             ttCGTGACCAGTTCAATAGCATGAAGAAAGACCTCGGTTACGGTCTGCTGCTGAAACGCTATACGCCAAACGTGGCTGATGCGACTGAAGCGCAGATTCAACAGGCAACCAAAGACTCCATCCCGCGTGTAGCGCCGCTGTACTTTGCGt                                                                                                                                          >  2:1170377/1‑150 (MQ=255)
               cGTGACCAGTTCAATAGCATGAAGAAAGACCTCGGTTACGGTCTGCTGCTGAAACGCTATACGCCAAACGTGGCTGATGCGACTGAAGCGCAGATTCAACAGGCAACCAAAGACTCCATCCCGCGTGTAGCGCCGCTGTACTTTGCGTTc                                                                                                                                        <  1:2763574/150‑1 (MQ=255)
               cGTGACCAGTTCAATAGCATGAAGAAAGACCTCGGTTACGGTCTGCTGCTGAAACGCTATACGCCAAACGTGGCTGATGCGACTGAAGCGCAGATTCAACAGGCAACCAAAGACTCCATCCCGCGTGTAGCGCCGCTGTACTTTGCGTTc                                                                                                                                        <  1:188970/150‑1 (MQ=255)
               cGTGACCAGTTCAATAGCATGAAGAAAGACCTCGGTTACGGTCTGCTGCTGAAACGCTATACGCCAAACGTGGCTGATGCGACTGAAGCGCAGATTCAACAGGCAACCAAAGACTCCATCCCGCGTGTAGCGCCGCTGTACTTTGCGTTc                                                                                                                                        >  2:1836044/1‑150 (MQ=255)
                gTGACCAGTTCAATAGCATGAAGAAAGACCTCGGTTACGGTCTGCTGCTGAAACGCTATACGCCAAACGTGGCTGATGCGACTGAAGCGCAGATTCAACAGGCAACCAAAGACTCCATCCCGCGTGTAGCGCCGCTGTACTTTGCGTTcc                                                                                                                                       <  1:2255714/150‑1 (MQ=255)
                gTGACCAGTTCAATAGCATGAAGAAAGACCTCGGTTACGGTCTGCTGCTGAAACGCTATACGCCAAACGTGGCTGATGCGACTGAAGCGCAGATTCAACAGGCAACCAAAGACTCCATCCCGCGTGTAGCGCCGCTGTACTTTGCGTTcc                                                                                                                                       >  1:2334936/1‑150 (MQ=255)
                gTGACCAGTTCAATAGCATGAAGAAAGACCTCGGTTACGGTCTGCTGCTGAAACGCTATACGCCAAACGTGGCTGATGCGACTGAAGCGCAGATACAACAGTCAACCAAAGACTCCATCCCGCGTGTAGCGCCGCTGTACTTTGCGTTcc                                                                                                                                       >  1:1912662/1‑150 (MQ=255)
                 tGACCAGTTCAATAGCATGAAGAAAGACCTCGGTTACGGTCTGCTGCTGAAACGCTATACGCCAAACGTGGCTGATGCGACTGAAGCGCAGATTCAACAGGCAACCAAAGACTCCATCCCGCGTGTAGCGCCGCTGTACTTTGCGTTCCg                                                                                                                                      >  1:4383948/1‑150 (MQ=255)
                 tGACCAGTTCAATAGCATGAAGAAAGACCTCGGTTACGGTCTGCTGCTGAAACGCTATACGCCAAACGTGGCTGATGCGACTGAAGCGCAGATTCAACAGGCAACCAAAGACTCCATCCCGCGTGTAGCGCCGCTGTACTTTGCGTTCCg                                                                                                                                      >  1:4382479/1‑150 (MQ=255)
                 tGACCAGTTCAATAGCATGAAGAAAGACCTCGGTTACGGTCTGCTGCTGAAACGCTATACGCCAAACGTGGCTGATGCGACTGAAGCGCAGATTCAACAGGCAACCAAAGACTCCATCCCGCGTGTAGCGCCGCTGTACTTTGCGTTCCg                                                                                                                                      <  1:3972575/150‑1 (MQ=255)
                   aCCAGTTCAATAGCATGAAGAAAGACCTCGGTTACGGTCTGCTGCTGAAACGCTATACGCCAAACGTGGCTGATGCGACTGAAGCGCAGATTCAACAGGCAACCAAAGACTCCATCCCGCGTGTAGCGCCGCTGTACTTTGCGTTCCGTa                                                                                                                                    >  1:4339953/1‑150 (MQ=255)
                       gTTCAATAGCATGAAGAAAGACCTCGGTTACGGTCTGCTGCTGAAACGCTATACGCCAAACGTGGCTGATGCGACTGAAGCGCAGATTCAACAGGCAACCAAAGACTCCATCCCGCGTGTAGCGCCGCTGTACTTTGCGTTCCGTATCAt                                                                                                                                <  2:1032327/150‑1 (MQ=255)
                        ttCAATAGCATGAAGAAAGACCTCGGTTACGGTCTGCTGCTGAAACGCTATACGCCAAACGTGGCTGATGCGACTGAAGCGCAGATTCAACAGGCAACCAAAGACTCCATCCCGCGTGTAGCGCCGCTGTACTTTGCGTTCCGTATCAtg                                                                                                                               >  1:3916416/1‑150 (MQ=255)
                           aaTAGCATGAAGAAAGACCTCGGTTACGGTCTGCTGCTGAAACGCTATACGCCAAACGTGGCTGATGCGACTGAAGCGCAGATTCAACAGGCAACCAAAGACTCCATCCCGCGTGTAGCGCCGCTGTACTTTGCGTTCCGTATCAtggtg                                                                                                                            >  1:987315/1‑150 (MQ=255)
                            aTAGCATGAAGAAAGACCTCGGTTACGGTCTGCTGCTGAAACGCTATACGCCAAACGTGGCTGATGCGACTGAAGCGCAGATTCAACAGGCAACCAAAGACTCCATCCCGCGTGTAGCGCCGCTGTACTTTGCGTTCCGTATCAtggtgg                                                                                                                           >  2:3238417/1‑150 (MQ=255)
                                  tGAAGAAAGACCTCGGTTACGGTCTGCTGCTGAAACGCTATACGCCAAACGTGGCTGATGCGACTGAAGCGCAGATTCAACAGGCAACCAAAGACTCCATCCCGCGTGTAGCGCCGCTGTACTTTGCGTTCCGTATCATGGTGGCGtgtg                                                                                                                     <  1:908415/150‑1 (MQ=255)
                                    aagaaagaCCTCGGTTACGGTCTGCTGCTGAAACGCTATACGCCAAACGTGGCTGATGCGACTGAAGCGCAGATTCAACAGGCAACCAAAGACTCCATCCCGCGTGTAGCGCCTCTGTACTTTGCGTTCCGTATCATGGTGGCGTGTGGc                                                                                                                   >  1:2345899/1‑150 (MQ=255)
                                      gaaagaCCTCGGTTACGGTCTGCTGCTGAAACGCTATACGCCAAACGTGGCTGATGCGACTGAAGCGCAGATTCAACAGGCAACCAAAGACTCCATCCCGCGTGTAGCGCCGCTGTACTTTGCGTTCCGTATCATGGTGGCGTGTGGCtt                                                                                                                 <  2:756465/150‑1 (MQ=255)
                                      gaaagaCCTCGGTTACGGTCTGCTGCTGAAACGCTATACGCCAAACGTGGCTGATGCGACTGAAGCGCAGATTCAACAGGCAACCAAAGACTCCATCCCGCGTGTAGCGCCGCTGTACTTTGCGTTCCGTATCATGGTGGCGTGTGGCtt                                                                                                                 <  2:4096641/150‑1 (MQ=255)
                                         agaCCTCGGTTACGGTCTGCTGCTGAAACGCTATACGCCAAACGTGGCTGATGCGACTGAAGCGCAGATTCAACAGGCAACCAAAGACTCCATCCCGCGTGTAGCGCCGCTGTACTTTGCGTTCCGTATCATGGTGGCGTGTGGCTTCCt                                                                                                              <  2:1444600/150‑1 (MQ=255)
                                            ccTCGGTTACGGTCTGCTGCTGAAACGCTATACGCCAAACGTGGCTGATGCGACTGAAGCGCAGATTCAACAGGCAACCAAAGACTCCATCCCGCGTGTAGCGCCGCTGTACTTTGCGTTCCGTATCATGGTGGCGTGTGGCTTCCTGCt                                                                                                           >  1:3601984/1‑150 (MQ=255)
                                            ccTCGGTTACGGTCTGCTGCTGAAACGCTATACGCCAAACGTGGCTGATGCGACTCAAGCGCAGATTCAACAGGCAACCAAAGACTCCATCCCGCGTGTAGCGCCGCTGTACTTTGCGTTCCGTATCATGGTGGCGTGTGGCTTCCTGCt                                                                                                           >  1:2939009/1‑150 (MQ=255)
                                             cTCGGTTACGGTCTGCTGCTGAAACGCTATACGCCAAACGTGGCTGATGCGACTGAAGCGCAGATTCAACAGGCAACCAAAGACTCCATCCCGCGTGTAGCGCCGCTGTACTTTGCGTTCCGTATCATGGTGGCGTGTGGCTTCCTGCtt                                                                                                          <  2:707169/150‑1 (MQ=255)
                                             cTCGGTTACGGTCTGCTGCTGAAACGCTATACGCCAAACGTGGCTGATGCGACTGAAGCGCAGATTCAACAGGCAACCAAAGACTCCATCCCGCGTGTAGCGCCGCTGTACTTTGCGTTCCGTATCATGGTGGCGTGTGGCTTCCTGCtt                                                                                                          <  2:1035687/150‑1 (MQ=255)
                                                 gTTACGGTCTGCTGCTGAAACGCTATACGCCAAACGTGGCTGATGCGACTGAAGCGCAGATTCAACAGGCAACCAAAGACTCCATCCCGCGTGTAGCGCCGCTGTACTTTGCGTTCCGTATCATGGTGGCGTGTGGCTTCCTGCTTCTgg                                                                                                      >  1:2850058/1‑150 (MQ=255)
                                                  ttACGGTCTGCTGCTGAAACGCTATACGCCAAACGTGGCTGATGCGACTGAAGCGCAGATTCAACAGGCAACCAAAGACTCCATCCCGCGTGTAGCGCCGCTGTACTTTGCGTTCCGTATCATGGTGGCGTGTGGCTTCCTGCTTCTGGc                                                                                                     >  1:1378723/1‑150 (MQ=255)
                                                   tACGGTCTGCTGCTGAAACGCTATACGCCAAACGTGGCTGATGCGACTGAAGCGCAGATTCAACAGGCAACCAAAGACTCCATCCCGCGTGTAGCGCCGCTGTACTTTGCGTTCCGTATCATGGTGGCGTGTGGCTTCCTGCTTCTGGCa                                                                                                    >  2:3107657/1‑150 (MQ=255)
                                                   tACGGTCTGCTGCTGAAACGCTATACGCCAAACGTGGCTGATGCGACTGAAGCGCAGATTCAACAGGCAACCAAAGACTCCATCCCGCGTGTAGCGCCGCTGTACTTTGCGTTCCGTATCATGGTGGCGTGTGGCTTCCTGCTTCTGGCa                                                                                                    >  1:2128106/1‑150 (MQ=255)
                                                    aCGGTCTGCTGCTGAAACGCTATACGCCAAACGTGGCTGTTGCGACTGAAGCGCAGATTCAACAGGCAACCAAAGACTCCATCCCGCGTGTAGCGCCGCTGTACTTTGCGTTCCGTATCATGGTGGCGTGTGGCTTCCTGCTTCTGGCaa                                                                                                   <  1:3674370/150‑1 (MQ=255)
                                                    aCGGTCTGCTGCTGAAACGCTATACGCCAAACGTGGCTGATGCGACTGAAGCGCAGATTCAACAGGCAACCAAAGACTCCATCCCGCGTGTAGCGCCGCTGTACTTTGCGTTCCGTATCATGGTGGCGTGTGGCTTCCTGCTTCTGGCaa                                                                                                   <  2:1127979/150‑1 (MQ=255)
                                                      ggTCTGCTGCTGAAACGCTATACGCCAAACGTGGCTGATGCGACTGAAGCGCAGATTCAACAGGCAACCAAAGACTCCATCCCGCGTGTAGCGCCGCTGTACTTTGCGTTCCGTATCATGGTGGCGTGTGGCTTCCTGCTTCTGGCAatc                                                                                                 >  1:1431292/1‑150 (MQ=255)
                                                       gTCTGCTGCTGAAACGCTATACGCCAAACGTGGCTGATGCGACTGAAGCGCAGATTCAACAGGCAACCAAAGACTCCATCCCGCGTGTAGCGCCGCTGTACTTTGCGTTCCGTATCATGGTGGCGTGTGGCTTCCTGCTTCTGGCAatca                                                                                                <  1:3248920/150‑1 (MQ=255)
                                                             tgctgAAACGCTATACGCCAAACGTGGCTGATGCGACTGAAGCGCAGATTCAACAGGCAACCAAAGACTCCATCCCGCGTGTAGCGCCGCTGTACTTTGCGTTCCGTATCATGGTGGCGTGTGGCTTCCTGCTTCTGGCAATCATCgcgc                                                                                          >  1:3335432/1‑150 (MQ=255)
                                                                tgAAACGCTATACGCCAAACGTGGCTGATGCGACTGAAGCGCAGATTCAACAGGCAACCAAAGACTCCATCCCGCGTGTAGCGCCGCTGTACTTTGCGTTCCGTATCATGGTGGCGTGTGGCTTCCTGCTTCTGGCAATCATCGCGctct                                                                                       <  2:233279/150‑1 (MQ=255)
                                                                 gAAACGCTATACGCCAAACGTGGCTGATGCGACTGAAGCGCAGATTCAACAGGCAACCAAAGACTCCATCCCGCGTGTAGCGCCGCTGTACTTTGCGTTCCGTATCATGGTGGCGTGTGGCTTCCTGCTTCTGGCAATCATCGCGctctc                                                                                      <  2:62953/150‑1 (MQ=255)
                                                                   aaCGCTATACGCCAAACGTGGCTGATGCGACTGAAGCGCAGATTCAACAGGCAACCAAAGACTCCATCCCGCGTGTAGCGCCGCTGTACTTTGCGTTCCGTATCATGGTGGCGTGTGGCTTCCTGCTTCTGGCAATCATCGCGCTCTCtt                                                                                    <  2:2567341/150‑1 (MQ=255)
                                                                       cTATACGCCAAACGTGGCTGATGCGACTGAAGCGCAGATTCAACAGGCAACCAAAGACTCCATCCCGCGTGTAGCGCCGCTGTACTTTGCGTTCCGTATCATGGTGGCGTGTGGCTTCCTGCTTCTGGCAATCATCGCGCTCTCTTTCTg                                                                                <  1:2428798/150‑1 (MQ=255)
                                                                       cTATACGCCAAACGTGGCTGATGCGACTGAAGCGCAGATTCAACAGGCAACCAAAGACTCCATCCCGCGTGTAGCGCCGCTGTACTTTGCGTTCCGTATCATGGTGGCGTGTGGCTTCCTGCTTCTGGCAATCATCGCGCTCTCTTTCTg                                                                                >  1:1737136/1‑150 (MQ=255)
                                                                       cTATACGCCAAACGTGGCTGATGCGACTGAAGCGCAGATTCAACAGGCAACCAAAGACTCCATCCCGCGTGTAGCGCCGCTGTACTTTGCGTTCCGTATCATGGTGGCGTGTGGCTTCCTGCTTCTGGCAATCATCGCGCTCTCTTTCTg                                                                                <  1:2948214/150‑1 (MQ=255)
                                                                       cTATACGCCAAACGTGGCTGATGCGACTGAAGCGCAGATTCAACAGGCAACCAAAGACTCCATCCCGCGTGTAGCGCCGCTGTACTTTGCGTTCCGTATCATGGTGGCGTGTGGCTTCCTGCTTCTGGCAATCATCGCGCTCTCTTTCTg                                                                                <  1:2843014/150‑1 (MQ=255)
                                                                        tataCGCCAAACGTGGCTGATGCGACTGAAGCGCAGATTCAACAGGCAACCAAAGACTCCATCCCGCGTGTAGCGCCTCTGTACTTTGCGTTCCGTATCATGGTGGCGTGTGGCTTCCTGCTTCTGGCAATCATCGCGCTCTCTTTCTgg                                                                               >  2:1375994/1‑150 (MQ=255)
                                                                             gCCAAACGTGGCTGATGCGACTGAAGCGCAGATTCAACAGGCAACCAAAGACTCCATCCCGCGTGTAGCGCCGCTGTACTTTGCGTTCCGTATCATGGTGGCGTGTGGCTTCCTGCTTCTGGCAATCATCGCGCTCTCTTTCTGGAgtgt                                                                          <  2:166110/150‑1 (MQ=255)
                                                                              ccAAACGTGGCTGATGCGACTGAAGCGCAGATTCAACAGGCAACCAAAGACTCCATCCCGCGTGTAGCGCCTCTGTACTTTGCGTTCCGTATCATGGTGGCGTGTGGCTTCCTGCTTCTGGCAATCATCGCGCTCTCTTTCTGGAGTGTc                                                                         <  1:2583382/150‑1 (MQ=255)
                                                                              ccAAACGTGGCTGATGCGACTGAAGCGCAGATTCAACAGGCAACCAAAGACTCCATCCCGCGTGTAGCGCCGCTGTACTTTGCGTTCCGTATCATGGTGGCGTGTGGCTTCCTGCTTCTGGCAATCATCGCGCTCTCTTTCTGGAGTGTc                                                                         >  1:1107936/1‑150 (MQ=255)
                                                                                aaaCGTGGCTGATGCGACTGAAGCGCAGATTCAACAGGCAACCAAAGACTCCATCCCGCGTGTAGCGCCGCTGTACTTTGCGTTCCGTATCATGGTGGCGTGTGGCTTCCTGCTTCTGGCAATCATCGCGCTCTCTTTCTGGAGTGTCAt                                                                       >  1:3155357/1‑150 (MQ=255)
                                                                                 aaCGTGGCTGATGCGACTGAAGCGCAGATTCAACAGGCAACCAAAGACTCCATCCCGCGTGTAGCGCCGCTGTACTTTGCGTTCCGTATCATGGTGGCGTGTGGCTTCCTGCTTCTGGCAATCATCGCGCTCTCTTTCTGGAGTGTCATc                                                                      >  2:2965891/1‑150 (MQ=255)
                                                                                     tGGCTGATGCGACTGAAGCGCAGATTCAACAGGCAACCAAAGACTCCATCCCGCGTGTAGCGCCGCTGTACTTTGCGTTCCGTATCATGGTGGCGTGTGGCTTCCTGCTTCTGGCAATCATCGCGCTCTCTTTCTGGAGTGTCATCCGCa                                                                  >  1:1467422/1‑150 (MQ=255)
                                                                                     tGGCTGATGCGACTGAAGCGCAGATTCAACAGGCAACCAAAGACTCCATCCCGCGTGTAGCGCCGCTGTACTTTGCGTTCCGTATCATGGTGGCGTGTGGCTTCCTGCTTCTGGCAATCATCGCGCTCTCTTTCTGGAGTGTCATCCGCa                                                                  <  2:458933/150‑1 (MQ=255)
                                                                                     tGGCTGATGCGACTGAAGCGCAGATTCAACAGGCAACCAAAGACTCCATCCCGCGTGTAGCGCCGCTGTACTTTGCGTTCCGTATCATGGTGGCGTGTGGCTTCCTGCTTCTGGCAATCATCGCGCTCTCTTTCTGGAGTGTCATCCGCa                                                                  <  2:2747130/150‑1 (MQ=255)
                                                                                        cTGATGCGACTGAAGCGCAGATTCAACAGGCAACCAAAGACTCCATCCCGCGTGTAGCGCCGCTGTACTTTGCGTTCCGTATCATGGTGGCGTGTGGCTTCCTGCTTCTGGCAATCATCGCGCTCTCTTTCTGGAGTGTCATCCGCAAcc                                                               <  2:3630519/150‑1 (MQ=255)
                                                                                           aTGCGACTGAAGCGCAGATTCAACAGGCAACCAAAGACTCCATCCCGCGTGTAGCGCCGCTGTACTTTGCGTTCCGTATCATGGTGGCGTGTGGCTTCCTGCTTCTGGCAATCATCGCGCTCTCTTGCTGGAGTGTCATCCGCAACCGCa                                                            <  1:4456373/150‑1 (MQ=255)
                                                                                              cGACTGAAGCGCAGATTCAACAGGCAACCAAAGACTCCATCCCGCGTGTAGCGCCGCTGTACTTTGCGTTCCGTATCATGGTGGCGTGTGGCTTCCTGCTTCTGGCAATCATCGCGCTCTCTTTCTGGAGTGTCATCCGCAACCGCATTg                                                         >  2:1467634/1‑150 (MQ=255)
                                                                                                aCTGAAGCGCAGATTCAACAGGCAACCAAAGACTCCATCCCGCGTGTAGCGCCGCTGTACTTTGCGTTCCGTATCATGGTGGCGTGTGGCTTCCTGCTTCTGGCAATCATCGCGCTCTCTTTCTGGAGTGTCATCCGCAACCGCATTGGc                                                       >  2:647945/1‑150 (MQ=255)
                                                                                                aCTGAAGCGCAGATTCAACAGGCAACCAAAGACTCCATCCCGCGTGTAGCGCCGCTGTACTTTGCGTTCCGTATCATGGTGGCGTGTGGCTTCCTGCTTCTGGCAATCATCGCGCTCTCTTTCTGGAGTGTCATCCGCAACCGCATTGGc                                                       <  2:847294/150‑1 (MQ=255)
                                                                                                aCTGAAGCGCAGATTCAACAGGCAACCAAAGACTCCATCCCGCGTGTAGCGCCGCTGTACTTTGCGTTCCGTATCATGGTGGCGTGTGGCTTCCTGCTTCTGGCAATCATCGCGCTCTCTTTCTGGAGTGTCATCCGCAACCGCATTGGc                                                       <  1:61828/150‑1 (MQ=255)
                                                                                                 cTGAAGCGCAGATTCAACAGGCAACCAAAGACTCCATCCCGCGTGTAGCGCCGCTGTACTTTGCGTTCCGTATCATGGTGGCGTGTGGCTTCCTGCTTCTGGCAATCATCGCGCTCTCTTTCTGGAGTGTCATCCGCAACCGCATTGGCg                                                      >  1:4262814/1‑150 (MQ=255)
                                                                                                   gAAGCGCAGATTCAACAGGCAACCAAAGACTCCATCCCGCGTGTAGCGCCGCTGTACTGTGCGTTCCGTATCATGGTGGCGTGTGGCTTCCTGCTTCTGGCAATCATCGCGCTCTCTTTCTGGAGTGTCATCCGCAACCGCATTGGCgag                                                    <  2:3624555/150‑1 (MQ=255)
                                                                                                    aaGCGCAGATTCAACAGGCAACCAAAGACTCCATCCCGCGTGTAGCGCCGCTGTACTTTGCGTTCCGTATCATGGTGGCGTGTGGCTTCCTGCTTCTGGCAATCATCGCGCTCTCTTTCTGGAGTGTCATCCGCAACCGCATTGGCgaga                                                   >  2:3449026/1‑150 (MQ=255)
                                                                                                    aaGCGCAGATTCAACAGGCAACCAAAGACTCCATCCCGCGTGTAGCGCCGCTGTACTTTGCGTTCCGTATCATGGTGGCGTGTGGCTTCCTGCTTCTGGCAATCATCGCGCTCTCTTTCTGGAGTGTCATCCGCAACCGCATTGGCgaga                                                   >  2:3449081/1‑150 (MQ=255)
                                                                                                    aaGCGCAGATTCAACAGGCAACCAAAGACTCCATCCCGCGTGTAGCGCCGCAGTACTTTGCGTTCCGTATCATGGTGGCGTGTGGCTTCCTGCTTCTGGCAATCAACGCGCTCTCTTTCTGGAGTGTCATCCGCAACCGCATTGGCgaga                                                   >  2:1943801/1‑150 (MQ=255)
                                                                                                     aGCGCAGATTCAACAGGCACCCAAAGACTCCATCCCGCGTGTAGCGCCGCTGTACTTTGCGTTCCGTATCATGGTGGCGTGTGGCTTCCTGCTTCTGGCAATCATCGCGCTCTCTTTCTGGAGTGTCATCCGCAACCGCATTGGCGAGaa                                                  >  2:3979098/1‑150 (MQ=255)
                                                                                                        gcAGATTCAACAGGCAACCAAAGACTCCATCCCGCGTGTAGCGCCGCTGTACTTTGCGTTCCGTATCATGGTGGCGTGTGGCTTCCTGCTTCTGGCAATCATCGCGCTCTCTTTCTGGAGTGTCATCCGCAACCGCATTGGCGAGaaaaa                                               <  1:574880/150‑1 (MQ=255)
                                                                                                           gATTCAACAGGCAACCAAAGACTCCATCCCGCGTGTAGCGCCGCTGTACTTTGCGTTCCGTATCATGGTGGCGTGTGGCTTCCTGCTTCTGGCAATCATCGCGCTCTCTTTCTGGAGTGTCATCCGCAACCGCATTGGCGAGAAAAAATg                                            <  2:212313/150‑1 (MQ=255)
                                                                                                           gATTCAACAGGCAACCAAAGACTCCATCCCGCGTGTAGCGCCGCTGTACTTTGCGTTCCGTATCATGGTGGCGTGTGGCTTCCTGCTTCTGGCAATCATCGCGCTCTCTTTCTGGAGTGTCATCCGCAACCGCATTGGCGAGAAAAAATg                                            <  2:212072/150‑1 (MQ=255)
                                                                                                            aTTCAACAGGCAACCAAAGACTCCATCCCGCGTGTAGCGCCGCTGTACTTTGCGTTCCGTATCATGGTGGCGTGTGGCTTCCTGCTTCTGGCAATCATCGCGCTCTCTTTCTGGAGTGTCATCCGCAACCGCATTGGCGAGAAAAAATgg                                           <  2:4267457/150‑1 (MQ=255)
                                                                                                             ttCAACAGGCAACCAAAGACTCCATCCCGCGTGTAGCGCCGCTGTACTTTGCGTTCCGTATCATGGTGGCGTGTGGCTTCCTGCTTCTGGCAATCATCGCGCTCTCTTTCTGGAGTGTCATCCGCAACCGCATTGGCGAGAAAAAATGGc                                          >  2:2483339/1‑150 (MQ=255)
                                                                                                              tCAACAGGCAACCAAAGACTCCATCCCGCGTGTAGCGCCGCTGTACTTTGCGTTCCGTATCATGGTGGCGTGTGGCTTCCTGCTTCTGGCAATCATCGCGCTCTCTTTCTGGAGTGTCATCCGCAACCGCATTGGCGAGAAAAAATGGCt                                         <  2:1902974/150‑1 (MQ=255)
                                                                                                                aaCAGGCAACCAAAGACTCCATCCCGCGTGTAGCGCCGCTGTACTTTGCGTTCCGTATCATGGTGGCGTGTGGCTTCCTGCTTCTGGCAATCATCGCGCTCTCTTTCTGGAGTGTCATCCGCAACCGCATTGGCGAGAAAAAATGGCTTc                                       >  2:3767409/1‑150 (MQ=255)
                                                                                                                aaCAGGCAACCAAAGACTCCATCCCGCGTGTAGCGCCGCTGTACTTTGCGTTCCGTATCATGGTGGCGTGTGGCTTCCTGCTTCTGGCAATCATCGCGCTCTCTTTCTGGAGTGTCATCCGCAACCGCATTGGCGAGAAAAAATGGCTTc                                       <  2:3380145/150‑1 (MQ=255)
                                                                                                                 aCAGGCAACCAAAGACTCCATCCCGCGTGTAGCGCCGCTGTACTTTGCGTTCCGTATCATGGTGGCGTGTGGCTTCCTGCTTCTGGCCATCATCGCGCTCTCTTTCTGGAGTGTCATCCGCAACCGCATTGGCGAGAAAAAATGGCTTCt                                      >  2:1781751/1‑150 (MQ=255)
                                                                                                                 aCAGGCAACCAAAGACTCCATCCCGCGTGTAGCGCCGCTGTACTTTGCGTTCCGTATCATGGTGGCGTGTGGCTTCCTGCTTCTGGCAATCATCGCGCTCTCTTTCTGGAGTGTCATCCGCAACCGCATTGGCGAGAAAAAATGGCTTCt                                      <  1:3418059/150‑1 (MQ=255)
                                                                                                                 aCAGGCAACCAAAGACTCCATCCCGCGTGTAGCGCCGCTGTACTTTGCGTTCCGTATCATGGTGGCGTGTGGCTTCCTGCTTCTGGCAATCATCGCGCTCTCTTTCTGGAGTGTCATCCGCAACCGCATTGGCGAGAAAAAATGGCTTCt                                      <  1:2093097/150‑1 (MQ=255)
                                                                                                                 aCAGGCAACCAAAGACTCCATCCCGCGTGTAGCGCCGCTGTACTTTGCGTTCCGTATCATGGTGGCGTGTGGCTTCCTGCTTCTGGCAATCATCGCGCTCTCTTTCTGGAGTGTCATCCGCAACCGCATTGGCGAGAAAAAATGGCTTCt                                      <  2:2193006/150‑1 (MQ=255)
                                                                                                                    ggCAACCAAAGACTCCATCCCGCGTGTAGCGCCGCTGTACTAAGCGTTCCGTATCATGGTGTCGTGTGGCTTCCTGCTTCTGGCAATCATCGCGCACACATTCTGGAGTGTCATCCGCAACCGCATAGGAGAGACAAAATGGCTTCTgcg                                   >  2:645669/1‑150 (MQ=255)
                                                                                                                      cAACCAAAGACTCCATCCCGCGTGTAGCGCCGCTGTACTTTGCGTTCCGTATCATGGTGGCGTGTGGCTTCCTGCTTCTGGCAATCATCGCGCTCTCTTTCTGGAGTGTCATCCGCAACCGCATTGGCGAGAAAAAATGGCTTCTgcgcg                                 >  2:1198708/1‑150 (MQ=255)
                                                                                                                        aCCAAAGACTCCATCCCGCGTGTAGCGCCGCTGTACTTTGCGTTCCGTATCATGGTGGCGTGTGGCTTCCTGCTTCTGGCAATCATCGCGCTCTCTTTCTGGAGTGTCATCCGCAACCGCATTGGCGAGAAAAAATGGCTTCTGcgcgcc                               <  1:4050366/150‑1 (MQ=255)
                                                                                                                             aGACTCCATCCCGCGTGTAGCGCCGCTGTACTTTGCGTTCCGTATCATGGTGGCGTGTGGCTTCCTGCTTCTGGCAATCATCGCGCTCTCTTTCTGGAGTGTCATCCGCAACCGCATTGGCGAGAAAAAATGGCTTCTGCGCGCCGCGCt                          >  2:13072/1‑150 (MQ=255)
                                                                                                                                  ccATCCCGCGTGTAGCGCCGCTGTACTTTGCGTTCCGTATCATGGTGGCGTGTGGCTTCCTGCTTCTGGCAATCATCGCGCTCTCTTTCTGGAGTGTCATCCGCAACCGCATTGGCGAGAAAAAATGGCTTCTGCGCGCCGCGCTGTACg                     >  2:3683785/1‑150 (MQ=255)
                                                                                                                                    aTCCCGCGTGTAGCGCCGCTGTACTTTGCGTTCCGTATCATGGTGGCGTGTGGCTTCCTGCTTCTGGCAATCATCGCGCTCTCTTTCTGGAGTGTCATCCGCAACCGCATTGGCGAGAAAAAATGGCTTCTGCGCGCCGCGCTGTACGGt                   <  1:4129220/150‑1 (MQ=255)
                                                                                                                                      cccGCGTGTAGCGCCGCTGTACTTTGCGTTCCGTATCATGGTGGCGTGTGGCTTCCTGCTTCTGGCAATCATCGCGCTCTCTTTCTGGAGTGTCATCCGCAACCGCATTGGCGAGAAAAAATGGCTTCTGCGCGCCGCGCTGTACGGTAt                 <  1:137085/150‑1 (MQ=255)
                                                                                                                                      cccGCGTGTAGCGCCGCTGTACTTTGCGTTCCGTATCATGGTGGCGTGTGGCTTCCTGCTTCTGGCAATCATCGCGCTCTCTTTCTGGAGTGTCATCCGCAACCGCATTGGCGAGAAAAAATGGCTTCTGCGCGCCGCGCTGTACGGTAt                 <  1:2874416/150‑1 (MQ=255)
                                                                                                                                      cccGCGTGTAGCGCCGCTGTACTTTGCGTTCCGTATCATGGTGGCGTGTGGCTTCCTGCTTCTGGCAATCATCGCGCTCTCTTTCTGGAGTGTCATCCGCAACCGCATTGGCGAGAAAAAATGGCTTCTGCGCGCCGCGCTGTACGGTAt                 <  1:1571613/150‑1 (MQ=255)
                                                                                                                                       ccGCGTGTAGCGCCGCTGTACTTTGCGTTCCGTATCATGGTGGCGTGTGGCTTCCTGCTTCTGGCAATCATCGCGCTCTCTTTCTGGAGTGTCATCCGCAACCGCATTGGCGAGAAAAAATGGCTTCTGCGCGCCGCGCTGTACGGTAtt                >  2:3738881/1‑150 (MQ=255)
                                                                                                                                       ccGCGTGTAGCGCCGCTGTACTTTGCGTTCCGTATCATGGTGGCGTGTGGCTTCCTGCTTCTGGCAATCATCGCGCTCTCTTTCTGGAGTGTCATCCGCAACCGCATTGGCGAGAAAAAATGGCTTCTGCGCGCCGCGCTGTACGGTAtt                <  1:531208/150‑1 (MQ=255)
                                                                                                                                         gcgTGTAGCGCCGCTGTACTTTGCGTTCCGTATCATGGTGGCGTGTGGCTTCCTGCTTCTGGCAATCATCGCGCTCTCTTTCTGGAGTGTCATCCGCAACCGCATTGGCGAGAAAAAATGGCTTCTGCGCGCCGCGCTGTACGGTATTcc              >  2:1599756/1‑150 (MQ=255)
                                                                                                                                         gcgTGTAGCGCCGCGGTACTTTGCGTTCCGTATCATGGTGGCGTGTGGCTTCCTGCTTCTGGCAATCATCGCGCTCTCTTTCTGGAGTGTCATCCGCAACCGCATTGGCGAGAAAAAATGGCTTCTGCGCGCCGCGCTGTACGGTATTcc              <  2:4272826/150‑1 (MQ=255)
                                                                                                                                            tgtAGCGCCGCTGTACTTTGCGTTCCGTATCATGGTGGCGTGTGGCTTCCTGCTTCTGGCAATCATCGCGCTCTCTTTCTGGAGTGTCATCCGCAACCGCATTGGCGAGAAAAAATGGCTTCTGCGCGCCGCGCTGTACGGTATTCCGCt           >  2:1609105/1‑150 (MQ=255)
                                                                                                                                               aGCGCCGCTGTACTTTGCGTTCCGTATCATGGTGGCGTGTGGCTTCCTGCTTCTGGCAATCATCGCGCTCTCTTTCTGGAGTGTCATCCGCAACCGCATTGGCGAGAAAAAATGGCTTCTGCGCGCCGCGCTGTACGGTATTCCGCTGcc        <  2:821902/150‑1 (MQ=255)
                                                                                                                                                 cgccgcTGTACTTTGCGTTCCGTATCATGGTGGCGTGTGGCTTCCTGCTTCTGGCAATCATCGCGCTCTCTTTCTGGAGTGTCATCCGCAACCGCATTGGCGAGAAAAAATGGCTTCTGCGCGCCGCGCTGTACGGTATTCCGCTGCCGt      >  2:4540026/1‑150 (MQ=255)
                                                                                                                                                  gcctcTGTACTTTGCGTTCCGTATCATGGTGGCGTGTGGCTTCCTGCTTCTGGCAATCATCGCGCTCTCTTTCTGGAGTGTCATCCGCAACCGCATTGGCGAGAAAAAATGGCTTCTGCGCGCCGCGCTGTACGGTATTCCGCTGCCGTg     >  1:2785270/1‑150 (MQ=255)
                                                                                                                                                     gcTGTACTTTGCGTTCCGTATCATGGTGGCGTGTGGCTTCCTGCTTCTGGCAATCATCGCGCTCTCTTTCTGGAGTGTCATCCGCAACCGCATTGGCGAGAAAAAATGGCTTCTGCGCGCCGCGCTGTACGGTATTCCGCTGCCGTGGAt  >  1:1169154/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
ACCGACCAGGCGGTTCGTGACCAGTTCAATAGCATGAAGAAAGACCTCGGTTACGGTCTGCTGCTGAAACGCTATACGCCAAACGTGGCTGATGCGACTGAAGCGCAGATTCAACAGGCAACCAAAGACTCCATCCCGCGTGTAGCGCCGCTGTACTTTGCGTTCCGTATCATGGTGGCGTGTGGCTTCCTGCTTCTGGCAATCATCGCGCTCTCTTTCTGGAGTGTCATCCGCAACCGCATTGGCGAGAAAAAATGGCTTCTGCGCGCCGCGCTGTACGGTATTCCGCTGCCGTGGAT  >  CP009273/767919‑768217

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: