Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,661,108 C→G 5.5% G33A (GGG→GCG)  hcaT ← putative 3‑phenylpropionic transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,661,1080CG 5.5% 214.3 / 6.0 73G33A (GGG→GCG) hcaTputative 3‑phenylpropionic transporter
Reads supporting (aligned to +/- strand):  ref base C (36/33);  new base G (2/2);  total (38/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.28e-01

GAGAAGTGTCAGCAGTGCCAGCACGCGCAAGGCGGAAATCAGGCGGGAAGGATCGCTGACGCGGGGCGCGATGAGCAAACTCCCGAGGAAACGGGCAACCAGACCTGCCCCCAATAACAGGCCGATGGTTTCTGGCGTTAAACCAATCCCTTTAAGCCAGACGCTCCAGAAAGGTAGAAAAATGCCGTAACTAAAAAAGTATGTGAAATAGCCGAGCGCCAACCAGCGCGTGGATTGCAAAACCATGAGTCCCTCCCGTCAGGAGGCGTTAGTCTGGGGTTAATCGAATGTTGG  >  CP009273/2660960‑2661253
                                                                                                                                                    |                                                                                                                                                 
gagaAGTGTCAGCAGTGCCAGCACGCGCAAGGCGGAAATCAGGCGGGAAGGATCGCTGACGCGGGGCGCGATGAGCAAACTCCCGGGGAAACGGGCAACCAGACCTGCCCCCAATAACAGGCCGATGGTTTCTGGCGTTAAACCAATcc                                                                                                                                                   <  1:725273/149‑1 (MQ=255)
 agaAGTGTCAGCAGTGCCAGCACGCGCAAGGCGGAAATCAGGCGGGGAGGATCGCTGACGCGGGGCGCGATGAGCAAACTCCCGAGGAAACGGGCAACCAGACCTGCCCCCAATAACAGGCCGATGGTTTCTGGCGTTAAACCAATccc                                                                                                                                                  >  2:201442/1‑149 (MQ=255)
        tCAGCAGTGCCAGCACGCGCAAGGCGGAAATCAGGCGGGAAGGATCGCTGACGCGGGGCGCGATGAGCAAACTCCCGAGGAAACGGGCAACCAGACCTGCCCCCAATAACAGGCCGCTGGTTTCTGGCGTTAAACCAATCCCTTTAAGc                                                                                                                                           >  1:1099358/1‑149 (MQ=255)
           gcagGGCCAGCACGCGCAAGGCGGAAATCAGGCGGGAAGGATCGCTGACGCGGGGCGCGATGAGCAAACTCCCGAGGAAACGGGCAACCAGACCTGCCCCCAATAACAGGCCGATGGTTTCTGGCGTTAAACCAATCCCTTTAAGCCAg                                                                                                                                        <  1:208752/149‑1 (MQ=255)
            cagTGCCAGCACGCGCAAGGCGGAAATCAGGCGGGGAGGATCGCTGACGCGGGGCGCGATGAGCAAACTCCCGAGGAAACGGGCAACCAGACCTGCCCCCAATAACAGGCCGATGGTTTCTGGCGGGAAACCAATCCCTTTAAGCCAGa                                                                                                                                       >  1:1130044/1‑149 (MQ=255)
            cagTGCCAGCACGCGCAAGGCGGAAATCAGGCGGGAAGGATCGCTGACGCGGGGCGCGATGAGCAAACTCCCGAGGAAACGGGCAACCAGACCTGCCCCCAATAACAGGCCGATGGTTTCTGGCGTTAAACCAATCCCTTTAAGCCAGa                                                                                                                                       <  1:554034/149‑1 (MQ=255)
               tGCCAGCACGCGCAAGGCGGAAATCAGGCGGGAAGGATCGCTGACGCGGGGCGCGATGAGCAAACTCCCGAGGAAACGGGCAACCAGACCTGCCCCCAATAACAGGCCGATGGTTTCTGGCGTTAAACCAATCCCTTTAAGCCAGACGc                                                                                                                                    >  1:963170/1‑149 (MQ=255)
                 ccAGCACGCGCAAGGCGGAAATCAGGCGGGAAGGATCGCTGACGCGGGGCGCGATGAGCAAACTCCCGAGGAAACGGGCAACCAGACCTGCCCCCAATAACAGGCCGATGGTTTCTGGCGTTAAACCAATCCCTTTAAGCCAGACGCTc                                                                                                                                  >  2:50455/1‑149 (MQ=255)
                 ccAGCACGCGCAAGGCGGAAATCAGGCGGGAAGGATCGCTGACGCGGGGCGCGATGAGCAAACTCCCGAGGAAACGGGCAACCAGACCTGCCCCCAATAACAGGCCGATGGTTTCTGGCGTTAAACCAATCCCTTTAAGCCAGACGCTc                                                                                                                                  >  1:322420/1‑149 (MQ=255)
                 ccAGCACGCGCAAGGCGGAAATCAGGCGGGAAGGATCGCTGACGCGGGGCGCGATGAGCAAACTCCCGAGGAAACGGGCAACCAGACCTGCCCCCAATAACAGGCCGATGGTTTCTGCCGTTAAACCAATCCCACTAAGGCGGAcccct                                                                                                                                  >  2:549601/1‑145 (MQ=255)
                   aGCCCGCGCAATGCGGAAATCAGGCGGGAAGGATCGCTGACGCGGGGCGCGATGAGCAAACTCCCGAGGAAACGGGCAACCAGACCTGCCCCCCATAACAGGCCGATGGTTTCTGGCGTTAAACCAACCCCTTTAAGCCAGACGCTCCa                                                                                                                                >  2:12898/1‑149 (MQ=255)
                   aGCACGCGCAAGGCGGAAATCAGGCGGGAAGGATCGCTGACGCGGGGCGCGATGAGCAAACTCCCGAGGAAACGGGCAACCAGACCTGCCCCCAATAACAGGCCGATGGTTTCTGGCGTTAAACCAATCCCTTTAAGCCAGACGCTCCa                                                                                                                                >  1:1192213/1‑149 (MQ=255)
                   aGCACGCGCAAGGCGGAAATCAGGCGGGAAGGATCGCTGACGCGGGGCGCGATGAGCAAACTCCCGAGGAAACGGGCAACCAGACCTGCCCCCAATAACAGGCCGATGGTTTCTGGCGTTAAACCAATCCCTTTAAGCCAGACGCTCCa                                                                                                                                >  1:445889/1‑149 (MQ=255)
                      aCGCGCAAGGCGGAAATCAGGCGGGAAGGATCGCTGACGCGGGGCGCGATGAGCAAACTCCCGAGGAAACGGGCAACCAGACCTGCCCCCAATAACAGGCCGATGGTTTCTGGCGTTAAACCAATCCCTTTAAGCCAGAAGCTCCCGaa                                                                                                                             >  1:500088/1‑149 (MQ=255)
                            gaGGCGGAAATCAGGCGGGAAGGATCGCTGACGGGGGGCGCGATGAGCAAACTCCCGGGGGAACGGGCAACCAGACCTGCCCCCAATAACAGGACGATGGTTTATGGCGTTAAGCCGCTACCTTTAAGCCACACGCTCCAGAAAGGTaa                                                                                                                       >  2:33299/2‑148 (MQ=255)
                            aaGGCGGAAATCAGGCGGGAAGGATCGCTGACGCGGGGCGCGATGAGCAAACTCCCGAGGAACCGGGCACCCCGACCTGCCCCCAATAACAGGCCGATGGGTTCTGGCGTTAAACCAATCCCTTTAAGCCGGAGACTCCGGAAAGGAAg                                                                                                                       >  2:1151378/1‑149 (MQ=255)
                            aaGGCGGAAATCAGGCGGGAAGGATCGCTGACGCGGGGCGCGATGAGCAAACTCCCGAGGAAACGGGCCACCAGACCTGCCCCCAATAACAGGCCGATGGTTTCTGGCGTTAAACCAATCCCTTTAAGCCAGACGCTCCAGAAAGGTAg                                                                                                                       >  2:979541/1‑149 (MQ=255)
                               gCGGAAATCAGGCGGGAAGGATCGCTGACGCGGGGCGCGATGAGCAAACTCCCGAGGAAACGGGCAACCAGACCTGCCCCCAATAACAGGCCGATGGTTTCTGGCGTTAAACCAATCCCTTTAAGCCAGACGCTCCAGAAAGGTAGaaa                                                                                                                    >  2:588932/1‑149 (MQ=255)
                                cGGAAATCAGGCGGGAAGGATCGCTGACGCGGGGCGCGATGAGCAAACTCCCGAGGAAACGGGCAACCAGACCTGCCCCCAATAACAGGCCGATGGTTTCTGGCGTTAAACCAATCGCTTTAAGCCAGACGCTCCAGAAAGGTAGaaaa                                                                                                                   >  1:1183866/1‑149 (MQ=255)
                                cGGAAATCAGGCGGGAAGGATCGCTGACGCGGGGCGCGATGAGCAAACTCCCGAGGAAACGGGCAACCAGACCTGCCCCCAATAACAGGCCGATGGTTTCTGGCGTTAAACCAATCGCTTTAAGCCAGACGCTCCAGAAAGGTAGaaaa                                                                                                                   >  2:708462/1‑149 (MQ=255)
                                cGGAAATCAGGCGGGAAGGATCGCTGACGCGGGGCGCGATGAGCAAACTCCCGAGGAAACGGGCAACCAGACCTGCCCCCAATAACAGGCCGATAGTTTCTGGCGTTAAACCAATCCCTTTAAGCCAGACGCGCCAGAAAGGGAGaaaa                                                                                                                   >  1:53792/1‑149 (MQ=255)
                                  gAAATCAGGCGGGAAGGATCGCGGACGCGGGGCGCGATGGGCAAACTCCCGAGGAAACGGGCAACCAGACCTGCCCCCAATAACAGGCCGATGGTTTCTGGCGTTAAACCAATCCCTTTAAGCCAGACGCTCCAGAAAGGTAGAAAAAt                                                                                                                 <  2:1007583/149‑1 (MQ=255)
                                             ggAAGGATCGCTGACGCGGGGCGCGATGAGCAAACTCCCGAGGAAACGGGCAACCAGACCTGCCCCCAATAACAGGCCGATGGTATCTGGCGTTAAACCAATCCCTTTAAGCCAGACGCTCCAGAAAGGTAGAAAAATGCCGTAACTaa                                                                                                      <  1:943746/149‑1 (MQ=255)
                                                      gCTGACGCGGGGCGCGATGAGCAAACTCCCGAGGAAACGGGCAACCAGACCTGCCCCCAATAACAGGCCGATGGTTTCTGGCGTTAAACCAATCCCTTTAAGCCAGACGCTCCAGAAAGGTAGAAAAATGCCGTAACTAAAAAAGTAtg                                                                                             <  2:338607/149‑1 (MQ=255)
                                                        tGACGCGGGGCGCGATGAGCAAACTCCCGAGGAAACGGGCAACCAGACCTGCCCCCAATCACAGGCCGATGGTTTCTGGCGTTAAACCAATCCCTTTAAGCCAGACGCTCCAGAAAGGTAGAAAAATGCCGTAACTCAAAAAGTAcgtg                                                                                           >  2:943567/1‑149 (MQ=255)
                                                              ggggCGCGATGAGCAAACTCCCGAGGAAACGGGCAACCAGACCTGCCCCCAATAACAGGCCGATGGTTTCTGGCGTTAAACCAATCCCTTTAAGCCAGACGCTCCAGAACGGTCGAAAAATGCCGTAACTAAAAAAGTATGTGCAATAg                                                                                     >  1:310784/1‑149 (MQ=255)
                                                              ggggCGCGATGAGCAAACTCCCGAGGAAACGGGCAACCAGACCTGCCCCCAATAACAGGCCGATGGTTTCTGGCGTTAAACCAATCCCTTTAAGCCAGACGCTCCAGAAAGGTAGAAAAATGCCGTAACTAAAAAAGTATGTGAAATAg                                                                                     >  1:743948/1‑149 (MQ=255)
                                                              ggggCGCGATGAGCAAACTCCCGAGGAAACGGGCAACCAGACCTGCCCCCAATAACAGGCCGATGGTTTCTGGCGTTAAACCAATCCCTTTAAGCCAGACGCTCCAGAAAGGTAGAAAAATGCCGTAACTAAAAAAGTATGTGAAATAg                                                                                     >  1:743960/1‑149 (MQ=255)
                                                                ggCGCGATGAGCAAACTCCCGAGGAAACGGGCAACCAGACCTGCCCCCAATAACAGGCCGATGGTTTCTGGCGTTAAACCAATCCCTTTAAGCCAGACGCTCCAGAAAGGTAGAAAAATGCCGTAACTAAAAAAGTATGTGAAATAGcc                                                                                   <  2:322420/149‑1 (MQ=255)
                                                                 gcgcGATGAGCAAACTCCCGAGGAAACGGGCAACCAGACCTGCCCCCAATAACAGGCCGATGGTTTCTGGCGTTAAACCAATCCCTTTAAGCCAGACGCTCCAGAAAGGTAGAAAAATGCCGTAACTAAAAAAGTATGTGAAATAGCCg                                                                                  >  2:1180274/1‑149 (MQ=255)
                                                                  cgcgATGAGCAAACTCCCGAGGAAACGGGCAACCAGACCTGCCCCCAATAACAGGCCGATGGTTTCTGGCGTTAAACCAATCCCTTTAAGCCAGACGCTCCAGAAAGGTAGAAAAATGCCGTAACTAAAAAAGTATGTGAAATAGCCGa                                                                                 >  1:709086/1‑149 (MQ=255)
                                                                  cgcgATGAGCAAACCGCCTAGGAAACGGGCAACCAGACATACCCCCAATAACAGGCCGATGGTTTCTGGCGTTAAGCCAATCCCTTTAAGCCAGACGCTCCAGAAAGGTAGAAAAATGCCGTATCTAAAAAAGTATGTGAAATAGCCgt                                                                                 >  1:912283/1‑148 (MQ=255)
                                                                   gcgATGAGCAAACTCCCGAGGAAACGGGCAACCAGACCTGCCCCCAATAACAGGCCGATGGTTTCTGGCGTTAAACCAATCCCTTTAAGCCAGACGCTCCAGAAAGGTAGAAAAATGCCGTAACTAAAAAAGTATGTGAAATAGCCGAg                                                                                >  2:865491/1‑149 (MQ=255)
                                                                    cgATGAGCAAACGCCCGAGGAAACGGGCAACCAGACCTGCCCCCAATAACAGGCCGATGGTTTCTGGCGTTAAACCAATCCCTTTAAGCCAGACGCTCCAGAAAGGTAGAAAAATGCCGTAACTAAAAAAGTATGTGAAATAGCCGAgc                                                                               <  1:837516/149‑1 (MQ=255)
                                                                      aTGAGCAAACTCCCGAGGAAACGGGCAACCAGACCTGCCCCCAATACCAGGCCGATGGTTTCTGGCGTTAAACCAATCCCTTGAAGCCAGACGCTCCAGAAAGGTAGAAAAATGCCGTAACTAAAAAAGTATGTGAAATAGCCGAgcgc                                                                             <  2:855972/149‑1 (MQ=255)
                                                                                  ccGGGGAACCGGGCAACCAGACCTGCCCCCAACAACAGGCCGATGGTTTCTGGCGTTAAACCAATCCCTTTAAGCCAGACGCTCCAGAAAGGTAGAAAAATGCCGTAACTAAAAAAGTATGTGAAATAGCCGAGCGCCAACCAGCGCGt                                                                 <  2:942900/149‑1 (MQ=255)
                                                                                  ccGAGGAAACGGGCAACCAGACCTGCCCCCAATAACAGGCCGATGGTTTCTGGCGTTAAACCAATCCCTTGAAGCCAGACGCTCCAGAAAGGTAGAAAAATGCCGTAACTAAAAAAGTATGTGAAATAGCCGAGCGCCAACCAGCGCGt                                                                 >  2:994582/1‑149 (MQ=255)
                                                                                       gAAGCGGGCACTCAGACCTGCTCCCAATAACCGGCCGATGGTTTCTGGCGTTAAACCAATCCCTTTAAGCCAGACGCTCCAGAAAGTTAGAAAAATGCCGTAACTAAAAAAGTATGTGAAATAGCCGAGCGCCAACCAGCGCGTGGAtt                                                            <  2:144139/149‑1 (MQ=255)
                                                                                       gAAACGGGCAACCAGACCTGCCCCCAATAACAGGCCGATGGTTTCTGGCGTTAAACCAATCCCTTTAAGCCAGACGCTCCAGAAAGGTAGAAAAATGCCGTAACTAAAAAAGTATGTGAAATAGCCGAGCGCCAACCAGCGCGTGGAtt                                                            <  1:1134962/149‑1 (MQ=255)
                                                                                       gAAACGGGCAACCAGACCTGCCCCCAATAACAGGCCGATGGTTTCTGGCGTTAAACCAATCCCTTTAAGCCAGACGCTCCAGAAAGGTAGAAAAATGCCGTAACTAAAAAAGTATGTGAAATAGCCGAGCGCCAACCAGCGCGTGGAtt                                                            <  1:612480/149‑1 (MQ=255)
                                                                                       gAAACGGGCAACCAGACCTGCCCCCAATAACAGGCCGATGGTTTCTGGCGTTAAACCAATCCCTTTAAGCCAGACGCTCCAGAAAGGTAGAAAAATGCCGTAACTAAAAAAGTATGTGAAATAGCCGAGCGCCAACCAGCGCGAGGAtt                                                            >  2:1094346/1‑149 (MQ=255)
                                                                                       gAAACGGGCAACCAGACCTGCCCCCAATAACAGGCCGATGGTTTCTGGCGTTAAACCAATCCCTTTAAGCCAGACGCTCCAGAAAGGTAGAAAAATGCCGTAACTAAAAAAGTATGCGAAATAGCCGAGCGCCAACCAGCGCGTGGAtt                                                            >  2:1094444/1‑149 (MQ=255)
                                                                                       gAAACGGGCAACCAGACCTGCCCCCAATAACAGGCCGATGGTTTCTGGCGTTAAACCAATCCCTTTAAGCCAGACGCTCCAGAAAGGTAGAAAAAAGCCATAACTAAAAAAGTATGTGAAATAGCCGAGCGCCAACCAGCGCGTGGAtt                                                            <  2:155021/149‑1 (MQ=255)
                                                                                        aaaCGGGCAACCAGACCTGCCCCCAATAACAGGCCGATGGTTTCTGGCGTTAAACCAATCCCTTTAAGCCAGACGCTCCAGAAAGGTAGAAAAATGCCGTAACTAAAAAAGTATGTGAAATAGCCGAGCGCCAACCAGCGCGTGGATTg                                                           <  1:875482/149‑1 (MQ=255)
                                                                                             ggCAACCAGACCTGCCCCCAATAACAGGCCGATGGTTTCTGGCGTTAAACCAATCCCTTTAAGCCAGACGCTCCAGAAAGGTAGAAAAATGCCGTAACTAAAAAAGTATGTGAAATAGCCGAGCGCCAACCAGCGCGTGGATTGCaaaa                                                      <  2:800835/149‑1 (MQ=255)
                                                                                             ggCAACCAGACCTGCCCCCAATAACAGGCCGATGGTTTCAGGCGTTAAACCAATCCCTTTAAGCCAGACGCTCCAGAAAGGTAGAAAAATGCCGTAACTAAAAAAGTATGTGAAATAGCCGAGCGCCTACCAGCGCGTGGATTGCaaaa                                                      <  1:878062/149‑1 (MQ=255)
                                                                                                   cAGACCTGCCTCCAATAACAGGCCGATGGTTTCTGGCGTTAAACCAAACCCCTTAAGCCAGACGCTCCAGAAAGGTAGAAAAATGCCGTAACTAAAAAAGTATGTGAAATAGCCGAGCGCCAACCAGCGCGTGGATTGCAAAACCATGa                                                <  2:367623/149‑1 (MQ=255)
                                                                                                    cgaactccccccAATAACAGGCCGATGGTTTCTGGCGTTCAACCAACCCCCTTAAGCCAGACGCTCCAGAAAGGTAGAAAAATGCCGTAACTAAACAAGTATGTGAAATAGCCGAGCGCCAACCAGCGCGTGGATTGCAAAACCATGAg                                               <  1:201442/142‑1 (MQ=255)
                                                                                                           ccccccACAACAGGCCGAGGGTTGCTGGCGTTAAACCAATCCCTTTAAGCCAGACGCTCCAGAAAGGTAGAAAAATGCCGTAACTAAAAAAGTATGTGAAATAGCCGAGCGCCAACCAGCGCGTGGATTGCAAAACCATGAGtccctcc                                        <  2:607508/149‑1 (MQ=255)
                                                                                                           cccccAATAACAGGCCGATGGTTTCTGGCGTTAAACCAATCGCTTTAAGCCAGACGCTCCAGAAAGGTAGAAAAATGCCGTAACTAAAAAAGTATGTGAAATAGCCGAGCGCCAACCAGCGCGTGGATTGCAAAACCATGAGtccctcc                                        <  2:1183866/149‑1 (MQ=255)
                                                                                                           cccccAATAACAGGCCGATGGTTTCTGGCGTTAAACCAATCGCTTTAAGCCAGACGCTCCAGAAAGGTAGAAAAATGCCGTAACTAAAAAAGTATGTGAAATAGCCGAGCGCCAACCAGCGCGTGGATTGCAAAACCATGAGtccctcc                                        <  1:708462/149‑1 (MQ=255)
                                                                                                           cccccAATAACAGGCCGATGGTTTCTGGCGTTAAACCAATCCCTGTAAGCCAGACGCTCCAGAAAGGTAGAAAAATGCCGTAACTAAAAAAGTATGTGAAATAGCCGAGCGCCAACCAGCGCGTGGATTGCAAAACCATGAGtccctcc                                        <  1:81402/149‑1 (MQ=255)
                                                                                                              ccAATAACAGGCCGATGGTTTCTGGCGTTAAACCAATCCCTTTAAGCCAGACGCTCCAGAAAGGTAGAAAAATGCCGTAACTAAAAAAGTATGTGAAATAGCCGACCGCCAACCAGCGCGTGGAGTGCAAAACCATGAGTCCCTCCCgc                                     >  2:1054302/1‑148 (MQ=255)
                                                                                                                attAACAGGCCGATGGTGTCTGGCGTTAAACCAATCCCTTTAAGCCAGACGCTCCAGAAAGGTAGAAAAATGCCGTAACTAAAAAAGTATGTGAAATAGCCGAGCGCCAACCAGCGCGTGGATTGCAAAACCATGAGTCCCTCCCGTCa                                   <  2:220689/147‑1 (MQ=255)
                                                                                                                aaTAACAGGCCGATGGTTTCTGGCGTTACAACGATGCCTGTAAGCCAGACGCTCCAGAAAGGTAGAAAAATGCCGTAACTAAAAAAGTATGTGAAATAGCCGAGCGCCAACCAGCGCGGGGATTGCAAAACCATGAGTCCCTCCCGTCa                                   >  1:1226440/1‑149 (MQ=255)
                                                                                                                 aTAACAGGCCGATGGTTTCTGGCGTTAAACCAATCCCTTTAAGCCAGACGCTCCAGAAAGGTAGAAAAATGCCGTAACTAAAAAAGTATGTGAAATAGCCGAGCGCCAACCAGCGCGTGGATTGCAAAACCATGAGTCCCTCCCGTCag                                  <  1:1180274/149‑1 (MQ=255)
                                                                                                                  tAACAGGCCGATGGTTTCTGGCGTTAAACCAATCCCTTTAAGCCAGACGCTCCAGAAAGGTAGAAAAATGCCGTAACTAAAAAAGTATGTGAAATAGCCGAGCGCCAACCAGCGCGTGGCTTGCAAAACCATGAGTCCCTCCCGTCagg                                 >  1:529868/1‑149 (MQ=255)
                                                                                                                   aaCCGGCCGATGGTTTCTGGCGTTAAACCAATCCCTTTAAGCCAGACGCTCCAGAAAGGTAGAAAAATGTCGTAACTAAAAAAGTATGTGAAATAGCCGAGCGCCAACCAGCGCGTGGATTGCAAAACCATGAGTCCCTCCCGTCagga                                <  1:670615/149‑1 (MQ=255)
                                                                                                                    aCAGGCCGATGGTTTCTGGCGTTAAACCAATCCCTTTAAGCCAGCCGCTCCAGAAAGGTAGAAAAATGCCGTAACTAAAAAAGTATGTGAAATAGCCGAGCGCCAACCAGCGCGTGGATTGCAAAACCATGAGTCCCTCCCGTCaggag                               <  1:1105944/149‑1 (MQ=255)
                                                                                                                    aCAGGCCGAGGGTGGCGGGCGGGAAACCAATCCCTTTAAGCCAGACGCTCCAGAAAGGTAGAAAAATGCCGTAACTAAAAAAGTATGTGAAATAGCCGAGCGCCAACCAGCGCGTGGATTGCAAAACCATGAGTCCCTCCCGTCaggag                               <  2:1117268/149‑1 (MQ=255)
                                                                                                                      aGGCCGATGGTTTCTGGCGTTAAACCAATCCCTTTAAGCCAGACGCTCCAGAAAGGTAGAAAAATGCCGTAACTAAAAAAGTATGTGAAATAGCCGAGCGCCAACCAGCGCGTGGATTGCAAAACCATGAGTCCCTCCCGTCAGGAGGc                             >  2:420635/1‑149 (MQ=255)
                                                                                                                       ggCCGATGGTTTCTGGCGTTAAACCAATCCCTTTAAGCCAGACGCTCCAGAAAGGTAG‑AAAATGCCGTAACTAAAAAAGTATGTGAAATAGCCGAGCGCCAACCAGCGCGTGGATTGCAAAACCATGAGTCCCTCCCGTCAGGAGGCGt                           <  1:313583/149‑1 (MQ=255)
                                                                                                                          cGATGGTGTACGGCTTCAAACCAAGCCCGTTGAGCGTGACGTTCCAGAAAGGTAGAAAAATGCCGTAACTAAAAAAGTATGTGAAATTGCCGAGCGCCAACCAGCGCGTGGATTGCAAAACCAGGAGTCCCTCCCGTCAGGAGGCGTTa                         <  2:590540/149‑1 (MQ=255)
                                                                                                                               gtgtCTGGCGTTAAACCAATCCCTTTAAGCCAGACGCTCCAGAAAGGTAGAAAAATGCCGTAACTAAAAAAGTATGTTAAATAGCCGAGCGCCAACCAGCGCGTGGATTGCAAAACCATGAGTCCCTCCCGTCAGGAGGCGTTAGTCTg                    <  2:123398/146‑1 (MQ=255)
                                                                                                                               gTTTCTGGCGTTAAACCAATCCCTTTAAGCCAGACGCTCCAGAAAGGTAGAAAAATGCCGTAACTAAAAAAGTATGTGAAATAGCCGGGCGCCAACCAGCGCGTGGATTGCAAAACCATGAGTCCCTCCCGTCAGGAGGCGTTAGTCTg                    <  2:152482/149‑1 (MQ=255)
                                                                                                                                        gTTAAACCAATCCCTTTAAGCCAGACGCTCCAGAAAGGTAGAAAAATGCCGTAACTAAAAAAGTATGTGAAATAGCCGAGCGCCAACCAGCGCGTGGATTGCAAAACCATGAGTCCCTCCCGTCAGGAGGCGTTAGTCTGGGGTTAAtg           >  2:453747/1‑148 (MQ=255)
                                                                                                                                        gTTAAACCAATCCCTTTAAGCCAGACGCTCCAGAAAGGTAGAAAAATGCCGTAACTAAAAAAGTATGTGAAATAGCCGAGCGCCAACCAGCGCGTGGATTGCAAAACCATGAGTCCCTCCCGTCAGGAGGCGTTAGTCTGGGGTTAATc           >  1:528290/1‑149 (MQ=255)
                                                                                                                                        gTTAAACCAATCCCTTTAAGCCAGACGCTCCAGAAAGGTAGAAAAATGCCGTAACTAAAAAAGTATGTGAAATAGCCGAGCGCCAACCAGCGCGTGGATTGCAAAACCATGAGTCCCTCCCGTCAGGAGGCGTTAGTCTGGGGTTAATc           >  1:553333/1‑149 (MQ=255)
                                                                                                                                        gTTAAACCAATCCCTTTAAGCCAGACGCTCCAGAAAGGTAGAAAAATGCCGTAACTAAAAAAGTATGGGAAATAGCCGAGCGCCAACCAGCGCGTGGATTGCAAAACCATGAGTCCCTCCCGCGAGGAGGCGCTAGTCTGGGGTTAATc           >  1:751741/1‑149 (MQ=255)
                                                                                                                                           aaaCCAATCCCTTTAAGCCAGGCGCGCCAGAAAGGTAGAAAAATGCCGTAACTAAAAAAGTATGTGAAATAGCCGAGCGCCAACCAGCGCGTGGATTGCAAAACCATGAGTCCCTCCCGTCAGGAGGCGTTAGTCTGGGGTTAATCGaa        <  1:25450/149‑1 (MQ=255)
                                                                                                                                           aaaCCAATCCCTTTAAGCCAGACGCTCCAGAAAGGTAGAAAAATGCCGTAACTAAAAAACTATGTGAAATAGCCGAGCGCCAACCAGCGCGTGGATTGCAAAACCATGAGTCCCTCCCATCAGTAGGCGTTAGTCTGGGGTTAATCGaa        <  1:643752/149‑1 (MQ=255)
                                                                                                                                            aaCCAATTCCTTTAAGCCAGACGCTCCAGAAAGGTAGGAAAGTGCCGTAACTAAGAAAGTATGTGAAATAGCCGAGCGCCAACCAGCGCGTGGATTGCAAAACCATGAGTCCTTCCCGTCAGGAGGCGTTAGTCTGGGGTTAATCGAAt       <  2:935598/149‑1 (MQ=255)
                                                                                                                                            aaCCAATCCCTTTAAGCCAGACGCTCCAGAAAGGTAGAAAAATGCCGTAACTAAAAAAGTATGTGAAATAGCCGAGCGCCAACCAGCGCGTGGATTGCAAAACCATGAGTCCCTCCCGTCAGGGGGCGTTAGTCTGGGGTTAATCGAAt       >  1:310337/1‑149 (MQ=255)
                                                                                                                                            aaCCAATCCCTTTAAGCCAGACGCTCCAGAAAGGTAGAAAAATGCCGTAACTAAAAAAGTATGTGAAATAGCCGAGCGCCAACCAGCGCGTGGATTGCAAAACCATGAGTCCCTCCCGTCAGGAGGCGTTAGTCTGGGGTTAATCGAAt       >  1:1095799/1‑149 (MQ=255)
                                                                                                                                                 aTCCCTTTAAGCCAGACGCTCCAGAAAGGTAGAAAAATGCCGTAACTAAAAAAGTATGTGAAATAGCCGAGCGCCAACCAGCGCGTGGATTGCAAAACCATGAGTCCCTCCCGTCAGGAGGCGTTAGTCTGGGGTTAATCGAATGTTgg  <  2:709086/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
GAGAAGTGTCAGCAGTGCCAGCACGCGCAAGGCGGAAATCAGGCGGGAAGGATCGCTGACGCGGGGCGCGATGAGCAAACTCCCGAGGAAACGGGCAACCAGACCTGCCCCCAATAACAGGCCGATGGTTTCTGGCGTTAAACCAATCCCTTTAAGCCAGACGCTCCAGAAAGGTAGAAAAATGCCGTAACTAAAAAAGTATGTGAAATAGCCGAGCGCCAACCAGCGCGTGGATTGCAAAACCATGAGTCCCTCCCGTCAGGAGGCGTTAGTCTGGGGTTAATCGAATGTTGG  >  CP009273/2660960‑2661253

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑