Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 310,198 C→A 45.2% pseudogene (131/872 nt) eaeH → pseudogene, attaching and effacing protein homology;factor; Not classified

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273310,1980CA45.2% 184.8 / 299.0 246pseudogene (131/872 nt)eaeHpseudogene, attaching and effacing protein homology;factor; Not classified
Reads supporting (aligned to +/- strand):  ref base C (68/67);  new base A (57/54);  total (125/121)
Fisher's exact test for biased strand distribution p-value = 8.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.46e-01

Only 100 of 250 total aligned reads displayed.

TCAGGAAATTATCTCAATGTCACATTATAAAACAGGTCATAAACAACCACGATTTCGTTATTCAGTTCTGGCCCGCTGCGTGGCGTGGGCAAATATCTCTGTTCAGGTTCTTTTTCCACTCGCTGTCACCTTTACCCCAGTAATGGCGGCACGTGCGCAGCATGCGGTTCAGCCACGGTTGAGCATGGGAAATACTACGGTAACTGCTGATAATAACGTGGAGAAAAATGTCGCGTCGTTTGCCGCAAATGCCGGGACATTTTTAAGCAGTCAGCCAGATAGCGATGCGACA  >  CP009273/310052‑310343
                                                                                                                                                  |                                                                                                                                                 
tcaGGAAATTATCTCAATGTCACATTATAAAACAGGTCATAAACAACCACGATTTCGTTATTCAGTTCTGGCCCGCTGCGTGGCGTGGGCAAATATCTCTGTTCAGGTTCTTTTTCCACTCGCTGTCACCTTTACCCCAGTAATggcg                                                                                                                                                  >  2:102572/1‑148 (MQ=255)
       aTTATCTCAATGTCACATTATAAAACAGGTCATAAACAACCACGATTTCGTTATTCAGTTCTGGCCCGCTGCGTGGCGTGGGCAAATATCTCTGTTCAGGTTCTTTTTCCACTCGCTGTCACCTTTACCCCAGTAATGGAGGCACGTg                                                                                                                                           <  1:2885101/148‑1 (MQ=255)
       aTTATCTCAATGTCACATTATAAAACAGGTCATAAACAACCACGATTTCGTTATTCAGTTCTGGCCCGCTGCGTGGCGTGGGCAAATATCTCTGTTCAGGTTCTTTTTCCACTCGCTGTCACCTTTACCCCAGTAATGGAGGCACGTg                                                                                                                                           <  1:1408324/148‑1 (MQ=255)
         tATCTCAATGTCACATTATAAAACAGGTCATAAACAACCACGATTTCGTTATTCAGTTCTGGCCCGCTGCGTGGCGTGGGCAAATATCTCTGTTCAGGTTCTTTTTCCACTCGCTGTCACCTTTACCCCAGTAATGGAGGCACGTgcg                                                                                                                                         <  2:3183445/148‑1 (MQ=255)
          aTCTCAATGTCACATTATAAAACAGGTCATAAACAACCACGATTTCGTTATTCAGTTCTGGCCCGCTGCGTGGCGTGGGCAAATATCTCTGTTCAGGTTCTTTTTCCACTCGCTGTCACCTTTACCCCAGTAATGGAGGCACGTgcgc                                                                                                                                        <  1:2507220/148‑1 (MQ=255)
              cAATGTCACATTATAAAACAGGTCATAAACAACCACGATTTCGTTATTCAGTTCTGGCCCGCTGCGTGGCGTGGGCAAATATCTCTGTTCAGGTTCTTTTTCCACTCGCTGTCACCTTTACCCCAGTAATGGAGGCACGTGCgcagca                                                                                                                                    <  2:2994677/148‑1 (MQ=255)
                aTGTCACATTATAAAACAGGTCATAAACAACCACGATTTCGTTATTCAGTTCTGGCCCGCTGCGTGGCGTGGGCAAATATCTCTGTTCAGGTTCTTTTTCCACTCGCTGTCACCTTTACCCCAGTAATGGAGGCACGTGCGCAGCATg                                                                                                                                  >  2:2224748/1‑148 (MQ=255)
                    cacaTTATAAAACAGGTCATAAACAACCACGATTTCGTTATTCAGTTCTGGCCCGCTGCGTGGCGTGGGCAAATATCTCTGTTCAGGTTCTTTTTCCACTCGCTGTCACCTTTACCCCAGTAATGGCGGCACGTGCGCAGCATGCGGt                                                                                                                              <  1:1992425/148‑1 (MQ=255)
                      caTTATAAAACAGGTCATAAACAACCACGATTTCGTTATTCAGTTCTGGCCCGCTGCGTGGCGTGGGCAAATATCTCTGTTCAGGTTCTTTTTCCACTCGCTGTCACCTTTACCCCAGTAATGGCGGCACGTGCGCAGCATGCGGTTc                                                                                                                            >  1:822756/1‑148 (MQ=255)
                      caTTATAAAACAGGTCATAAACAACCACGATTTCGTTATTCAGTTCTGGCCCGCTGCGTGGCGTGGGCAAATATCTCTGTTCAGGTTCTTTTTCCACTCGCTGTCACCTTTACCCCAGTAATGGAGGCACGTGCGCAGCATGCGGTTc                                                                                                                            <  1:757736/148‑1 (MQ=255)
                      caTTATAAAACAGGTCATAAACAACCACGATTTCGTTATTCAGTTCTGGCCCGCTGCGTGGCGTGGGCAAATATCTCTGTTCAGGTTCTTTTTCCACTCGCTGTCACCTTTACCCCAGTAATGGAGGCACGTGCGCAGCATGCGGTTc                                                                                                                            <  1:1552370/148‑1 (MQ=255)
                       aTTATAAAACAGGTCATAAACAACCACGATTTCGTTATTCAGTTCTGGCCCGCTGCGTGGCGTGGGCAAATATCTCTGTTCAGGTTCTTTTTCCACTCGCTGTCACCTTTACCCCAGTAATGGAGGCACGTGCGCAGCATGCGGTTCa                                                                                                                           >  2:186472/1‑148 (MQ=255)
                             aaaCAGGTCATAAACAACCACGATTTCGTTATTCAGTTCTGGCCCGCTGCGTGGCGTGGGCAAATATCTCTGTTCAGGTTCTTTTTCCACTCGCTGTCACCTTTACCCCAGTAATGGAGGCACGTGCGCAGCATGCGGTTCAGCCACg                                                                                                                     >  1:3063986/1‑148 (MQ=255)
                               aCAGGTCATAAACAACCACGATTTCGTTATTCAGTTCTCGCCCGCTTCGTGGCGTGGGCAAATATCTCTGTTCAGGTTCTTTTTCCACTCGCTGTCACCTTTACCCCAGTAATGGCGGCACGTGCGCAGCATGCGGTTCAGCCACGGt                                                                                                                   >  2:1546582/1‑148 (MQ=255)
                                cAGGTCATAAACAACCACGATTTCGTTATTCAGTTCTGGCCCGCTGCGTGGCGTGGGCAAATATCTCTGTTCAGGTTCTTTTTCCACTCGCTGTCACCTTTACCCCAGTAATGGCGGCACGTGCGCAGCATGCGGTTCAGCCACGGtt                                                                                                                  <  1:2935167/148‑1 (MQ=255)
                                  ggTCATAAACAACCACGATTTCGTTATTCAGTTCTGGCCCGCTGCGTGGCGTGGGCAAATATCTCTGTTCAGGTTCTTTTTCCACTCGCTGTCACCTTTACCCCAGTAATGGCGGCACGTGCGCAGCATGCGGTTCAGCCACGGTTGa                                                                                                                >  2:2729656/1‑148 (MQ=255)
                                  ggTCATAAACAACCACGATTTCGTTATTCAGTTCTGGCCCGCTGCGTGGCGTGGGCAAATATCTCTGTTCAGGTTCTTTTTCCACTCGCTGTCACCTTTACCCCAGTAATGGCGGCACGTGCGCAGCATGCGGTTCAGCCACGGTTGa                                                                                                                >  1:561799/1‑148 (MQ=255)
                                    tCATAAACAACCACGATTTCGTTATTCAGTTCTGGCCCGCTGCGTGGCGTGGGCAAATATCTCTGTTCAGGTTCTTTTTCCACTCGCTGTCACCTTTACCCCAGTAATGGAGGCACGTGCGCAGCATGCGGTTCAGCCACGGTTGAGc                                                                                                              >  2:2626841/1‑148 (MQ=255)
                                        aaaCAACCACGATTTCGTTATTCAGTTCTGGCCCGCTGCGTGGCGTGGGCAAATATCTCTGTTCAGGTTCTTTTTCCACTCGCTGTCACCTTTACCCCAGTAATGGAGGCACGTGCGCAGCATGCGGTTCAGCCACGGTTGAGCATgg                                                                                                          >  1:2530983/1‑148 (MQ=255)
                                          acaacCACGATTTCGTTATTCAGTTCTGGCCCGCTGCGTGGCGTGGGCAAATATCTCTGTTCAGGTTCTTTTTCCACTCGCTGTCACCTTTACCCCAGTAATGGCGGCACGTGCGCAGCATGCGGTTCAGCCACGGTTGAGCATGGGa                                                                                                        >  1:2873936/1‑148 (MQ=255)
                                          acaacCACGATTTCGTTATTCAGTTCTGGCCCGCTGCGTGGCGTGGGCAAATATCTCTGTTCAGGTTCTTTTTCCACTCGCTGTCACCTTTACCCCAGTAATGGCGGCACGTGCGCAGCATGCGGTTCAGCCACGGTTGAGCATGGGa                                                                                                        >  1:2128035/1‑148 (MQ=255)
                                             acCACGATTTCGTTATTCAGTTCTGGCCCGCTGCGTGGCGTGGGCAAATATCTCTGTTCAGGTTCTTTTTCCACTCGCTGTCACCTTTACCCCAGTAATGGCGGCACGTGCGCAGCATGCGGTTCAGCCACGGTTGAGCATGGGAAAt                                                                                                     >  2:2406333/1‑148 (MQ=255)
                                             acCACGATTTCGTTATTCAGTTCTGGCCCGCTGCGTGGCGTGGGCAAATATCTCTGTTCAGGTTCTTTTTCCACTCGCTGTCACCTTTACCCCAGTAATGGCGGCACGTGCGCAGCATGCGGTTCAGCCACGGTTGAGCATGGGAAAt                                                                                                     >  1:351270/1‑148 (MQ=255)
                                               cACGATTTCGTTATTCAGTTCTGGCCCGCTGCGTGGCGTGGGCAAATATCTCTGTTCAGGTTCTTTTTCCACTCGCTGTCACCTTTACCCCAGTAATGGAGGCACGTGCGCAGCATGCGGTTCAGCCACGGTTGAGCATGGGAAAtac                                                                                                   >  2:339253/1‑148 (MQ=255)
                                                aCGATTTCGTTATTCAGTTCTGGCCCGCTGCGTGGCGTGGGCAAATATCTCTGTTCAGGTTCTTTTTCCACTCGCTGTCACCTTTACCCCAGTAATGGAGGCACGTGCGCAGCATGCGGTTCAGCCACGGTTGAGCATGGGAAAtact                                                                                                  <  2:2277603/148‑1 (MQ=255)
                                                 cGATTTCGTTATTCAGTTCTGGCCCGCTGCGTGGCGTGGGCAAATATCTCTGTTCAGGTTCTTTTTCCACTCGCTGTCACCTTTACCCCAGTAATGGAGGCACGTGCGCAGCATGCGGTTCAGCCACGGTTGAGCATGGGAAAtacta                                                                                                 >  2:870865/1‑148 (MQ=255)
                                                    tttCGTTATTCAGTTCTGGCCCGCTGCGTGGCGTGGGCAAATATCTCTGTTCAGGTTCTTTTTCCACTCGCTGTCACCTTTACCCCAGTAATGGAGGCACGTGCGCAGCATGCGGTTCAGCCACGGTTGAGCATGGGAAATACTACgg                                                                                              >  1:937622/1‑148 (MQ=255)
                                                    tttCGTTATTCAGTTCTGGCCCGCTGCGTGGCGTGGGCAAATATCTCTGTTCAGGTTCTTTTTCCACTCGCTGTCACCTTTACCCCAGTAATGGAGGCACGTGCGCAGCATGCGGTTCAGCCACGGTTGAGCATGGGAAATACTACgg                                                                                              >  1:1572114/1‑148 (MQ=255)
                                                      tCGTTATTCAGTTCTGGCCCGCTGCGTGGCGTGGGCAAATATCTCTGTTCAGGTTCTTTTTCCACTCGCTGTCACCTTTACCCCAGTAATGGCGGCACGTGCGCAGCATGCGGTTCAGCCACGGTTGAGCATGGGAAATACTACGGTa                                                                                            >  1:392222/1‑148 (MQ=255)
                                                      tCGTTATTCAGTTCTGGCCCGCTGCGTGGCGTGGGCAAATATCTCTGTTCAGGTTCTTTTTCCACTCGCTGTCACCTTTACCCCAGTAATGGAGGCACGTGCGCAGCATGCGGTTCAGCCACGGTTGAGCATGGGAAATACTACGGTa                                                                                            >  1:282289/1‑148 (MQ=255)
                                                       cGTTATTCAGTTCTGGCCCGCTGCGTGGCGTGGGCAAATATCTCTGTTCAGGTTCTTTTTCCACTCGCTGTCACCTTTACCCCAGTAATGGCGGCACGTGCGCAGCATGCGGTTCAGCCACGGTTGAGCATGGGAAATACTACGGTaa                                                                                           <  1:1301515/148‑1 (MQ=255)
                                                          tATTCAGTTCTGGCCCGCTGCGTGGCGTGGGCAAATATCTCTGTTCAGGTTCTTTTTCCACTCGCTGTCACCTTTACCCCAGTAATGGCGGCACGTGCGCAGCATGCGGTTCAGCCACGGTTGAGCATGGGAAATACTACGGTAActg                                                                                        <  2:2846503/148‑1 (MQ=255)
                                                                 ttCTGGCCCGCTGCGTGGCGTGGGCAAATATCTCTGTTCAGGTTCTTTTTCCACTCGCTGTCACCTTTACCCCAGTAATGGCGGCACGTGCGCAGCATGCGGTTCAGCCACGGTTGAGCATGGGAAATACTACGGTAACTGCTGataa                                                                                 >  2:471356/1‑148 (MQ=255)
                                                                 ttCTGGCCCGCTGCGTGGCGTGGGCAAATATCTCTGTTCAGGTTCTTTTTCCACTCGCTGTCACCTTTACCCCAGTAATGGCGGCACGTGCGCAGCATGCGGTTCAGCCACGGTTGAGCATGGGAAATACTACGGTAACTGCTGataa                                                                                 <  2:976149/148‑1 (MQ=255)
                                                                 ttCTGGCCCGCTGCGTGGCGTGGGCAAATATCTCTGTTCAGGTTCTTTTTCCACTCGCTGTCACCTTTACCCCAGTAATGGAGGCACGTGCGCAGCATGCGGTTCAGCCACGGTTGAGCATGGGAAATACTACGGTAACTGCTGataa                                                                                 >  1:2403848/1‑148 (MQ=255)
                                                                     ggCCCGCTGCGTGGCGTGGGCAAATATCTCTGTTCAGGTTCTTTTTCCACTCGCTGTCACCTTTACCCCAGTAATGGCGGCACGTGCGCAGCATGCGGTTCAGCCACGGTTGAGCATGGGAAATACTACGGTAACTGCTGATAATAAc                                                                             >  2:1357396/1‑148 (MQ=255)
                                                                     ggCCCGCTGCGTGGCGTGGGCAAATATCTCTGTTCAGGTTCTTTTTCCACTCGCTGTCACCTTTACCCCAGTAATGGCGGCACGTGCGCAGCATGCGGTTCAGCCACGGTTGAGCATGGGAAATACTACGGTAACTGCTGATAATAAc                                                                             >  2:2477965/1‑148 (MQ=255)
                                                                      gCCCGCTGCGTGGCGTGGGCAAATATCTCTGTTCAGGTTCTTTTTCCACTCGCTGTCACCTTTACCCCAGTAATGGCGGCACGTGCGCAGCATGCGGTTCAGCCACGGTTGAGCATGGGAAATACTACGGTAACTGCTGATAATAACg                                                                            <  1:1843191/148‑1 (MQ=255)
                                                                       cccGCTGCGTGGCGTGGGCAAATATCTCTGTTCAGGTTCTTTTTCCACTCGCTGTCACCTTTACCCCAGTAATGGAGGCACGTGCGCAGCATGCGGTTCAGCCACGGTTGAGCATGGGAAATACTACGGTAACTGCTGATAATAACGt                                                                           <  1:3127698/148‑1 (MQ=255)
                                                                       cccGCTGCGTGGCGTGGGCAAATATCTCTGTTCAGGTTCTTTTTCCACTCGCTGTCACCTTTACCCCAGTAATGGAGGCACGTGCGCAGCATGCGGTTCAGCCACGGTTGAGCATGGGAAATACTACGGTAACTGCTGATAATAACGt                                                                           <  1:3152463/148‑1 (MQ=255)
                                                                         cGCTGCGTGGCGTGGGCAAATATCTCTGTTCAGGTTCTTTTTCCACTCGCTGTCACCTTTACCCCAGTAATGGCGGCACGTGCGCAGCATGCGGTTCAGCCACGGTTGAGCATGGGAAATACTACGGTAACTGCTGATAATAACGTgg                                                                         <  2:2957421/148‑1 (MQ=255)
                                                                           cTGCGTGGCGTGGGCAAATATCTCTGTTCAGGTTCTTTTTCCACTCGCTGTCACCTTTACCCCAGTAATGGCGGCACGTGCGCAGCATGCGGTTCAGCCACGGTTGAGCATGGGAAATACTACGGTAACTGCTGATAATAACGTGgag                                                                       >  1:382044/1‑148 (MQ=255)
                                                                           cTGCGTGGCGTGGGCAAATATCTCTGTTCAGGTTCTTTTTCCACTCGCTGTCACCTTTACCCCAGTAATGGCGGCACGTGCGCAGCATGCGGTTCAGCCACGGTTGAGCATGGGAAATACTACGGTAACTGCTGATAATAACGTGgag                                                                       >  1:2205399/1‑148 (MQ=255)
                                                                            tGCGTGGCGTGGGCAAATATCTCTGTTCAGGTTCTTTTTCCACTCGCTGTCACCTTTACCCCAGTAATGGCGGCACGTGCGCAGCATGCGGTTCAGCCACGGTTGAGCATGGGAAATACTACGGTAACTGCTGATAATAACGTGgaga                                                                      >  1:1891407/1‑148 (MQ=255)
                                                                             gcgtggcgtgGGCAAATATCTCTGTTCAGGTTCTTTTTCCACTCGCTGTCACCTTTACCCCAGTAATGGCGGCACGTGCGCAGCATGCGGTTCAGCCACGGTTGAGCATGGGAAATACTACGGTAACTGCTGATAATAACGTGGAGaa                                                                     >  1:279433/1‑148 (MQ=255)
                                                                             gcgtggcgtgGGCAAATATCTCTGTTCAGGTTCTTTTTCCACTCGCTGTCACCTTTACCCCAGTAATGGCGGCACGTGCGCAGCATGCGGTTCAGCCACGGTTGAGCATGGGAAATACTACGGTAACTGCTGATAATAACGTGGAGaa                                                                     >  1:2496895/1‑148 (MQ=255)
                                                                                 ggcgtggGCAAATATCTCTGTTCAGGTTCTTTTTCCACTCGCTGTCACCTTTACCCCAGTAATGGCGGCACGTGCGCAGCATGCGGTTCAGCCACGGTTGAGCATGGGAAATACTACGGTAACTGCTGATAATAACGTGGAGAAAAAt                                                                 <  2:775934/148‑1 (MQ=255)
                                                                                   cgtggGCAAATATCTCTGTTCAGGTTCTTTTTCCACTCGCTGTCACCTTTACCCCAGTAATGGAGGCACGTGCGCAGCATGCGGTTCAGCCACGGTTGAGCATGGGAAATACTACGGTAACTGCTGATAATAACGTGGAGAAAAATGt                                                               <  2:2040116/148‑1 (MQ=255)
                                                                                     tggGCAAATATCTCTGTTCAGGTTCTTTTTCCACTCGCTGTCACCTTTACCCCAGTAATGGCGGCACGTGCGCAGCATGCGGTTCAGCCACGGTTGAGCATGGGAAATACTACGGTAACTGCTGATAATAACGTGGAGAAAAATGTcg                                                             <  1:3431759/148‑1 (MQ=255)
                                                                                       ggCAAATATCTCTGTTCAGGTTCTTTTTCCACTCGCTGTCACCTTTACCCCAGTAATGGCGGCACGTGCGCAGCATGCGGTTCAGCCACGGTTGAGCATGGGAAATACTACGGTAACTGCTGATAATAACGTGGAGAAAAATGTcgcg                                                           <  1:1084119/148‑1 (MQ=255)
                                                                                         cAAATATCTCTGTTCAGGTTCTTTTTCCACTCGCTGTCACCTTTACCCCAGTAATGGCGGCACGTGCGCAGCATGCGGTTCAGCCACGGTTGAGCATGGGAAATACTACGGTAACTGCTGATAATAACGTGGAGAAAAATGTCGcgtc                                                         <  2:1891755/148‑1 (MQ=255)
                                                                                           aaTATCTCTGTTCAGGTTCTTTTTCCACTCGCTGTCACCTTTACCCCAGTAATGGCGGCACGTGCGCAGCATGCGGTTCAGCCACGGTTGAGCATGGGAAATACTACGGTAACTGCTGATAATAACGTGGAGAAAAATGTCGcgtcgt                                                       <  1:3290031/148‑1 (MQ=255)
                                                                                            atatCTCTGTTCAGGTTCTTTTTCCACTCGCTGTCACCTTTACCCCAGTAATGGCGGCACGTGCGCAGCATGCGGTTCAGCCACGGTTGAGCATGGGAAATACTACGGTAACTGCTGATAATAACGTGGAGAAAAATGTCGCGTCGtt                                                      >  2:1609144/1‑148 (MQ=255)
                                                                                             tatCTCTGTTCAGGTTCTTTTTCCACTCGCTGTCACCTTTACCCCAGTAATGGAGGCACGTGCGCAGCATGCGGTTCAGCCACGGTTGAGCATGGGAAATACTACGGTAACTGCTGATAATAACGTGGAGAAAAATGTCGCGTCGttt                                                     >  1:1430545/1‑148 (MQ=255)
                                                                                              atCTCTGTTCAGGTTCTTTTTCCACTCGCTGTCACCTTTACCCCAGTAATGGCGGCACGTGCGCAGCATGCGGTTCAGCCACGGTTGAGCATGGGAAATACTACGGTAACTGCTGATAATAACGTGGAGAAAAATGTCGCGTCGTTTg                                                    <  1:919847/148‑1 (MQ=255)
                                                                                               tctcTGTTCAGGTTCTTTTTCCACTCGCTGTCACCTTTACCCCAGTAATGGAGGCACGTGCGCAGCATGCGGTTCAGCCACGGTTGAGCATGGGAAATACTACGGTAACTGCTGATAATAACGTGGAGAAAAATGTCGCGTCGTTTGc                                                   <  2:1763707/148‑1 (MQ=255)
                                                                                               tctcTGTTCAGGTTCTTTTTCCACTCGCTGACACCTTTACCCCAGTAATGGCGGCACGTGCGCAGCATGCGGTTCAGCCACGGATGAGCATGGGAAATACTACGGTAACTGCTGATAATAACGTGGAGAAAAATGTCGCGTCGTTTGc                                                   >  2:653631/1‑148 (MQ=255)
                                                                                                  ctGTTCAGGTTCTTTTTCCACTCGCTGTCACCTTTACCCCAGTAATGGAGGCACGTGCGCAGCATGCGGTTCAGCCACGGTTGAGCATGGGAAATACTACGGTAACTGCTGATAATAACGTGGAGAAAAATGTCGCGTCGTTTGCCGc                                                >  2:1659883/1‑148 (MQ=255)
                                                                                                  ctGTTCAGGTTCTTTTTCCACTCGCTGTCACCTTTACCCCAGTAATGGAGGCACGTGCGCAGCATGCGGTTCAGCCACGGTTGAGCATGGGAAATACTACGGTAACTGCTGATAATAACGTGGAGAAAAATGTCGCGTCGTTTGCCGc                                                >  2:1060620/1‑148 (MQ=255)
                                                                                                     ttCAGGTTCTTTTTCCACTCGCTGTCACCTTTACCCCAGTAATGGCGGCACGTGCGCAGCATGCGGTTCAGCCACGGTTGAGCATGGGAAATACTACGGTAACTGCTGATAATAACGTGGAGAAAAATGTCGCGTCGTTTGCCGCaaa                                             <  2:966797/148‑1 (MQ=255)
                                                                                                     ttCAGGTTCTTTTTCCACTCGCTGTCACCTTTACCCCAGTAATGGCGGCACGTGCGCAGCATGCGGTTCAGCCACGGTTGAGCATGGGAAATACTACGGTAACTGCTGATAATAACGTGGAGAAAAATGTCGCGTCGTTTGCCGCaaa                                             >  1:3441806/1‑148 (MQ=255)
                                                                                                     ttCAGGTTCTTTTTCCACTCGCTGTCACCTTTACCCCAGTAATGGAGGCACGTGCGCAGCATGCGGTTCAGCCACGGTTGAGCATGGGAAATACTACGGTAACTGCTGATAATAACGTGGAGAAAAATGTCGCGTCGTTTGCCGCaaa                                             <  2:1232603/148‑1 (MQ=255)
                                                                                                         ggTTCTTTTTCCACTCGCTGTCACCTTTACCCCAGTAATGGCGGCACGTGCGCAGCATGCGGTTCAGCCACGGTTGAGCATGGGAAATACTACGGTAACTGCTGATAATAACGTGGAGAAAAATGTCGCGTCGTTTGCCGCAAATGcc                                         <  1:1627599/148‑1 (MQ=255)
                                                                                                          gTTCTTTTTCCACTCGCTGTCACCTTTACCCCAGTAATGGCGGCACGTGCGCAGCATGCGGTTCAGCCACGGTTGAGCATGGGAAATACTACGGTAACTGCTGATAATAACGTGGAGAAAAATGTCGCGTCGTTTGCCGCAAATGCCg                                        >  2:444994/1‑148 (MQ=255)
                                                                                                           ttCTTTTTCCACTCGCTGTCACCTTTACCCCAGTAATGGCGGCACGTGCGCAGCATGCGGTTCAGCCACGGTTGAGCATGGGAAATACTACGGTCACTTCTGATCATAACGTGGAGAAAAATGTCTCGACGTGTGCCGCAAATgaccg                                       >  2:2663645/1‑144 (MQ=255)
                                                                                                           ttCTTTTTCCACTCGCTGTCACCTTTACCCCAGTAATGGAGGCACGTGCGCAGCATGCGGTTCAGCCACGGTTGAGCATGGGAAATACTACGGTAACTGCTGATAATAACGTGGAGAAAAATGTCGCGTCGTTTGCCGCAAATGCCgg                                       >  1:721447/1‑148 (MQ=255)
                                                                                                              tttttCCACTCGCTGTCACCTTTACCCCAGTAATGGAGGCACGTGCGCAGCATGCGGTTCAGCCACGGTTGAGCATGGGAAATACTACGGTAACTGCTGATAATAACGTGGAGAAAAATGTCGCGTCGTTTGCCGCAAATGCCGGGAc                                    >  1:2714358/1‑148 (MQ=255)
                                                                                                              tttttCCACTCGCTGTCACCTGGACCCCAGTAAGGGCGGCACGGGCGCAGCATGCGGTTCAGCCACGGTTGAGCAGGGGAAATACTACGGTACCTGCTGATAATAACGTGGAGAACAATGTCGCGTCGTGTGCCGCAAATGCCGGGAc                                    <  2:1417560/148‑1 (MQ=255)
                                                                                                                 ttCCACTCGCTGTCACCTTTACCCCAGTAATGGAGGCACGTGCGCAGCATGCGGTTCAGCCACGGTTGAGCATGGGAAATACTACGGTAACTGCTGATAATAACGTGGAGAAAAATGTCGCGTCGTTTGCCGCAAATGCCGGGACAtt                                 <  1:1955107/148‑1 (MQ=255)
                                                                                                                  tCCACTCGCTGTCACCTTTACCCCAGTAATGGCGGCACGTGCGCAGCATGCGGTTCAGCCACGGTTGAGCATGGGAAATACTACGGTAACTGCTGATAATAACGTGGAGAAAAATGTCGCGTCGTTTGCCGCAAATGCCGGGACAttt                                <  2:2128035/148‑1 (MQ=255)
                                                                                                                  tCCACTCGCTGTCACCTTTACCCCAGTAATGGCGGCACGTGCGCAGCATGCGGTTCAGCCACGGTTGAGCATGGGAAATACTACGGTAACTGCTGATAATAACGTGGAGAAAAATGTCGCGTCGTTTGCCGCAAATGCCGGGACAttt                                <  1:1057255/148‑1 (MQ=255)
                                                                                                                      cTCGCTGTCACCTTTACCCCAGTAATGGCGGCACGTGCGCAGCATGCGGTTCAGCCACGGTTGAGCATGGGAAATACTACGGTAACTGCTGATAATAACGTGGAGAAAAATGTCGCGTCGTTTGCCGCAAATGCCGGGACATTTTTaa                            <  2:450265/148‑1 (MQ=255)
                                                                                                                      cTCGCTGTCACCTTTACCCCAGTAATGGAGGCACGTGCGCAGCATGCGGTTCAGCCACGGTTGAGCATGGGAAATACTACGGTAACTGCTGATAATAACGTGGAGAAAAATGTCGCGTCGTTTGCCGCAAATGCCGGGACATTTTTaa                            >  2:2888438/1‑148 (MQ=255)
                                                                                                                        cGCTGTCACCTTTACCCCAGTAATGGAGGCACGTGCGCAGCATGCGGTTCAGCCACGGTTGAGCATGGGAAATACTACGGTAACTGCTGATAATAACGTGGAGAAAAATGTCGCGTCGTTTGCCGCAAATGCCGGGACATTTTTAAGc                          >  2:1943031/1‑148 (MQ=255)
                                                                                                                        cGCTGTCACCTTTACCCCAGTAATGGAGGCACGTGCGCAGCATGCGGTTCAGCCACGGTTGAGCATGGGAAATACTACGGTAACTGCTGATAATAACGTGGAGAAAAATGTCGCGTCGTTTGCCGCAAATGCCGGGACATTTTTAAGc                          >  2:366481/1‑148 (MQ=255)
                                                                                                                         gCTGTCACCTTTACCCCAGTAATGGAGGCACGTGCGCAGCATGCGGTTCAGCCACGGTTGAGCATGGGAAATACTACGGTAACTGCTGATAATAACGTGGAGAAAAATGTCGCGTCGTTTGCCGCAAATGCCGGGACATTTTTAAGCa                         <  1:1943031/148‑1 (MQ=255)
                                                                                                                           tGTCACCTTTACCCCAGTAATGGCGGCACGTGCGCAGCATGCGGTTCAGCCACGGTTGAGCATGGGAAATACTACGGTAACTGCTGATAATAACGTGGAGAAAAATGTCGCGTCGTTTGCCGCAAATGCCGGGACATTTTTAAGCAGt                       <  1:1858633/148‑1 (MQ=255)
                                                                                                                           tGTCACCTTTACCCCAGTAATGGAGGCACGTGCGCAGCATGCGGTTCAGCCACGGTTGAGCATGGGAAATACTACGGTAACTGCTGATAATAACGTGGAGAAAAATGTCGCGTCGTTTGCCGCAAATGCCGGGACATTTTTAAGCAGt                       >  2:3310322/1‑148 (MQ=255)
                                                                                                                            gTCACCTTTACCCCAGTAATGGCGGCACGTGCGCAGCATGCGGTTCAGCCACGGTTGAGCATGGGAAATACTACGGTAACTGCTGATAATAACGTGGAGAAAAATGTCGCGTCGTTTGCCGCAAATGCCGGGACATTTTTAAGCAGTc                      >  1:2108958/1‑148 (MQ=255)
                                                                                                                            gTCACCTTTACCCCAGTAATGGAGGCACGTGCGCAGCATGCGGTTCAGCCACGGTTGAGCATGGGAAATACTACGGTAACTGCTGATAATAACGTGGAGAAAAATGTCGCGTCGTTTGCCGCAAATGCCGGGACATTTTTAAGCAGTc                      >  1:2386820/1‑148 (MQ=255)
                                                                                                                            gTCACCTTTACCCCAGTAATGGAGGCACGTGCGCAGCATGCGGTTCAGCCACGGTTGAGCATGGGAAATACTACGGTAACTGCTGATAATAACGTGGAGAAAAATGTCGCGTCGTTTGCCGCAAATGCCGGGACATTTTTAAGCAGTc                      >  1:3078620/1‑148 (MQ=255)
                                                                                                                            gTCACCTTTACCCCAGTAATGGAGGCACGTGCGCAGCATGCGGTTCAGCCACGGTTGAGCATGGGAAATACTACGGTAACTGCTGATAATAACGTGGAGAAAAATGTCGCGTCGTTTGCCGCAAATGCCGGGACATTTTTAAGCAGTc                      >  1:1364677/1‑148 (MQ=255)
                                                                                                                            gTCACCTTTACCCCAGTAATGGAGGCACGTGCGCAGCATGCGGTTCAGCCACGGTTGAGCATGGGAAATACTACGGTAACTGCTGATAATAACGTGGAGAAAAATGTCGCGTCGTTTGCCGCAAATGCCGGGACATTTTTAAGCAGTc                      >  1:2272926/1‑148 (MQ=255)
                                                                                                                                 cTTTACCCCAGTAATGGCGGCACGTGCGCAGCATGCGGTTCAGCCACGGTTGAGCATGGGAAATACTACGGTAACTGCTGATAATAACGTGGAGAAAAATGTCGCGTCGTTTGCCGCAAATGCCGGGACATTTTTAAGCAGTCAGCCa                 >  2:3505185/1‑148 (MQ=255)
                                                                                                                                 cTTTACCCCAGTAATGGAGGCACGTGCGCAGCATGCGGTTCAGCCACGGTTGAGCATGGGAAATACTACGGTAACTGCTGATAATAACGTGGAGAAAAATGTCGCGTCGTTTGCCGCAAATGCCGGGACATTTTTAAGCAGTCAGCCa                 >  2:253467/1‑148 (MQ=255)
                                                                                                                                  tttACCCCAGTAATGGCGGCACGTGCGCAGCATGCGGTTCAGCCACGGTTGAGCATGGGAAATACTACGGTAACTGCTGATAATAACGTGGAGAAAAATGTCGCGTCGTTTGCCGCAAATGCCGGGACATTTTTAAGCAGTCAGCCAg                <  2:3111169/148‑1 (MQ=255)
                                                                                                                                  tttACCCCAGTAATGGCGGCACGTGCGCAGCATGCGGTTCAGCCACGGTTGAGCATGGGAAATACTACGGTAACTGCTGATAATAACGTGGAGAAAAATGTCGCGTCGTTTGCCGCAAATGCCGGGACATTTTTAAGCAGTCAGCCAg                <  2:1194220/148‑1 (MQ=255)
                                                                                                                                     aCCCCAGTAATGGCGGCACGTGCGCAGCATGCGGTTCAGCCACGGTTGAGCATGGGAAATACTACGGTAACTGCTGATAATAACGTGGAGAAAAATGTCGCGTCGTTTGCCGCAAATGCCGGGACATTTTTAAGCAGTCAGCCCGATa             >  1:2601739/1‑148 (MQ=255)
                                                                                                                                     aCCCCAGTAATGGCGGCACGTGCGCAGCATGCGGTTCAGCCACGGTTGAGCATGGGAAATACTACGGTAACTGCTGATAATAACGTGGAGAAAAATGTCGCGTCGTTTGCCGCAAATGCCGGGACATTTTTAAGCAGTCAGCCAGATa             >  1:3495962/1‑148 (MQ=255)
                                                                                                                                      ccccAGTAATGGCGGCACGTGCGCAGCATGCGGTTCAGCCACGGTTGAGCATGGGAAATACTACGGTAACTGCTGATAATAACGTGGAGAAAAATGTCGCGTCGTTTGCCGCAAATGCCGGGACATTTTTAAGCAGTCAGCCAGATAg            <  1:1681980/148‑1 (MQ=255)
                                                                                                                                        ccAGTAATGGCGGCACGTGCGCAGCATGCGGTTCAGCCACGGTTGAGCATGGGAAATACTACGGTAACTGCTGATAATAACGTGGAGAAAAATGTCGCGTCGTTTGCCGCAAATGCCGGGACATTTTTAAGCAGTCAGCCAGATAGCg          <  1:298304/148‑1 (MQ=255)
                                                                                                                                         cAGTAATGGCGGCACGTGCGCAGCATGCGGTTCAGCCACGGTTGAGCATGGGAAATACTACGGTAACTGCTGATAATAACGTGGAGAAAAATGTCGCGTCGTTTGCCGCAAATGCCGGGACATTTTTAAGCAGTCAGCCAGATAGCGa         <  1:1759531/148‑1 (MQ=255)
                                                                                                                                         cAGTAATGGCGGCACGTGCGCAGCATGCGGTTCAGCCACGGTTGAGCATGGGAAATACTACGGTAACTGCTGATAATAACGTGGAGAAAAATGTCGCGTCGTTTGCCGCAAATGCCGGGACATTTTTAAGCAGTCAGCCAGATAGCGa         >  1:2647156/1‑148 (MQ=255)
                                                                                                                                         cAGTAATGGCGGCACGTGCGCAGCATGCGGTTCAGCCACGGTTGAGCATGGGAAATACTACGGTAACTGCTGATAATAACGTGGAGAAAAATGTCGCGTCGTTTGCCGCAAATGCCGGGACATTTTTAAGCAGTCAGCCAGATAGCGa         <  1:3298705/148‑1 (MQ=255)
                                                                                                                                          aGTAATGGCGGCACGTGCGCAGCATGCGGTTCAGCCACGGTTGAGCATGGGAAATACTACGGTAACTGCTGATAATAACGTGGAGAAAAATGTCGCGTCGTTTGCCGCAAATGCCGGGACATTTTTAAGCAGTCAGCCAGATAGCGAt        <  1:2056046/148‑1 (MQ=255)
                                                                                                                                           gTAATGGAGGCACGTGCGCAGCATGCGGTTCAGCCACGGTTGAGCATGGGAAATACTACGGTAACTGCTGATAATAACGTGGAGAAAAATGTCGCGTCGTTTGCCGCAAATGCCGGGACATTTTTAAGCAGTCAGCCAGATAGCGATg       <  1:3530824/148‑1 (MQ=255)
                                                                                                                                            tAATGGCGGCACGTGCGCAGCATGCGGTTCAGCCACGGTTGAGCATGGGAAATACTACGGTAACTGCTGATAATAACGTGGAGAAAAATGTCGCGTCGTTTGCCGCAAATGCCGGGACATTTTTAAGCAGTCAGCCAGATAGCGATGc      >  2:1342802/1‑148 (MQ=255)
                                                                                                                                             aaTGGCGGCACGTGCGCAGCATGCGGTTCAGCCACGGTTGAGCATGGGAAATACTACGGTAACTGCTGATAATAACGTGGAGAAAAATGTCGCGTCGTTTGCCGCAAATGCCGGGACATTTTTAAGCAGTCAGCCAGATAGCGATGCg     <  1:471356/148‑1 (MQ=255)
                                                                                                                                             aaTGGAGGCACGTGCGCAGCATGCGGTTCAGCCACGGTTGAGCATGGGAAATACTACGGTAACTGCTGATAATAACGTGGAGAAAAATGTCGCGTCGTTTGCCGCAAATGCCGGGACATTTTTAAGCAGTCAGCCAGATAGCGATGCg     >  2:3369486/1‑148 (MQ=255)
                                                                                                                                              aTGGAGGCACGTGCGCAGCATGCGGTTCAGCCACGGTTGAGCATGGGAAATACTACGGTAACTGCTGATAATAACGTGGAGAAAAATGTCGCGTCGTTTGCCGCAAATGCCGGGACATTTTTAAGCAGTCAGCCAGATAGCGATGCGa    <  2:3073043/148‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
TCAGGAAATTATCTCAATGTCACATTATAAAACAGGTCATAAACAACCACGATTTCGTTATTCAGTTCTGGCCCGCTGCGTGGCGTGGGCAAATATCTCTGTTCAGGTTCTTTTTCCACTCGCTGTCACCTTTACCCCAGTAATGGCGGCACGTGCGCAGCATGCGGTTCAGCCACGGTTGAGCATGGGAAATACTACGGTAACTGCTGATAATAACGTGGAGAAAAATGTCGCGTCGTTTGCCGCAAATGCCGGGACATTTTTAAGCAGTCAGCCAGATAGCGATGCGACA  >  CP009273/310052‑310343

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: