Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,646,372 C→A 100% 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,3720CA97.5% 739.5 / 8.3 235pseudogene (565/1158 nt)intQpseudogene, Qin prophage; phage integrase family;Phage or Prophage Related
Reads supporting (aligned to +/- strand):  ref base C (3/3);  new base A (134/95);  total (137/98)
Fisher's exact test for biased strand distribution p-value = 6.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 235 total aligned reads displayed.

TAATCCGTTTAAGGGGATCTCGCCATTAAAGGAATCAAGAACCATTCCGGATCCTCTTTCGCGGGAAGAATTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTCGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAG  >  CP009273/1646227‑1646515
                                                                                                                                                 |                                                                                                                                               
tAATCCGTTTAAGGGGATCTCGCCATTAAAGGAATCAAGAACCATTCCGGATCCTCTTTCGCGGGAAGAATTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGc                                                                                                                                               >  2:3523355/1‑148 (MQ=255)
tAATCCGTTTAAGGGGATCTCGCCATTAAAGGAATCAAGAACCATTCCGGATCCTCTTTCGCGGGAAGAATTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGc                                                                                                                                               >  2:1960901/1‑148 (MQ=255)
tAATCCGTTTAAGGGGATCTCGCCATTAAAGGAATCAAGAACCATTCCGGATCCTCTTTCGCGGGAAGAATTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGc                                                                                                                                               >  2:1759985/1‑148 (MQ=255)
 aaTCCGTTTAAGGGGATCTCGCCATTAAAGGAATCAAGAACCATTCCGGATCCTCTTTCGCGGGAAGAATTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGcc                                                                                                                                              <  2:2567848/148‑1 (MQ=255)
   tCCGTTTAAGGGGATCTCGCCATTAAAGGAATCAAGAACCATTCCGGATCCTCTTTCGCGGGAAGAATTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTCGCCCt                                                                                                                                            >  2:1471509/1‑148 (MQ=255)
       tttAAGGGGATCTCGCCATTAAAGGAATCAAGAACCATTCCGGATCCTCTTTCGCGGGAAGAATTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTg                                                                                                                                        >  2:2283762/1‑148 (MQ=255)
                   tcGCCATTAAAGGAATCAAGAACCATTCCGGATCCTCTTTCGCGGGAAGAATTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCAc                                                                                                                            <  2:3446786/148‑1 (MQ=255)
                   tcGCCATTAAAGGAATCAAGAACCATTCCGGATCCTCTTTCGCGGGAAGAATTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCAc                                                                                                                            <  1:3328280/148‑1 (MQ=255)
                    cGCCATTAAAGGAATCAAGAACCATTCCGGATCCTCTTTCGCGGGAAGAATTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACt                                                                                                                           >  1:808462/1‑148 (MQ=255)
                     gCCATTAAAGGAATCAAGAACCATTCCGGATCCTCTTTCGCGGGAAGAATTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACtt                                                                                                                          <  2:2298445/148‑1 (MQ=255)
                        aTTAAAGGAATCAAGAACCATTCCGGATCCTCTTTCGCGGGAAGAATTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGa                                                                                                                       <  2:937739/148‑1 (MQ=255)
                         ttAAAGGAATCAAGAACCATTCCGGATCCTCTTTCGCGGGAAGAATTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGAt                                                                                                                      <  1:1653214/148‑1 (MQ=255)
                          tAAAGGAATCAAGAACCATTCCGGATCCTCTTTCGCGGGAAGAATTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATg                                                                                                                     >  1:1594118/1‑148 (MQ=255)
                           aaaGGAATCAAGAACCATTCCGGATCCTCTTTCGCGGGAAGAATTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATgg                                                                                                                    <  2:541186/148‑1 (MQ=255)
                             aGGAATCAAGAACCATTCCGGATCCTCTTTCGCGGGAAGAATTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGgga                                                                                                                  <  1:350276/148‑1 (MQ=255)
                              ggAATCAAGAACCAATCCGGATCCTCTTTCGCGGGAAGAATTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTAATTGTGTGTTTCTGTTTATACTGGATTTAGCCCTGGTGAGCTGTGTGCACTTGGATGggag                                                                                                                 >  1:1139427/1‑148 (MQ=255)
                               gAATCAAGAACCATTCCGGATCCTCTTTCGCGGGAAGAATTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGggagg                                                                                                                >  2:3205151/1‑148 (MQ=255)
                                 aTCAAGAACCATTCCGGATCCTCTTTCGCGGGAAGAATTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGAc                                                                                                              <  2:491773/148‑1 (MQ=255)
                                   cAAGAACCATTCCGGATCCTCTTTCGCGGGAAGAATTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACAt                                                                                                            <  1:3170718/148‑1 (MQ=255)
                                   cAAGAACCATTCCGGATCCTCTTTCGCGGGAAGAATTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACAt                                                                                                            <  2:3124778/148‑1 (MQ=255)
                                     aGAACCATTCCGGATCCTCTTTCGCGGGAAGAATTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAg                                                                                                          <  1:3097359/148‑1 (MQ=255)
                                      gAACCATTCCGGATCCTCTTTCGCGGGAAGAATTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGa                                                                                                         <  2:1125012/148‑1 (MQ=255)
                                      gAACCATTCCGGATCCTCTTTCGCGGGAAGAATTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGa                                                                                                         >  2:1231037/1‑148 (MQ=255)
                                         ccATTCCGGATCCTCTTTCGCGGGAAGAATTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCt                                                                                                      >  2:376322/1‑148 (MQ=255)
                                           agtCCGGATCCTCTTTCGCGGGAAGAATTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGa                                                                                                    <  2:2361469/146‑1 (MQ=255)
                                           aTTCCGGATCCTCTTTCGCGGGAAGAATTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGa                                                                                                    >  2:1505857/1‑148 (MQ=255)
                                            ttCCGGATCCTCTTTCGCGGGAAGAATTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGaa                                                                                                   >  2:2141926/1‑148 (MQ=255)
                                               cGGATCCTCTTTCGCGGGAAGAATTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGaaaaa                                                                                                >  1:1510140/1‑148 (MQ=255)
                                                 gATCCTCTTTCGCGGGAAGAATTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATg                                                                                              >  2:1380711/1‑148 (MQ=255)
                                                  aTCCTCTTTCGCGGGAAGAATTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATgg                                                                                             >  2:1251343/1‑148 (MQ=255)
                                                      tctTTCGCGGGAAGAATTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACa                                                                                         >  1:584302/1‑148 (MQ=255)
                                                       ctTTCGCGGGAAGAATTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACaa                                                                                        <  1:150128/148‑1 (MQ=255)
                                                          tCGCGGGAAGAATTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAatga                                                                                     <  2:1352043/148‑1 (MQ=255)
                                                           cgcgGGAAGAATTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAatgat                                                                                    <  2:2680299/148‑1 (MQ=255)
                                                            gcgGGAAGAATTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAatgatg                                                                                   <  1:900883/148‑1 (MQ=255)
                                                             cgGGAAGAATTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAAtgatga                                                                                  <  2:2966960/148‑1 (MQ=255)
                                                               ggAAGAATTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAatg                                                                                      >  2:16023/1‑142 (MQ=255)
                                                               ggAAGAATTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATc                                                                                <  1:1231024/148‑1 (MQ=255)
                                                                gaagaaTTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCa                                                                               >  2:2416589/1‑148 (MQ=255)
                                                                gaagaaTTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCa                                                                               >  2:2996504/1‑148 (MQ=255)
                                                                 aagaaTTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAg                                                                              >  1:2018674/1‑148 (MQ=255)
                                                                 aagaaTTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAg                                                                              >  2:1560074/1‑148 (MQ=255)
                                                                   gaaTTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGga                                                                            <  1:1683651/148‑1 (MQ=255)
                                                                     aTTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGgaga                                                                          >  2:402367/1‑148 (MQ=255)
                                                                      tttATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGaa                                                                         >  1:270737/1‑148 (MQ=255)
                                                                        tATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAAt                                                                       <  1:873960/148‑1 (MQ=255)
                                                                         aTCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAAtt                                                                      <  2:329031/148‑1 (MQ=255)
                                                                          tCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAAttt                                                                     >  2:1501805/1‑148 (MQ=255)
                                                                           ccGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTa                                                                    <  1:1049050/148‑1 (MQ=255)
                                                                             gTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGc                                                                  >  1:1294938/1‑148 (MQ=255)
                                                                              tCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCa                                                                 >  1:2338919/1‑148 (MQ=255)
                                                                              tCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCa                                                                 >  1:3428540/1‑148 (MQ=255)
                                                                               cTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCaa                                                                >  1:2659979/1‑148 (MQ=255)
                                                                                 tATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCaaaa                                                              >  1:2329248/1‑148 (MQ=255)
                                                                                   tCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGa                                                            <  2:399994/148‑1 (MQ=255)
                                                                                      aTGCGTGCAGAAATCAGCAAGCAAATAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCg                                                         <  2:1066551/148‑1 (MQ=255)
                                                                                      aTGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCg                                                         >  2:1574430/1‑148 (MQ=255)
                                                                                      aTGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCg                                                         >  1:260637/1‑148 (MQ=255)
                                                                                      aTGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGACCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCg                                                         >  2:126052/1‑148 (MQ=255)
                                                                                       tGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGt                                                        >  1:682731/1‑148 (MQ=255)
                                                                                       tGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGt                                                        >  1:132662/1‑148 (MQ=255)
                                                                                          gTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTc                                                     <  2:2469972/148‑1 (MQ=255)
                                                                                           tGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCa                                                    >  1:441660/1‑148 (MQ=255)
                                                                                            gCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCAc                                                   <  1:587644/148‑1 (MQ=255)
                                                                                            gCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCAc                                                   >  1:1598182/1‑148 (MQ=255)
                                                                                                  aTCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTAcc                                             >  2:2521023/1‑148 (MQ=255)
                                                                                                   tCAGCAATCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCa                                            <  2:3096699/148‑1 (MQ=255)
                                                                                                   tCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCa                                            <  1:3072452/148‑1 (MQ=255)
                                                                                                   tCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCa                                            >  1:2128307/1‑148 (MQ=255)
                                                                                                     agcaagcaAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCaaa                                          >  1:1662615/1‑148 (MQ=255)
                                                                                                     agcaagcaAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCaaa                                          >  2:2202552/1‑148 (MQ=255)
                                                                                                      gcaagcaaAAATTTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCaaaa                                         <  2:3312658/148‑1 (MQ=255)
                                                                                                       caagcaaAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCaaaaa                                        >  2:2825237/1‑148 (MQ=255)
                                                                                                         agcaaAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAaca                                      <  2:927574/148‑1 (MQ=255)
                                                                                                          gcaaAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAATAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAacac                                     <  2:891060/148‑1 (MQ=255)
                                                                                                          gcaaAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAacac                                     <  1:258301/148‑1 (MQ=255)
                                                                                                              aaaaTTTATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGc                                 >  1:2543954/1‑148 (MQ=255)
                                                                                                                  tttATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGa                             >  2:2167677/1‑148 (MQ=255)
                                                                                                                  tttATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGa                             <  1:2335881/148‑1 (MQ=255)
                                                                                                                  tttATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGa                             >  1:2201726/1‑148 (MQ=255)
                                                                                                                   ttATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGTAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGaa                            >  2:2078677/1‑148 (MQ=255)
                                                                                                                   ttATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGac                            >  2:3239434/1‑147 (MQ=255)
                                                                                                                   ttATGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGaa                            >  2:3242752/1‑148 (MQ=255)
                                                                                                                     aTGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAAcc                          >  2:2623891/1‑148 (MQ=255)
                                                                                                                      tGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCa                         <  1:2467967/148‑1 (MQ=255)
                                                                                                                      tGGTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCa                         <  2:2435258/148‑1 (MQ=255)
                                                                                                                       ggTGTGTTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCaa                        <  1:883260/148‑1 (MQ=255)
                                                                                                                          gtgtTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCg                     <  1:1960901/148‑1 (MQ=255)
                                                                                                                          gtgtTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCg                     <  2:859958/148‑1 (MQ=255)
                                                                                                                            gtTTCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCggg                   >  2:562835/1‑148 (MQ=255)
                                                                                                                              ttCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGtc                 >  1:123795/1‑148 (MQ=255)
                                                                                                                               tCTGTTTATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGtca                >  1:1454030/1‑148 (MQ=255)
                                                                                                                                   tttATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGtcattca            >  1:2976770/1‑148 (MQ=255)
                                                                                                                                    ttATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGtcattcat           <  2:1672821/148‑1 (MQ=255)
                                                                                                                                    ttATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGtcattcat           >  1:2166923/1‑148 (MQ=255)
                                                                                                                                    ttATACTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGtcattcat           <  1:1193371/148‑1 (MQ=255)
                                                                                                                                     tataCTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCACAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATc          <  1:1529714/148‑1 (MQ=255)
                                                                                                                                       taCTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCtt        <  1:1560074/148‑1 (MQ=255)
                                                                                                                                        aCTGGAGTTAGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCtta       <  2:705764/148‑1 (MQ=255)
                                                                                                                                         cTGGAGTTCGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCttat      <  1:1471509/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                               
TAATCCGTTTAAGGGGATCTCGCCATTAAAGGAATCAAGAACCATTCCGGATCCTCTTTCGCGGGAAGAATTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTCGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAG  >  CP009273/1646227‑1646515

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: