Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,646,372 C→A 92.9% pseudogene (565/1158 nt) intQ → pseudogene, Qin prophage; phage integrase family;Phage or Prophage Related

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,646,3720CA92.9% 675.6 / 42.0 239pseudogene (565/1158 nt)intQpseudogene, Qin prophage; phage integrase family;Phage or Prophage Related
Reads supporting (aligned to +/- strand):  ref base C (12/5);  new base A (104/118);  total (116/123)
Fisher's exact test for biased strand distribution p-value = 7.78e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

Only 100 of 239 total aligned reads displayed.

TAATCCGTTTAAGGGGATCTCGCCATTAAAGGAATCAAGAACCATTCCGGATCCTCTTTCGCGGGAAGAATTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTCGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCC  >  CP009273/1646227‑1646517
                                                                                                                                                 |                                                                                                                                                 
 aaTCCGTTTAAGGGGATCTCGCCATTAAAGGAATCAAGAACCATTCCGGATCCTCTTTCGCGGGAAGAATTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGcc                                                                                                                                                >  2:1601348/1‑148 (MQ=255)
 aaTCCGTTTAAGGGGATCTCGCCATTAAAGGAATCAAGAACCATTCCGGATCCTCTTTCGCGGGAAGAATTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGcc                                                                                                                                                <  1:219863/148‑1 (MQ=255)
 aaTCCGTTTAAGGGGATCTCGCCATTAAAGGAATCAAGAACCATTCCGGATCCTCTTTCGCGGGAAGAATTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGcc                                                                                                                                                <  2:6057/148‑1 (MQ=255)
   tCCGTTTAAGGGGATCTCGCCATTAAAGGAATCAAGAACCATTCCGGATCCTCTTTCGCGGGAAGAATTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCt                                                                                                                                              <  2:2469789/148‑1 (MQ=255)
      gTTTAAGGGGATCTCGCCATTAAAGGAATCAAGCACCATTCCGGATCCTCTTTCGCGGGAAGAATTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGt                                                                                                                                           <  2:1804701/148‑1 (MQ=255)
      gTTTAAGGGGATCTCGCCATTAAAGGAATCAAGAACCATTCCGGATCCTCTTTCGCGGGAAGAATTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTCGCCCTGGt                                                                                                                                           <  2:531083/148‑1 (MQ=255)
      gTTTAAGGGGATCTCGCCATTAAAGGAATCAAGAACCATTCCGGATCCTCTTTCGCGGGAAGAATTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGt                                                                                                                                           <  1:2500934/148‑1 (MQ=255)
      gTTTAAGGGGATCTCGCCATTAAAGGAATCAAGAACCATTCCGGATCCTCTTTCGCGGGAAGAATTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGt                                                                                                                                           <  2:1804960/148‑1 (MQ=255)
          aaGGGGATCTCGCCATTAAAGGAATCAAGAACCATTCCGGATCCTCTTTCGCGGGAAGAATTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGc                                                                                                                                       <  1:708407/148‑1 (MQ=255)
          aaGGGGATCTCGCCATTAAAGGAATCAAGAACCATTCCGGATCCTCTTTCGCGGGAAGAATTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGc                                                                                                                                       >  2:1939860/1‑148 (MQ=255)
             gggATCTCGCCATTAAAGGAATCAAGAACCATTCCGGATCCTCTTTCGCGGGAAGAATTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTCGCCCTGGTGAGCtgt                                                                                                                                    >  2:2900807/1‑148 (MQ=255)
             gggATCTCGCCATTAAAGGAATCAAGAACCATTCCGGATCCTCTTTCGCGGGAAGAATTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTCGCCCTGGTGAGCtgt                                                                                                                                    >  2:1594963/1‑148 (MQ=255)
                aTCTCGCCATTAAAGGAATCAAGAACCATTCCGGATCCTCTTTCGCGGGAAGAATTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCtgtgtg                                                                                                                                 >  1:1324100/1‑148 (MQ=255)
                 tctcGCCATTAAAGGAATCAAGAACCATTCCGGATCCTCTTTCGCGGGAAGAATTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGc                                                                                                                                >  1:2653028/1‑148 (MQ=255)
                      ccATTAAAGGAATCAAGAACCATTCCGGATCCTCTTTCGCGGGAAGAATTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTg                                                                                                                           <  1:2834447/148‑1 (MQ=255)
                       cATTAAAGGAATCAAGAACCATTCCGGATCCTCTTTCGCGGGAAGAATTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTgg                                                                                                                          >  1:10763/1‑148 (MQ=255)
                       cATTAAAGGAATCAAGAACCATTCCGGATCCTCTTTCGCGGGAAGAATTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTgg                                                                                                                          >  1:2208270/1‑148 (MQ=255)
                       cATTAAAGGAATCAAGAACCATTCCGGATCCTCTTTCGCGGGAAGAATTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTgg                                                                                                                          <  1:653159/148‑1 (MQ=255)
                         ttAAAGGAATCAAGAACCATTCCGGATCCTCTTTCGCGGGAAGAATTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGAt                                                                                                                        <  2:1124565/148‑1 (MQ=255)
                          tAAAGGAATCAAGAACCATTCCGGATCCTCTTTCGCGGGAAGAATTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATg                                                                                                                       >  1:1913176/1‑148 (MQ=255)
                          tAAAGGAATCAAGAACCATTCCGGATCCTCTTTCGCGGGAAGAATTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATg                                                                                                                       >  2:114685/1‑148 (MQ=255)
                            aaGGAATCAAGAACCATTCCGGATCCTCTTTCGCGGGAAGAATTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATggg                                                                                                                     >  2:2995804/1‑148 (MQ=255)
                            aaGGAATCAAGAACCATTCCGGATCCTCTTTCGCGGGAAGAATTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATggg                                                                                                                     <  2:1398256/148‑1 (MQ=255)
                               gAATCAAGAACCATTCCGGATCCTCTTTCGCGGGAAGAATTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGggagg                                                                                                                  <  1:348883/148‑1 (MQ=255)
                                  tCAAGAACCATTCCGGATCCTCTTTCGCGGGAAGAATTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACa                                                                                                               >  2:3084327/1‑148 (MQ=255)
                                   cAAGAACCATTCCGGATCCTCTTTCGCGGGAAGAATTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACAt                                                                                                              <  1:1061227/148‑1 (MQ=255)
                                   cAAGAACCATTCCGGATCCTCTTTCGCGGGAAGAATTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACAt                                                                                                              <  1:1931502/148‑1 (MQ=255)
                                     aGAACCATTCCGGATCCTCTTTCGCGGGAAGAATTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAg                                                                                                            >  1:2826877/1‑148 (MQ=255)
                                       aaCCATTCCGGATCCTCTTTCGCGGGAAGAATTTATCAGTCTTATCGATGAGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGAt                                                                                                          <  2:2653226/148‑1 (MQ=255)
                                        aCCATTCCGGATCCTCTTTCGCGGGAAGAATTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATc                                                                                                         >  1:9006/1‑148 (MQ=255)
                                         ccATTCCGGATCCTCTTTCGCGGGAAGAATTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCt                                                                                                        >  1:1619075/1‑148 (MQ=255)
                                            ttCCGGATCCTCTTTCGCGGGAAGAATTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGaa                                                                                                     <  1:2167866/148‑1 (MQ=255)
                                            ttCCGGATCCTCTTTCGCGGGAAGAATTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGaa                                                                                                     <  1:760686/148‑1 (MQ=255)
                                              ccGGATCCTCTTTCGCGGGAAGAATTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGaaaa                                                                                                   <  2:2035857/148‑1 (MQ=255)
                                               cGGATCCTCTTTCGCGGGAAGAATTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGaaaaa                                                                                                  <  2:1267084/148‑1 (MQ=255)
                                                   tCCTCTTTCGCGGGAAGAATTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTCGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGa                                                                                              <  1:1007655/148‑1 (MQ=255)
                                                    ccTCTTTCGCGGGAAGAATTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGaa                                                                                             >  2:1743696/1‑148 (MQ=255)
                                                      tctTTCGCGGGAAGAATTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACa                                                                                           <  1:1485594/148‑1 (MQ=255)
                                                      tctTTCGCGGGAAGAATTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACa                                                                                           <  1:620473/148‑1 (MQ=255)
                                                       ctTTCGCGGGAAGAATTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACaa                                                                                          >  1:1148184/1‑148 (MQ=255)
                                                       ctTTCGCGGGAAGAATTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACaa                                                                                          >  1:1118004/1‑148 (MQ=255)
                                                         ttCGCGGGAAGAATTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAatg                                                                                        >  1:1377933/1‑148 (MQ=255)
                                                          tCGCGGGAAGAATTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAatga                                                                                       <  1:2395642/148‑1 (MQ=255)
                                                          tCGCGGGAAGAATTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAatga                                                                                       <  1:3071300/148‑1 (MQ=255)
                                                          tCGCGGGAAGAATTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAatga                                                                                       <  1:2027458/148‑1 (MQ=255)
                                                               ggAAGAATTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATc                                                                                  <  2:2161983/148‑1 (MQ=255)
                                                               ggAAGAATTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATc                                                                                  >  1:2882533/1‑148 (MQ=255)
                                                                 aagaaTTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTATATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAATTGGAACAATGATGATCAg                                                                                >  2:457035/1‑148 (MQ=255)
                                                                  agaaTTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAgg                                                                               <  1:2247180/148‑1 (MQ=255)
                                                                   gaaTTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGga                                                                              <  1:1601348/148‑1 (MQ=255)
                                                                    aaTTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGgag                                                                             <  1:433105/148‑1 (MQ=255)
                                                                       ttATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGaaa                                                                          <  1:1234139/148‑1 (MQ=255)
                                                                          tCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAAttt                                                                       <  2:2815615/148‑1 (MQ=255)
                                                                          tCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAAttt                                                                       <  1:88343/148‑1 (MQ=255)
                                                                          tCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTAATACTGGAGTTAGCCCTGGTGAGCTGTGTGAACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAAttt                                                                       >  2:1880772/1‑148 (MQ=255)
                                                                           ccGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTa                                                                      <  1:1977815/148‑1 (MQ=255)
                                                                               cTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCaa                                                                  <  2:1703454/148‑1 (MQ=255)
                                                                               cTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCaa                                                                  <  2:2458006/148‑1 (MQ=255)
                                                                               cTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCaa                                                                  >  2:2703960/1‑148 (MQ=255)
                                                                               cTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCaa                                                                  >  2:329759/1‑148 (MQ=255)
                                                                                 tATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCaaaa                                                                >  2:561957/1‑148 (MQ=255)
                                                                                   tCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGa                                                              <  2:706431/148‑1 (MQ=255)
                                                                                       tGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGt                                                          >  2:1430394/1‑148 (MQ=255)
                                                                                       tGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGt                                                          <  2:3171931/148‑1 (MQ=255)
                                                                                           tGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCa                                                      >  2:920555/1‑148 (MQ=255)
                                                                                           tGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCa                                                      <  1:277871/148‑1 (MQ=255)
                                                                                           tGCAGAAATCAGAAAGCAATAAATTTATGGTGTGTTTCTGTTTATTCTTGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACAGTTTCa                                                      <  1:139195/148‑1 (MQ=255)
                                                                                             cAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACg                                                    >  1:1297350/1‑148 (MQ=255)
                                                                                                  aTCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCTTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTAcc                                               >  1:710583/1‑148 (MQ=255)
                                                                                                  aTCAGCAAGCAAAAAATTTATGGTGTGATACTGTTTATACTGGAGTAAGCCTTGGTGAGCTGTGTGCACTTGGATTGGAGGACATAGATCTGAAAAATGTAACAATGATGATCAGGAGAAATTTAGCAAAAGACAGTTTCACGGTAcc                                               >  1:2474814/1‑148 (MQ=255)
                                                                                                   tCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCa                                              >  2:823906/1‑148 (MQ=255)
                                                                                                    cAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCaa                                             >  2:2288954/1‑148 (MQ=255)
                                                                                                    cAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCaa                                             >  2:2289192/1‑148 (MQ=255)
                                                                                                    cAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCaa                                             >  2:578035/1‑148 (MQ=255)
                                                                                                        aagcaaAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAac                                         <  1:2130498/148‑1 (MQ=255)
                                                                                                          gcaaAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAacac                                       <  2:2410826/148‑1 (MQ=255)
                                                                                                           caaAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAAcaca                                      <  2:2422435/148‑1 (MQ=255)
                                                                                                           caaAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAAcaca                                      <  2:762461/148‑1 (MQ=255)
                                                                                                             aaaaaTTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAgg                                    >  2:662828/1‑148 (MQ=255)
                                                                                                              aaaaTTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGc                                   >  1:2393381/1‑148 (MQ=255)
                                                                                                               aaaTTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCg                                  <  1:2979878/148‑1 (MQ=255)
                                                                                                                aaTTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCgg                                 <  1:1729109/148‑1 (MQ=255)
                                                                                                                aaTTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCgg                                 <  1:54123/148‑1 (MQ=255)
                                                                                                                aaTTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCgg                                 <  1:2322882/148‑1 (MQ=255)
                                                                                                                 aTTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCggg                                >  1:533554/1‑148 (MQ=255)
                                                                                                                  tttATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGa                               >  1:2719524/1‑148 (MQ=255)
                                                                                                                   ttATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGaa                              >  1:3138662/1‑148 (MQ=255)
                                                                                                                         tgtgtTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATc                        >  1:2660002/1‑148 (MQ=255)
                                                                                                                           tgtTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCgg                      <  1:2947139/148‑1 (MQ=255)
                                                                                                                               tCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGtca                  >  1:7917/1‑148 (MQ=255)
                                                                                                                                cTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGtcat                 <  1:2392884/148‑1 (MQ=255)
                                                                                                                                 tGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGtcatt                <  2:3220346/148‑1 (MQ=255)
                                                                                                                                  gTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGtcattc               <  2:2565844/148‑1 (MQ=255)
                                                                                                                                    ttATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGtcattcat             <  1:1698030/148‑1 (MQ=255)
                                                                                                                                    ttATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGtcattcat             >  2:1268056/1‑148 (MQ=255)
                                                                                                                                    ctataCCGGAGTGAGCCCTGGCGAGCTGTGTGCACTGGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGtcattcat             <  1:1835378/147‑1 (MQ=255)
                                                                                                                                       taCTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCtt          <  2:3279362/148‑1 (MQ=255)
                                                                                                                                       taCTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCtt          >  1:608700/1‑148 (MQ=255)
                                                                                                                                            gAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTaa     >  2:167151/1‑148 (MQ=255)
                                                                                                                                               ttCGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGcc  >  2:1466130/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                                 
TAATCCGTTTAAGGGGATCTCGCCATTAAAGGAATCAAGAACCATTCCGGATCCTCTTTCGCGGGAAGAATTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTCGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCC  >  CP009273/1646227‑1646517

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: