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,8760CA98.8% 850.0 / 0.3 264P76P (CCG→CCTyecMputative metal‑binding enzyme
Reads supporting (aligned to +/- strand):  ref base C (0/4);  new base A (133/127);  total (133/131)
Fisher's exact test for biased strand distribution p-value = 5.92e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 264 total aligned reads displayed.

AGCGTTTCTGGATCACCCGGCAGAACGATTTCAATATGTTCCCAACCTTCGTGTGGGTAACGTTTTTCCCCGGGCCACGGTAGCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATCGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGCCAGCGTTCAGCAGTGGCGTTTTGATGGCAACGCAAAGAAATGTGATCGGCAGTCAGCGGAGTGATATTCAACCCCAGACGGCGGG  >  CP009273/1953730‑1954021
                                                                                                                                                  |                                                                                                                                                 
aGCGTTTCTGGATCACCCGGCAGAACGATTTCAATATGTTCCCAACCTTCGTGTGGGTAACGTTTTTCCCCGGGCCACGGTAGCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAgg                                                                                                                                                 >  1:3665341/1‑149 (MQ=255)
 gCGTTTCTGGATCACCCGGCAGAACGATTTCAATATGTTCCCAACCTTCGTGTGGGTAACGTTTTTCCCCGGGCCACGGTAGCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGt                                                                                                                                                <  2:3517941/149‑1 (MQ=255)
 gCGTTTCTGGATCACCCGGCAGAACGATTTCAATATGTTCCCAACCTTCGTGTGGGTAACGTTTTTCCCCGGGCCACGGTAGCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGt                                                                                                                                                <  1:2098763/149‑1 (MQ=255)
  cGTTTCTGGATCACCCGGCAGAACGATTTCAATATGGTCCCAACCTTCGTGTGGGTAACGTTTTTCCCCGGGCCACGGTAGCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTc                                                                                                                                               <  1:404132/149‑1 (MQ=255)
   gTTTCTGGATCACCCGTCAGAACGATTTCAATATGTTCCCAACCTTCGTGTGGGTAACGTTTTTCCCCGGGCCACGGTAGCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCt                                                                                                                                              <  2:368943/149‑1 (MQ=255)
       cTGGATCACCCGGCAGAACGATTTCAATATGTTCCCAACCTTCGTGTGGGTAACGTTTTTCCCCGGGCCACGGTAGCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCTGCCa                                                                                                                                          <  1:2048498/149‑1 (MQ=255)
        tGGATCACCCGGCAGAACGATTTCAATATGTTCCCAACCTTCGTGTGGGTAACGTTTTTCCCCGGGCCACGGTAGCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCTGCCAt                                                                                                                                         <  1:820011/149‑1 (MQ=255)
         ggATCACCCGGCAGAACGATTTCAATATGTTCCCAACCTTCGTGTGGGTAACGTTTTTCCCCGGGCCACGGTAGCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCTGCCAtt                                                                                                                                        >  2:2886847/1‑149 (MQ=255)
           aTCACCCGGCAGAACGATTTCAATATGTTCCCAACCTTCGTGTGGGTAACGTTTTTCCCCGGGCCACGGTAGCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCTGCCATTGa                                                                                                                                      <  2:737620/149‑1 (MQ=255)
           aTCACCCGGCAGAACGATTTCAATATGTTCCCAACCTTCGTGTGGGTAACGTTTTTCCCCGGGCCACGGTAGCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCTGCCATTGa                                                                                                                                      <  2:4105414/149‑1 (MQ=255)
            tCACCCGGCAGAACGATTTCAATATGTTCCCAACCTTCGAGTGGGTAACGTTTTTCCCCGGGCCACGGTAGCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCTGCCATTGAt                                                                                                                                     >  1:1784191/1‑149 (MQ=255)
             cACCCGGCAGAACGATTTCAATATGTTCCCAACCTTCGTGTGGGTAACGTTTTTCCCCGGGCCACGGTAGCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCTGCCATTGATc                                                                                                                                    <  2:1698918/149‑1 (MQ=255)
                  ggCAGAACGATTTCAATATGTTCCCAACCTTCGTGTGGGTAACGTTTTTCCCCGGGCCACGGTAGCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCTGCCATTGATCATAtt                                                                                                                               <  1:1991/149‑1 (MQ=255)
                   gCAGAACGATTTCAATATGTTCCCAACCTTCGTGTGGGTAACGTTTTTCCCCGGGCCACGGTAGCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCTGCCATTGATCATAttt                                                                                                                              >  2:356674/1‑149 (MQ=255)
                   gCAGAACGATTTCAATATGTTCCCAACCTTCGTGTGGGTAACGTTTTTCCCCGGGCCACGGTAGCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCTGCCATTGATCATAttt                                                                                                                              <  2:308615/149‑1 (MQ=255)
                      gAACGATTTCAATATGTTCCCAACCTTCGTGTGGGTAACGTTTTTCCCCGGGCCACGGTAGCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCTGCCATTGATCATATTTTCt                                                                                                                           <  1:155950/149‑1 (MQ=255)
                       aaCGATTTCAATATGTTCCCAACCTTCGTGTGGGTAACGTTTTTCCCCGGGCCACGGTAGCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCTGCCATTGATCATATTTTCTg                                                                                                                          >  2:2736743/1‑149 (MQ=255)
                       aaCGATTTCAATATGTTCCCAACCTTCGTGTGGGTAACGTTTTTCCCCGGGCCACGGTAGCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCTGCCATTGATCATATTTTCTg                                                                                                                          >  2:561086/1‑149 (MQ=255)
                           aTTTCAATATGTTCCCAACCTTCGTGTGGGTAACGTTTTTCCCCGGGCCACGGTAGCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCTGCCATTGATCATATTTTCTGACaa                                                                                                                      >  2:406413/1‑149 (MQ=255)
                           aTTTCAATATGTTCCCAACCTTCGTGTGGGTAACGTTTTTCCCCGGGCCACGGTAGCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCTGCCATTGATCATATTTTCTGACaa                                                                                                                      >  2:2473460/1‑149 (MQ=255)
                            tttCAATATGTTCCCAACCTTCGTGTGGGTAACGTTTTTCCCCGGGCCACGGTAGCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCTGCCATTGATCATATTTTCTGACaaa                                                                                                                     <  2:204465/149‑1 (MQ=255)
                             ttCAATATGTTCCCAACCTTCGTGTGGGTAACGTTTTTCCCCGGGCCACGGTAGCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCTGCCATTGATCATATTTTCTGACaaaa                                                                                                                    >  2:3293961/1‑149 (MQ=255)
                               cAATATGTTCCCAACCTTCGTGTGGGTAACGTTTTTCCCCGGGCCACGGTAGCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCTGCCATTGATCATATTTTCTGACAAAAGc                                                                                                                  >  2:298001/1‑149 (MQ=255)
                                aaTATGTTCCCAACCTTCGTGTGGGTAACGTTTTTCCCCGGGCCACGGTAGCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCTGCCATTGATCATATTTTCTGACAAAAGCt                                                                                                                 >  2:3464875/1‑149 (MQ=255)
                                aaTATGTTCCCAACCTTCGTGTGGGTAACGTTTTTCCCCGGGCCACGGTAGCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCTGCCATTGATCATATTTTCTGACAAAAGCt                                                                                                                 <  1:300296/149‑1 (MQ=255)
                                 atatGTTCCCAACCTTCGTGTGGGTAACGTTTTTCCCCGGGCCACGGTAGCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCTGCCATTGATCATATTTTCTGACAAAAGCTc                                                                                                                >  2:2065288/1‑149 (MQ=255)
                                   atGTTCCCAACCTTCGTGTGGGTAACGTTTTTCCCCGGGCCACGGTAGCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGc                                                                                                              >  2:2829016/1‑149 (MQ=255)
                                     gTTCCCAACCTTCGTGTGGGTAACGTTTTTCCCCGGGCCACGGTATAGACACAATAGAACAATGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCca                                                                                                            >  1:4034526/1‑149 (MQ=255)
                                            aCCTTCGTGTGGGTAACGTTTTTCCCCGGGCCACGGTAGCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGtt                                                                                                     >  1:1353904/1‑149 (MQ=255)
                                            aCCTTCGTGTGGGTAACGTTTTTCCCCGGGCCACGGTAGCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGtt                                                                                                     >  1:2789563/1‑149 (MQ=255)
                                            aCCTTCGTGTGGGTAACGTTTTTCCCCGGGCCACGGTAGCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGtt                                                                                                     <  1:1505620/149‑1 (MQ=255)
                                            aCCTTCGTGTGGGTAACGTTTTTCCCCGGGCCACGGTAGCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGtt                                                                                                     <  1:1921011/149‑1 (MQ=255)
                                             ccTTCGTGTGGGTAACGTTTTTCCCCGGGCCACGGTAGCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTc                                                                                                    >  2:653749/1‑149 (MQ=255)
                                               ttCGTGTGGGTAACGTTTTTCCCCGGGCCACGGTAGCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCaa                                                                                                  <  1:3153089/149‑1 (MQ=255)
                                                 cGTGTGGGTAACGTTTTTCCCCGGGCCACGGTAGCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAAc                                                                                                <  2:1181816/149‑1 (MQ=255)
                                                 cGTGTGGGTAACGTTTTTCCCCGGGCCACGGTAGCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAAc                                                                                                <  2:3276774/149‑1 (MQ=255)
                                                   tgtgGGTAACGTTTTTCCCCGGGCCACGGTAGCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAAccc                                                                                              >  1:2586088/1‑149 (MQ=255)
                                                   tgtgGGTAACGTTTTTCCCCGGGCCACGGTAGCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAAccc                                                                                              <  2:3602926/149‑1 (MQ=255)
                                                     tgGGTAACGTTTTTCCCCGGGCCACGGTAGCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCcgc                                                                                            <  2:4014553/149‑1 (MQ=255)
                                                      gggTAACGTTTTTCCCCGGGCCACGGTAGCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCcgcg                                                                                           >  1:3871780/1‑149 (MQ=255)
                                                        gTAACGTTTTTCCCCGGGCCACGGTAGCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGAc                                                                                         >  1:458311/1‑149 (MQ=255)
                                                           aCGTTTTTCCCCGGGCCACGGTAGCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGcc                                                                                      <  2:2607359/149‑1 (MQ=255)
                                                            cGTTTTTCCCCGGGCCACGGTAGCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGCCa                                                                                     <  1:1435216/149‑1 (MQ=255)
                                                              tttttCCCCGGGCCACGGTAGCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGCCAGc                                                                                   >  2:1970680/1‑149 (MQ=255)
                                                              tttttCCCCGGGCCACGGTAGCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGCCAGc                                                                                   >  1:1060692/1‑149 (MQ=255)
                                                              tttttCCCCGGGCCACGGTAGCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGCCAGc                                                                                   >  1:599960/1‑149 (MQ=255)
                                                                    cccGGGCCACGGTAGCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGCCAGCGTTcag                                                                             <  1:2828874/149‑1 (MQ=255)
                                                                     ccGGGCCACGGTAGCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGCCAGCGTTcagc                                                                            >  1:4101716/1‑149 (MQ=255)
                                                                      cGGGCCACGGTAGCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGCCAGCGTTcagca                                                                           <  2:1986778/149‑1 (MQ=255)
                                                                        ggCCACGGTAGCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGCCAGCGTTCAGCAGt                                                                         <  2:1685927/149‑1 (MQ=255)
                                                                         gCCACGGTAGCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGCCAGCGTTCAGCAGTg                                                                        >  1:1478650/1‑149 (MQ=255)
                                                                           cACGGTAGCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCTGCCATTGATCATTTTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGCCAGCGTTCAGCAGTGGc                                                                      <  1:3617440/149‑1 (MQ=255)
                                                                           cACGGTAGCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGCCAGCGTTCAGCAGTGGc                                                                      <  1:2680079/149‑1 (MQ=255)
                                                                             cGGTAGCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGCCAGCGTTCAGCAGTGGCGt                                                                    >  1:827379/1‑149 (MQ=255)
                                                                               gTAGCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGCCAGCGTTCAGCAGTGGCGttt                                                                  <  1:2088684/149‑1 (MQ=255)
                                                                               gTAGCTACACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCTGCCATTGATCTTATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGCCAGCGTTCAGCAGTGGCGttt                                                                  <  1:259174/149‑1 (MQ=255)
                                                                                tAGCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGCCAGCGTTCAGCAGTGGCGtttt                                                                 >  1:187712/1‑149 (MQ=255)
                                                                                tAGCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGCCAGCGTTCAGCAGTGGCGtttt                                                                 <  1:2282196/149‑1 (MQ=255)
                                                                                tAGCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGCCAGCGTTCAGCAGTGGCGtttt                                                                 >  1:2263091/1‑149 (MQ=255)
                                                                                  gCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGCCAGCGTTCAGCAGTGGCGTTTTGa                                                               <  2:1379716/149‑1 (MQ=255)
                                                                                     ccACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGCCAGCGTTCAGCAGTGGCGTTTTGATgg                                                            >  1:883021/1‑149 (MQ=255)
                                                                                       acaATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGCCAGCGTTCAGCAGTGGCGTTTTGATGGCa                                                          >  1:1303114/1‑149 (MQ=255)
                                                                                         aaTAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGCCAGCGTTCAGCAGTGGCGTTTTGATGGCAAc                                                        >  1:32019/1‑149 (MQ=255)
                                                                                            aGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGCCAGCGTTCAGCAGTGGCGTTTTGATGGCAACGCa                                                     >  1:2943910/1‑149 (MQ=255)
                                                                                             gAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGCCAGCGTTCAGCAGTGGCGTTTTGATGGCAACGCaa                                                    >  2:2257739/1‑149 (MQ=255)
                                                                                                aaCTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGCCAGCGTTCAGCAGTGGCGTTTTGATGGCAACGCAAAGa                                                 <  2:678244/149‑1 (MQ=255)
                                                                                                 aCTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGCCAGCGTTCAGCAGTGGCGTTTTGATGGCAACGCAAAGaa                                                >  1:3792575/1‑149 (MQ=255)
                                                                                                 aCTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGCCAGCGTTCAGCAGTGGCGTTTTGATGGCAACGCAAAGaa                                                >  2:3856909/1‑149 (MQ=255)
                                                                                                  cTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAAGAGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGCCAGCGTTCAGCAGTGGCGTTTTGATGGCAACGCAAAGaaa                                               <  1:2153638/149‑1 (MQ=255)
                                                                                                    gCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAACTAGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGCCAGCGTTCAGCAGTGGCGTTTTGATGGCAACGCAAAGAAAtg                                             <  2:626595/149‑1 (MQ=255)
                                                                                                      cAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGCCAGCGTTCAGCAGTGGCGTTTTGATGGCAACGCAAAGAAAtgtg                                           <  1:3672717/149‑1 (MQ=255)
                                                                                                       aGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGCCAGCGTTCAGCAGTGGCGTTTTGATGGCAACGCAAAGAAATGTGa                                          <  1:544969/149‑1 (MQ=255)
                                                                                                       aGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGCCAGCGTTCAGCAGTGGCGTTTTGATGGCAACGCAAAGAAATGTGa                                          <  1:2886847/149‑1 (MQ=255)
                                                                                                        gTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGCCAGCGTTCAGCAGTGGCGTTTTGATGGCAACGCAAAGAAATGTGAt                                         <  1:1839227/149‑1 (MQ=255)
                                                                                                         tGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGCCAGCGTTCAGCAGTGGCGTTTTGATGGCAACGCAAAGAAATGTGATc                                        <  2:474992/149‑1 (MQ=255)
                                                                                                             cAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGCCAGCGTTCAGCAGTGGCGTTTTGATGGCAACGCAAAGAAATGTGATCGGCa                                    <  1:812368/149‑1 (MQ=255)
                                                                                                             cAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGCCAGCGTTCAGCAGTGGCGTTTTGATGGCAACGCAAAGAAATGTGATCGGCa                                    <  1:1698110/149‑1 (MQ=255)
                                                                                                              aaCCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGCCAGCGTTCAGCAGTGGCGTTTTGATGGCAACGCAAAGAAATGTGATCGGCAg                                   <  2:1835887/149‑1 (MQ=255)
                                                                                                                ccTGTACCGGTTCATGCAATTTAAACAGACAAATCGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGCCAGCGTTCAGCAGTGGCGTTTTGATGGCAACGCAAAGAAATGTGATCGGCAGTc                                 <  1:1004642/149‑1 (MQ=255)
                                                                                                                ccTGTACCGGTTCATGCAATTTAAACAGACAAATCGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGCCAGCGTTCAGCAGTGGCGTTTTGATGGCAACGCAAAGAAATGTGATCGGCAGTc                                 <  2:2328709/149‑1 (MQ=255)
                                                                                                                ccTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGCCAGCGTTCAGCAGTGGCGTTTTGATGGCAACGCAAAGAAATGTGATCGGCAGTc                                 >  2:918256/1‑149 (MQ=255)
                                                                                                                ccTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGCCAGCGTTCAGCAGTGGCGTTTTGATGGCAACGCAAAGAAATGTGATCGGCAGTc                                 <  2:1560522/149‑1 (MQ=255)
                                                                                                                ccTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGCCAGCGTTCAGCAGTGGCGTTTTGATGGCAACGCAAAGAAATGTGATCGGCAGTc                                 <  1:3293961/149‑1 (MQ=255)
                                                                                                                 cTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGCCAGCGTTCAGCAGTGGCGTTTTGATGGCAACGCAAAGAAATGTGATCGGCAGTCa                                >  2:2351057/1‑149 (MQ=255)
                                                                                                                  tGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGCCAGCGTTCAGCAGTGGCGTTTTGATGGCAACGCAAAGAAATGTGATCGGCAGTCAg                               >  1:690665/1‑149 (MQ=255)
                                                                                                                      ccGGTTCATGCAATTTAAACAGACAAATCGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGCCAGCGTTCAGCAGTGGCGTTTTGATGGCAACGCAAAGAAATGTGATCGGCAGTCAGCGGa                           <  2:2953646/149‑1 (MQ=255)
                                                                                                                        ggTTCATGCAATTTAAACAGACAAATAGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGCCAGCGTTCAGCAGTGGCGTTTTGATGGCAACGCAAAGAAATGTGATCGGCAGTCAGCGGAGt                         <  2:2784999/149‑1 (MQ=255)
                                                                                                                        ggTTCATGCAATTTAAACAGACAAATAGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGCCAGCGTTCAGCAGTGGCGTTTTGATGGCAACGCAAAGAAATGTGATCGGCAGTCAGCGGAGt                         >  2:1510452/1‑149 (MQ=255)
                                                                                                                        ggTTCATGCAATTTAAACAGACAAATAGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGCCAGCGTTCAGCAGTGGCGTGTTGATGGCAACGCAAAGAAATGTGATCGGCAGTCAGCGGAGt                         >  2:1078879/1‑149 (MQ=255)
                                                                                                                         gTTCATGCAATTTAAACAGACAAATAGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGCCAGCGTTCAGCAGTGGCGTTTTGATGGCAACGCAAAGAAATGTGATCGGCAGTCAGCGGAGTg                        >  2:2658772/1‑149 (MQ=255)
                                                                                                                           tCATGCAATTTAAACAGACAAATAGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGCCAGCGTTCAGCAGTGGCGTTTTGATGGCAACGCAAAGAAATGTGATCGGCAGTCAGCGGAGTGat                      <  1:2351503/149‑1 (MQ=255)
                                                                                                                             aTGCAATTTAAACAGACAAATAGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGCCAGCGTTCAGCAGTGGCGTTTTGATGGCAACGCAAAGAAATGTGATCGGCAGTCAGCGGAGTGatat                    >  1:725131/1‑149 (MQ=255)
                                                                                                                              tGCAATTTAAACAGACAAATAGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGCCAGCGTTCAGCAGTGGCGTTTTGATGGCAACGCAAAGAAATGTGATCGGCAGTCAGCGGAGTGATAtt                   >  2:3223741/1‑149 (MQ=255)
                                                                                                                                  aTTTAAACAGACAAATAGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGCCAGCGTTCAGCAGTGGCGTTTTGATGGCAACGCAAAGAAATGTGATCGGCAGTCAGCGGAGTGATATTCAAc               >  1:2005008/1‑149 (MQ=255)
                                                                                                                                   tttAAACAGACAAATAGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGCCTGCGTTCAGCAGTGGCGTTTTGATGGCAACGCACAGAAATGTGATCGGCAGTCAGCGGAGTGATATTCAAcc              >  1:3087071/1‑149 (MQ=255)
                                                                                                                                   tttAAACAGACAAATAGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGCCAGCGTTCAGCAGTGGCGTTTTGATGGCAACGCAAAGAAATGTGATCGGCAGTCAGCGGAGTGATATTCAAcc              >  1:3508559/1‑149 (MQ=255)
                                                                                                                                   tttAAACAGACAAATAGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGCCAGCGTTCAGCAGTGGCGTTTTGATGGCAACGCAAAGAAATGTGATCGGCAGTCAGCGGAGTGATATTCAAcc              >  2:816744/1‑149 (MQ=255)
                                                                                                                                    ttAAACAGACAAATAGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGCCAGCGTTCAGCAGTGGCGTTTTGATGGCAACGCAAAGAAATGTGATCGGCAGTCAGCGGAGTGATATTCAAccc             <  1:3841159/149‑1 (MQ=255)
                                                                                                                                       aaCAGACAAATAGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGCCAGCGTTCAGCAGTGGCGTTTTGATGGCAACGCAAAGAAATGTGATCGGCAGTCAGCGGAGTGATATTCAACCCCAg          >  2:310706/1‑149 (MQ=255)
                                                                                                                                               aaTAGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGCCAGCGTTCAGCAGTGGCGTTTTGATGGCAACGCAAAGAAATGTGATCGGCAGTCAGCGGAGTGATATTCAACCCCAGACGGCggg  >  2:3739679/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
AGCGTTTCTGGATCACCCGGCAGAACGATTTCAATATGTTCCCAACCTTCGTGTGGGTAACGTTTTTCCCCGGGCCACGGTAGCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATCGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGCCAGCGTTCAGCAGTGGCGTTTTGATGGCAACGCAAAGAAATGTGATCGGCAGTCAGCGGAGTGATATTCAACCCCAGACGGCGGG  >  CP009273/1953730‑1954021

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: