Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,183,054 C→A 16.4% A151S (GCG→TCG)  thiF ← adenylyltransferase, modifies ThiS C‑terminus

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,183,0540CA16.4% 1056.9 / 142.7 345A151S (GCG→TCG) thiFadenylyltransferase, modifies ThiS C‑terminus
Reads supporting (aligned to +/- strand):  ref base C (148/140);  new base A (28/29);  total (176/169)
Fisher's exact test for biased strand distribution p-value = 7.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 357 total aligned reads displayed.

CAAAGTGCCCATAACCCCGACCACCGGGCCAACCACGCCCGCCGTGCGGCAGTTGCGTTCTGGCTCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTCCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCT  >  CP009273/4182905‑4183202
                                                                                                                                                     |                                                                                                                                                    
cAAAGTGCCCATAACCCCGACCACCGGGCCAACCACGCCCGCCGTGCGGCAGTTGCGTTCTGGCTCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACcgc                                                                                                                                                      >  2:4294375/1‑150 (MQ=255)
  aaGTGCCCATAACCCCGACCACCGGGCCAACCACGCCCGCCGTGCGGCAGTTGCGTTCTGGCTCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCgcgc                                                                                                                                                    >  2:3469475/1‑150 (MQ=255)
   aGTGCCCATAACCCCGACCACCGGGCCAACCACGCCCGCCGTGCGGCAGTTGCGTTCTGGCTCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCt                                                                                                                                                   >  2:1898003/1‑150 (MQ=255)
          aTAACCCCGACCACCGGGCCAACCACGCCCGCCGTGCGGCAGTTGCGTTCTGGCTCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTg                                                                                                                                            >  2:3222734/1‑150 (MQ=255)
          aTAACCCCGACCACCGGGCCAACCACGCCCGCCGTGCGGCAGTTGCGTTCTGGCTCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTg                                                                                                                                            >  1:380778/1‑150 (MQ=255)
           tAACCCCGACCACCGGGCCAACCACGCCCGCCGTGCGGCAGTTGCGTTCTGGCTCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGAGCTGGCGGTGa                                                                                                                                           >  1:2219166/1‑150 (MQ=255)
            aaCCCCGACCACCGGGCCAACCACGCCCGCCGTGCGGCAGTTGCGTTCTGGCTCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGAGCTGGCGGTGAt                                                                                                                                          <  1:4065112/150‑1 (MQ=255)
             aCCCCGACCACCGGGCCAACCACGCCCGCCGTGCGGCAGTTGCGTTCTGGCTCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCGGGCGGTGATa                                                                                                                                         >  2:3989095/1‑150 (MQ=255)
               cccGACCACCGGGCCAACCACGCCCGCCGTGCGGCAGTTGCGTTCTGGCTCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAg                                                                                                                                       <  2:3726283/150‑1 (MQ=255)
                  gACCACCGGGCCAACCACGCCCGCCGTGCGGCAGTTGCGTTCTGGCTCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAgcgg                                                                                                                                    >  2:2708536/1‑150 (MQ=255)
                    ccaccGGGCCAACCACGCCCGCCGTGCGGCAGTTGCGTTCTGGCTCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAgcggcg                                                                                                                                  >  1:2993483/1‑150 (MQ=255)
                     caccGGGCCAACCACGCCCGCCGTGCGGCAGTTGCGTTCTGGCTCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCgt                                                                                                                                 <  1:2680569/150‑1 (MQ=255)
                     caccGGGCCAACCACGCCCGCCGTGCGGCAGTTGCGTTCTGGCTCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCgt                                                                                                                                 <  1:2425301/150‑1 (MQ=255)
                     caccGGGCCAACCACGCCCGCCGTGCGGCAGTTGCGTTCTGGCTCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCgt                                                                                                                                 <  2:100330/150‑1 (MQ=255)
                      accGGGCCAACCACGCCCGCCGTGCGGCAGTTGCGTTCTGGCTCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCgtg                                                                                                                                <  1:3342503/150‑1 (MQ=255)
                       ccGGGCCAACCACGCCCGCCGTGCGGCAGTTGCGTTCTGGCTCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACGGACCGCCAAATCCGACCGCGCGGGCGGTGATAAGCGGCgtgt                                                                                                                               >  2:585070/1‑150 (MQ=255)
                          ggCCAACCACGCCCGCCGTGCGGCAGTTGCGTTCTGGCTCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGa                                                                                                                            >  2:3591619/1‑150 (MQ=255)
                          ggCCAACCACGCCCGCCGTGCGGCAGTTGCGTTCTGGCTCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGa                                                                                                                            <  1:924493/150‑1 (MQ=255)
                           gCCAACCACGCCCGCCGTGCGGCAGTTGCGTTCTGGCTCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAg                                                                                                                           <  1:1245662/150‑1 (MQ=255)
                            ccAACCACGCCCGCCGTGCGGCAGTTGCGTTCTGGCTCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGt                                                                                                                          <  2:3651134/150‑1 (MQ=255)
                             cAACCACGCCCGCCGTGCGGCAGTTGCGTTCTGGCTCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTg                                                                                                                         <  1:4087077/150‑1 (MQ=255)
                                 cACGCCCGCCGTGCGGCAGTTGCGTTCTGGCTCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCAc                                                                                                                     >  1:1948359/1‑150 (MQ=255)
                                    gcccgccGTGCGGCAGTTGCGTTCTGGCTCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCa                                                                                                                  >  2:4089748/1‑150 (MQ=255)
                                     cccgccGTGCGGCAGTTGCGTTCTGGCTCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGAGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAg                                                                                                                 >  1:1445487/1‑150 (MQ=255)
                                       cgccgTGCGGCAGTTGCTATCTGGCTCCTGGTTATCAGGCCACAGGCATCGGAAACACCCATGCTCACAGGGCGGCGTCACTACCATCAACTGACCGCCAAATCAGACCGAGCTGGCGGTCATAAGCGGCGTGTTGAGTGCCACGCAggc                                                                                                               >  2:2391707/1‑150 (MQ=255)
                                        gccgTGCGGCAGTTGCGTTCTGGCTCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAggcg                                                                                                              <  1:1640172/150‑1 (MQ=255)
                                            tGCGGCAGTTGCGTTCTGGCTCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCAt                                                                                                          <  2:1786853/150‑1 (MQ=255)
                                               ggCAGGTGCGTTCTGGCTCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTaa                                                                                                       <  2:2819705/150‑1 (MQ=255)
                                                  aGTTGCGTTCTGGCTCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAAtct                                                                                                    <  2:2499688/150‑1 (MQ=255)
                                                        gTTCTGGCTCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGAGCCACGCAGGCGGCATTAATCTCCTGGc                                                                                              <  2:2261136/150‑1 (MQ=255)
                                                          tCTGGCTCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTCCTGGCGa                                                                                            <  2:4333601/150‑1 (MQ=255)
                                                           cTGGCTCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGAGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTCCTGGCGAg                                                                                           <  1:2981563/150‑1 (MQ=255)
                                                            tGGCTCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTCCTGGCGAGt                                                                                          >  2:3893151/1‑150 (MQ=255)
                                                               cTCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTCCTGGCGAGTCGc                                                                                       >  2:894049/1‑150 (MQ=255)
                                                                     gTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTCCTGGCGAGTCGCCATAtt                                                                                 <  1:2110028/150‑1 (MQ=255)
                                                                      ttATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTCCTGGCGAGTCGCCATATTg                                                                                >  2:3097708/1‑150 (MQ=255)
                                                                       tATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTCCTGGCGAGTCGCCATATTGt                                                                               >  2:2014170/1‑150 (MQ=255)
                                                                       tATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTCCTGGCGAGTCGCCATATTGt                                                                               <  1:2734385/150‑1 (MQ=255)
                                                                        aTCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTCCTGGCGAGTCGCCATATTGTc                                                                              >  2:2532551/1‑150 (MQ=255)
                                                                           tGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTCCTGGCGAGTCGCCATATTGTCGGt                                                                           >  1:2325069/1‑150 (MQ=255)
                                                                             gCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTCCTGGCGAGTCGCCATATTGTCGGTAc                                                                         <  1:94300/150‑1 (MQ=255)
                                                                              ccACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTCCTGGCGAGTCGCCATATTGTCGGTACa                                                                        <  2:1079045/150‑1 (MQ=255)
                                                                              ccACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGAGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTCCTGGCGAGTCGCCATATTGTCGGTACa                                                                        >  2:765776/1‑150 (MQ=255)
                                                                                 caGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTCCTGGCGAGTCGCCATATTGTCGGTACAGTc                                                                     >  2:2048152/1‑150 (MQ=255)
                                                                                   ggCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTCCTGGCGAGTCGCCATATTGTCGGTACAGTCGa                                                                   >  1:3942591/1‑150 (MQ=255)
                                                                                     cAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTCCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGc                                                                 <  1:567910/150‑1 (MQ=255)
                                                                                      aGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTCCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGca                                                                <  1:2545209/150‑1 (MQ=255)
                                                                                      aGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTCCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGca                                                                >  1:3529350/1‑150 (MQ=255)
                                                                                      aGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGAGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTCCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGca                                                                >  1:3471624/1‑150 (MQ=255)
                                                                                        cGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTCCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGcacc                                                              <  1:4164964/150‑1 (MQ=255)
                                                                                        cGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTCCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGcacc                                                              <  2:217457/150‑1 (MQ=255)
                                                                                          gTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGAGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTCCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGcaccac                                                            >  1:1967488/1‑150 (MQ=255)
                                                                                            aaCACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATCAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTCCTGGCGAGTCGCCCTAGTGTCGGTACAGTCGAGCACCACAt                                                          >  2:1685404/1‑150 (MQ=255)
                                                                                            aaCACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTCCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACAt                                                          >  2:149459/1‑150 (MQ=255)
                                                                                              cacCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTCCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCg                                                        <  2:2985753/150‑1 (MQ=255)
                                                                                              cacCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTCCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCg                                                        >  2:1605381/1‑150 (MQ=255)
                                                                                                  ccTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTCCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGccc                                                    <  2:4152184/150‑1 (MQ=255)
                                                                                                   cTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTCCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCg                                                   <  1:3985174/150‑1 (MQ=255)
                                                                                                    tGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTCCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGt                                                  >  1:449158/1‑150 (MQ=255)
                                                                                                      cTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTCCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGc                                                >  2:699094/1‑150 (MQ=255)
                                                                                                        cccAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTCCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCaa                                              >  2:2342981/1‑150 (MQ=255)
                                                                                                          cAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTCCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAAcc                                            >  2:4228572/1‑150 (MQ=255)
                                                                                                          cAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTCCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAAcc                                            >  1:510180/1‑150 (MQ=255)
                                                                                                             ggcggcGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTCCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCa                                         >  1:137901/1‑150 (MQ=255)
                                                                                                               cggcgTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTCCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATc                                       >  1:3703232/1‑150 (MQ=255)
                                                                                                                ggcgTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTCCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCt                                      >  1:3214663/1‑150 (MQ=255)
                                                                                                                ggcgTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTCCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCt                                      <  2:1093135/150‑1 (MQ=255)
                                                                                                                ggcgTCAGTACCATCAACTGACCGCCAAATCCGACCGAGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTCCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCt                                      >  2:191824/1‑150 (MQ=255)
                                                                                                                 gcgTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTCCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCtt                                     <  1:684600/150‑1 (MQ=255)
                                                                                                                  cgTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTCCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCttt                                    <  1:4309371/150‑1 (MQ=255)
                                                                                                                   gTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTCCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCtttt                                   <  1:4262286/150‑1 (MQ=255)
                                                                                                                   gTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTCCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCtttt                                   >  2:3411389/1‑150 (MQ=255)
                                                                                                                   gTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTCCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCtttt                                   >  2:3022144/1‑150 (MQ=255)
                                                                                                                   gTCAGTACCATCAACTGACCGCCAAATCCGACCGAGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTCCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCtttt                                   >  1:1606820/1‑150 (MQ=255)
                                                                                                                   gTCAGTACCATCAACTGACCGCCAAATCCGACCGAGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTCCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCtttt                                   <  1:984459/150‑1 (MQ=255)
                                                                                                                      aGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTCCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAAc                                >  1:1416512/1‑150 (MQ=255)
                                                                                                                       gTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTCCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACg                               <  2:2286441/150‑1 (MQ=255)
                                                                                                                        tACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTCCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGc                              >  2:1308168/1‑150 (MQ=255)
                                                                                                                           cATCAACTGACCGCCAAATCCGACCGAGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTCCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTc                           >  1:2005678/1‑150 (MQ=255)
                                                                                                                             tCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTCCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCAc                         >  2:212045/1‑150 (MQ=255)
                                                                                                                             tCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTCCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCAc                         >  2:182936/1‑150 (MQ=255)
                                                                                                                               aaCTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTCCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCAccc                       >  1:3897386/1‑150 (MQ=255)
                                                                                                                                aCTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTCCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCg                      <  2:2651115/150‑1 (MQ=255)
                                                                                                                                 cTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTCCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGt                     >  2:3367220/1‑150 (MQ=255)
                                                                                                                                  tGACCGCCAAATCCGACCGAGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTCCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGtt                    >  2:613537/1‑150 (MQ=255)
                                                                                                                                   gACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTCCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTa                   <  1:3097708/150‑1 (MQ=255)
                                                                                                                                   gACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTCCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTa                   <  2:1496952/150‑1 (MQ=255)
                                                                                                                                     ccGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTCCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAAc                 <  1:1013723/150‑1 (MQ=255)
                                                                                                                                       gCCAAATCCGACCGAGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTCCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCg               >  1:4245616/1‑150 (MQ=255)
                                                                                                                                       gCCAAATCCGACCGAGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTCCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCg               <  2:828389/150‑1 (MQ=255)
                                                                                                                                        ccAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTCCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCgt              >  1:3552936/1‑150 (MQ=255)
                                                                                                                                         cAAATCCGACCGAGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTCCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCgtt             <  1:28494/150‑1 (MQ=255)
                                                                                                                                          aaaTCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTCCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCgttg            >  2:447665/1‑150 (MQ=255)
                                                                                                                                           aaTCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGATTGCCACGCAGGCGGCATTAATATCCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCgttgt           >  1:1109136/1‑150 (MQ=255)
                                                                                                                                             tCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTCCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGttgttg         <  2:3833375/150‑1 (MQ=255)
                                                                                                                                              ccgaccgaGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTCCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTtgttgt        <  2:1385392/150‑1 (MQ=255)
                                                                                                                                                  ccGCGGTGGCGGTGATAAGCGGCGTTTTGAGTGCCACGCAGGCGGCATTAATCTCCTGGCGAGTCGCAATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATg    <  1:3469345/150‑1 (MQ=255)
                                                                                                                                                   cgcgCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTCCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGc   <  2:1200968/150‑1 (MQ=255)
                                                                                                                                                   cgcgCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTCCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGc   <  1:1816561/150‑1 (MQ=255)
                                                                                                                                                    gcgcTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTCCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCt  <  2:3044868/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CAAAGTGCCCATAACCCCGACCACCGGGCCAACCACGCCCGCCGTGCGGCAGTTGCGTTCTGGCTCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTCCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCT  >  CP009273/4182905‑4183202

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: