Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,646,993 C→T 18.5% intergenic (+28/+19) intQ → / ← ydfJ pseudogene, Qin prophage; phage integrase family;Phage or Prophage Related/pseudogene, MFS transporter family; interrupted by Qin prophage;Phage or Prophage Related; putative transport protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,646,9930CT18.5% 661.3 / 121.5 254intergenic (+28/+19)intQ/ydfJpseudogene, Qin prophage; phage integrase family;Phage or Prophage Related/pseudogene, MFS transporter family; interrupted by Qin prophage;Phage or Prophage Related; putative transport protein
Reads supporting (aligned to +/- strand):  ref base C (117/89);  new base T (32/15);  total (149/105)
Fisher's exact test for biased strand distribution p-value = 1.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 255 total aligned reads displayed.

GGTATTTCAGGTATACGGAAAATGGATGTCTGAAAACAATAATGCACAGGTAGCTTTGTTAAATACACAGTTAAGCGAGTTTGCCCCAACCATGCCCCATAACGAAGCAATGAAAAATTAATTTAATATTTATCAAATAGTTAACACGCATGACTCTTGAAATCCATAAATTCAAGCGCAGTGCCCAGCCATCCCGATACTGCTGCTTTCACCAAATCCTTAGTGCTTCTTTCGTGTTTTTCTATTGTCATAATGGTTATCTCTAAAAAAGAGGTAAGATGCGTACTACTTACTCGCC  >  CP009273/1646845‑1647142
                                                                                                                                                    |                                                                                                                                                     
ggTATTTCAGGTATACGGAAAATGGATGTCTGAAAACAATAATGCACAGGTAGCTTTGTTAAATACACAGTTAAGCGAGTTTGCCCCAACCATGCCCCATAACGAAGCAATGAAAAATTAATTTAATATTTATCAAATAGTTAACACGCa                                                                                                                                                      >  1:3905391/1‑150 (MQ=255)
     ttCAGGTATACGGAAAATGGATGTCTGAAAACAATAATGCACAGGTAGCTTTGTTAAATACACAGTTAAGCGAGTTTGCCCCAACCATGCCCCATAACGAAGCAATGAAAAATTAATTTAATATTTATCAAATAGTTAACACGCATGAct                                                                                                                                                 <  2:934468/150‑1 (MQ=255)
       cAGGTATACGGAAAATGGATGTCTGAAAACAATAATGCACAGGTAGCTTTGTTAAATACACAGTTAAGCGAGTTTGCCCCAACCATGCCCCATAACGAAGCAATGAAAAATTAATTTAATATTTATCAAATAGTTAACACGCATGActct                                                                                                                                               >  1:2934400/1‑150 (MQ=255)
         ggTATACGGAAAATGGATGTCTGAAAACAATAATGCACAGGTAGCTTTGTTAAATACACAGTTAAGCGAGTTTGCCCCAACCATGCCCCATAACGAAGCAATGAAAAATTAATTTAATATTTATCAAATAGTTAACACGCATGACTCTTg                                                                                                                                             >  1:1171993/1‑150 (MQ=255)
              aCGGAAAATGGATGTCTGAAAACAATAATGCACAGGTAGCTTTGTTAAATACACAGTTAAGCGAGTTTGCCCCAACCATGCCCCATAACGAAGCAATGAAAAATTAATTTAATATTTATCAAATAGTTAACACGCATGACTCTTGAAATc                                                                                                                                        <  1:33985/150‑1 (MQ=255)
                 gAAAATGGATGTCTGAAAACAATAATGCACAGGTAGCTTTGTTAAATACACAGTTAAGCGAGTTTGCCCCAACCATGCCCCATAACGAAGCAATGAAAAATTAATTTAATATTTATCAAATAGTTAACACGCATGACTCTTGAAATCCAt                                                                                                                                     <  2:4587952/150‑1 (MQ=255)
                  aaaaTGGATGTCTGAAAACAATAATGCACAGGTAGCTTTGTTAAATACACAGTTAAGCGAGTTTGCCCCAACCATGCCCCATAACGAAGCAATGAAAAATTAATTTAATATTTATCAAATAGTTAACACGTATGACTCTTGAAATCCATa                                                                                                                                    <  2:1884781/150‑1 (MQ=255)
                      tGGATGTCTGAAAACAATAATGCACAGGTAGCTTTGTTAAATACACAGTTAAGCGAGTTTGCCCCAACCATGCCCCATAACGAAGCAATGAAAAATTAATTTAATATTTATCAAATAGTTAACACGCATGACTCTTGAAATCCATAAAtt                                                                                                                                <  1:2067894/150‑1 (MQ=255)
                       ggATGTCTGAAAACAATAATGCACAGGTAGCTTTGTTAAATACACAGTTAAGCGAGTTTGCCCCAACCTTGCCCCATAACGAAGCAATGAAAAATTAATTTAATATTTATCAAATAGTTAACACGCATGACTCTTGAAATCCATAAATTc                                                                                                                               <  2:3948283/150‑1 (MQ=255)
                            tCTGAAAACAATAATGCACAGGTAGCTTTGTTAAATACACAGTTAAGCGAGTTTGCCCCAACCATGCCCCATAACGAAGCAATGAAAAATTAATTTAATATTTATCAAATAGTTAACACGCATGACTCTTGAAATCCATAAATTCAAgcg                                                                                                                          >  1:821893/1‑150 (MQ=255)
                               gAAAACAATAATGCACAGGTAGCTTTGTTAAATACACAGTTAAGCGAGTTTGCCCCAACCATGCCCCATAACGAAGCAATGAAAAATTAATTTAATATTTATCAAATAGTTAACACGCATGACTCTTGAAATCCATAAATTCAAGCGCAg                                                                                                                       >  1:546299/1‑150 (MQ=255)
                                aaaaCAATAATGCACAGGTAGCTTTGTTAAATACACAGTTAAGCGAGTTTGCCCCAACCATGCCCCATAACGAAGCAATGAAAAATTAATTTAATATTTATCAAATAGTTAACACGCATGACTCTTGAAATCCATAAATTCAAGCGCAGt                                                                                                                      >  2:610700/1‑150 (MQ=255)
                                aaaaCAATAATGCACAGGTAGCTTTGTTAAATACACAGATAAGAGAGTTTGCCCCAACCATGACCCATAACGAAGCAATGAAAAATTAATATAATATTTATCAAAAAGTTAACACGCATGACTCTTGAAATCCATAAATTCAATCGCatt                                                                                                                      >  1:2866338/1‑148 (MQ=255)
                                  aaCAATAATGCACAGGTAGCTTTGTTAAATACACAGTTAAGCGAGTTTGCCCCAACCATGCCCCATAACGAAGCAATGAAAAATTAATTTAATATTTATCAAATAGTTAACACGTATGACTCTTGAAATCCATAAATTCAAGCGCAGTGc                                                                                                                    <  2:3008720/150‑1 (MQ=255)
                                  aaCAATAATGCACAGGTAGCTTTGTTAAATACACAGTTAAGCGAGTTTGCCCCAACCATGCCCCATAACGAAGCAATGAAAAATTAATTTAATATTTATCAAATAGTTAACACGTATGACTCTTGAAATCCATAAATTCAAGCGCAGTGc                                                                                                                    <  2:3008854/150‑1 (MQ=255)
                                  aaCAATAATGCACAGGTAGCTTTGTTAAATACACAGTTAAGCGAGTTTGCCCCAACCATGCCCCATAACGAAGCAATGAAAAATTAATTTAATATTTATCAAATAGTTAACACGCATGACTCTTGAAATCCATAAATTCAAGCGCAGTGc                                                                                                                    >  1:4068851/1‑150 (MQ=255)
                                    cAATAATGCACAGGTAGCTTTGTTAAATACACAGTTAAGCGAGTTTGCCCCAACCATGCCCCATAACGAAGCAATGAAAAATTAATTTAATATTTATCAAATAGTTAACACGTATGACTCTTGAAATCCATAAATTCAAGCGCAGTGccc                                                                                                                  >  2:3344106/1‑150 (MQ=255)
                                     aataatGCACAGGTAGCTTTGTTAAATACACAGTTAAGCGAGTTTGCCCCAACCATGCCCCATAACGAAGCAATGAAAAATTAATTTAATATTTATCAAATAGTTAACACGCATGACTCTTGAAATCCATAAATTCAAGCGCAGTGCCCa                                                                                                                 >  2:2164440/1‑150 (MQ=255)
                                       taatGCACAGGTAGCTTTGTTAAATACACAGTTAAGCGAGTTTGCCCCAACCATGCCCCATAACGAAGCAATGAAAAATTAATTTAATATTTATCAAATAGTTAACACGCATGACTCTTGAAATCCATAAATTCAAGCGCAGTGCCCAGc                                                                                                               >  2:3811553/1‑150 (MQ=255)
                                        aatGCACAGGTAGCTTTGTTAAATACACAGTTAAGCGAGTTTGCCCCAACCATGCCCCATAACGAAGCAATGAAAAATTAATTTAATATTTATCAAATAGTTAACACGCATGACTCTTGAAATCCATAAATTCAAGCGCAGTGCCCAGcc                                                                                                              <  1:1461600/150‑1 (MQ=255)
                                        aatGCACAGGTAGCTTTGTTAAATACACAGTTAAGCGAGTTTGCCCCAACCATGCCCCATAACGAAGCAATGAAAAATTAATTTAATATTTATCAAATAGTTAACACGCATGACTCTTGAAATCCATAAATTCAAGCGCAGTGCCCAGcc                                                                                                              <  1:3358732/150‑1 (MQ=255)
                                          tGCACAGGTAGCTTTGTTAAATACACAGTTAAGCGAGTTTGCCCCAACCATGCCCCATAACGAAGCAATGAAAAATTAATTTAATATTTATCAAATAGTTAACACGTATGACTCTTGAAATCCATAAATTCAAGCGCAGTGCCCAGCCAt                                                                                                            >  2:1185000/1‑150 (MQ=255)
                                           gCACAGGTAGCTTTGTTAAATACACAGTTAAGCGAGTTTGCCCCAACCATGCCCCATAACGAAGCAATGAAAAATTAATTTAATATTTATCAAATAGTTAACACGCATGACTCTTGAAATCCATAAATTCAAGCGCAGTGCCCAGCCATc                                                                                                           <  1:2818468/150‑1 (MQ=255)
                                           gCACAGGTAGCTTTGTTAAATACACAGTTAAGCGAGTTTGCCCCAACCATGCCCCATAACGAAACAATGAAAAATTAATTTAATATTTATCAAATAGTTAACACGCATGACTCTTGAAATCCATAAATTCAAGCGCAGTGCCCAGCCATc                                                                                                           >  2:1168111/1‑150 (MQ=255)
                                            cacaGGTAGCTTTGTTAAATACACAGTTAAGCGAGTTTGCCCCAACCATGCCCCATAACGAAGCAATGAAAAATTAATTTAATATTTATCAAATAGTTAACACGCATGACTCTTGAAATCCATAAATTCAAGCGCAGTGCCCAGCCATcc                                                                                                          <  1:2593262/150‑1 (MQ=255)
                                               aGGTAGCTTTGTTAAATACACAGTTAAGCGAGTTTGCCCCAACCATGCCCCATAACGAAGCAATGAAAAATTAATTTAATATTTATCAAATAGTTAACACGCATGACTCTTGAAATCCATAAATTCAAGCGCAGTGCCCAGCCATCCCGa                                                                                                       <  2:2308562/150‑1 (MQ=255)
                                                 gTAGCTTTGTTAAATACACAGTTAAGCTAGTTTGCCCCAACCATGCCCCATACCGAAGCAATGAAAAATTAATTTAATATTTATCAAATAGTTAACACGCATGACTCTTGAAATCCATAAATTCAAGCGCAGTGCCCAGCCATCCCGATa                                                                                                     <  2:2583809/150‑1 (MQ=255)
                                                 gTAGCTTTGTTAAATACACAGTTAAGCGAGTTTGCCCCAACCATGCCCCATAACGAAGCAATGAAAAATTAATTTAATATTTATCAAATAGTTAACACGCATGACTCTTGAAATCCATAAATTCAAGCGCAGTGCCCAGCCATCCCGATa                                                                                                     <  2:2016632/150‑1 (MQ=255)
                                                 gTAGCTTTGTTAAATACACAGTTAAGCGAGTTTGCCCCAACCATGCCCCATAACGAAGCAATGAAAAATTAATTTAATATTTATCAAATAGTTAACACGCATGACTCTTGAAATCCATAAATTCAAGCGCAGTGCCCAGCCATCCCGATa                                                                                                     >  2:2717807/1‑150 (MQ=255)
                                                   aGCTTTGTTAAATACACAGTTAAGCGAGTTTGCCCCAACCATGCCCCATAACGAAGCAATGAAAAATTAATTTAATATTTATCAAATAGTTAACACGCATGACTCTTGAAATCCATAAATTCAAGCGCAGTGCCCAGCCATCCCGATAct                                                                                                   <  1:3013239/150‑1 (MQ=255)
                                                   aGCTTTGTTAAATACACAGTTAAGCGAGTTTGCCCCAACCATGCCCCATAACGAAGCAATGAAAAATTAATTTAATATTTATCAAATAGTTAACACGCATGACTCTTGAAATCCATAAATTCAAGCGCAGTGCCCAGCCATCCCGATAct                                                                                                   >  2:4228532/1‑150 (MQ=255)
                                                      tttGTTAAATACACAGTTAAGCGAGTTTGCCCCAACCATGCCCCATAACGAAGCAATGAAAAATTAATTTAATATTTATCAAATAGTTAACACGCATGACTCTTGAAATCCATAAATTCAAGCGCAGTGCCCAGCCATCCCGATActgct                                                                                                <  2:1103272/150‑1 (MQ=255)
                                                      tttGTTAAATACACAGTTAAGCGAGTTTGCCCCAACCATGCCCCATAACGAAGCAATGAAAAATTAATTTAATATTTATCAAATAGTTAACACGCATGACTCTTGAAATCCATAAATTCAAGCGCAGTGCCCAGCCATCCCGATActgct                                                                                                <  2:3794661/150‑1 (MQ=255)
                                                           tAAATACACAGTTAAGCGAGTTTGCCCCAACCATGCCCCATAACGAAGCAATGAAAAATTAATTTAATATTTATCAAATAGTTAACACGTATGACTCTTGAAATCCATAAATTCAAGCGCAGTGCCCAGCCATCCCGATACTGCTGCttt                                                                                           >  1:496665/1‑150 (MQ=255)
                                                                  acaGTTAAGCGAGTTTGCCCCAACCATGCCCCATAACGAAGCAATGAAAAATTAATTTAATATTTATCAAATAGTTAACACGCATGACTCTTGAAATCCATAAATTNAAGCGCAGTNCCCAGCCATCCCGATACTGCTGCTTTCACCaaa                                                                                    >  1:3522346/1‑150 (MQ=255)
                                                                   caGTTAAGCGAGTTTGCCCCAACCATGCCCCATAACGAAGCAATGAAAAATTAATTTAATATTTATCAAATAGTTAACACGCATGACTCTTGAAATCCATAAATTCAAGCGCAGTGCCCAGCCATCCCGATACTGCTGCTTTCACCAAAt                                                                                   <  1:314493/150‑1 (MQ=255)
                                                                    aGTTAAGCGAGTTTGCCCCAACCATGCCCCATAACGAAGCAATGAAAAATTAATTTAATATTTATCAAATAGTTAACACGCATGACTCTTGAAATCCATAAATTCAAGCGCAGTGCCCAGCCATCCCGATACTGCTGCTTTCACCAAATc                                                                                  >  2:1549157/1‑150 (MQ=255)
                                                                     gTTAAGCGAGTTTGCCCCAACCATGCCCCATAACGAAGCAATGAAAAATTAATTTAATATTTATCAAATAGTTAACACGCATGACTCTTGAAATCCATAAATTCAAGCGCAGTGCCCAGCCATCCCGATACTGCTGCTTTCACCAAATcc                                                                                 >  1:4378906/1‑150 (MQ=255)
                                                                        aaGCGAGTTTGCCCCAACCATGCCCCATAACGAAGCAATGAAAAATTAATTTAATATTTATCAAATAGTTAACACGCATGACTCTTGAAATCCATAAATTCAAGCGCAGTGCCCAGCCATCCCGATACTGCTGCTTTCACCAAATCCTTa                                                                              <  1:2364174/150‑1 (MQ=255)
                                                                          gCGAGTTTGCCCCAACCATGCCCCATAACGAAGCAATGAAAAATTAATTTAATATTTATCAAATAGTTAACACGCATGACTCTTGAAATCCATAAATTCAAGCGCAGTGCCCAGCCATCCCGATACTGCTGCTTTCACCAAATCCTTAGt                                                                            <  1:4531499/150‑1 (MQ=255)
                                                                           cGAGTTTGCCCCAACCATGCCCCATAACGAAGCAATGAAAAATTAATTTAATATTTATCAAATAGTTAACACGCATGACTCTTGAAATCCATAAATTCAAGCGCAGTGCCCAGCCATCCCGATACTGCTGCTTTCACCAAATCCTTAGTg                                                                           >  1:2146837/1‑150 (MQ=255)
                                                                            gAGTTTGCCCCAACCATGCCCCATAACGAAGCAATGAAAAATTAATTTAATATTTATCAAATAGTTAACACGCATGACTCTTGAAATCCATAAATTCAAGCGCAGTGCCCAGCCATCCCGATACTGCTGCTTTCACCAAATCCTTAGTGc                                                                          >  1:3479397/1‑150 (MQ=255)
                                                                            gAGTTTGCCCCAACCATGCCCCATAACGAAGCAATGAAAAATTAATTTAATATTTATCAAATAGTTAACACGCATGACTCTTGAAATCCATAAATTCAAGCGCAGTGCCCAGCCATCCCGATACTGCTGCTTTCACCAAATCCTTAGTGc                                                                          >  1:952000/1‑150 (MQ=255)
                                                                                  gCCCCAACCATGCCCCATAACGAAGCAATGAAAAATTAATTTAATATTTATCAAATAGTTAACACGCATGACTCTTGAAATCCATAAATTCAAGCGCAGTGCCCAGCCATCCCGATACTGCTGCTTTCACCAAATCCTTAGTGCTTCttt                                                                    <  2:102022/150‑1 (MQ=255)
                                                                                   ccccAACCATGCCCCATAACGAAGCAATGAAAAATTAATTTAATATTTATCAAATAGTTAACACGCATGACTCTTGAAATCCATAAATTCAAGCGCAGTGCCCAGCCATCCCGATACTGCTGCTTTCACCAAATCCTTAGTGCTTCTTTc                                                                   >  1:2887404/1‑150 (MQ=255)
                                                                                    cccAACCATGCCCCATAACGAAGCAATGAAAAATTAATTTAATATTTATCAAATAGTTAACACGCATGACTCTTGAAATCCATAAATTCAAGCGCAGTGCCCAGCCATCCCGATACTGCTGCTTTCACCAAATCCTTAGTGCTTCTTTCg                                                                  >  1:3096293/1‑150 (MQ=255)
                                                                                      cAACCATGCCCCATAACGAAGCAATGAAAAATTAATTTAATATTTATCAAATAGTTAACACGCATGACTCTTGAAATCCATAAATTCAAGCGCAGTGCCCAGCCATCCCGATACTGCTGCTTTCACCAAATCCTTAGTGCTTCTTTCgtg                                                                <  2:3649863/150‑1 (MQ=255)
                                                                                            tGCCCCATAACGAAGCAATGAAAAATTAATTTAATATTTATCAAATAGTTAACACGCATGACTCTTGAAATCCATAAATTCAAGCGCAGTGCCCAGCCATCCCGATACTGCTGCTTTCACCAAATCCTTAGTGCTTCTTTCGTGTTTTTc                                                          >  1:1473703/1‑150 (MQ=255)
                                                                                            tGCCCCATAACGAAGCAATGAAAAATTAATTTAATATTTATCAAATAGTTAACACGCATGACTCTTGAAATCCATAAATTCAAGCGCAGTGCCCAGCCATCCCGATACTGCTGCTTTCACCAAATCCTTAGTGCTTCTTTCGTGTTTTTc                                                          >  1:3442964/1‑150 (MQ=255)
                                                                                               cccATAACGAAGCAATGAAAAATTAATTTAATATTTATCAAATAGTTAACACGTATGACTCTTGAAATCCATAAATTCAAGCGCAGTGCCCAGCCATCCCGATACTGCTGCTTTCACCAAATCCTTAGTGCTTCTTTCGTGTTTTTCTAt                                                       >  2:3206333/1‑150 (MQ=255)
                                                                                               cccATAACGAAGCAATGAAAAATTAATTTAATATTTATCAAATAGTTAACACGCATGACTCTTGAAATCCATAAATTCAAGCGCAGTGCCCAGCCATCCCGATACTGCTGCTTTCACCAAATCCTTAGTGCTTCTTTCGTGTTTTTCTAt                                                       <  2:631538/150‑1 (MQ=255)
                                                                                                  aTAACGAAGCAATGAAAAATTAATTTAATATTTATCAAATAGTTAACACGCATGACTCTTGAAATCCATAAATTCAAGCGCAGTGCCCAGCCATCCCGATACTGCTGCTTTCACCAAATCCTTAGTGCTTCTTTCGTGTTTTTCTATTGt                                                    <  2:4368740/150‑1 (MQ=255)
                                                                                                      cGAAGCAATGAAAAATTAATTTAATATTTATCAAATAGTTAACACGCATGACTCTTGAAATCCATAAATTCAAGCGCAGTGCCCAGCCATCCCGATACTGCTGCTTTCACCAAATCCTTAGTGCTTCTTTCGTGTTTTTCTATTGTCATa                                                <  1:2039569/150‑1 (MQ=255)
                                                                                                      cGAAGCAATGAAAAATTAATTTAATATTTATCAAATAGTTAACACGCATGACTCTTGAAATCCATAAATTCAAGCGCAGTGCCCAGCCATCCCGATACTGCTGCTTTCACCAAATCCTTAGTGCTTCTTTCGTGTTTTTCTATTGTCATa                                                >  1:1849513/1‑150 (MQ=255)
                                                                                                      cGAAGCAATGAAAAATTAATTTAATATTTATCAAATAGTTAACACGCATGACTCTTGAAATCCATAAATTCAAGCGCAGTGCCCAGCCATCCCGATACTGCTGCTTTCACCAAATCCTTAGTGCTTCTTTCGTGTTTTTCTATTGTCATa                                                >  1:4112479/1‑150 (MQ=255)
                                                                                                         aGCAATGAAAAATTAATTTAATATTTATCAAATAGTTAACACGTATGACTCTTGAAATCCATAAATTCAAGCGCAGTGCCCAGCCATCCCGATACTGCTGCTTTCACCAAATCCTTAGTGCTTCTTTCGTGTTTTTCTATTGTCATAATg                                             <  1:353915/150‑1 (MQ=255)
                                                                                                         aGCAATGAAAAATTAATTTAATATTTATCAAATAGTTAACACGCATGACTCTTGCAATCCATAAATTCAAGCGCAGTGCCCAGCCATCCCGATACTGCTGCTTTCACCAAATCCTTAGTGCTTCTTTCGTGTTTTTCTATTGTCATAATg                                             <  2:2929804/150‑1 (MQ=255)
                                                                                                         aGCAATGAAAAATTAATTTAATATTTATCAAATAGTTAACACGCATGACTCTTGAAATCCATAAATTCAAGCGCAGTGCCCAGCCATCCCGATACTGCTGCTTTCACCAAATCCTTAGTGCTTCTTTCGTGTTTTTCTATTGTCATAATg                                             <  2:1171927/150‑1 (MQ=255)
                                                                                                         aGCAATGAAAAATTAATTTAATATTTATCAAATAGTTAACACGCATGACTCTTGAAATCCATAAATTCAAGCGCAGTGCCCAGCCATCCCGATACTGCTGCTTTCACCAAATCCTTAGTGCTTCTTTCGTGTTTTTCTATTGTCATAATg                                             >  2:3270483/1‑150 (MQ=255)
                                                                                                          gCAATGAAAAATTAATTTAATATTTATCAAATAGTTAACACGCATGACTCTTGAAATCCATAAATTCAAGCGCAGTGCCCAGCCATCCCGATACTGCTGCTTTCACCAAATCCTTAGTGCTTCTTTCGTGTTTTTCTATTGTCATAATgg                                            >  2:31388/1‑150 (MQ=255)
                                                                                                          gCAATGAAAAATTAATTTAATATTTATCAAATAGTTAACACGCATGACTCTTGAAATCCATAAATTCAAGCGCAGTGCCCAGCCATCCCGATACTGCTGCTTTCACCAAATCCTTAGTGCTTCTTTCGTGTTTTTCTATTGTCATAATgg                                            >  1:1130785/1‑150 (MQ=255)
                                                                                                            aaTGAAAAATTAATTTAATATTTATCAAATAGTTAACACGCATGACTCTTGAAATCCATAAATTCAAGCGCAGTGCCCAGCCATCCCGATACTGCTGCTTTCACCAAATCCTTAGTGCTTCTTTCGTGTTTTTCTATTGTCATAATGGtt                                          >  2:1763190/1‑150 (MQ=255)
                                                                                                             aTGAAAAATTAATTTAATATTTATCAAATAGTTAACACGTATGACTCTTGAAATCCATAAATTCAAGCGCAGTGCCCAGCCATCCCGATACTGCTGCTTTCACCAAATCCTTAGTGCTTCTTTCGTGTTTTTCTATTGTCATAATGGTTa                                         >  2:3481360/1‑150 (MQ=255)
                                                                                                             aTGAAAAATTAATTTAATATTTATCAAATAGTTAACACGTATGACTCTTGAAATCCATAAATTCAAGCGCAGTGCCCAGCCATCCCGATACTGCTGCTTTCACCAAATCCTTAGTGCTTCTTTCGTGTTTTTCTATTGTCATAATGGTTa                                         >  2:3576561/1‑150 (MQ=255)
                                                                                                             aTGAAAAATTAATTTAATATTTATCAAATAGTTAACACGCATGACTCTTGAAATCCATAAATTCAAGCGCAGTGCCCAGCCATCCCGATACTGCTGCTTTCACCAAATCCTTAGTGCTTCTTTCGTGTTTTTCTATTGTCATAATGGTTa                                         <  2:2413095/150‑1 (MQ=255)
                                                                                                             aTGAAAAATTAATTTAATATTTATCAAATAGTTAACACGCATGACTCTTGAAATCCATAAATTCAAGCGCAGTGCCCAGCCATCCCGATACTGCTGCTTTCACCAAATCCTTAGTGCTTCTTTCGTGTTTTTCTATTGTCATAATGGTTa                                         <  2:1493174/150‑1 (MQ=255)
                                                                                                               gAAAAATTAATTTAATATTTATCAAATAGTTAACACGTATGACTCTTGAAATCCATAAATTCAAGCGCAGTGCCCAGCCATCCCGATACTGCTGCTTTCACCAAATCCTTAGTGCTTCTTTCGTGTTTTTCTATTGTCATAATGGTTAtc                                       <  1:3662802/150‑1 (MQ=255)
                                                                                                                aaaaaTTAATTTAATATTTATCAAATAGTTAACACGTATGACTCTTGAAATCCATAAATTCAAGCGCAGTGCCCAGCCATCCCGATACTGCTGCTTTCACCAAATCCTTAGTGCTTCTTTCGTGTTTTTCTATTGTCATAATGGTTAtct                                      >  1:3958250/1‑150 (MQ=255)
                                                                                                                 aaaaTTAATTTAATATTTATCAAATAGTTAACACGCATGACTCTTGAAATCCATAAATTCAAGCGCAGTGCCCAGCCATCCCGATACTGCTGCTTTCACCAAATCCTTAGTGCTTCTTTCGTGTTTTTCTATTGTCATAATGGTTAtctc                                     >  2:1820937/1‑150 (MQ=255)
                                                                                                                 aaaaTTAATTTAATATTTATCAAATAGTTAACACGCATGACTCTTGAAATCCATAAATTCAAGCGCAGTGCCCAGCCATCCCGATACTGCTGCTTTCACCAAATCCTTAGTGCTTCTTTCGTGTTTTTCTATTGTCATAATGGTTAtctc                                     >  2:2966798/1‑150 (MQ=255)
                                                                                                                  aaaTTAATTTAATATTTATCAAATAGTTAACACGCATGACTCTTGAAATCCATAAATTCAAGCGCAGTGCCCAGCCATCCCGATACTGCTGCTTTCACCAAATCCTTAGTGCTTCTTTCGTGTTTTTCTATTGTCATAATGGTTATctct                                    <  2:3103051/150‑1 (MQ=255)
                                                                                                                    attaattTAATATTTATCAAATAGTTAACACGCATGACTCTTGAAATCCATAAATTCAAGCGCAGTGCCCAGCCATCCCGATACTGCTGCTTTCACCAAATCCTTAGTGCTTCTTTCGTGTTTTTCTATTGTCATAATGGTTATCTCTaa                                  >  1:1925379/1‑150 (MQ=255)
                                                                                                                       aatttaatATTTATCAAATAGTTAACACGCATGACTCTTGAAATCCATAAATTCAAGCGCAGTGCCCAGCCATCCCGATACTGCTGCTTTCACCAAATCCTTAGTGCTTCTTTCGTGTTTTTCTATTGTCATAATGGTTATCTCTaaaaa                               >  2:3260327/1‑150 (MQ=255)
                                                                                                                        atttaatATTTATCAAATAGTTAACACGTATGACTCTTGAAATCCATAAATTCAAGCGCAGTGCCCAGCCATCCCGATACTGCTGCTTTCACCAAATCCTTAGTGCTTCTTTCGTGTTTTTCTATTGTCATAATGGTTATCTCTaaaaaa                              >  2:2152172/1‑150 (MQ=255)
                                                                                                                         tttaatATTTATCAAATAGTTAACACGTATGACTCTTGAAATCCATAAATTCAAGCGCAGTGCCCAGCCATCCCGATACTGCTGCTTTCACCAAATCCTTAGTGCTTCTTTCGTGTTTTTCTATTGTCATAATGGTTATCTCTAAAAAag                             >  1:2946338/1‑150 (MQ=255)
                                                                                                                             atatTTATCAAATAGTTAACACGCATGACTCTTGAAATCCATAAATTCAAGCGCAGTGCCCAGCCATCCCGATACTGCTGCTTTCACCAAATCCTTAGTGCTTCTTTCGTGTTTTTCTATTGTCATAATGGTTATCTCTAAAAAAGAGGt                         <  2:2866090/150‑1 (MQ=255)
                                                                                                                             atatTTATCAAATAGTTAACACGCATGACTCTTGAAATCCATAAATTCAAGCGCAGTGCCCAGCCATCCCGATACTGCTGCTTTCACCAAATCCTTAGTGCTTCTTTCGTGTTTTTCTATTGTCATAATGGTTATCTCTAAAAAAGAGGt                         >  1:1592902/1‑150 (MQ=255)
                                                                                                                             atatTTATCAAATAGTTAACACGCATGACTCTTGAAATCCATAAATTCAAGCGCAGTGCCCAGCCATCCCGATACTGCTGCTTTCACCAAATCCTTAGTGCTTCTTTCGTGTTTTTCTATTGTCATAATGGTTATCTCTAAAAAAGAGGt                         <  2:2763611/150‑1 (MQ=255)
                                                                                                                             atatTTATCAAATAGTTAACACGCATGACTCTTGAAATCCATAAATTCAAGCGCAGTGCCCAGCCATCCCGATACTGCTGCTTTCACCAAATCCTTAGTGCTTCTTTCGTGTTTTTCTATTGTCATAATGGTTATCTCTAAAAAAGAGGt                         <  1:2883943/150‑1 (MQ=255)
                                                                                                                             atatTTATCAAATAGTTAACACGCATGACTCTTGAAATCCATAAATTCAAGCGCAGTGCCCAGCCATCCCGATACTGCTGCTTTCACCAAATCCTTAGTGCTTCTTTCGTGTTTTTCTATTGTCATAATGGTTATCTCTAAAAAAGAGGt                         >  1:3325734/1‑150 (MQ=255)
                                                                                                                              tatTTATCAAATAGTTAACACGCATGACTCTTGAAATCCATAAATTCAAGCGCAGTGCCCAGCCATCCCGATACTGCTGCTTTCACCAAATCCTTAGTGCTTCTTTCGTGTTTTTCTATTGTCATAATGGTTATCTCTAAAAAAGAGGTa                        >  2:4111654/1‑150 (MQ=255)
                                                                                                                               atTTATCAAATAGTTAACACGCATGACTCTTGAAATCCATAAATTCAAGCGCAGTGCCCAGCCATCCCGATACTGCTGCTTTCACCAAATCCTTAGTGCTTCTTTCGTGTTTTTCTATTGTCATAATGGTTATCTCTAAAAAAGAGGTaa                       <  1:3990661/150‑1 (MQ=255)
                                                                                                                                 ttATCAAATAGTTAACACGTATGACTCTTGAAATCCATAAATTCAAGCGCAGTGCCCAGCCATCCCGATACTGCTGCTTTCACCAAATCCTTAGTGCTTCTTTCGTGTTTTTCTATTGTCATAATGGTTATCTCTAAAAAAGAGGTAAGa                     >  1:2202273/1‑150 (MQ=255)
                                                                                                                                  tATCAAATAGTTAACACGCATGACTCTTGAAATCCATAAATTCAAGCGCAGTGCCCAGCCATCCCGATACTGCTGCTTTCACCAAATCCTTAGTGCTTCTTTCGTGTTTTTCTATTGTCATAATGGTTATCTCTAAAAAAGAGGTAAGAt                    >  1:3168525/1‑150 (MQ=255)
                                                                                                                                   aTCAAATAGTTAACACGCATGACTCTTGAAATCCATAAATTCAAGCGCAGTGCCCAGCCATCCCGATACTGCTGCTTTCACCAAATCCTTAGTGCTTCTTTCGTGTTTTTCTATTGTCATAATGGTTATCTCTAAAAAAGAGGTAAGATg                   >  1:4489190/1‑150 (MQ=255)
                                                                                                                                   aTCAAATAGTTAACACGCATGACTCTTGAAATCCATAAATTCAAGCGCAGTGCCCAGCCATCCCGATACTGCTGCTTTCACCAAATCCTTAGTGCTTCTTTCGTGTTTTTCTATTGTCATAATGGTTATCTCTAAAAAAGAGGTAAGATg                   >  2:2323688/1‑150 (MQ=255)
                                                                                                                                    tCAAATAGTTAACACGCATGACTCTTGAAATCCATAAATTCAAGCGCAGTGCCCAGCCATCCCGATACTGCTGCTTTCACCAAATCCTTAGTGCTTCTTTCGTGTTTTTCTATTGTCATAATGGTTATCTCTAAAAAAGAGGTAAGATGc                  <  2:3610114/150‑1 (MQ=255)
                                                                                                                                        aTAGTTAACACGCATGACTCTTGAAATCCATAAATTCAAGCGCAGTGCCCAGCCATCCCGATACTGCTGCTTTCACCAAATCCTTAGTGCTTCTTTCGTGTTTTTCTATTGTCATAATGGTTATCTCTAAAAAAGAGGTAAGATGCGtac              >  2:4398277/1‑150 (MQ=255)
                                                                                                                                         tAGTTAACACGTATGACTCTTGAAATCCATAAATTCAAGCGCAGTGCCCAGCCATCCCGATACTGCTGCTTTCACCAAATCCTTAGTGCTTCTTTCGTGTTTTTCTATTGTCATAATGGTTATCTCTAAAAAAGAGGTAAGATGCGtact             >  2:4499742/1‑150 (MQ=255)
                                                                                                                                          aGTTAACACGCATGACTCTTGAAATCCATAAATTCAAGCGCAGTGCCCAGCCATCCCGATACTGCTGCTTTCACCAAATCCTTAGTGCTTCTTTCGTGTTTTTCTATTGTCATAATGGTTATCTCTAAAAAAGAGGTAAGATGCGtacta            >  2:1430364/1‑150 (MQ=255)
                                                                                                                                          aGTTAACACGCATGACTCTTGAAATCCATAAATTCAAGCGCAGTGCCCAGCCATCCCGATACTGCTGCTTTCACCAAATCCTTAGTGCTTCTTTCGTGTTTTTCTATTGTCATAATGGTTATCTCTAAAAAAGAGGTAAGATGCGtacta            >  2:1623167/1‑150 (MQ=255)
                                                                                                                                          aGTTAACACGCATGACTCTTGAAATCCATAAATTCAAGCGCAGTGCCCAGCCATCCCGATACTGCTGCTTTCACCAAATCCTTAGTGCTTCTTTCGTGTTTTTCTATTGTCATAATGGTTATCTCTAAAAAAGAGGTAAGATGCGtacta            >  2:2228042/1‑150 (MQ=255)
                                                                                                                                           gTTAACACGCATGACTCTTGAAATCCATAAATTCAAGCGCAGTGCCCAGCCATCCCGATACTGCTGCTTTCACCAAATCCTTAGTGCTTCTTTCGTGTTTTTCTATTGTCATAATGGTTATCTCTAAAAAAGAGGTAAGATGCGtactac           >  1:3957279/1‑150 (MQ=255)
                                                                                                                                            ttAACACGTATGACTCTTGAAATCCATAAATTCAAGCGCAGTGCCCAGCCATCCCGATACTGCTGCTTTCACCAAATCCTTAGTGCTTCTTTCGTGTTTTTCTATTGTCATAATGGTTATCTCTAAAAAAGAGGTAAGATGCGTactact          <  2:496665/150‑1 (MQ=255)
                                                                                                                                            ttAACACGCATGACTCTTGAAATCCATAAATTCAAGCGCAGTGCCCAGCCATCCCGATACTGCTGCTTTCACCAAATCCTTAGTGCTTCTTTCGTGTTTTTCTATTGTCATAATGGTTATCTCTAAAAAAGAGGTAAGATGCGTactact          >  2:2440314/1‑150 (MQ=255)
                                                                                                                                            ttAACACGCATGACTCTTGAAATCCATAAATTCAAGCGCAGTGCCCAGCCATCCCGATACTGCTGCTTTCACCAAATCCTTAGTGCTTCTTTCGTGTTTTTCTATTGTCATAATGGTTATCTCTAAAAAAGAGGTAAGATGCGTactact          <  1:4382860/150‑1 (MQ=255)
                                                                                                                                             tAACACGCATGACTCTTGAAATCCATAAATTCAAGCGCAGTGCCCAGCCATCCCGATACTGCTGCTTTCACCAAATCCTTAGTGCTTCTTTCGTGTTTTTCTATTGTCATAATGGTTATCTCTAAAAAAGAGGTAAGATGCGTACtactt         >  2:2769737/1‑150 (MQ=255)
                                                                                                                                               acacGCATGACTCTTGAAATCCATAAATTCAAGCGCAGTGCCCAGCCATCCCGATACTGCTGCTTTCACCAAATCCTTAGTGCTTCTTTCGTGTTTTTCTATTGTCATAATGGTTATCTCTAAAAAAGAGGTAAGATGCGTACtacttac       >  1:1770059/1‑150 (MQ=255)
                                                                                                                                                 acGCATGACTCTTGAAATCCATAAATTCAAGCGCAGTGCCCAGCCATCCCGATACTGCTGCTTTCACCAAATCCTTAGTGCTTCTTTCGTGTTTTTCTATTGTCATAATGGTTATCTCTAAAAAAGAGGTAAGATGCGTACTACTTACTc     <  1:1168740/150‑1 (MQ=255)
                                                                                                                                                    cATGACTCTTGAAATCCATAAATTCAAGCGCAGTGCCCAGCCATCCCGATACTGCTGCTTTCACCAAATCCTTAGTGCTTCTTTCGTGTTTTTCTATTGTCATAATGGTTATCTCTAAAAAAGAGGTAAGATGCGTACTACTTACTCGcc  <  1:1020407/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GGTATTTCAGGTATACGGAAAATGGATGTCTGAAAACAATAATGCACAGGTAGCTTTGTTAAATACACAGTTAAGCGAGTTTGCCCCAACCATGCCCCATAACGAAGCAATGAAAAATTAATTTAATATTTATCAAATAGTTAACACGCATGACTCTTGAAATCCATAAATTCAAGCGCAGTGCCCAGCCATCCCGATACTGCTGCTTTCACCAAATCCTTAGTGCTTCTTTCGTGTTTTTCTATTGTCATAATGGTTATCTCTAAAAAAGAGGTAAGATGCGTACTACTTACTCGCC  >  CP009273/1646845‑1647142

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: