Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,953,876 C→A 100% P76P (CCG→CCT yecM ← putative metal‑binding enzyme

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,953,8760CA97.8% 691.5 / 6.1 225P76P (CCG→CCTyecMputative metal‑binding enzyme
Reads supporting (aligned to +/- strand):  ref base C (3/2);  new base A (112/108);  total (115/110)
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.43e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 225 total aligned reads displayed.

CGTTTCTGGATCACCCGGCAGAACGATTTCAATATGTTCCCAACCTTCGTGTGGGTAACGTTTTTCCCCGGGCCACGGTAGCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATCGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGCCAGCGTTCAGCAGTGGCGTTTTGATGGCAACGCAAAGAAATGTGATCGGCAGTCAGCGGAGTGATATTCAACCCCAGACGGCGG  >  CP009273/1953732‑1954020
                                                                                                                                                |                                                                                                                                                
cGTTTCTGGATCACCCGGCAGAACGATTTCAATATGTTCCCAACCTTCGTGTGGGTAACGTTTTTCCCCGGGCCACGGTAGCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGt                                                                                                                                               >  2:2302830/1‑148 (MQ=255)
cGTTTCTGGATCACCCGGCAGAACGATTTCAATATGTTCCCAACCTTCGTGAGGGTAACGTTTTTCCCCGGGCCACGGTAGCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGt                                                                                                                                               <  2:3513782/148‑1 (MQ=255)
    tCTGGATCACCCGGCAGAACGATTTCAATATGTTCCCAACCTTCGTGTGGGTAACGTTTTTCCCCGGGCCACGGTAGCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCTGc                                                                                                                                           <  1:368825/148‑1 (MQ=255)
       ggATCACCCGGCAGAACGATTTCAATATGTTCCCAACCTTCGTGTGGGTAACGTTTTTCCCCGGGCCACGGTAGCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCTGCCAt                                                                                                                                        <  1:471226/148‑1 (MQ=255)
           cACCCGGCAGAACGATTTCAATATGTTCCCAACCTTCGTGTGGGTAACGTTTTTCCCCGGGCCACGGTAGCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCTGCCATTGAt                                                                                                                                    <  1:1399813/148‑1 (MQ=255)
           cACCCGGCAGAACGATTTCAATATGTTCCCAACCTTCGTGTGGGTAACGTTTTTCCCCGGGCCACGGTAGCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCTGCCATTGAt                                                                                                                                    >  2:2943292/1‑148 (MQ=255)
           cACCCGGCAGAACGATTTCAATATGTTCCCAACCTTCGTGTGGGTAACGTTTTTCCCCGGGCCACGGTAGCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCTGCCATTGAt                                                                                                                                    <  1:885811/148‑1 (MQ=255)
            aCCCGGCAGAACGATTTCAATATGTTCCCAACCTTCGTGTGGGTAACGTTTTTCCCCGGGCCACGGTAGCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATCGGTCTGCCATTGATc                                                                                                                                   <  2:3340551/148‑1 (MQ=255)
            aCCCGGCAGAACGATTTCAATATGTTCCCAACCTTCGTGTGGGTAACGTTTTTCCCCGGGCCACGGTAGCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCTGCCATTGATc                                                                                                                                   <  1:2194621/148‑1 (MQ=255)
            aCCCGGCAGAACGATTTCAATATGTTCCCAACCTTCGTGTGGGTAACGTTTTTCCCCGGGCCACGGTAGCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCTGCCATTGATc                                                                                                                                   <  2:1894417/148‑1 (MQ=255)
             cccGGCAGAACGATTTCAATATGTTCCCAACCTTCGTGTGGGTAACGTTTTTCCCCGGGCCACGGTAGCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCTGCCATTGATCa                                                                                                                                  >  2:438045/1‑148 (MQ=255)
              ccGGCAGAACGATTTCAATATGTTCCCAACCTTCGTGTGGGTAACGTTTTTCCCCGGGCCACGGTAGCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCACTTTAAACAGACAAATAGGTCTGCCATTGATCat                                                                                                                                 >  2:841940/1‑148 (MQ=255)
              ccGGCAGAACGATTTCAATATGTTCCCAACCTTCGTGTGGGTAACGTTTTTCCCCGGGCCACGGTAGCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCTGCCATTGATCat                                                                                                                                 <  1:3026340/148‑1 (MQ=255)
                gtcAGAACGATTTCAATCTGTTCCCAACCTTCGTGTGGGTAACGTTTTTCCCCGGGCCACGGTAGCTCCACAATAGAAAACTGCCAGTGCTCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCTGCCATTGATCatat                                                                                                                               <  2:3497191/146‑1 (MQ=255)
                ggCAGAACGATTTCAATATGTTCCCAACCTTCGTGTGGGTAACGTTTTTCCCCGGGCCACGGTAGCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCTGCCATTGATCatat                                                                                                                               >  1:2123312/1‑148 (MQ=255)
                ggCAGAACGATTTCAATATGTTCCCAACCTTCGTGTGGGTAACGTTTTTCCCCGGGCCACGGTAGCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCTGCCATTGATCatat                                                                                                                               <  1:423226/148‑1 (MQ=255)
                 gCAGAACGATTTCAATATGTTCCCAACCTTCGTGTGGGTAACGTTTTTCCCCGGGCCACGGTAGCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCTGCCATTGATCATAtt                                                                                                                              <  2:3189115/148‑1 (MQ=255)
                 gCAGAACGATTTCAATATGTTCCCAACCTTCGTGTGGGTAACGTTTTTCCCCGGGCCACGGTAGCTCCACAATAGAAAACTGCCAGTGCGCAAACTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCTGCCATTGATCATAtt                                                                                                                              >  2:3383385/1‑148 (MQ=255)
                  cAGAACGATTTCAATATGTTCCCAACCTTCGTGTGGGTAACGTTTTTCCCCGGGCCACGGTAGCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCTGCCATTGATCATAttt                                                                                                                             >  1:534174/1‑148 (MQ=255)
                    gAACGATTTCAATATGTTCCCAACCTTCGTGTGGGTAACGTTTTTCCCCGGGCCACGGTAGCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCTGCCATTGATCATATTTTc                                                                                                                           <  2:2268234/148‑1 (MQ=255)
                    gAACGATTTCAATATGTTCCCAACCTTCGTGTGGGTAACGTTTTTCCCCGGGCCACGGTAGCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCTGCCATTGATCATATTTTc                                                                                                                           >  2:2426623/1‑148 (MQ=255)
                     aaCGATTTCAATATGTTCCCAACCTTCGTGTGGGTAACGTTTTTCCCCGGGCCACGGTAGCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCTGCCATTGATCATATTTTCt                                                                                                                          >  2:2039752/1‑148 (MQ=255)
                       cGATTTCAATATGTTCCCAACCTTCGTGTGGGTAACGTTTTTCCCCGGGCCACGGTAGCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCTGCCATTGATCATATTTTCTGa                                                                                                                        >  1:2938073/1‑148 (MQ=255)
                          tttCAATATGTTCCCAACCTTCGTGTGGGTAACGTTTTTCCCCGGGCCACGGTAGCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCTTGCAATTTAAACAGACAAATAGGTCTGCCATTGATCATATTTTCTGACaa                                                                                                                     >  2:818739/1‑148 (MQ=255)
                          tttCAATATGTTCCCAACCTTCGTGTGGGTAACGTTTTTCCCCGGGCCACGGTAGCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCACTTTAAACAGACAAATAGGTCTGCCATTGATCATATTTTCTGACaa                                                                                                                     >  2:1842528/1‑148 (MQ=255)
                          tttCAATATGTTCCCAACCTTCGTGTGGGTAACGTTTTTCCCCGGGCCACGGTAGCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCTGCCATTGATCATATTTTCTGACaa                                                                                                                     <  1:2002084/148‑1 (MQ=255)
                            tCAATATGTTCCCAACCTTCGTGTGGGTAACGTTTTTCCCCGGGCCACGGTAGCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATCGGTCTGCCATTGATCATATTTTCTGACaaaa                                                                                                                   >  1:3019786/1‑148 (MQ=255)
                              aaTATGTTCCCAACCTTCGTGTGGGTAACGTTTTTCCCCGGGCCACGGTAGCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCTGCCATTGATCATATTTTCTGACAAAAGc                                                                                                                 <  2:1634821/148‑1 (MQ=255)
                              aaTATGTTCCCAACCTTCGTGTGGGTAACGTTTTTCCCCGGGCCACGGTAGCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCTGCCATTGATCATATTTTCTGACAAAAGc                                                                                                                 <  2:2003920/148‑1 (MQ=255)
                               atatGTTCCCAACCTTCGTGTGGGTAACGTTTTTCCCCGGGCCACGGTAGCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCTGCCATTGATCATATTTTCTGACAAAAGCt                                                                                                                >  2:709299/1‑148 (MQ=255)
                                tatGTTCCCAACCTTCGTGTGGGTAACGTTTTTCCCCGGGCCACGGTAGCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCTGCCATTGATCATATTTTCTGACAAAAGCTc                                                                                                               >  1:2500844/1‑148 (MQ=255)
                                 atGTTCCCAACCTTCGTGTGGGTAACGTTTTTCCCCGGGCCACGGTAGCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCg                                                                                                              <  2:1236925/148‑1 (MQ=255)
                                  tGTTCCCAACCTTCGTGTGGGTAACGTTTTTCCCCGGGCCACGGTAGCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGc                                                                                                             <  2:1128988/148‑1 (MQ=255)
                                  tGTTCCCAACCTTCGTGTGGGTAACGTTTTTCCCCGGGCCACGGTAGCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGc                                                                                                             >  1:2052829/1‑148 (MQ=255)
                                   gTTCCCAACCTTCGTGTGGGTAACGTTTTTCCCCGGGCCACGGTAGCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCACTTTAAACAGACAAATAGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGcc                                                                                                            >  2:2602190/1‑148 (MQ=255)
                                   gTTCCCAACCTTCGTGTGGGTAACGTTTTTCCCCGGGCCACGGTAGCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGcc                                                                                                            >  2:265622/1‑148 (MQ=255)
                                    ttCCCAACCTTCGTGTGGGTAACGTTTTTCCCCGGGCCACGGTAGCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCca                                                                                                           >  1:1625623/1‑148 (MQ=255)
                                       ccAACCTTCGTGTGGGTAACGTTTTTCCCCGGGCCACGGTAGCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCacac                                                                                                        >  1:1339043/1‑148 (MQ=255)
                                       ccAACCTTCGTGTGGGTAACGTTTTTCCCCGGGCCACGGTAGCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCacac                                                                                                        <  1:1340945/148‑1 (MQ=255)
                                        cAACCTTCGTGTGGGTAACGTTTTTCCCCGGGCCACGGTAGCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACt                                                                                                       >  1:2531774/1‑148 (MQ=255)
                                          aCCTTCGTGTGGGTAACGTTTTTCCCCGGGCCACGGTAGCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGt                                                                                                     <  1:673793/148‑1 (MQ=255)
                                            cTTCGTGTGGGTAACGTTTTTCCCCGGGCCACGGTAGCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTc                                                                                                   >  1:3145925/1‑148 (MQ=255)
                                               cGTGTGGGTAACGTTTTTCCCCGGGCCACGGTAGCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCaaa                                                                                                <  1:632411/148‑1 (MQ=255)
                                                    gggTAACGTTTTTCCCCGGGCCACGGTAGCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCcgc                                                                                           >  2:2560567/1‑148 (MQ=255)
                                                    gggTAACGTTTTTCCCCGGGCCACGGTAGCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCcgc                                                                                           >  2:479886/1‑148 (MQ=255)
                                                     ggTAACGTTTTTCCCTGGGCCACGGTAGCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCcgcg                                                                                          <  1:1685293/148‑1 (MQ=255)
                                                      gTAACGTTTTTCCCCGGGCCACGGTAGCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGa                                                                                         <  2:2920541/148‑1 (MQ=255)
                                                          cGTTTTTCCCCGGGCCACGGTAGCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGcc                                                                                     >  2:40956/1‑148 (MQ=255)
                                                             ttttCCCCGGGCCACGGTAGCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGCCAGc                                                                                  >  2:2697778/1‑148 (MQ=255)
                                                                 ccccGGGCCACGGTAGCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGCCAGCGTTc                                                                              >  1:1796402/1‑148 (MQ=255)
                                                                  cccGGGCCACGGTAGCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGCCAGCGTTca                                                                             >  1:1957476/1‑148 (MQ=255)
                                                                   ccGGGCCACGGTAGCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGCCAGCg                                                                                 <  2:531219/143‑1 (MQ=255)
                                                                    cGGGCCACGGTAGCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGCCAGCGTTcagc                                                                           >  1:1088398/1‑148 (MQ=255)
                                                                    cGGGCCACGGTAGCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGCCAGCGTTcagc                                                                           >  2:1455709/1‑148 (MQ=255)
                                                                      ggCCACGGTAGCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGCCAGCGTTcagcag                                                                         <  2:157895/148‑1 (MQ=255)
                                                                            ggTAGCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGCCAGCGTTCAGCAGTGGCGt                                                                   <  1:769514/148‑1 (MQ=255)
                                                                             gTAGCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGCCAGCGTTCAGCAGTGGCGtt                                                                  >  1:3255289/1‑148 (MQ=255)
                                                                             gTAGCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGCCAGCGTTCAGCAGTGGCGtt                                                                  >  1:2195110/1‑148 (MQ=255)
                                                                              tAGCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGCCAGCGTTCAGCAGTGGCGttt                                                                 >  2:1289275/1‑148 (MQ=255)
                                                                                gCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGCCAGCGTTCAGCAGTGGCGTTTTg                                                               <  1:3351637/148‑1 (MQ=255)
                                                                                  tCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGCCAGCGTTCAGCAGTGGCGTTTTGAt                                                             <  1:1358096/148‑1 (MQ=255)
                                                                                  tCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGCCAGCGTTCAGCAGTGGCGTTTTGAt                                                             <  1:691636/148‑1 (MQ=255)
                                                                                   ccACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGCCAGCGTTCAGCAGTGGCGTTTTGATg                                                            >  2:934208/1‑148 (MQ=255)
                                                                                   ccACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGCCAGCGTTCAGCAGTGGCGTTTTGATg                                                            <  2:2153451/148‑1 (MQ=255)
                                                                                          aGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGCCAGCGTTCAGCAGTGGCGTGTTGATGGCAACGc                                                     >  1:2368419/1‑148 (MQ=255)
                                                                                           gAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGCCAGCGTTCAGCAGTGGCGTTTTGATGGCAACGCa                                                    >  2:607479/1‑148 (MQ=255)
                                                                                            aaaaCTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGCCAGCGTTCAGCAGTGGCGTTTTGATGGCAACGCaa                                                   >  1:3337943/1‑148 (MQ=255)
                                                                                             aaaCTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGCCAGCGTTCAGCAGTGGCGTTTTGATGGCAACGCaaa                                                  >  2:1339044/1‑148 (MQ=255)
                                                                                             aaaCTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGCCAGCGTTCAGCAGTGGCGTTTTGATGGCAACGCaaa                                                  <  2:151339/148‑1 (MQ=255)
                                                                                             aaaCTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGCCAGCGTTCAGCAGTGGCGTTTTGATGGCAACGCaaa                                                  <  2:151303/148‑1 (MQ=255)
                                                                                                cTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGCCAGCGTTCAGCAGTGGCGTTTTGATGGCAACGCAAAGaa                                               <  1:2608366/148‑1 (MQ=255)
                                                                                                cTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGCCAGCGTTCAGCAGTGGCGTTTTGATGGCAACGCAAAGaa                                               <  1:2330681/148‑1 (MQ=255)
                                                                                                cTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGCCAGCGTTCAGCAGTGGCGTTTTGATGGCAACGCAAAGaa                                               <  1:2330775/148‑1 (MQ=255)
                                                                                                cTGACAGTGCGCAACCTGTACCGGTTCATGAAATTTAAACAGACAAATAGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGCCAGCGTTCAGCAGTGGCGTTTTGATGGCAACGCAAAGaa                                               <  1:1276888/148‑1 (MQ=255)
                                                                                                 tGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATCGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGCCAGCGTTCAGCAGTGGCGTTTTGATGGCAACGCAAAGaaa                                              <  2:1523224/148‑1 (MQ=255)
                                                                                                   ccAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGCCAGCGTTCAGCAGTGGCGTTTTGATGGCAACGCAAAGAAAtg                                            >  1:3300439/1‑148 (MQ=255)
                                                                                                    cAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGCCAGCGTTCAGCAGTGGCGTTTTGATGGCAACGCAAAGAAAtgt                                           >  1:2839381/1‑148 (MQ=255)
                                                                                                     aGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGCCAGCGTTCAGCAGTGGCGTTTTGATGGCAACGCAAAGAAAtgtg                                          >  1:3467765/1‑148 (MQ=255)
                                                                                                      gTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGCCAGCGTTCAGCAGTGGCGTTTTGATGGCAACGCAAAGAAATGTGa                                         <  2:2013159/148‑1 (MQ=255)
                                                                                                       tGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGCCAGCGTTCAGCAGTGGCGTTTTGATGGCAACGCAAAGAAATGTGAt                                        >  1:1974558/1‑148 (MQ=255)
                                                                                                         cgcAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGCCAGCGTTCAGCAGTGGCGTTTTGATGGCAACGCAAAGAAATGTGATCg                                      <  1:2139743/148‑1 (MQ=255)
                                                                                                          gcAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTGAATGCCGCGACGCCTTCTTTCAGCAGTGGCGTTTTGATGGCAACGCAAAGAAATGTGATCgg                                     <  1:3448121/148‑1 (MQ=255)
                                                                                                              ccTGTACCGGTTCATGCAATTTAAACAGACAAATCGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGCCAGCGTTCAGCAGTGGCGTTTTGATGGCAACGCAAAGAAATGTGATCGGCAGt                                 >  1:100429/1‑148 (MQ=255)
                                                                                                              ccTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGCCAGCGTTCAGCAGTGGCGTTTTGATGGCAACGCAAAGAAATGTGATCGGCAGt                                 <  1:2697778/148‑1 (MQ=255)
                                                                                                                 gTACCGGTTCATGCAATTTAAACAGACAAATAGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGCCAGCGTTCAGCAGTGGCGTTTTGATGGCAACGCAAAGAAATGTGATCGGCAGTCAg                              <  2:3418818/148‑1 (MQ=255)
                                                                                                                    ccGGTTCATGCAATTTAAACAGACAAATAGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGCCAGCGTTCAGCAGTGGCGTTTTGATGGCAACGCAAAGAAATGTGATCGGCAGTCAGCgg                           >  2:1198839/1‑148 (MQ=255)
                                                                                                                     cGGTTCATGCAATTTAAACAGACAAATAGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGCCAGCGTTCAGCAGTGGCGTTTTGATGGCAACGCAAAGAAATGTGATCGGCAGTCAGCGGa                          <  2:3025098/148‑1 (MQ=255)
                                                                                                                     cGGTTCATGCAATTTAAACAGACAAATAGGTCTGCCATTGAGCATATTTTCTGACAAAAGCTCGCCACCCTGTTCAAACCCGCGACGCCAGCGTTCAGCAGTGGCGTTTTGATGGCAACGCAAAGAAATGTGATCGGCAGTCAGCGGa                          <  2:180185/148‑1 (MQ=255)
                                                                                                                      ggTTCATGCAATTTAAACAGACAAATAGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGCCAGCGTTCAGCAGTGGCGTTTTGATGGCAACGCAAAGAAATGTGATCGGCAGTCAGCGGAg                         <  1:3371693/148‑1 (MQ=255)
                                                                                                                        ttCATGCAATTTAAACAGACAAATAGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGCCAGCGTTCAGCAGTGGCGTTTTGATGGCAACGCAAAGAAATGTGATCGGCAGTCAGCGGAGTg                       >  1:417005/1‑148 (MQ=255)
                                                                                                                          cATGCAATTTAAACAGACAAATAGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGCCAGCGTTCAGCAGTGGCGTTTTGATGGCAACGCAAAGAAATGTGATCGGCAGTCAGCGGAGTGat                     <  1:1509833/148‑1 (MQ=255)
                                                                                                                           aTGCAATTTAAACAGACAAATAGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACTCTGTTCAAACCCGCGACGCCAGCGTTCAGCAGTGGCGTTGTGATGGCAACGCAAAGAAATGTGATCGGCAGTCAGCGGAGTGata                    <  1:1018859/148‑1 (MQ=255)
                                                                                                                              cAATTTAAACAGACAAATAGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGCCAGCGTTCAGCAGTGGCGTTTTGATGGCAACGCAAAGAAATGTGATCGGCAGTCAGCGGAGTGATATTc                 >  2:2885246/1‑148 (MQ=255)
                                                                                                                               aaTTTAAACAGACAAATAGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGCCAGCGTTCAGCAGTGGCGTTTTGATGGCAACGCAAAGAAATGTGATCGGCAGTCAGCGGAGTGATATTCa                <  1:818739/148‑1 (MQ=255)
                                                                                                                                 tttAAACAGACAAATAGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGCCAGCGTTCAGCAGTGGCGTTTTGATGGCAACGCAAAGAAATGTGATCGGCAGTCAGCGGAGTGATATTCAAc              >  1:575860/1‑148 (MQ=255)
                                                                                                                                       cAGACAAATAGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGCCAGCGTTCAGCAGTGGCGTTTTGATGGCAACGCAAAGAAATGTGATCGGCAGTCAGCGGAGTGATATTCAACCCCAGa        >  1:2624038/1‑148 (MQ=255)
                                                                                                                                       cAGACAAATAGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGCCAGCGTTCAGCAGTGGCGTTTTGATGGCAACGCAAAGAAATGTGATCGGCAGTCAGCGGAGTGATATTCAACCCCAGa        <  2:957451/148‑1 (MQ=255)
                                                                                                                                       cAGACAAATAGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGCCAGCGTTCAGCAGTGGCGTGTTGATGGCCACGCAAAGAAATCTGATCGGCACTCAGAGGAGTGATATTCAAACCCAGa        >  1:981909/1‑148 (MQ=255)
                                                                                                                                        aGACAAATAGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGCCAGCGTTCAGCAGTGGCGTTTTGATGGCAACGCAAAGAAATGTGATCGGCAGTCAGCGGAGTGATATTCAACCCCAGAc       <  1:2302830/148‑1 (MQ=255)
                                                                                                                                           cAAATAGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGCCAGCGTTCAGCAGTGGCGTTTTGATGGCAACGCAAAGAAATGTGATCGGCAGTCAGCGGAGTGATATTCAACCCCAGAcggc    >  2:122219/1‑148 (MQ=255)
                                                                                                                                                |                                                                                                                                                
CGTTTCTGGATCACCCGGCAGAACGATTTCAATATGTTCCCAACCTTCGTGTGGGTAACGTTTTTCCCCGGGCCACGGTAGCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATCGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGCCAGCGTTCAGCAGTGGCGTTTTGATGGCAACGCAAAGAAATGTGATCGGCAGTCAGCGGAGTGATATTCAACCCCAGACGGCGG  >  CP009273/1953732‑1954020

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: