Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,139,910 G→T 92.1% R1055L (CGA→CTA)  yegE → putative diguanylate cyclase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,139,9100GT92.1% 989.9 / 68.5 305R1055L (CGA→CTA) yegEputative diguanylate cyclase
Reads supporting (aligned to +/- strand):  ref base G (9/15);  new base T (137/144);  total (146/159)
Fisher's exact test for biased strand distribution p-value = 3.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.34e-01

Only 100 of 306 total aligned reads displayed.

GCCAGGTGGGCCGCGATCTGCAAATCTTCAACTCGCTGAAAGCGAATATGGCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCGACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTCGATCTGATTTATGGTGAGGTGATTGCCGATGCCCAACCGCTGGATTTGCTGGTGAATAGTAGTTATTTCGCGATT  >  CP009273/2139763‑2140058
                                                                                                                                                   |                                                                                                                                                    
 ccAGGTGGGCCGCGATCTGCAAATCTTCAACTCGCTGAAAGCGAATATGGCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTAc                                                                                                                                                    >  1:4289392/1‑149 (MQ=255)
    ggTGGGCCGCGATCTGCAAATCTTCAACTCGCTGAAAGCGAATATGGCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCg                                                                                                                                                 <  1:5062988/149‑1 (MQ=255)
    ggTGGGCCGCGATCTGCAAATCTTCAACTCGCTGAAAGCGAATATGGCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCg                                                                                                                                                 <  1:3113711/149‑1 (MQ=255)
          ccGCGATCTGCAAATCTTCAACTCGCTGAAAGCGAATATGGCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATgg                                                                                                                                           >  1:528710/1‑148 (MQ=255)
          ccGCGATCTGCAAATCTTCAACTCGCTGAAAGCGAATATGGCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGa                                                                                                                                           <  2:4256993/149‑1 (MQ=255)
           cgcgATCTGCAAATCTTCAACTCGCTGAAAGCGAATATGGCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGaa                                                                                                                                          <  1:1669451/149‑1 (MQ=255)
            gcgATCTGCAAATCTTCAACTCGCTGAAAGCGAATATGGCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGaaa                                                                                                                                         >  2:1836328/1‑149 (MQ=255)
                tCTGCAAATCTTCAACTCGCTGAAAGCGAATATGGCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCGACTCGGGATGAAAACCa                                                                                                                                     >  1:769368/1‑149 (MQ=255)
                  tGCAAATCTTCAACTCGCTGAAAGCGAATATGGCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATc                                                                                                                                   >  2:2054545/1‑149 (MQ=255)
                  tGCAAATCTTCAACTCGCTGAAAGCGAATATGGCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATc                                                                                                                                   <  1:5006335/149‑1 (MQ=255)
                  tGCAAATCTTCAACTCGCTGAAAGCGAATATGGCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATc                                                                                                                                   >  2:4670291/1‑149 (MQ=255)
                   gCAAATCTTCAACTCGCTGAAAGCGAATATGGCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCg                                                                                                                                  >  1:2782015/1‑149 (MQ=255)
                     aaaTCTTCAACTCGCTGAAAGCGAATATGGCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGcc                                                                                                                                >  2:597629/1‑149 (MQ=255)
                      aaTCTTCAACTCGCTGAAAGCGAATATGGCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCg                                                                                                                               <  2:3033136/149‑1 (MQ=255)
                          ttCAACTCGCTGAAAGCGAATATGGCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGcc                                                                                                                           <  1:17991/149‑1 (MQ=255)
                           tCAACTCGCTGAAAGCGAATATGGCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCa                                                                                                                          >  2:270071/1‑149 (MQ=255)
                            cAACTCGCTGAAAGCGAATATGGCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAg                                                                                                                         >  2:3125699/1‑149 (MQ=255)
                             aaCTCGCTGAAAGCGAATATGGCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGt                                                                                                                        >  1:803134/1‑149 (MQ=255)
                              aCTCGCTGAAAGCGAATATGGCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTc                                                                                                                       <  1:2447496/149‑1 (MQ=255)
                               cTCGCTGAAAGCGAATATGGCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCg                                                                                                                      >  2:4014920/1‑149 (MQ=255)
                                 cGCCGAAAGCGAATATGGCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGtt                                                                                                                    <  1:4121543/149‑1 (MQ=255)
                                    tGAAAGCGAATATGGCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTa                                                                                                                 >  2:3769076/1‑149 (MQ=255)
                                     gAAAGCGAATATGGCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTAc                                                                                                                >  1:2329043/1‑149 (MQ=255)
                                       aaGCGAATATGGCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTAccc                                                                                                              >  1:5019109/1‑149 (MQ=255)
                                          cGAATATGGCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTGTTACCCGTa                                                                                                           >  1:596099/1‑149 (MQ=255)
                                            aaTATGGCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGt                                                                                                         <  2:2049352/149‑1 (MQ=255)
                                             atatGGCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTg                                                                                                        >  2:4354380/1‑149 (MQ=255)
                                               atGGCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTATGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGAt                                                                                                      >  2:1612474/1‑149 (MQ=255)
                                               atGGCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGTACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCATAGTGAt                                                                                                      >  1:883787/1‑149 (MQ=255)
                                               atGGCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGAt                                                                                                      <  1:3766750/149‑1 (MQ=255)
                                               atGGCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGAt                                                                                                      >  1:2933357/1‑149 (MQ=255)
                                                tGGCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATg                                                                                                     <  2:3856604/149‑1 (MQ=255)
                                                tGGCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTGTTACCCTTAGTGATg                                                                                                     >  1:1313276/1‑149 (MQ=255)
                                                  gCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGt                                                                                                                        <  2:2890233/128‑1 (MQ=255)
                                                  gCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGa                                                                                                   >  2:3935109/1‑149 (MQ=255)
                                                  gCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGa                                                                                                   <  1:3043224/149‑1 (MQ=255)
                                                  gCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGa                                                                                                   >  2:804821/1‑149 (MQ=255)
                                                  gCAGATTACCTGCTACTTGATGATGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGa                                                                                                   >  2:4691041/1‑149 (MQ=255)
                                                     gggttccTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATAc                                                                                                <  1:4475761/144‑1 (MQ=255)
                                                     gATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTATGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATAc                                                                                                >  2:1493940/1‑149 (MQ=255)
                                                       ttACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGc                                                                                              <  2:941859/149‑1 (MQ=255)
                                                        tACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCt                                                                                             <  2:125117/149‑1 (MQ=255)
                                                        tACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCGACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCt                                                                                             >  1:4245477/1‑149 (MQ=255)
                                                          ccTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCttt                                                                                           <  2:3162176/149‑1 (MQ=255)
                                                                  ttGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATc                                                                                   <  2:3364041/149‑1 (MQ=255)
                                                                  ttGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATc                                                                                   <  1:1172150/149‑1 (MQ=255)
                                                                      tGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTATGGCATCGGCg                                                                               >  1:202411/1‑149 (MQ=255)
                                                                             ttATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTCGATCt                                                                        <  1:2121057/149‑1 (MQ=255)
                                                                               aTGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTCGATCTGa                                                                      <  2:4977503/149‑1 (MQ=255)
                                                                                tGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTCGATCTGAt                                                                     <  2:2314897/149‑1 (MQ=255)
                                                                                tGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTCGATCTGAt                                                                     <  1:3917126/149‑1 (MQ=255)
                                                                                tGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCGACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTCGATCTGAt                                                                     <  2:3093883/149‑1 (MQ=255)
                                                                                 gcgcCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTCGATCTGAtt                                                                    <  2:5068844/149‑1 (MQ=255)
                                                                                      aaCGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTCGATCTGATTTATgg                                                               <  2:1050340/149‑1 (MQ=255)
                                                                                      aaCGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTCGATCTGATTTATgg                                                               <  1:830246/149‑1 (MQ=255)
                                                                                       aCGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTCGATCTGATTTATggt                                                              >  1:4689262/1‑149 (MQ=255)
                                                                                       aCGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTCGATCTGATTTATggt                                                              <  2:4074951/149‑1 (MQ=255)
                                                                                          tGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCGTAGTGATGGATACGCGTTCTGGCATCGGCGTCGCTCGGATTTATggtgag                                                           >  1:1410513/1‑149 (MQ=255)
                                                                                           gCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTCGATCTGATTTATggtgagg                                                          <  1:3004429/149‑1 (MQ=255)
                                                                                            cAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTCGATCTGATTTATggtgaggt                                                         <  2:1980084/149‑1 (MQ=255)
                                                                                               ggTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTCGATCTGATTTATGGTGAGGTGAt                                                      <  2:70085/149‑1 (MQ=255)
                                                                                                gTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTCGATCTGATTTATGGTGAGGTGAtt                                                     <  1:1744334/149‑1 (MQ=255)
                                                                                                gTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCGTAGTGATGGATACGCTTTCTGGCATCGGCGTCGATCTGATTTATGGTGAGGTGAtt                                                     >  2:4506624/1‑149 (MQ=255)
                                                                                                 tAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTCGATCTGATTTATGGTGAGGTGATTg                                                    <  1:2404408/149‑1 (MQ=255)
                                                                                                  aaTTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTCGATCTGATTTATGGTGAGGTGATTGc                                                   <  2:428220/149‑1 (MQ=255)
                                                                                                   aTTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTCGATCTGATTTATGGTGAGGTGATTGcc                                                  <  1:4840888/149‑1 (MQ=255)
                                                                                                   aTTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATAAGCTTTATTGCATCGGCGTCGATCTGATTTATGGTGAGGTGATTGcc                                                  >  1:1119855/1‑149 (MQ=255)
                                                                                                   aTTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCGACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTCGATCTGATTTATGGTGAGGTGATTGcc                                                  <  2:189912/149‑1 (MQ=255)
                                                                                                      tGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTCGATCTGATTTATGGTGAGGTGATTGCCGAt                                               <  2:1166368/149‑1 (MQ=255)
                                                                                                      tGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTCGATCTGATTTATGGTGAGGTGATTGCCGAt                                               <  2:3075977/149‑1 (MQ=255)
                                                                                                       gatggatgAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTCGATCTGATTTATGGTGAGGTGATTTCCGATg                                              >  1:5165717/1‑149 (MQ=255)
                                                                                                         tggatgAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTATGGCAACGGCGTCGATCGGATTTAAGGTGAGGTGATTGCCGATGcc                                            >  1:1645097/1‑149 (MQ=255)
                                                                                                               agaTGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTCGATCTGATTTATGGTGAGGTGATTTCCGATGCCCAACCg                                      >  2:3801012/1‑149 (MQ=255)
                                                                                                               agaTGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTAACCTTAGTGTTTGATAAGCTTTCTGGCATCGGCGTCGATCTGTTTTATTGTGTGGTGATTTCCGATTCCCAAccc                                      >  1:2706038/1‑148 (MQ=255)
                                                                                                               agaTGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTGTTAACCTTAGTGATGGATACGCTTTCTGGCATCGGCGTCGATCGGATTTATGGTGAGGTGATTTCCGATGCCCCAccc                                      >  1:1008968/1‑148 (MQ=255)
                                                                                                                gaTGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTGTTACCCTTAGTGATGGATACGCTTTCATGCATCGGCGTCGATCTGATTTATTGTGAGGTGATTGCCGATGCCCAACCGc                                     >  1:1056922/1‑149 (MQ=255)
                                                                                                                    cTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTCGATCTGATTTATGGTGAGGTGATTGCCGATGCCCAACCGCTGGa                                 <  1:912747/149‑1 (MQ=255)
                                                                                                                     tGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTCGATCTGATTTATGGTGAGGTGATTGCCGATGCCCAACCGCTGGAt                                <  2:4842463/149‑1 (MQ=255)
                                                                                                                     tGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTCGATCTGATTTATGGTGAGGTGATTGCCGATGCCCAACCGCTGGAt                                <  1:66564/149‑1 (MQ=255)
                                                                                                                     tGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTCGATCTGATTTATGGTGAGGTGATTGCCGATGCCCAACCGCTGGAt                                <  1:66643/149‑1 (MQ=255)
                                                                                                                     tGATTACGATTATTCAGGGGCACGCTCAGCGACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTCGATCTGATTTATGGTGAGGTGATTGCCGATGCCCAACCGCTGGAt                                <  2:261070/149‑1 (MQ=255)
                                                                                                                      gATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTCGATCTGATTTATGGTGAGGTGATTGCCGATGCCCAACCGCTGGAtt                               <  1:3975599/149‑1 (MQ=255)
                                                                                                                        ttACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTCGATCTGATTTATGGTGAGGTGATTGCCGATGCCCAACCGCTGGATTTg                             >  2:547292/1‑149 (MQ=255)
                                                                                                                         tACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTGTCTGGCATCGGCGTCGATCTGATTTATGGTGAGGTGATTGCCGATGCCCAACCGCTGGATTTGc                            <  1:3643719/149‑1 (MQ=255)
                                                                                                                          aCGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTCGATCTGATTTATGGTGAGGTGATTGCCGATGCCCAACCGCTGGATTTGCt                           <  1:4919656/149‑1 (MQ=255)
                                                                                                                          aCGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTCGATCTGATTTATGGTGAGGTGATTGCCGATGCCCAACCGCTGGATTTGCt                           <  1:4972243/149‑1 (MQ=255)
                                                                                                                           cGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTCGATCTGATTTATGGTGAGGTGATTGCCGATGCCCAACCGCTGGATTTGCTg                          <  2:3652072/149‑1 (MQ=255)
                                                                                                                           cGATTATTCAGGGGCACGCTCAGCGACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTCGATCTGATTTATGGGGAGGTGATTGCCGATTCCCAACCGCAGGATTTGCTg                          >  2:8489/1‑149 (MQ=255)
                                                                                                                              ttattCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCTGGCCAGTCGTTTAAACATTAGTGTTGGATAAGATTTATTGAATAGGCGTCGGAGTGTTTTATTGGTAGGTTATTTGCGATGACCAACCCATTGATTTGCTtgtg                       >  1:1497340/1‑149 (MQ=255)
                                                                                                                                attCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGCTACGCTTTCTGGCATCGGCGTCGATCTGGTTTATGGTGTGGTGATTGCCGATGCCCAACCGCTGGATTTGCTGGTgga                     >  2:1366156/1‑147 (MQ=255)
                                                                                                                                 ttCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTCGATCTGATTTATGGTGAGGTGATTGCCGATGCCCAACCGCTGGATTTGCTGGTGAAt                    >  2:1207992/1‑149 (MQ=255)
                                                                                                                                       ggCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTCGATCTGATTTATGGTGAGGTGATTGCCGATGCCCAACCGCTGGATTTGCTGGTGAATagtagt              <  1:2209219/149‑1 (MQ=255)
                                                                                                                                        gCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTCGATCTGATTTATGGTGAGGTGATTGCCGATGCCCAACCGCTGGATTTGCTGGTGAATAGTAGtt             <  2:1347501/149‑1 (MQ=255)
                                                                                                                                          aCGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTCGATCTGAGTTATGGTGAGGTGATTGCCGATGCCCAACCGCTGGATTTGCTGGTGAATAGTAGTTAt           >  2:237954/1‑149 (MQ=255)
                                                                                                                                            gCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTCGATCTGATTTATGGTGAGGTGATTGCCGATGCCCAACCGCTGGATTTGCTGGTGAATAGTAGTTAttt         <  1:340987/149‑1 (MQ=255)
                                                                                                                                             cTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTCGATCTGATTTATGGTGAGGTGATTGCCGATGCCCAACCGCTGGATTTGCTGGTGAATAGTAGTTATTTc        <  2:3592072/149‑1 (MQ=255)
                                                                                                                                             cTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTCGATCTGATTTATGGTGAGGTGATTGCCGATGCCCAACCGCTGGATTTGCTGGTGAATAGTAGTTATTTc        <  2:530406/149‑1 (MQ=255)
                                                                                                                                             cTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTCGATCTGATTTATGGTGAGGTGATTGCCGATGCCCAACCGCTGGATTTGCTGGTGAATAGTAGTTATTTc        <  2:1688534/149‑1 (MQ=255)
                                                                                                                                               cAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTTACCTTAGTTATGGATAAGCTTTCTGGCATCGGGGGCGATGGGTTTTTTTGTGGGGGGGTTTGCGGTGCCCAACCCCGGGAGTTGCTGGTGGATAGTATTTATTTagc      >  1:2638310/1‑149 (MQ=255)
                                                                                                                                                   gACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGCTGGATACGCGTTCTGGCATCGGCGTCGCTCTGATTTATGGGGCGGTGCTTGCCGATGCCCAACCCCGGGCTTTGCTGGGGAATAGTAGTTATTTCGCGagt  >  2:2251358/1‑147 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
GCCAGGTGGGCCGCGATCTGCAAATCTTCAACTCGCTGAAAGCGAATATGGCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCGACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTCGATCTGATTTATGGTGAGGTGATTGCCGATGCCCAACCGCTGGATTTGCTGGTGAATAGTAGTTATTTCGCGATT  >  CP009273/2139763‑2140058

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: