Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,628,207 G→A 100% G110E (GGA→GAA)  pinQ → Qin prophage; putative site‑specific recombinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,628,2070GA99.1% 497.5 / ‑2.7 216G110E (GGA→GAA) pinQQin prophage; putative site‑specific recombinase
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base A (106/108);  total (107/109)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-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 222 total aligned reads displayed.

AATGTGGTGATCAATTGATTGTGACAAAACTGGATCGCCTTGGTTGTAATGCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGTAAGAGCCCGCGGCGCAGGGAAACGTTTTGGTCGACCACCTGTGTTAAATGAAGAACAGAAACAG  >  CP009273/1628060‑1628352
                                                                                                                                                   |                                                                                                                                                 
aaTGTGGTGATCAATTGATTGTGACAAAACTGGATCGCCTTGGTTGTAATGCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGaaa                                                                                                                                                 >  2:2935655/1‑150 (MQ=17)
   gtggtgATCAATTGATTGTGACAAAACTGGATCGCCTTGGTTGTAATGCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGaaaaaa                                                                                                                                              >  2:4538768/1‑150 (MQ=17)
     ggtgATCAATTGATTGTGACAAAACTGGATCGCCTTGGTTGTAATGCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAatg                                                                                                                                            >  1:2321805/1‑150 (MQ=17)
           cAATTGATTGTGACAAAACTGGATCGCCTTGGTTGTAATGCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAatgatgatg                                                                                                                                      >  2:1943052/1‑150 (MQ=17)
           cAATTGATTGTGACAAAACTGGATCGCCTTGGTTGTAATGCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAatgatgatg                                                                                                                                      >  2:1162047/1‑150 (MQ=17)
           cAATTGATTGTGACAAAACTGGATCGCCTTGGTTGTAATGCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAatgatgatg                                                                                                                                      >  1:4805649/1‑150 (MQ=17)
               tgattgTGACAAAACTGGATCGCCTTGGTTGTAATGCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAg                                                                                                                                  <  2:3302000/150‑1 (MQ=17)
                 attgTGACAAAACTGGATCGCCTTGGTTGTAATGCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTa                                                                                                                                <  2:1308860/150‑1 (MQ=17)
                   tgtgACAAAACTGGATCGCCTTGGTTGTAATGCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAAt                                                                                                                              <  2:3241783/150‑1 (MQ=17)
                    gtgACAAAACTGGATCGCCTTGGTTGTAATGCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAAtt                                                                                                                             <  2:1912310/150‑1 (MQ=17)
                    gtgACAAAACTGGATCGCCTTGGTTGTAATGCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAAtt                                                                                                                             >  2:3659291/1‑150 (MQ=17)
                       aCAAAACTGGATCGCCTTGGTTGTAATGCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAATTTca                                                                                                                          <  2:382061/150‑1 (MQ=17)
                        cAAAACTGGATCGCCTTGTTTGTAATGCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTTCTTAGCATTGTGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAATTTcag                                                                                                                         >  2:1686102/1‑150 (MQ=255)
                          aaaCTGGATCGCCTTGGTTGTAATGCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAATTTcagca                                                                                                                       >  2:994373/1‑150 (MQ=17)
                          aaaCTGGATCGCCTTGGTTGTAATGCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAATTTcagca                                                                                                                       <  1:3629547/150‑1 (MQ=17)
                          aaaCTGGATCGCCTTGGTTGTAATGCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAATTTcagca                                                                                                                       <  2:2242111/150‑1 (MQ=17)
                          aaaCTGGATCGCCTTGGTTGTAATGCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAATTTcagca                                                                                                                       >  1:4651088/1‑150 (MQ=17)
                          aaaCTGGATCGCCTTGGTTGTAATGCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAATTTcagca                                                                                                                       >  1:4648331/1‑150 (MQ=17)
                          aaaCTGGATCGCCTTGGTTGTAATGCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAATTTcagca                                                                                                                       >  2:763549/1‑150 (MQ=17)
                           aaCTGGATCGCCTTGGTTGTAATGCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAATTTcagcag                                                                                                                      >  2:3526790/1‑150 (MQ=17)
                            aCTGGATCGCCTTGGTTGTAATGCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGt                                                                                                                     <  2:905064/150‑1 (MQ=17)
                               ggATCGCCTTGGTTGTAATGCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCGc                                                                                                                  >  2:3927262/1‑150 (MQ=17)
                                   cGCCTTGGTTGTAATGCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGaa                                                                                                              >  1:1677227/1‑150 (MQ=17)
                                     ccTTGGTTGTAATGCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAAtt                                                                                                            >  1:2596982/1‑150 (MQ=17)
                                     ccTTGGTTGTAATGCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAAtt                                                                                                            >  1:2406584/1‑150 (MQ=17)
                                       ttGGTTGTAATGCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTg                                                                                                          <  2:2791993/150‑1 (MQ=32)
                                       ttGGTTGTAATGCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTg                                                                                                          >  1:2947813/1‑150 (MQ=17)
                                         ggTTGTAATGCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGaa                                                                                                        <  2:3149603/150‑1 (MQ=17)
                                          gTTGTAATGCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCAGTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAAc                                                                                                       >  1:311043/1‑150 (MQ=17)
                                            tGTAATGCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACga                                                                                                     >  2:5075598/1‑150 (MQ=17)
                                              tAATGCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACgaga                                                                                                   <  1:2117806/150‑1 (MQ=17)
                                              tAATGCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACgaga                                                                                                   >  1:4519598/1‑150 (MQ=17)
                                               aaTGCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGTCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGAc                                                                                                  >  2:2075156/1‑150 (MQ=17)
                                                aTGCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGAcc                                                                                                 <  1:3915346/150‑1 (MQ=17)
                                                aTGCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGAcc                                                                                                 <  2:235211/150‑1 (MQ=17)
                                                   cAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCttt                                                                                              <  1:315846/150‑1 (MQ=17)
                                                    aaTGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCtttt                                                                                             >  2:1815853/1‑150 (MQ=17)
                                                    aaTGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCtttt                                                                                             <  1:2715741/150‑1 (MQ=17)
                                                        gATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGCTGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACtt                                                                                         <  1:4010144/150‑1 (MQ=17)
                                                          tatatggAAAACATTTGTACAACTGACCGAAACAGGTATCAGAGTTCATTGCTTAGCATTGGGGTGCATTGACCTGACCAGTCCAACAGAAAAAATGATTATGAAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGa                                                                                       <  1:715535/145‑1 (MQ=11)
                                                          tatCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGa                                                                                       <  1:1029577/150‑1 (MQ=17)
                                                          tatCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGa                                                                                       >  2:4714778/1‑150 (MQ=17)
                                                            tCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAAc                                                                                     >  1:3564285/1‑150 (MQ=17)
                                                             cAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACg                                                                                    >  2:3214062/1‑150 (MQ=17)
                                                              aGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGc                                                                                   >  1:3967961/1‑150 (MQ=17)
                                                              aGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGc                                                                                   >  1:4127918/1‑150 (MQ=17)
                                                                gAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCAc                                                                                 <  1:4998959/150‑1 (MQ=17)
                                                                gAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCAc                                                                                 >  1:4065773/1‑150 (MQ=17)
                                                                 aaaaCAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAACTGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACt                                                                                >  2:814113/1‑150 (MQ=17)
                                                                 aaaaCAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACt                                                                                <  2:1219824/150‑1 (MQ=17)
                                                                    aCAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCAt                                                                             >  2:4137761/1‑150 (MQ=17)
                                                                      aGTGGAACATCTGACCGAAACAGGTATCTGAGTGAATTGCTTAGCATTGGGGGGCATTGTCCTGACCAGTCCAACAGATAAAATGTTGATGAAAGTAATTTCAGAAGTCGCTGAATTTGAACGATAACTTTTACTTGAACGCACTCATTc                                                                           <  1:2935655/150‑1 (MQ=11)
                                                                      aGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTc                                                                           <  2:403526/150‑1 (MQ=17)
                                                                      aGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTc                                                                           >  1:57183/1‑150 (MQ=17)
                                                                      aGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTc                                                                           <  2:228691/150‑1 (MQ=17)
                                                                       gTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATTCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTcc                                                                          >  2:1182514/1‑150 (MQ=17)
                                                                         ggAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCgg                                                                        >  1:1885395/1‑150 (MQ=17)
                                                                                cTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGTa                                                                 <  2:2406584/150‑1 (MQ=17)
                                                                                 tGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGTaa                                                                >  2:3049481/1‑150 (MQ=17)
                                                                                  gACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGTAag                                                               >  2:4965058/1‑150 (MQ=17)
                                                                                  gACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGTAag                                                               >  1:2783649/1‑150 (MQ=17)
                                                                                  gACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGTAag                                                               >  1:5144971/1‑150 (MQ=17)
                                                                                   aCCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGTAaga                                                              <  1:1203229/150‑1 (MQ=17)
                                                                                      gAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGAGAGTAAGAGcc                                                           >  2:364757/1‑150 (MQ=17)
                                                                                       aaaCAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGTAAGAGccc                                                          <  1:2466961/150‑1 (MQ=17)
                                                                                       aaaCAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGTAAGAGccc                                                          <  1:4540835/150‑1 (MQ=17)
                                                                                       aaaCAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGTAAGAGccc                                                          <  2:1954489/150‑1 (MQ=17)
                                                                                       aaaCAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGTAAGAGccc                                                          >  2:5036589/1‑150 (MQ=17)
                                                                                        aaCAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGTAAGAGCCcg                                                         >  1:2711484/1‑150 (MQ=17)
                                                                                        aaCAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGTAAGAGCCcg                                                         <  1:115815/150‑1 (MQ=17)
                                                                                           aGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGTAAGAGCCCgcgg                                                      <  1:4398686/150‑1 (MQ=17)
                                                                                              tATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGTAAGAGCCCGCGgcgc                                                   >  1:4846540/1‑150 (MQ=17)
                                                                                                 cAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGTAAGAGCCCGCGGCGCAgg                                                <  1:5067618/150‑1 (MQ=17)
                                                                                                      tGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGTAAGAGCCCGCGGCGCAGGGAAAc                                           <  2:4558467/150‑1 (MQ=17)
                                                                                                           tGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGTAAGAGCCCGCGGCGCAGGGAAACGtttt                                      <  1:3741950/150‑1 (MQ=17)
                                                                                                            gCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGTAAGAGCCCGCGGCGCAGGGAAACGTTTTg                                     >  1:4375943/1‑150 (MQ=17)
                                                                                                            gCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGTAAGAGCCCGCGGCGCAGGGAAACGTTTTg                                     <  1:2625506/150‑1 (MQ=17)
                                                                                                            gCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGTAAGAGCCCGCGGCGCAGGGAAACGTTTTg                                     <  1:2018183/150‑1 (MQ=17)
                                                                                                                 gCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTGACTTGAACGCACTCATTCCGGGATAGTAAGAGCCCGCGGCGCAGGGAAACGTTTTGGTCGa                                >  2:2938621/1‑150 (MQ=17)
                                                                                                                  caggggggggCATTGACCTGACCAGTCCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGTAAGAGCCCGCGGCGCTGGGAAACGTTTTGGTCGac                               >  2:4322364/5‑150 (MQ=17)
                                                                                                                   aTTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGTAAGAGCCCGCGGCGCAGGGAAACGTTTTGGTCGacc                              >  1:2968279/1‑150 (MQ=17)
                                                                                                                    ttGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGTAAGAGCCCGCGGCGCAGGGAAACGTTTTGGTCGacna                             <  1:4221046/150‑1 (MQ=17)
                                                                                                                     tGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGTAAGAGCCCGCGGCGCAGGGAAACGTTTTGGTCGaccac                            >  1:1086972/1‑150 (MQ=18)
                                                                                                                     tGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGTAAGAGCCCGCGGCGCAGGGAAACGTTTTGGTCGaccac                            >  1:1000129/1‑150 (MQ=18)
                                                                                                                     tGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGTAAGAGCCCGCGGCGCAGGGAAACGTTTTGGTCGaccac                            <  1:4164433/150‑1 (MQ=18)
                                                                                                                        ggggCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGTAAGAGCCCGCGGCGCAGGGAAACGTTTTGGTCGACCACCtg                         <  1:3052137/150‑1 (MQ=18)
                                                                                                                          ggCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGTAAGAGCCCGCGGCGCAGGGAAACGTTTTGGTCGACCACCtgtg                       <  1:2098877/150‑1 (MQ=17)
                                                                                                                            cATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGTAAGAGCCCGCGGCGCAGGGAAACGTTTTGGTCGACCACCTGTGtt                     <  2:1052704/150‑1 (MQ=17)
                                                                                                                              ttGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGTAAGAGCCCGCGGCGCAGGGAAACGTTTTGGTCGACCACCTGTGTTaa                   <  1:2804845/150‑1 (MQ=17)
                                                                                                                                gacctgaccAGTCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGTAAGAGCCCGCGGCGCAGGGAAACGTTTTGGTCGACCACCTGTGTTAAAt                 >  1:825330/1‑150 (MQ=17)
                                                                                                                                  cctgaccAGTCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGTAAGAGCCCGCGGCGCAGGGAAACGTTTTGGTCGACCACCTGTGTTAAATga               >  2:2147397/1‑150 (MQ=17)
                                                                                                                                    tgaccAGTCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGTAAGAGCCCGCGGCGCAGGGAAACGTTTTGGTCGACCACCTGTGTTAAATgaag             >  2:3386484/1‑150 (MQ=17)
                                                                                                                                      accAGTCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGTAAGAGCCCGCGGCGCAGGGAAACGTTTTGGTCGACCACCTGTGTTAAATgaagaa           <  1:595927/150‑1 (MQ=17)
                                                                                                                                        cAGTCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGTAAGAGCCCGCGGCGCAGGGAAACGTTTTGGTCGACCACCTGTGTTAAATGAAGAACa         <  2:4648331/150‑1 (MQ=17)
                                                                                                                                         aGTCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGTAAGAGCCCGCGGCGCAGGGAAACGTTTTGGTCGACCACCTGTGTTAAATGAAGAACAg        <  2:3960033/150‑1 (MQ=17)
                                                                                                                                          gTCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGTAAGAGCCCGCGGCGCAGGGAAACGTTTTGGTCGACCACCTGTGTTAAATGAAGAACAGa       <  1:4288827/150‑1 (MQ=17)
                                                                                                                                           tCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGCTAGTAAGAGCCCGCGGCGCAGGGAAACGTTTTGGTCGACCACCTGTGTTCAATCAAGCACAGaa      >  2:4834101/1‑150 (MQ=11)
                                                                                                                                           tCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGTAAGAGCCCGCGGCGCAGGGAAACGTTTTGGTCGACCACCTGTGTTAAATGAAGAACAGaa      >  2:1609636/1‑150 (MQ=17)
                                                                                                                                           tCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGTAAGAGCCCGCGGCGCAGGGAAACGTTTTGGTCGACCACCTGTGTTAAATGAAGAACAGaa      <  1:3480126/150‑1 (MQ=17)
                                                                                                                                               aCAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTGTACTTGAACGCACTCATTCCGGGATAGTAAGAGCCCGCGGCGCAGGGAAACGTTTTGGTCGACCACCTGTGTTAAATGAAGAACAGAAACAg  <  2:5131713/150‑1 (MQ=17)
                                                                                                                                                   |                                                                                                                                                 
AATGTGGTGATCAATTGATTGTGACAAAACTGGATCGCCTTGGTTGTAATGCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGTAAGAGCCCGCGGCGCAGGGAAACGTTTTGGTCGACCACCTGTGTTAAATGAAGAACAGAAACAG  >  CP009273/1628060‑1628352

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: