Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,646,993 C→T 46.8% 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,9930CT46.8% 76.5 / inf 233intergenic (+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 (62/62);  new base T (71/38);  total (133/100)
Fisher's exact test for biased strand distribution p-value = 2.42e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.74e-01

Only 100 of 233 total aligned reads displayed.

GGTATTTCAGGTATACGGAAAATGGATGTCTGAAAACAATAATGCACAGGTAGCTTTGTTAAATACACAGTTAAGCGAGTTTGCCCCAACCATGCCCCATAACGAAGCAATGAAAAATTAATTTAATATTTATCAAATAGTTAACACGCATGACTCTTGAAATCCATAAATTCAAGCGCAGTGCCCAGCCATCCCGATACTGCTGCTTTCACCAAATCCTTAGTGCTTCTTTCGTGTTTTTCTATTGTCATAATGGTTATCTCTAAAAAAGAGGTAAGATGCGTACTACTTACTC  >  CP009273/1646845‑1647139
                                                                                                                                                    |                                                                                                                                                  
ggTATTTCAGGTATACGGAAAATGGATGTCTGAAAACAATAATGCACAGGTAGCTTTGTTAAATACACAGTTAAGCGAGTTTGCCCCAACCATGCCCCATAACGAAGCAATGAAAAATTAATTTAATATTTATCAAATAGTTAACACGCa                                                                                                                                                   <  1:2743990/150‑1 (MQ=255)
 gTATTTCAGGTATACGGAAAATGGATGTCTGAAAACAATACTGCACAGGTAGCTTTGTTAAATACACAGTTAAGCGAGTTTGCCCCAACCATGCCCCATAACGAAGCAATGAAAAATTAATTTAATATTTATCAAATAGTTAACACGCAt                                                                                                                                                  <  2:1578839/150‑1 (MQ=255)
  tATTTCAGGTATACGGAAAATGGATGTCTGAAAACAATAATGCACAGGTAGCTTTGTTAAATACACAGTTAAGCGAGTTTGCCCCAACCATGCCCCATAACGAAGCAATGAAAAATTAATTTAATATTTATCAAATAGTTAACACGTATg                                                                                                                                                 >  1:3329499/1‑150 (MQ=255)
  tATTTCAGGTATACGGAAAATGGATGTCTGAAAACAATAATGCACAGGTAGCTTTGTTAAATACACAGTTAAGCGAGTTTGCCCCAACCATGCCCCATAACGAAGCAATGAAAAATTAATTTAATATTTATCAAATAGTTAACACGTATg                                                                                                                                                 <  1:457119/150‑1 (MQ=255)
  tATTTCAGGTATACGGAAAATGGATGTCTGAAAACAATAATGCACAGGTAGCTTTGTTAAATACACAGTTAAGCGAGTTTGCCCCAACCATGCCCCATAACGAAGCAATGAAAAATTAATTTAATATTTATCAAATAGTTAACACGCATg                                                                                                                                                 <  1:2470557/150‑1 (MQ=255)
   aTTTCAGGTATACGGAAAATGGATGTCTGAAAACAATAATGCACAGGTAGCTTTGTTAAATACACAGTTAAGCGAGTTTGCCCCAACCATGCCCCATAACGAAGCAATGAAAAATTAATTTAATATTTATCAAATAGTTAACACGCATGa                                                                                                                                                <  2:3344287/150‑1 (MQ=255)
           tataCGGAAAATGGATGTCTGAAAACAATAATGCACAGGTAGCTTTGTTAAATACACAGTTAAGCGAGTTTGCCCCAACCATGCCCCATAACGAAGCAATGAAAAATTAATTTAATATTTATCAAATAGTTAACACGCATGACTCTTGaa                                                                                                                                        >  2:3231474/1‑150 (MQ=255)
             taCGGAAAATGGATGTCTGAAAACAATAATGCACAGGTAGCTTTGTTAAATACACAGTTAAGCGAGTTTGCCCCAACCATGCCCCATAACGAAGCAATGAAAAATTAATTTAATATTTATCAAATAGTTAACACGTATGACTCTTGAAAt                                                                                                                                      <  1:1226175/150‑1 (MQ=255)
                ggAAAATGGATGTCTGAAAACAATAATGCACAGGTAGCTTTGTTAAATACACAGTTAAGCGAGTTTGCCCCAACCATGCCCCATAACGAAGCAATGAAAAATTAATTTAATATTTATCAAATAGTTAACACGCATGACTCTTGAAATCCa                                                                                                                                   >  2:943479/1‑150 (MQ=255)
                 gAAAATGGATGTCTGAAAACAATAATGCACAGGTAGCTTTGTTAAATACACAGTTAAGCGAGTTTGCCCCAACCATGCCCCATAACGAAGCAATGAAAAATTAATTTAATATTTATCAAATAGTTAACACGTATGACTCTTGAAATCCAt                                                                                                                                  >  1:28574/1‑150 (MQ=255)
                  aaaaTGGATGTCTGAAAACAATAATGCACAGGTAGCTTTGTTAAATACACAGTTAAGCGAGTTTGCCCCAACCATGCCCCATAACGAAGCAATGAAAAATTAATTTAATATTTATCAAATAGTTAACACGTATGACTCTTGAAATCCATa                                                                                                                                 <  1:1677193/150‑1 (MQ=255)
                     aTGGATGTCTGAAAACAATAATGCACAGGTAGCTTTGTTAAATACACAGTTAAGCGAGTTTGCCCCAACCATGCCCCATAACGAAGCAATGAAAAATTAATTTAATATTTATCAAATAGTTAACACGTATGACTCTTGAAATCCATAAAt                                                                                                                              >  2:3122603/1‑150 (MQ=255)
                     aTGGATGTCTGAAAACAATAATGCACAGGTAGCTTTGTTAAATACACAGTTAAGCGAGTTTGCCCCAACCATGCCCCATAACGAAGCAATGAAAAATTAATTTAATATTTATCAAATAGTTAACACGCATGACTCTTGAAATCCATAAAt                                                                                                                              >  1:1760995/1‑150 (MQ=255)
                      tGGATGTCTGAAAACAATAATGCACAGGTAGCTTTGTTAAATACACAGTTAAGCGAGTTTGCCCCAACCATGCCCCATAACGAAGCAATGAAAAATTAATTTAATATTTATCAAATAGTTAACACGCATGACTCTTGAAATCCATAAAtt                                                                                                                             <  2:1263972/150‑1 (MQ=255)
                      tGGATGTCTGAAAACAATAATGCACAGGTAGCTTTGTTAAATACACAGTTAAGCGAGTTTGCCCCAACCATGCCCCATAACGAAGCAATGAAAAATTAATTTAATATTTATCAAATAGTTAACACGCATGACTCTTGAAATCCATAAAtt                                                                                                                             <  1:806936/150‑1 (MQ=255)
                       ggATGTCTGAAAACAATAATGCACAGGTAGCTTTGTTAAATACACAGTTAAGCGAGTTTGCCCCAACCATGCCCCATAACGAAGCAATGAAAAATTAATTTAATATTTATCAAATATTTAACACGTATGACTCTTGAAATCCCAAAATTc                                                                                                                            >  1:214125/1‑150 (MQ=255)
                           gTCTGAAAACAATAATGCACAGGTAGCTTTGTTAAATACACAGTTAAGCGAGTTTGCCCCAACCATGCCCCATAACGAAGCAATGAAAAATTAATTTAATATTTATCAAATAGTTAACACGCATGACTCTTGAAATCCATAAATTCAAgc                                                                                                                        <  1:3705447/150‑1 (MQ=255)
                            tCTGAAAACAATAATGCACAGGTAGCTTTGTTAAATACACAGTTAAGCGAGTTTGCCCCAACCATGCCCCATAACGAAGCAATGAAAAATTAATTTAATATTTATCAAATAGTTAACACGTATGACTCTTGAAATCCATAAATTCAAgcg                                                                                                                       <  1:33600/150‑1 (MQ=255)
                                 aaaCAATAATGCACAGGTAGCTTTGTTAAATACACAGTTAAGCGAGTTTGCCCCAACCATGCCCCATAACGAAGCAATGAAAAATTAATTTAATATTTATCAAATAGTTAACACGCATGACTCTTGAAATCCATAAATTCAAGCGCAGTg                                                                                                                  >  1:1069954/1‑150 (MQ=255)
                                     aataatGCACAGGTAGCTTTGTTAAATACACAGTTAAGCGAGTTTGCCCCAACCATGCCCCATAACGAAGCAATGAAAAATTAATTTAATATTTATCAAATAGTTAACACGTATGACTCTTGAAATCCATAAATTCAAGCGCAGTGCCCa                                                                                                              >  2:568717/1‑150 (MQ=255)
                                     aataatGCACAGGTAGCTTTGTTAAATACACAGTTAAGCGAGTTTGCCCCAACCATGCCCCATAACGAAGCAATGAAAAATTAATTTAATATTTATCAAATAGTTAACACGCATGACTCTTGAAATCCATAAATTCAAGCGCAGTGCCCa                                                                                                              <  1:2935794/150‑1 (MQ=255)
                                      ataatGCACAGGTAGCTTTGTTAAATACACAGTTAAGCGAGTTTGCCCCAACCATGCCCCATAACGAAGCAATGAAAAATTAATTTAATATTTATCAAATAGTTAACACGCATGACTCTTGAAATCCATAAATTCAAGCGCAGTGCCCAg                                                                                                             >  2:1357788/1‑150 (MQ=255)
                                        aatGCACAGGTAGCTTTGTTAAATACACAGTTAAGCGAGTTTGCCCCAACCATGCCCCATAACGAAGCAATGAAAAATTAATTTAATATTTATCAAATAGTTAACACGTATGACTCTTGAAATCCATAAATTCAAGCGCAGTGCCCAGcc                                                                                                           <  2:3512453/150‑1 (MQ=255)
                                          tGCACAGGTAGCTTTGTTAAATACACAGTTAAGCGAGTTTGCCCCAACCATGCCCCATAACGAAGCAATGAAAAATTAATTTAATATTTATCAAATAGTTAACACGCATGACTCTTGAAATCCATAAATTCAAGCGCAGTGCCCAGCCAt                                                                                                         <  2:1948570/150‑1 (MQ=255)
                                           gCACAGGTAGCTTTGTTAAATACACAGTTAAGCGAGTTTGCCCCAACCATGCCCCATAACGAAGCAATGAAAAATTAATTTAATATTTATCAAATAGTTAACACGTATGACTCTTGAAATCCATAAATTCAAGCGCAGTGCCCAGCCATc                                                                                                        <  2:1116029/150‑1 (MQ=255)
                                           gCACAGGTAGCTTTGTTAAATACACAGTTAAGCGAGTTTGCCCCAACCATGCCCCATAACGAAGCAATGAAAAATTAATTTAATATTTATCAAATAGTTAACACGCATGACTCTTGAAATCCATAAATTCAAGCGCAGTGCCCAGCCATc                                                                                                        <  1:903160/150‑1 (MQ=255)
                                              caGGTAGCTTTGTTAAATACACAGTTAAGCGAGTTTGCCCCAACCATGCCCCATAACGAAGCAATGAAAAATTAATTTAATATTTATCAAATAGTTAACACGCATGACTCTTGAAATCCATAAATTCAAGCGCAGTGCCCAGCCATCCCg                                                                                                     <  1:1973411/150‑1 (MQ=255)
                                              caGGTAGCTTTGTTAAATACACAGTTAAGCGAGTTTGCCCCAACCATGCCCCATAACGAAGCAATGAAAAATTAATTTAATATTTATCAAATAGTTAACACGCATGACTCTTGAAATCCATAAATTCAAGCGCAGTGCCCAGCCATCCCg                                                                                                     <  1:3181931/150‑1 (MQ=255)
                                                 gTAGCTTTGTTAAATACACAGTTAAGCGAGTTTGCCCCAACCATGCCCCATAACGAAGCAATGAAAAATTAATTTAATATTTATCAAATAGTTAACACGTATGACTCTTGAAATCCATAAATTCAAGCGCAGTGCCCAGCCATCCCGATa                                                                                                  >  2:2949045/1‑150 (MQ=255)
                                                 gTAGCTTTGTTAAATACACAGTTAAGCGAGTTTGCCCCAACCATGCCCCATAACGAAGCAATGAAAAATTAATTTAATATTTATCAAATAGTTAACACGCATGACTCTTGAAATCCATAAATTCAAGCGCAGTGCCCAGCCATCCCGATa                                                                                                  >  1:1501542/1‑150 (MQ=255)
                                                  tAGCTTTGTTAAATACACAGTTAAGCGAGTTTGCCCCAACCATGCCCCATAACGAAGCAATGAAAAATTAATTTAATATTTATCAAATAGTTAACACGTATGACTCTTGAAATCCATAAATTCAAGCGCAGTGCCCAGCCATCCCGATAc                                                                                                 >  1:663066/1‑150 (MQ=255)
                                                   aGCTTTGTTAAATACACAGTTAAGCGAGTTTGCCCCAACCATGCCCCATAACGAAGCAATGAAAAATTAATTTAATATTTATCAAATAGTTAACACGCATGACTCTTGAAATCCATAAATTCAAGCGCAGTGCCCAGCCATCCCGATAct                                                                                                <  1:218392/150‑1 (MQ=255)
                                                   aGCTTTGTTAAATACACAGTTAAGCGAGTTTGCCCCAACCATGCCCCATAACGAAGCAATGAAAAATTAATTTAATATTTATCAAATAGTTAACACGCATGACTCTTGAAATCCATAAATTCAAGCGCAGTGCCCAGCCATCCCGATAct                                                                                                >  1:3873378/1‑150 (MQ=255)
                                                       ttGTTAAATACACAGTTAAGCGAGTTTGCCCCAACCATGCCCCATAACGAAGCAATGAAAAATTAATTTAATATTTATCAAATAGTTAACACGCATGACTCTTGAAATCCATAAATTCAAGCGCAGTGCCCAGCCATCCCGATActgctg                                                                                            <  2:278301/150‑1 (MQ=255)
                                                         gTTAAATACAAAGTTAAGCGAGTTTGCCCCAACCATGCCCCATAACGAATCAATGAAAAATTAATTTAATAAATATCAAATAGTTAACACGCATGACTATTTAAATCCAAAAATTCAAGCCCAGAGCCCAGCCATCCCGATACTgctgct                                                                                          >  2:1834168/1‑150 (MQ=255)
                                                          ttAAATACACAGTTAAGCGAGTTTGCCCCAACCATGCCCCATAACGAAGCAATGAAAAATTAATTTAATATTTATCAAATAGTTAACACGCATGACTCTTGAAATCCATAAATTCAAGCGCAGTGCCCAGCCATCCCGATACTGCTGCtt                                                                                         >  2:2622261/1‑150 (MQ=255)
                                                           tAAATACACAGTTAAGCGAGTTTGCCCCAACCATGCCCCATAACGAAGCAATGAAAAATTAATTTAATATTTATCAAATAGTTAACACGTATGACTCTTGAAATCCATAAATACAAGCGCAGTGCCCCGCCATCCCGATACTGCTGCttt                                                                                        >  2:342609/1‑150 (MQ=255)
                                                           tAAATACACAGTTAAGCGAGTTTGCCCCAACCATGCCCCATAACGAAGCAATGAAAAATTAATTTAATATTTATCAAATAGTTAACACGCATGACTCTTGAAATCCATAAATTCAAGCGCAGTGCCCAGCCATCCCGATACTGCTGCttt                                                                                        <  2:1641680/150‑1 (MQ=255)
                                                            aaaTACACAGTTAAGCGAGTTTGCCCCAACCATGCCCCATAACGAAGCAATGAAAAATTAATTTAATATTTATCAAATAGTTAACACGTATGACTCTTGAAATCCATAAATTCAAGCGCAGTGCCCAGCCATCCCGATACTGCTGCTTTc                                                                                       >  1:413118/1‑150 (MQ=255)
                                                             aaTACACAGTTAAGCGAGTTTGCCCCAACCATGCCCCATAACGAAGCAATGAAAAATTAATTTAATATTTATCAAATAGTTAACACGTATGACTCTTGAAATCCATAAATTCAAGCGCAGTGCCCAGCCATCCCGATACTGCTGCTTTCa                                                                                      >  2:3384456/1‑150 (MQ=255)
                                                             aaTACACAGTTAAGCGAGTTTGCCCCAACCATGCCCCATAACGAAGCAATGAAAAATTAATTTAATATTTATCAAATAGTTAACACGCATGACTCTTGAAATCCATAAATTCAAGCGCAGTGCCCAGCCATCCCGATACTGCTGCTTTCa                                                                                      <  1:3247457/150‑1 (MQ=255)
                                                                     gTTAAGCGAGTTTGCCCCAACCATGCCCCATAACGAAGCAATGAAAAATTAATTTAATATTTATCAAATAGTTAACACGTATGACTCTTGAAATCCATAAATTCAAGCGCAGTGCCCAGCCATCCCGATACTGCTGCTTTCACCAAATcc                                                                              <  1:2798960/150‑1 (MQ=255)
                                                                     gTTAAGCGAGTTTGCCCCAACCATGCCCCATAACGAAGCAATGAAAAATTAATTTAATATTTATCAAATAGTTAACACGCATGACTCTTGAAATCCATAAATTCAAGCGCAGTGCCCAGCCATCCCGATACTGCTGCTTTCACCAAATcc                                                                              <  1:1806259/150‑1 (MQ=255)
                                                                      ttAAGCGAGTTTGCCCCAACCATGCCCCATAACGAAGCAATGAAAAATTAATTTAATATTTATCAAATAGTTAACACGTATGACTCTTGAAATCCATAAATTCAAGCGCAGTGCCCAGCCATCCCGATACTGCTGCTTTCACCAAATCCt                                                                             <  2:2939340/150‑1 (MQ=255)
                                                                        aaGCGAGTTTGCCCCAACCATGCCCCATAACGAAGCAATGAAAAATTAATTTAATATTTATCAAATAGTTAACACGTATGACTCTTGAAATCCATAAATTCAAGCGCAGTGCCCAGCCATCCCGATACTGCTGCTTTCACCAAATCCTTa                                                                           >  2:1784215/1‑150 (MQ=255)
                                                                        aaGCGAGTTTGCCCCAACCATGCCCCATAACGAAGCAATGAAAAATTAATTTAATATTTATCAAATAGTTAACACGTATGACTCTTGAAATCCATAAATTCAAGCGCAGTGCCCAGCCATCCCGATACTGCTGCTTTCACCAAATCCTTa                                                                           >  2:3189196/1‑150 (MQ=255)
                                                                        aaGCGAGTTTGCCCCAACCATGCCCCATAACGAAGCAATGAAAAATTAATTTAATATTTATCAAATAGTTAACACGTATGACTCTTGAAATCCATAAATTCAAGCGCAGTGCCCAGCCATCCCGATACTGCTGCTTTCACCAAATCCTTa                                                                           >  2:974850/1‑150 (MQ=255)
                                                                               tttGCCCCAACCATGCCCCATAACGAAGCAATGAAAAATTAATTTAATATTTATCAAATAGTTAACACGTATGACTCTTGAAATCCATAAATTCAAGCGCAGTGCCCAGCCATCCCGATACTGCTGCTTTCACCAAATCCTTAGTGcttc                                                                    >  1:600530/1‑150 (MQ=255)
                                                                                ttGCCCCAACCATGCCCCATAACGAAGCAATGAAAAATTAATTTAATATTTATCAAATAGTTAACACGTATGACTCTTGAAATCCATAAATTCAAGCGCAGTGCCCAGCCATCCCGATACTGCTGCTTTCACCAAATCCTTAGTGcttct                                                                   <  1:1035002/150‑1 (MQ=255)
                                                                                  gCCCCAACCATGCCCCATAACGAAGCAATGAAAAATTAATTTAATATTTATCAAATAGTTAACACGTATGACTCTTGAAATCCATAAATTCAAGCGCAGTGCCCAGCCATCCCGATACTGCTGCTTTCACCAAATCCTTAGTGCTTCttt                                                                 >  1:2631878/1‑150 (MQ=255)
                                                                                    cccAACCATGCCCCATAACGAAGCAATGAAAAATTAATTTAATATTTATCAAATAGTTAACACGCATGACTCTTGAAATCCATAAATTCAAGCGCAGTGCCCAGCCATCCCGATACTGCTGCTTTCACCAAATCCTTAGTGCTTCTTTCg                                                               <  2:3802595/150‑1 (MQ=255)
                                                                                      cAACCATGCCCCATAACGAAGCAATGAAAAATTAATTTAATATTTATCAAATAGTTAACACGCATGACTCTTGAAATCCATAAATTCAAGCGCAGTGCCCAGCCATCCCGATACTGCTGCTTTCACCAAATCCTTAGTGCTTCTTTCgtg                                                             >  2:3152880/1‑150 (MQ=255)
                                                                                        aCCATGCCCCATAACGAAGCAATGAAAAATTAATTTAATATTTATCAAATAGTTAACACGTATGACTCTTGAAATCCATAAATTCAAGCGCAGTGCCCAGCCATCCCGATACTGCTGCTTTCACCAAATCCTTAGTGCTTCTTTCGTGtt                                                           >  1:3598433/1‑150 (MQ=255)
                                                                                        aCCATGCCCCATAACGAAGCAATGAAAAATTAATTTAATATTTATCAAATAGTTAACACGCATGACTCTTGAAATCCATAAATTCAAGCGCAGTGCCCAGCCATCCCGATACTGCTGCTTTCACCAAATCCTTAGTGCTTCTTTCGTGtt                                                           >  2:1735729/1‑150 (MQ=255)
                                                                                               cccATAACGAAGCAATGAAAAATTAATTTAATATTTATCAAATAGTTAACACGCATGACTCTTGAAATCCATAAATTCAAGCGCAGTGCCCAGCCATCCCGATACTGCTGCTTTCACCAAATCCTTAGTGCTTCTTTCGTGTTTTTCTAt                                                    <  2:2162111/150‑1 (MQ=255)
                                                                                                ccATAACGAAGCAATGAAAAATTAATTTAATATTTATCAAATAGTTAACACGTATGACTCTTGAAATCCATAAATTCAAGCGCAGTGCCCAGCCATCCCGATACTGCTGCTTTCACCAAATCCTTAGTGCTTCTTTCGTGTTTTTCTAtt                                                   >  1:2536725/1‑150 (MQ=255)
                                                                                                 cATAACGAAGCAATGAAAAATTAATTTAATATTTATCAAATAGTTAACACGTATGACTCTTGAAATCCATAAATTCAAGCGCAGTGCCCAGCCATCCCGATACTGCTGCTTTCACCAAATCCTTAGTGCTTCTTTCGTGTTTTTCTATTg                                                  >  1:1618482/1‑150 (MQ=255)
                                                                                                  aTAACGAAGCAATGAAAAATTAATTTAATATTTATCAAATAGTTAACACGTATGACTCTTGAAATCCATAAATTCAAGCGCAGTGCCCAGCCATCCCGATACTGCTGCTTTCACCAAATCCTTAGTGCTTCTTTCGTGTTTTTCTATTGt                                                 <  2:1801229/150‑1 (MQ=255)
                                                                                                   tAACGAAGCAATGAAAAATTAATTTAATATTTATCAAATAGTTAACACGCATGACTCTTGAAATCCATAAATTCAAGCGCAGTGCCCAGCCATCCCGATACTGCTGCTTTCACCAAATCCTTAGTGCTTCTTTCGTGTTTTTCTATTGTc                                                >  1:2471103/1‑150 (MQ=255)
                                                                                                   tAACGAAGCAATGAAAAATTAATTTAATATTTATCAAATAGTTAACACGCATGACTCTTGAAATCCATAAATTCAAGCGCAGTGCCCAGCCATCCCGATACTGCTGCTTTCACCAAATCCTTAGTGCTTCTTTCGTGTTTTTCTATTGTc                                                <  1:1932459/150‑1 (MQ=255)
                                                                                                    aaCGAAGCAATGAAAAATTAATTTAATATTTATCAAATAGTTAACACGTATGACTCTTGAAATCCATAAATTCAAGCGCAGTGCCCAGCCATCCCGATACTGCTGCTTTCACCAAATCCTTAGTGCTTCTTTCGTGTTTTTCTATTGTCa                                               >  1:2166392/1‑150 (MQ=255)
                                                                                                     aCGAAGCAATGAAAAATTAATTTAATATTTATCAAATAGTTAACACGCATGACTCTTGAAATCCATAAATTCAAGCGCAGTGCCCAGCCATCCCGATACTGCTGCTTTCACCAAATCCTTAGTGCTTCTTTCGTGTTTTTCTATTGTCAt                                              <  1:3401843/150‑1 (MQ=255)
                                                                                                      cGAAGCAATGAAAAATTAATTTAATATTTATCAAATAGTTAACACGCATGACTCTTGAACTCCATAAATTCAAGCGCAGTGCCCAGCCATCCCGATACTGCTGCTTTCACCAAATCCTTAGTGCTTCTTTCGTGTTTTTCTATTGTCATa                                             >  2:1672476/1‑150 (MQ=255)
                                                                                                      cGAAGCAATGAAAAATTAATTTAATATTTATCAAATAGTTAACACGCATGACTCTTGAAATCCATAAATTCAAGCGCAGTGCCCAGCCATCCCGATACTGCTGCTTTCACCAAATCCTTAGTGCTTCTTTCGTGTTTTTCTATTGTCATa                                             >  2:1672544/1‑150 (MQ=255)
                                                                                                      cGAAGCAATGAAAAATTAATTTAATATTTATCAAATAGTTAACACGCATGACTCTTGAAATCCATAAATTCAAGCGCAGTGCCCAGCCATCCCGATACTGCTGCTTTCACCAAATCCTTAGTGCTTCTTTCGTGTTTTTCTATTGTCATa                                             <  2:3645185/150‑1 (MQ=255)
                                                                                                       gAAGCAATGAAAAATTAATTTAATATTTATCAAATAGTTAACACGTATGACTCTTGAAATCCATAAATTCAAGCGCAGTGCCCAGCCATCCCGATACTGCTGCTTTCACCAAATCCTTAGTGCTTCTTTCGTGTTTTTCTATTGTCATaa                                            >  1:2023682/1‑150 (MQ=255)
                                                                                                       gAAGCAATGAAAAATTAATTTAATATTTATCAAATAGTTAACACGTATGACTCTTGAAATCCATAAATTCAAGCGCAGTGCCCAGCCATCCCGATACTGCTGCTTTCACCAAATCCTTAGTGCTTCTTTCGTGTTTTTCTATTGTCATaa                                            <  2:1395694/150‑1 (MQ=255)
                                                                                                       gAAGCAATGAAAAATTAATTTAATATTTATCAAATAGTTAACACGTATGACTCTTGAAATCCATAAATTCAAGCGCAGAGCCCAGCCATCCCGATACTGCTGCTTTCACCAAATCCTTAGTGCTTCTTTCGTGTTTTTCTATTGTCATaa                                            >  1:2828298/1‑150 (MQ=255)
                                                                                                        aaGCAATGAAAAATTAATTTAATATTTATCAAATAGTTAACACGTATGACTCTTGAAATCCATAAATTCAAGCGCAGTGCCCAGCCATCCCGATACTGCTGCTTTCACCAAATCCTTAGTGCTTCTTTCGTGTTTTTCTATTGTCATAAt                                           >  2:1634981/1‑150 (MQ=255)
                                                                                                         aGCAATGAAAAATTAATTTAATATTTATCAAATAGTTAACACGCATGACTCTTGAAATCCATAAATTCAAGCGCAGTGCCCAGCCATCCCGATACTGCTGCTTTCACCAAATCCTTAGTGCTTCTTTCGTGTTTTTCTATTGTCATAATg                                          <  2:654823/150‑1 (MQ=255)
                                                                                                           cAATGAAAAATTAATTTAATATTTATCAAATAGTTAACACGTATGACTCTTGAAATCCATAAATTCAAGCGCAGTGCCCAGCCATCCCGATACTGCTGCTTTCACCAAATCCTTAGTGCTTCTTTCGTGTTTTTCTATTGTCATAATGGt                                        <  2:680072/150‑1 (MQ=255)
                                                                                                            aaTGAAAAATTAATTTAATATTTATCAAATAGTTAACACGCATGACTCTTGAAATCCATAAATTCAAGCGCAGTGCCCAGCCATCCCGATACTGCTGCTTTCACCAAATCCTTAGTGCTTCTTTCGTGTTTTTCTATTGTCATAATGGtt                                       >  2:1114235/1‑150 (MQ=255)
                                                                                                              tGAAAAATTAATTTAATATTTATCAAATAGTTAACACGTATGACTCTTGAAATCCATAAATTCAAGCGCAGTGCCCAGCCATCCCGATACTGCTGCTTTCACCAAATCCTTAGTGCTTCTTTCGTGTTTTTCTATTGTCATAATGGTTAt                                     >  1:3641870/1‑150 (MQ=255)
                                                                                                              tGAAAAATTAATTTAATATTTATCAAATAGTTAACACGCATGACTCTTGAAATCCATAAATTCAAGCGCAGTGCCCAGCCATCCCGATACTGCTGCTTTCACCAAATCCTTAGTGCTTCTTTCGTGTTTTTCTATTGTCATAATGGTTAt                                     >  2:460425/1‑150 (MQ=255)
                                                                                                              tGAAAAATTAATTTAATATTTATCAAATAGTTAACACGCATGACTCTTGAAATCCATAAATTCAAGCGCAGTGCCCAGCCATCCCGATACTGCTGCTTTCACCAAATCCTTAGTGCTTCTTTCGTGTTTTTCTATTGTCATAATGGTTAt                                     >  1:3842176/1‑150 (MQ=255)
                                                                                                               gAAAAATTAATTTAATATTTATCAAATAGTTAACACGTATGACTCTTGAAATCCATAAATTCAAGCGCAGTGCCCAGCCATCCCGATACTGCTGCTTTCACCAAATCCTTAGTGCTTCTTTCGTGTTTTTCTATTGTCATAATGGTTAtc                                    >  1:960886/1‑150 (MQ=255)
                                                                                                               gAAAAATTAATTTAATATTTATCAAATAGTTAACACGCATGACTCTTGAAATCCATAAATTCAAGCGCAGTGCCCAGCCATCCCGATACTGCTGCTTTCACCAAATCCTTAGTGCTTCTTTCGTGTTTTTCTATTGTCATAATGGTTAtc                                    >  2:1434206/1‑150 (MQ=255)
                                                                                                                aaaaaTTAATTTAATATTTATCAAATAGTTAACACGTATGACTCTTGAAATCCATAAATTCAAGCGCAGTGCCCAGCCATCCCGATACTGCTGCTTTCACCAAATCCTTAGTGCTTCTTTCGTGTTTTTCTATTGTCATAATGGTTAtct                                   <  1:69038/150‑1 (MQ=255)
                                                                                                                aaaaaTTAATTTAATATTTATCAAATAGTTAACACGCATGACTCTTGAAATCCATAAATTCAAGCGCAGTGCCCAGCCATCCCGATACTGCTGCTTTCACCAAATCCTTAGTGCTTCTTTCGTGTTTTTCTATTGTCATAATGGTTAtct                                   >  2:1810978/1‑150 (MQ=255)
                                                                                                                 aaaaTTAATTTAATATTTATCAAATAGTTAACACGTATGACTCTTGAAATCCATAAATTCAAGCGCAGTGCCCAGCCATCCCGATACTGCTGCTTTCACCAAATCCTTAGTGCTTCTTTCGTGTTTTTCTATTGTCATAATGGTTAtctc                                  >  1:3658180/1‑150 (MQ=255)
                                                                                                                  aaaTTAATTTAATATTTATCAAATAGTTAACACGCATGACTCTTGAAATCCATAAATTCAAGCGCAGTGCCCAGCCATCCCGATACTGCTGCTTTCACCAAATCCTTAGTGCTTCTTTCGTGTTTTTCTATTGTCATAATGGTTATctct                                 >  2:874436/1‑150 (MQ=255)
                                                                                                                   aattaattTAATATTTATCAAATAGTTAACACGTATGACTCTTGAAATCCATAAATTCAAGCGCAGTGCCCAGCCATCCCGATACTGCTGCTTTCACCAAATCCTTAGTGCTTCTTTCGTGTTTTTCTATTGTCATAATGGTTATCTCTa                                >  1:633446/1‑150 (MQ=255)
                                                                                                                   aattaattTAATATTTATCAAATAGTTAACACGCATGACTCTTGAAATCCATAAATTCAAGCGCAGTGCCCAGCCATCCCGATACTGCTGCTTTCACCAAATCCTTAGTGCTTCTTTCGTGTTTTTCTATTGTCATAATGGTTATCTCTa                                <  1:3687137/150‑1 (MQ=255)
                                                                                                                   aattaattTAATATTTATCAAATAGTTAACACGCATGACTCTTGAAATCCATAAATTCAAGCGCAGTGCCCAGCCATCCCGATACTGCTGCTTTCACCAAATCCTTAGTGCTTCTTTCGTGTTTTTCTATTGTCATAATGGTTATCTCTa                                <  2:4079616/150‑1 (MQ=255)
                                                                                                                    attaattTAATATTTATCAAATAGTTAACACGTATGACTCTTGAAATCCATAAATTCAAGCGCAGTGCCCAGCCATCCCGATACTGCTGCTTTCACCAAATCCTTAGTGCTTCTTTCGTGTTTTTCTATTGTCATAATGGTTATCTCTaa                               <  2:1156211/150‑1 (MQ=255)
                                                                                                                       aatttaatATTTATCAAATAGTTAACACGCATGACTCTTGAAATCCATAAATTCAAGCGCAGTGCCCAGCCATCCCGATACTGCTGCTTTCACCAAATCCTTAGTGCTTCTTTCGTGTTTTTCTATTGTCATAATGGTTATCTCTaaaaa                            >  2:3250214/1‑150 (MQ=255)
                                                                                                                           taatATTTATCAAATAGTTAACACGTATGACTCTTGAAATCCATAAATTCAAGCGCAGTGCCCAGCCATCCCGATACTGCTGCTTTCACCAAATCCTTAGTGCTTCTTTCGTGTTTTTCTATTGTCATAATGGTTATCTCTAAAAAagag                        <  1:2787089/150‑1 (MQ=255)
                                                                                                                           taatATTTATCAAATAGTTAACACGTATGACTCTTGAAATCCATAAATTCAAGCGCAGTGCCCAGCCATCCCGATACTGCTGCTTTCACCAAATCCTTAGTGCTTCTTTCGTGTTTTTCTATTGTCATAATGGTTATCTCTAAAAAagag                        >  2:3826282/1‑150 (MQ=255)
                                                                                                                                   aTCAAATAGTTAACACGTATGACTCTTGAAATCCATAAATTCAAGCGCAGTGCCCAGCCATCCCGATACTGCTGCTTTCACCAAATCCTTAGTGCTTCTTTCGTGTTTTTCTATTGTCATAATGGTTATCTCTAAAAAAGAGGTAAGATg                <  1:2895127/150‑1 (MQ=255)
                                                                                                                                   aTCAAATAGTTAACACGCATGACTCTTGAAATCCATAAATTCAAGCGCAGTGCCCAGCCATCCCGATACTGCTGCTTTCACCAAATCCTTAGTGCTTCTTTCGTGTTTTTCTATTGTCATAATGGTTATCTCTAAAAAAGAGGTAAGATg                >  2:671981/1‑150 (MQ=255)
                                                                                                                                    tCAAATAGTTAACACGTATGACTCTTGAAATCCATAAATTCAAGCGCAGTGCCCAGCCATCCCGATACTGCTGCTTTCACCAAATCCTTAGTGCTTCTTTCGTGTTTTTCTATTGTCCTAATGGTTATCTCTAAAAAAGAGGTAAGATGc               >  1:3483374/1‑150 (MQ=255)
                                                                                                                                    tCAAATAGTTAACACGCATGACTCTTGAAATCCATAAATTCAAGCGCAGTGCCCAGCCATCCCGATACTGCTGCTTTCACCAAATCCTTAGTGCTTCTTTCGTGTTTTTCTATTGTCATAATGGTTATCTCTAAAAAAGAGGTAAGATGc               <  2:435556/150‑1 (MQ=255)
                                                                                                                                      aaaTAGTTAACACGTATGACTCTTGAAATCCATAAATTCAAGCGCAGTGCCCAGCCATCCCGATACTGCTGCTTTCACCAAATCCTTAGTGCTTCTTTCGTGTTTTTCTATTGTCATAATGGTTATCTCTAAAAAAGAGGTAAGATGCGt             <  1:2237752/150‑1 (MQ=255)
                                                                                                                                      aaaTAGTTAACACGCATGACTCTTGAAATCCATAAATTCAAGCGCAGTGCCCAGCCATCCCGATACTGCTGCTTTCACCAAATCCTTAGTGCTTCTTTCGTGTTTTTCTATTGTCATAATGGTTATCTCTAAAAAAGAGGTAAGATGCGt             >  1:1868143/1‑150 (MQ=255)
                                                                                                                                       aaTAGTTAACACGTATGACTCTTGAAATCCATAAATTCCAGCGCAGTGCCCAGCCATCCCGATACTGCTGCTTTCACCAAATCCTTAGTGCTTCTTTCGTGTTTTTCTATTGTCATAATGGTTATCTCTAAAAAAGAGGTAAGATGCGta            >  2:2864770/1‑150 (MQ=255)
                                                                                                                                       aaTAGTTAACACGCATGACTCTTGAAATCCATAAATTCAAGCGCAGTGCCCAGCCATCCCGATACTGCTGCTTTCACCAAATCCTTAGTGCTTCTTTCGTGTTTTTCTATTGTCATAATGGTTATCTCTAAAAAAGAGGTAAGATGCGta            >  1:3423881/1‑150 (MQ=255)
                                                                                                                                        aTAGTTAACACGCATGACTCTTGAAATCCATAAATTCAAGCGCAGTGCCCAGCCATCCCGATACTGCTGCTTTCACCAAATCCTTAGTGCTTCTTTCGTGTTTTTCTATTGTCATAATGGTTATCTCTAAAAAAGAGGTAAGATGCGtac           <  1:1304145/150‑1 (MQ=255)
                                                                                                                                         tAGTTAACACGCATGACTCTTGAAATCCATAAATTCAAGCGCAGTGCCCAGCCATCCCGATACTGCTGCTTTCACCAAATCCTTAGTGCTTCTTTCGTGTTTTTCTATTGTCATAATGGTTATCTCTAAAAAAGAGGTAAGATGCGtact          <  1:3572965/150‑1 (MQ=255)
                                                                                                                                           gTTAACACGTATGACTCTTGAAATCCATAAATTCAAGCGCAGTGCCCAGCCATCCCGATACTGCTGCTTTCACCAAATCCTTAGTGCTTCTTTCGTGTTTTTCTATTGTCATAATGGTTATCTCTAAAAAAGAGGTAAGATGCGtactac        >  1:3787190/1‑150 (MQ=255)
                                                                                                                                             tAACACGTATGACTCTTGAAATCCATAAATTCAAGCGCAGTGCCCAGCCATCCCGATACTGCTGCTTTCACCAAATCCTTAGTGCTTCTTTCGTGTTTTTCTATTGTCATAATGGTTATCTCTAAAAAAGAGGTAAGATGCGTACtactt      >  1:476864/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
GGTATTTCAGGTATACGGAAAATGGATGTCTGAAAACAATAATGCACAGGTAGCTTTGTTAAATACACAGTTAAGCGAGTTTGCCCCAACCATGCCCCATAACGAAGCAATGAAAAATTAATTTAATATTTATCAAATAGTTAACACGCATGACTCTTGAAATCCATAAATTCAAGCGCAGTGCCCAGCCATCCCGATACTGCTGCTTTCACCAAATCCTTAGTGCTTCTTTCGTGTTTTTCTATTGTCATAATGGTTATCTCTAAAAAAGAGGTAAGATGCGTACTACTTACTC  >  CP009273/1646845‑1647139

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: