Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,265,089 T→A 100% intergenic (‑120/‑310) yjcB ← / → yjcC putative inner membrane protein/putative membrane‑anchored cyclic‑di‑GMP phosphodiesterase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,265,0890TA99.1% 716.6 / ‑2.6 216intergenic (‑120/‑310)yjcB/yjcCputative inner membrane protein/putative membrane‑anchored cyclic‑di‑GMP phosphodiesterase
Reads supporting (aligned to +/- strand):  ref base T (2/0);  new base A (92/122);  total (94/122)
Fisher's exact test for biased strand distribution p-value = 1.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 223 total aligned reads displayed.

AAGAACCACGCCAGTGGTGAGGGTAGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGTGGCATTAAGGGCATATCTTCCCGCCGCCTCTGCATTCCTGTAGGAAATTAATT  >  CP009273/4264940‑4265235
                                                                                                                                                     |                                                                                                                                                  
   aaCCACGCCAGTGGTGAGGGTAGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCTGTTTGTATCa                                                                                                                                                 >  2:3971879/1‑150 (MQ=255)
       aCGCCAGTGGTGAGGGTAGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCTGTTTGTATCATTaa                                                                                                                                             >  2:417036/1‑150 (MQ=255)
         gCCAGTGGTGAGGGTAGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCTGTTTGTATCATTAAtt                                                                                                                                           >  2:3985860/1‑150 (MQ=255)
          ccAGTGGTGAGGGTAGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCTGTTTGTATCATTAAttt                                                                                                                                          <  2:773729/150‑1 (MQ=255)
            aGTGGTGAGGGTAGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCTGTTTGTATCATTAATTTtg                                                                                                                                        >  1:3167009/1‑150 (MQ=255)
             gtggtgAGGGTAGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCTGTTTGTATCATTAATTTtgt                                                                                                                                       <  2:4057710/150‑1 (MQ=255)
                gtgAGGGTAGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCTGTTTGTATCATTAATTTTGTGaa                                                                                                                                    >  2:2972108/1‑150 (MQ=255)
                    gggTAGCCATCAAACAATTCCTCTGAAACGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCTGTTTGTATCATTAATTTTGTGAACtat                                                                                                                                >  1:1754958/1‑150 (MQ=255)
                     ggTAGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTCCACTCTGGCTATCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCTGTTTGTATCATTAATTTTGTGAACtata                                                                                                                               >  1:2893626/1‑150 (MQ=255)
                      gTAGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCTGTTTGTATCATTAATTTTGTGAACTatat                                                                                                                              <  1:2374306/150‑1 (MQ=255)
                        aGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCTGTTTGTATCATTAATTTTGTGAACTATATca                                                                                                                            >  1:920948/1‑150 (MQ=255)
                        aGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCTGTTTGTATCATTAATTTTGTGAACTATATca                                                                                                                            <  2:1385857/150‑1 (MQ=255)
                         gCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCTGTTTGTATCATTAATTTTGTGAACTATATcac                                                                                                                           >  1:618499/1‑150 (MQ=255)
                          ccATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCTGTTTGTATCATTAATTTTGTGAACTATATcaca                                                                                                                          <  1:4125892/150‑1 (MQ=255)
                          ccATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCTGTTTGTATCATTAATTTTGTGAACTATATcaca                                                                                                                          <  1:2833794/150‑1 (MQ=255)
                          ccATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTACAGAAAACAAAAATGAAACACTCTGTTTGTATCATTAATTTTGTGAACTATATcaca                                                                                                                          >  2:3685184/1‑150 (MQ=255)
                              cAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCTGTTTGTATCATTAATTTTGTGAACTATATCACAattg                                                                                                                      >  2:2989219/1‑150 (MQ=255)
                              cAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCTGTTTGTATCATTAATTTTGTGAACTATATCACAattg                                                                                                                      >  2:2882212/1‑150 (MQ=255)
                              cAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCTGTTTGTATCATTAATTTTGTGAACTATATCACAattg                                                                                                                      >  2:5050626/1‑150 (MQ=255)
                                 aCAATTCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCTGTTTGTATCATTAATTTTGTGAACTATATCACAattgatt                                                                                                                   >  1:1493556/1‑150 (MQ=255)
                                  cAATTCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCTGTTTGTATCATTAATTTTGTGAACTATATCACAattgattg                                                                                                                  >  1:4309931/1‑150 (MQ=255)
                                     ttCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCTGTTTGTATCATTAATTTTGTGAACTATATCACAATTGAttgttt                                                                                                               >  2:5118801/1‑150 (MQ=255)
                                      tCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCTGTTTGTATCATTAATTTTGTGAACTATATCACAATTGAttgtttg                                                                                                              >  2:4678778/1‑150 (MQ=255)
                                       ccTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCTGTTTGTATCATTAATTTTGTGAACTATATCACAATTGAttgtttgt                                                                                                             <  1:4947881/150‑1 (MQ=255)
                                       ccTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCTGTTTGTATCATTAATTTTGTGAACTATATCACAATTGAttgtttgt                                                                                                             <  2:2487044/150‑1 (MQ=255)
                                        ctctGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCTGTTTGTATCATTAATTTTGTGAACTATATCACAATTGATtgtttgtt                                                                                                            <  1:2937535/150‑1 (MQ=255)
                                          ctGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCTGTTTGTATCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAg                                                                                                          >  2:4759941/1‑150 (MQ=255)
                                            taaaaGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCTGTTTGTATCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGcc                                                                                                        >  2:4709027/2‑150 (MQ=255)
                                             aaaaGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCTGTTTGTATCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCa                                                                                                       >  1:3251687/1‑150 (MQ=255)
                                              aaaGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCTGTTTGTATCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCat                                                                                                      <  1:3070645/150‑1 (MQ=255)
                                               aaGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCTGTTTGTATCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCata                                                                                                     >  1:2105135/1‑150 (MQ=255)
                                                 gagCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCTGTTTGTATCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATAtt                                                                                                   <  2:1391801/150‑1 (MQ=255)
                                                 gagCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCTGTTTGTATCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATAtt                                                                                                   >  1:1191037/1‑150 (MQ=255)
                                                   gCCGATGCCCTGGACAGCGCGGCGTTGTTCTCGGGAGGCGTGGTACACTCTGGCTATCGCGGTGCTTGCAGAAAACAAAAATGAAACACTCTGTTTGTATCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTTGCCATATTAg                                                                                                 <  1:3428240/150‑1 (MQ=255)
                                                   gCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCTGTTTGTATCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAg                                                                                                 <  1:134024/150‑1 (MQ=255)
                                                   gCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCTGTTTGTATCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAg                                                                                                 <  1:134089/150‑1 (MQ=255)
                                                   gCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCTGTTTGTATCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAg                                                                                                 <  1:736295/150‑1 (MQ=255)
                                                     cGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCTGTTTGTATCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGc                                                                                               >  2:24657/1‑150 (MQ=255)
                                                      gATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCTGTTTGTATCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGcc                                                                                              >  2:1679070/1‑150 (MQ=255)
                                                         gCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCTGTTTGTATCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTg                                                                                           <  1:805815/150‑1 (MQ=255)
                                                           ccTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCTGTTTGTATCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGAc                                                                                         >  2:2591677/1‑150 (MQ=255)
                                                           ccTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCTGTTTGTATCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGAc                                                                                         >  1:4513332/1‑150 (MQ=255)
                                                             tGGACAGCGCGGCGTTGTTCACGGTATGCGTGGTACAATCTGGCTATCGCGGGGCTTGCAGAAAACAAAAATGAAACACTATGTTTGTATGATTAATTTTGTTAACTATAAAACAAGTGATTGTTTGTTATACATATAATGCAGTTACtt                                                                                       >  2:4844456/1‑150 (MQ=255)
                                                             tGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCTGTTTGTATTATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACtt                                                                                       >  1:4494265/1‑150 (MQ=255)
                                                             tGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCTGTTTGTATCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACtt                                                                                       <  2:2177379/150‑1 (MQ=255)
                                                              ggACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCTGTTTGTATCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACttt                                                                                      <  1:2562737/150‑1 (MQ=255)
                                                                  aGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCTGTTTGTATCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTAtt                                                                                  <  1:650764/150‑1 (MQ=255)
                                                                  aGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCTGTTTGTATCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTAtt                                                                                  >  1:3666039/1‑150 (MQ=255)
                                                                  aGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCTGTTTGTATCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTAtt                                                                                  <  2:2566874/150‑1 (MQ=255)
                                                                   gcgcGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCTGTTTGTATCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTg                                                                                 <  1:4560449/150‑1 (MQ=255)
                                                                   gcgcGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCTGTTTGTATCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTg                                                                                 <  2:803245/150‑1 (MQ=255)
                                                                     gcggcgTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCTGTTTGTATCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCt                                                                               <  2:4771149/150‑1 (MQ=255)
                                                                           ttgttCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCTGTTTGTATCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAg                                                                         <  1:2391404/150‑1 (MQ=255)
                                                                           ttgttCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCTGTTTGTATCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAg                                                                         <  1:2208469/150‑1 (MQ=255)
                                                                              ttCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCTGTTTGTATCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGAtt                                                                      <  1:1272131/150‑1 (MQ=255)
                                                                                cACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCTGTTTGTATCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTGGGCCGTGACTTTTATTGCTGTACAGATTAt                                                                    <  2:271945/150‑1 (MQ=255)
                                                                                  cGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCTGTTTGTATCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTAtgt                                                                  >  2:3243554/1‑150 (MQ=255)
                                                                                    ggAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCTGTTTGTATCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTgg                                                                <  2:1580567/150‑1 (MQ=255)
                                                                                      aGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCTGTTTGTATCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGtt                                                              <  2:3833502/150‑1 (MQ=255)
                                                                                         cGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCTGTTTGTATCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGttttt                                                           >  1:2454805/1‑150 (MQ=255)
                                                                                             gTACACTCTGGCTATCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCTGTTTGTATCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTCCAGATTATGTGGTTTTTCAGt                                                       >  1:3951989/1‑150 (MQ=255)
                                                                                             gTACACTCTGGCTATCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCTGTTTGTATCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGt                                                       <  2:3675680/150‑1 (MQ=255)
                                                                                               acacTCTGGCTATCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCTGTTTGTATCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGTgg                                                     >  2:1450009/1‑150 (MQ=255)
                                                                                                cacTCTGGCTATCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCTGTTTGTATCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGTGGc                                                    >  2:2197560/1‑150 (MQ=255)
                                                                                                   tctGGCTATCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCTGTTTGTATCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGTGGCAtt                                                 <  1:3350096/150‑1 (MQ=255)
                                                                                                    ctGGCTATCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCTGTTTGTATCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGTGGCATTa                                                <  2:4097282/150‑1 (MQ=255)
                                                                                                    ctGGCTATCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCTGTTTGTATCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGTGGCATTa                                                >  2:3462291/1‑150 (MQ=255)
                                                                                                      ggCTATCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCTGTTTGTATCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGTGGCATTAAg                                              <  2:2148505/150‑1 (MQ=255)
                                                                                                       gCTATCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCTGTTTGTATCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGTGGCATTAAgg                                             <  1:848720/150‑1 (MQ=255)
                                                                                                         tATCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCTGTTTGTATCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGTGGCATTAAGGGc                                           <  1:1944597/150‑1 (MQ=255)
                                                                                                          aTCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCTGTTTGTATCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGTGGCATTAAGGGCa                                          >  2:385304/1‑150 (MQ=255)
                                                                                                            cgcgGGGCTTGCAGCAAACAAAAATGAAACACTCTGTTAGTATCATTAATATTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGGACAGATTATGTGGTTTTTCAGTGGCATTAAGGGCacc                                        <  1:2210346/150‑3 (MQ=255)
                                                                                                            cgcgGGGCTTGCAGAAAACAAAAATGAAACACTCTGTTTGTATCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGTGGCATTAAGGGCata                                        <  2:3341460/150‑1 (MQ=255)
                                                                                                                 ggCTTGCAGAAAACAAAAATGAAACACTCTGTTTGTATCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGTGGCATTAAGGGCATATCTTc                                   >  2:1314019/1‑150 (MQ=255)
                                                                                                                      gCAGAAAACAAAAATGAAACACTCTGTTTGTATCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGTGGCATTAAGGGCATATCTTCccgcc                              <  2:3069670/150‑1 (MQ=255)
                                                                                                                      gCAGAAAACAAAAATGAAACACTCTGTTTGTATCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGTGGCATTAAGGGCATATCTTCccgcc                              <  1:5050626/150‑1 (MQ=255)
                                                                                                                      gCAGAAAACAAAAATGAAACACTCTGTTTGTATCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGTGGCATTAAGGGCATATCTTCccgcc                              <  1:3975362/150‑1 (MQ=255)
                                                                                                                      gCAGAAAACAAAAATGAAACACTCTGTTTGTATCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTATTCAGTGGCATTAAGGGCATATCTTCccgcc                              >  1:1083692/1‑150 (MQ=255)
                                                                                                                       cAGAAAACAAAAATGAAACACTCTGTTTGTATCATTAATATTTTGAACTATATCACAATTGATAGTTTGTTAGCCATATAAGGCAGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGTGGCATTAAGGGAATATCTTCccgccg                             >  2:493972/1‑150 (MQ=255)
                                                                                                                         gAAAACAAAAATGAAACACTCTGTTTGTATCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGTGGCATTAAGGGCATATCTTCCCgccgcc                           <  1:5027245/150‑1 (MQ=255)
                                                                                                                         gAAAACAAAAATGAAACACTCTGTTTGTATCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGTGGCATTAAGGGCATATCTTCCCgccgcc                           <  2:2231606/150‑1 (MQ=255)
                                                                                                                         gAAAACAAAAATGAAACACTCTGTTTGTATCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGTGGCATTAAGGGCATATCTTCCCgccgcc                           <  1:3612639/150‑1 (MQ=255)
                                                                                                                           aaacaAAAATGAAACACTCTGTTTGTATCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGTGGCATTAAGGGCATATCTTCCCGCCGCctc                         <  1:328367/148‑1 (MQ=255)
                                                                                                                              caAAAATGAAACACTCTGTTTGTATCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGTGGCATTAAGGGCATATCTTCCCGCCGCCTCTGc                      >  1:3462230/1‑150 (MQ=255)
                                                                                                                              caAAAATGAAACACTCTGTTTGTATCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGTGGCATTAAGGGCATATCTTCCCGCCGCCTCTGc                      >  1:4899290/1‑150 (MQ=255)
                                                                                                                               aaaaaTGAAACACTCTGTTTGTATCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGTGGCATTAAGGGCATATCTTCCCGCCGCCTCTGCa                     >  1:2838559/1‑150 (MQ=255)
                                                                                                                               aaaaaTGAAACACTCTGTTTGTATCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGTGGCATTAAGGGCATATCTTCCCGCCGCCTCTGCa                     >  2:4451043/1‑150 (MQ=255)
                                                                                                                               aaaaaTGAAACACTCTGTTTGTATCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGTGGCATTAAGGGCATATCTTCCCGCCGCCTCTGCa                     <  2:1493556/150‑1 (MQ=255)
                                                                                                                                aaaaTGAAACACTCTGTTTGTATCATTAATTTTGTGAACTATATCAAAATTGATTGTTTGTTAGCCATATTAGGCCTTGACTTTTATTGCTGTACAGATTATGTGGTTATTCAGTGGCATTAAGGGCATATCTTCCCGCCGCCTCTGCAt                    >  2:860015/1‑150 (MQ=255)
                                                                                                                                  aaTGAAACACTCTGTTTGTATCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGTGGCATTAAGGGCATATCTTCCCGCCGCCTCTGCATTc                  <  2:3368723/150‑1 (MQ=255)
                                                                                                                                   aTGAAACACTCTGTTTGTATCATTAATTTTGTGACCTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGTGGCATTAAGGGCATATCTTCCCGCCGCCTCTGCATTcc                 <  2:600135/150‑1 (MQ=255)
                                                                                                                                   aTGAAACACTCTGTTTGTATCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGTGGCATTAAGGGCATATCTTCCCGCCGCCTCTGCATTcc                 >  2:4207062/1‑150 (MQ=255)
                                                                                                                                   aTGAAACACTCTGTTTGTATCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAACCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGTGGCATTAAGGGCATATCTTCCCGCCGCCTCTGCATTcc                 <  1:4440159/150‑1 (MQ=255)
                                                                                                                                     gAAACACTCTGTTTGTATCATTAATTTTGTGAACTAGATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGTGGCATTAAGGGCATATCTTCCCGCCGCCTCTGCATTCCTg               <  2:1834728/150‑1 (MQ=255)
                                                                                                                                      aaaCACTCTGTTTGTATCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGTGGCATTAAGGGCATATCTTCCCGCCGCCTCTGCATTCCTGt              <  2:919726/150‑1 (MQ=255)
                                                                                                                                       aaCACTCTGTTTGTATCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGTGGCATTAAGGGCATATCTTCCCGCCGCCTCTGCATTCCTGTa             <  2:4705535/150‑1 (MQ=255)
                                                                                                                                       aaCACTCTGTTTGTATCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGTGGCATTAAGGGCATATCTTCCCGCCGCCTCTGCATTCCTGTa             <  1:1679065/150‑1 (MQ=255)
                                                                                                                                        acacTCTGTTTGTATCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGTGGCATTAAGGGCATATCTTCCCGCCGCCTCTGCATTCCTGTAg            <  2:4583109/150‑1 (MQ=255)
                                                                                                                                                tttgtatCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGTGGCATTAAGGGCATATCTTCCCGCCGCCTCTGCATTCCTGTAGGAaattaa    <  1:213928/150‑1 (MQ=255)
                                                                                                                                                  tgtttCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGTGGCATTAAGGGCATATCTTCCCGCCGCCTCTGCATTCCTGTAGGAaattaatt  >  1:2027758/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
AAGAACCACGCCAGTGGTGAGGGTAGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGTGGCATTAAGGGCATATCTTCCCGCCGCCTCTGCATTCCTGTAGGAAATTAATT  >  CP009273/4264940‑4265235

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: