Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,461,872 G→T 45.2% pseudogene (2224/2513 nt) ydbA → pseudogene, autotransporter homolog; interrupted by IS2 and IS30

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,461,8720GT45.2% 177.9 / 274.1 228pseudogene (2224/2513 nt)ydbApseudogene, autotransporter homolog; interrupted by IS2 and IS30
Reads supporting (aligned to +/- strand):  ref base G (51/74);  new base T (49/54);  total (100/128)
Fisher's exact test for biased strand distribution p-value = 3.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.99e-01

Only 100 of 229 total aligned reads displayed.

TCCGTGGGTATCAACGTTTACGGCAGTGATAATAACGTGACGCTGGATGGAAAGTTAACTGTTGTATCAGACAGTGAGGTTACTTCTCGTCAGAGTAATTTATTTGATGGCAGCGCAGAGAAAACGTCAGGTCTGGTTGTGATTGGCGATGGCAATACCGTTAATATGAATGGTGGACTTGAACTGATTGGAGAGAAAAACGCGCTTGCAGATGGGTCGCAGGTTACTTCCTTGCGCACAGGATATAGTTATACCAGCGTTATTGTCGTTAGTGGTGAGTCGTCGGTATAT  >  CP009273/1461725‑1462015
                                                                                                                                                   |                                                                                                                                               
 ccGTGGGTATCAACGTTTACGGCAGTGATAATAACGTGACGCTGGATGGAAAGTTAACTGTTGTATCAGACAGTGAGGTTACTTCTCGTCAGAGTAATTTATTTGATGGCAGCGCAGAGAAAACGTCAGGTCTGGTTGTGATTGGCGa                                                                                                                                                >  1:3427549/1‑148 (MQ=255)
  cGTGGGTATCAACGTTTACGGCAGTGATAATAACGTGACGCTGGATGGAAAGTTAACTGTTGTATCAGACAGTGAGGTTACTTCTCGTCAGAGTAATTTATTTGATGGCAGCGCAGAGAAAACGTCAGGTCTGGTTGTGATTGGCGAt                                                                                                                                               <  2:1154983/148‑1 (MQ=255)
      ggTATCAACGTTTACGGCAGTGATAATAACGTGACGCTGGATGGAAAGTTAACTGTTGTATCAGACAGTGAGGTTACTTCTCGTCAGAGTAATTTATTTGATGGCAGCGCAGAGAAAACGTCAGGTCTGGTTGTGATTGGCTATGGCa                                                                                                                                           <  2:199310/148‑1 (MQ=255)
      ggTATCAACGTTTACGGCAGTGATAATAACGTGACGCTGGATGGAAAGTTAACTGTTGTATCAGACAGTGAGGTTACTTCTCGTCAGAGTAATTTATTTGATGGCAGCGCAGAGAAAACGTCAGGTCTGGTTGTGATTGGCTATGGCa                                                                                                                                           <  2:2544515/148‑1 (MQ=255)
      ggTATCAACGTTTACGGCAGTGATAATAACGTGACGCTGGATGGAAAGTTAACTGTTGTATCAGACAGTGAGGTTACTTCTCGTCAGAGTAATTTATTTGATGGCAGCGCAGAGAAAACGTCAGGTCTGGTTGTGATTGGCGATGGCa                                                                                                                                           <  2:1020688/148‑1 (MQ=255)
        tATCAACGTTTACGGCAGTGATAATAACGTGACGCTGGATGGAAAGTTAACTGTTGTATCAGACAGTGAGGTTACTTCTCGTCAGAGTAATTTATTTGATGGCAGCGCAGAGAAAACGTCAGGTCTGGTTGTGATTGGCGATGGCAAt                                                                                                                                         <  1:2214084/148‑1 (MQ=255)
        tATCAACGTTTACGGCAGTGATAATAACGTGACGCTGGATGGAAAGTTAACTGTTGTATCAGACAGTGAGGTTACTTCTCGTCAGAGTAATTTATTTGATGGCAGCGCAGAGAAAACGTCAGGTCTGGTTGTGATTGGCGATGGCAAt                                                                                                                                         <  1:1899799/148‑1 (MQ=255)
          tCAACGTTTACGGCAGTGATAATAACGTGACGCTGGATGGAAAGTTAACTGTTGTATCAGACAGTGAGGTTACTTCTCGTCAGAGTAATTTATTTGATGGCAGCGCAGAGAAAACGTCAGGTCTGGTTGTGATTGGCTATGGCAATAc                                                                                                                                       >  2:3407299/1‑148 (MQ=255)
                tttACGGCAGTGATAATAACGTGACGCTGGATGGAAAGTTAACTGTTGTATCAGACAGTGAGGTTACTTCTCGTCAGAGTAATTTATTTGATGGCAGCGCAGAGAAAACGTCAGGTCTGGTTGTGATTGGCGATGGCAATACCGTTaa                                                                                                                                 >  1:249349/1‑148 (MQ=255)
                 ttACGGCAGTGATAATAACGTGACGCTGGATGGAAAGTTAACTGTTGTATCAGACAGTGAGGTTACTTCTCGTCAGAGTAATTTATTTGATGGCAGCGCAGAGAAAACGTCAGGTCTGGTTGTGATTGGCTATGGCAATACCGTTAat                                                                                                                                >  1:380304/1‑148 (MQ=255)
                  tACGGCAGTGATAATAACGTGACGCTGGATGGAAAGTTAACTGTTGTATCAGACAGTGAGGTTACTTCTCGTCAGAGTAATTTATTTGATGGCAGCGCAGAGAAAACGTCAGGTCTGGTTGTGATTGGCGATGGCAATACCGTTAata                                                                                                                               <  2:2065608/148‑1 (MQ=255)
                   aCGGCAGTGATAATAACGTGACGCTGGATGGAAAGTTAACTGTTGTATCAGACAGTGAGGTTACTTCTCGTCAGAGTAATTTATTTGATGGCAGCGCAGAGAAAACGTCAGGTCTGGTTGTGATTGGCTATGGCAATACCGTTAatat                                                                                                                              >  2:164813/1‑148 (MQ=255)
                      gCAGTGATAATAACGTGACGCTGGATGGAAAGTTAACTGTTGTATCAGACAGTGAGGTTACTTCTCGTCAGAGTAATTTATTTGATGGCAGCGCAGAGAAAACGTCAGGTCTGGTTGTGATTGGCGATGGCAATACCGTTAATATGaa                                                                                                                           <  1:1787349/148‑1 (MQ=255)
                      gCAGTGATAATAACGTGACGCTGGATGGAAAGTTAACTGTTGTATCAGACAGTGAGGTTACTTCTCGTCAGAGTAATTTATTTGATGGCAGCGCAGAGAAAACGTCAGGTCTGGTTGTGATTGGCGATGGCAATACCGTTAATATGaa                                                                                                                           <  1:2952485/148‑1 (MQ=255)
                       cAGTGATAATAACGTGACGCTGGATGGAAAGTTAACTGTTGTATCAGACAGTGAGGTTACTTCTCGTCAGAGTAATTTATTTGATGGCAGCGCAGAGAAAACGTCAGGTCTGGTTGTGATTGGCTATGGCAATACCGTTAATATGAAt                                                                                                                          >  2:1427571/1‑148 (MQ=255)
                       cAGTGATAATAACGTGACGCTGGATGGAAAGTTAACTGTTGTATCAGACAGTGAGGTTACTTCTCGTCAGAGTAATTTATTTGATGGCAGCGCAGAGAAAACGTCAGGTCTGGTTGTGATTGGCGATGGCAATACCGTTAATATGAAt                                                                                                                          <  2:631367/148‑1 (MQ=255)
                          tGATAATAACGTGACGCTGGATGGAAAGTTAACTGTTGTATCAGACAGTGAGGTTACTTCTCGTCAGAGTAATTTATTTGATGGCAGCGCAGAGAAAACGTCAGGTCTGGTTGTGATTGGCTATGGCAATACCGTTAATATGAAtggt                                                                                                                       <  1:576604/148‑1 (MQ=255)
                          tGATAATAACGTGACGCTGGATGGAAAGTTAACTGTTGTATCAGACAGTGAGGTTACTTCTCGTCAGAGTAATTTATTTGATGGCAGCGCAGAGAAAACGTCAGGTCTGGTTGTGATTGGCGATGGCAATACCGTTAATATGAAtggt                                                                                                                       <  2:1429750/148‑1 (MQ=255)
                           gATAATAACGTGACGCTGGATGGAAAGTTAACTGTTGTATCAGACAGTGAGGTTACTTCTCGTCAGAGTAATTTATTTGATGGCAGCGCAGAGAAAACGTCAGGTCTGGTTGTGATTGGCGATGGCAATACCGTTAATATGAAtggtg                                                                                                                      >  2:3479248/1‑148 (MQ=255)
                             taataaCGTGACGCTGGATGGAAAGTTAACTGTTGTATCAGACAGTGAGGTTACTTCTCGTCAGAGTAATTTATTTGATGGCAGCGCAGAGAAAACGTCAGGTCTGGTTGTGATTGGCGATGGCAATACCGTTAATATGAATGGTGGa                                                                                                                    <  2:2814635/148‑1 (MQ=255)
                             taataaCGTGACGCTGGATAGAAAGTTAACTGTTGTATCAGACAGTGAGGTTACTTCTCGTCAGAGTAATTTATTTGATGGCAGCGCAGAGAAAACGTCAGGTCTGGTTGTGATTGGCTATGGCAATACCGTTAATATGAATGGTGGa                                                                                                                    >  1:2937302/1‑148 (MQ=255)
                              aataaCGTGACGCTGGATGGAAAGTTAACTGTTGTATCAGACAGTGAGGTTACTTCTCGTCAGATTAATTTATTTGATGGCAGCGCAGAGAAAACGTCAGGTCTGGTTGTGATTGGCGATGGCAATACCGTTAATATGAATGGTGGAc                                                                                                                   <  2:228261/148‑1 (MQ=255)
                                taaCGTGACGCTGGATGGAAAGTTAACTGTTGTATCAGACAGTGAGGTTACTTCTCGTCAGAGTAATTTATTTGATGGCAGCGCAGAGAAAACGTCAGGTCTGGTTGTGATTGGCTATGGCAATACCGTTAATATGAATGGTGGACtt                                                                                                                 <  2:732396/148‑1 (MQ=255)
                                  aCGTGACGCTGGATGGAAAGTTAACTGTTGTATCAGACAGTGAGGTTACTTCTCGTCAGAGTAATTTATTTGATGGCAGCGCAGAGAAAACGTCAGGTCTGGTTGTGATTGGCTATGGCAATACCGTTAATATGAATGGTGGACTTGa                                                                                                               >  1:269573/1‑148 (MQ=255)
                                   cGTGACGCTGGATGGAAAGTTAACTGTTGTATCAGACAGTGAGGTTACTTCTCGTCAGAGTAATTTATTTGATGGCAGCGCAGAGAAAACGTCAGGTCTGGTTGTGATTGGCTATGGCAATACCGTTAATATGAATGGTGGACTTGaa                                                                                                              <  1:2570124/148‑1 (MQ=255)
                                     tGACGCTGGATGGAAAGTTAACTGTTGTATCAGACAGTGAGGTTACTTCTCGTCAGAGTAATTTATTTGATGGCAGCGCAGAGAAAACGTCAGGTCTGGTTGTGATTGGCTATGGCAATACCGTTAATATGAATGGTGGACTTGAACt                                                                                                            >  2:2404643/1‑148 (MQ=255)
                                     tGACGCTGGATGGAAAGTTAACTGTTGTATCAGACAGTGAGGTTACTTCTCGTCAGAGTAATTTATTTGATGGCAGCGCAGAGAAAACGTCAGGTCTGGTTGTGATTGGCGATGGCAATACCGTTAATATGAATGGTGGACTTGAACt                                                                                                            <  2:3455693/148‑1 (MQ=255)
                                      gACGCTGGATGGAAAGTTAACTGTTGTATCAGACAGTGAGGTTACTTCTCGTCAGAGTAATTTATTTGATGGCAGCGCAGAGAAAACGTCAGGTCTGGTTGTGATTGGCGATGGCAATACCGTTAATATGAATGGTGGACTTGAACTg                                                                                                           <  1:434022/148‑1 (MQ=255)
                                       aCGCTGGATGGAAAGTTAACTGTTGTATCAGACAGTGAGGTTACTTCTCGTCAGAGTAATTTATTTGATGGCAGCGCAGAGAAAACGTCAGGTCTGGTTGTGATTGGCTATGGCAATACCGTTAATATGAATGGTGGACTTGAACTGa                                                                                                          <  1:1930323/148‑1 (MQ=255)
                                        cGCTGGATGGAAAGTTAACTGTTGTATCAGACAGTGAGGTTACTTCTCGTCAGAGTAATTTATTTGATGGCAGCGCAGAGAAAACGTCAGGTCTGGTTGTGATTGGCTATGGCAATACCGTTAATATGAATGGTGGACTTGAACTGAt                                                                                                         <  2:3329752/148‑1 (MQ=255)
                                           tggatggaAAGTTAACTGTTGTATCAGACAGTGAGGTTACTTCTCGTCAGAGTAATTTATTTGATGGCAGCGCAGAGAAAACGTCAGGTCTGGTTGTGATTGGCGATGGCAATACCGTTAATATGAATGGTGGACTTGAACTGATTgg                                                                                                      <  2:1308357/148‑1 (MQ=255)
                                             gatggaAAGTTAACTGTTGTATCAGACAGTGAGGTTACTTCTCGTCAGAGTAATTTATTTGATGGCAGCGCAGAGAAAACGTCAGGTCTGGTTGTGATTGGCTATGGCAATACCGTTAATATGAATGGTGGACTTGAACTGATTGgag                                                                                                    >  2:2452085/1‑148 (MQ=255)
                                             gatggaAAGTTAACTGTTGTATCAGACAGTGAGGTTACTTCTCGTCAGAGTAATTTATTTGATGGCAGCGCAGAGAAAACGTCAGGTCTGGTTGTGATTGGCTATGGCAATACCGTTAATATGAATGGTGGACTTGAACTGATTGgag                                                                                                    >  2:2563129/1‑148 (MQ=255)
                                              atggaAAGTTAACTGTTGTATCAGACAGTGAGGTTACTTCTCGTCAGAGTAATTTATTTGATGGCAGCGCAGAGAAAACGTCAGGTCTGGTTGTGATTGGCTATGGCAATACCGTTAATATGAATGGTGGACTTGAACTGATTGgaga                                                                                                   <  2:1892690/148‑1 (MQ=255)
                                              atggaAAGTTAACTGTTGTATCAGACAGTGAGGTTACTTCTCGTCAGAGTAATTTATTTGATGGCAGCGCAGAGAAAACGTCAGGTCTGGTTGTGATTGGCTATGGCAATACCGTTAATATGAATGGTGGACTTGAACTGATTGgaga                                                                                                   >  2:1909556/1‑148 (MQ=255)
                                                ggaAAGTTAACTGTTGTATCAGACAGTGAGGTTACTTCTCGTCAGAGTAATTTATTTGATGGCAGCGCAGAGAAAACGTCAGGTCTGGTTGTGATTGGCTATGGCAATACCGTTAATATGAATGGTGGACTTGAACTGATTGgagaga                                                                                                 >  1:1866009/1‑148 (MQ=255)
                                                ggaAAGTTAACTGTTGTATCAGACAGTGAGGTTACTTCTCGTCAGAGTAATTTATTTGATGGCAGCGCAGAGAAAACGTCAGGTCTGGTTGTGATTGGCTATGGCAATACCGTTAATATGAATGGTGGACTTGAACTGATTGgagaga                                                                                                 >  2:1801171/1‑148 (MQ=255)
                                                  aaaGTTAACTGTTGTATCAGACAGTGAGGTTACTTCTCGTCAGAGTAATTTATTTGATGGCAGCGCAGAGAAAACGTCAGGTCTGGTTGTGATTGGCTATGGCAATACCGTTAATATGAATGGTGGACTTGAACTGATTGGAGAGaaa                                                                                               <  1:3202577/148‑1 (MQ=255)
                                                  aaaGTTAACTGTTGTATCAGACAGTGAGGTTACTTCTCGTCAGAGTAATTTATTTGATGGCAGCGCAGAGAAAACGTCAGGTCTGGTTGTGATTGGCGATGGCAATACCGTTAATATGAATGGTGGACTTGAACTGATTGGAGAGaaa                                                                                               <  1:11492/148‑1 (MQ=255)
                                                    aGTTAACTGTTGTATCAGACAGTGAGGTTACTTCTCGTCAGAGTAATTTATTTGATGGCAGCGCAGAGAAAACGTCAGGTCTGGTTGTGATTGGCTATGGCAATACCGTTAATATGAATGGTGGACTTGAACTGATTGGAGAGaaaaa                                                                                             <  2:2713201/148‑1 (MQ=255)
                                                    aGTTAACTGTTGTATCAGACAGTGAGGTTACTTCTCGTCAGAGTAATTTATTTGATGGCAGCGCAGAGAAAACGTCAGGTCTGGTTGTGATTGGCGATGGCAATACCGTTAATATGAATGGTGGACTTGAACTGATTGGAGAGaaaaa                                                                                             <  1:2821655/148‑1 (MQ=255)
                                                      ttAACTGTTGTATCAGACAGTGAGGTTACTTCTCGTCAGAGTAATTTATTTGATGGCAGCGCAGAGAAAACGTCAGGTCTGGTTGTGATTGGCGATGGCAATACCGTTAATATGAATGGTGGACTTGAACTGATTGGAGAGAAAAAcg                                                                                           >  2:793597/1‑148 (MQ=255)
                                                           tgttgtATCAGACAGTGAGGTTACTTCTCGTCAGAGTAATTTATTTGATGGCAGCGCAGAGAAAACGTCAGGTCTGGTTGTGATTGGCGATGGCAATACCGTTAATATGAATGGTGGACTTGAACTGATTGGAGAGAAAAACGCGCtt                                                                                      >  1:1716220/1‑148 (MQ=255)
                                                             ttgtATCAGACAGTGAGGTTACTTCTCGTCAGAGTAATTTATTTGATGGCAGCGCAGAGAAAACGTCAGGTCTGGTTGTGATTGGCTATGGCAATACCGTTAATATGAATGGTGGACTTGAACTGATTGGAGAGAAAAACGCGCTTGc                                                                                    <  2:2701121/148‑1 (MQ=255)
                                                              tgtATCAGACAGTGAGGTTACTTCTCGTCAGAGTAATTTATTTGATGGCAGCGCAGAGAAAACGTCAGGTCTGGTTGTGATTGGCTATGGCAATACCGTTAATATGAATGGTGGACTTGAACTGATTGGAGAGAAAAACGCGCTTGCa                                                                                   >  1:2324224/1‑148 (MQ=255)
                                                              tgtATCAGACAGTGAGGTTACTTCTCGTCAGAGTAATTTATTTGATGGCAGCGCAGAGAAAACGTCAGGTCTGGTTGTGATTGGCTATGGCAATACCGTTAATATGAATGGTGGACTTGAACTGATTGGAGAGAAAAACGCGCTTGCa                                                                                   >  1:2715063/1‑148 (MQ=255)
                                                               gtATCAGACAGTGAGGTTACTTCTCGTCAGAGTAATTTATTTGATGGCAGCGCAGAGAAAACGTCAGGTCTGGTTGTGATTGGCTATGGCAATACCGTTAATATGAATGGTGGACTTGAACTGATTGGAGAGAAAAACGCGCTTGCAg                                                                                  <  2:214817/148‑1 (MQ=255)
                                                               gtATCAGACAGTGAGGTTACTTCTCGTCAGAGTAATTTATTTGATGGCAGCGCAGAGAAAACGTCAGGTCTGGTTGTGATTGGCGATGGCAATACCGTTAATATGAATGGTGGACTTGAACTGATTGGAGAGAAAAACGCGCTTGCAg                                                                                  <  2:312786/148‑1 (MQ=255)
                                                                   cAGACAGTGAGGTTACTTCTCGTCAGAGTAATTTATTTGATGGCAGCGCAGAGAAAACGTCAGGTCTGGTTGTGATTGGCGATGGCAATACCGTTAATATGAATGGTGGACTTGAACTGATTGGAGAGAAAAACGCGCTTGCAGATgg                                                                              <  2:1323535/148‑1 (MQ=255)
                                                                     gACAGTGAGGTTACTTCTCGTCAGAGTAATTTATTTGATGGCAGCGCAGAGAAAACGTCAGGTCTGGTTGTGATTGGCTATGGCAATACCGTTAATATGAATGGTGGACTTGAACTGATTGGAGAGAAAAACGCGCTTGCAGATGGGt                                                                            <  1:490519/148‑1 (MQ=255)
                                                                     gACAGTGAGGTTACTTCTCGTCAGAGTAATTTATTTGATGGCAGCGCAGAGAAAACGTCAGGTCTGGTTGTGATTGGCTATGGCAATACCGTTAATATGAATGGTGGACTTGAACTGATTGGAGAGAAAAACGCGCTTGCAGATGGGt                                                                            >  2:45988/1‑148 (MQ=255)
                                                                        aGTGAGGTTACTTCTCGTCAGAGTAATTTATTTGATGGCAGCGCAGAGAAAACGTCAGGTCTGGTTGTGATTGGCGATGGCAATACCGTTAATATGAATGGTGGACTTGAACTGATTGGAGAGAAAAACGCGCTTGCAGATGGGTCGc                                                                         <  2:1484371/148‑1 (MQ=255)
                                                                         gTGAGGTTACTTCTCGTCAGAGTAATTTATTTGATGGCAGCGCAGAGAAAACGTCAGGTCTGGTTGTGATTGGCGATGGCAATACCGTTAATATGAATGGTGGACTTGAACTGATTGGAGAGAAAAACGCGCTTGCAGATGGGTCGCa                                                                        >  2:1875344/1‑148 (MQ=255)
                                                                           gtggTTAATTCTCGTCAGATTAATTTATTTGATGGCAGCGCAGAGAAAACGTCAGGTCTGGTTGTGATTGGCGATGGCAATACCGTTAATATGATTGGTGGACTTGAACTGATTGGAGAGAAAAACGCGCTTGCAGATGGGTCGCAgg                                                                      <  1:2929488/146‑1 (MQ=255)
                                                                           gAGGTTACTTCTCGTCAGAGTAATTTATTTGATGGCAGCGCAGAGAAAACGTCAGGTCTGGTTGTGATTGGCTATGGCAATACCGTTAATATGAATGGTGGACTTGAACTGATTGGAGAGAAAAACGCGCTTGCAGATGGGTCGCAgg                                                                      <  2:2006173/148‑1 (MQ=255)
                                                                           gAGGTTACTTCTCGTCAGAGTAATTTATTTGATGGCAGCGCAGAGAAAACGTCAGGTCTGGTTGTGATTGGCGATGGCAATACCGTTAATATGAATGGTGGACTTGAACTGATTGGAGAGAAAAACGCGCTTGCAGATGGGTCGCAgg                                                                      <  2:1272601/148‑1 (MQ=255)
                                                                                 aCTTCTCGTCAGAGTAATTTATTTGATGGCAGCGCAGAGAAAACGTCAGGTCTGGTTGTGATTGGCTATGGCAATACCGTTAATATGAATGGTGGACTTGAACTGATTGGAGAGAAAAACGCGCTTGCAGATGGGTCGCAGGTTACtt                                                                <  1:1600405/148‑1 (MQ=255)
                                                                                 aCTTCTCGTAAGAGTAATTTATTTGAGGGCAGCGCATAGAAAACGTCAGGTCTGGTTGTGATTGGCGATGGCAAGACCGTTAATATGACTGGTGGACTTGAACTGATTGGAGAGAAAAACGCGCTTGCAGATGGGTCGCAGGTTACtt                                                                <  1:2799558/148‑1 (MQ=255)
                                                                                   ttCTCGTCAGAGTAATTTATTTGATGGCAGCGCAGAGAAAACGTCAGGTCTGGTTGTGATTGGCTATGGCAATACCGTTAATATGAATGGTGGACTTGAACTGATTGGAGAGAAAAACGCGCGTGCAGATGGGTCGCAGGTTACTTcc                                                              <  1:1250304/148‑1 (MQ=255)
                                                                                   ttCTCGTCAGAGTAATTTATTTGATGGCAGCGCAGAGAAAACGTCAGGTCTGGTTGTGATTGGCGATGGCAATACCGTTAATATGAATGGTGGACTTGAACTGATTGGAGAGAAAAACGCGCTTGCAGATGGGTCGCAGGTTACTTcc                                                              <  2:1777372/148‑1 (MQ=255)
                                                                                   ttCTCGTCAGAGTAATTTATTTGATGGCAGCGCAGAGAAAACGTCAGGTCTGGTTGTGATTGGCGATGGCAATACCGTTAATATGAATGGTGGACTTGAACTGATTGGAGAGAAAAACGCGCTTGCAGATGGGTCGCAGGTTACTTcc                                                              <  2:3539177/148‑1 (MQ=255)
                                                                                    tctcGTCAGAGTAATTTATTTGATGGCAGCGCAGAGAAAACGTCAGGTCTGGTTGTGATTGGCGATGGCAATACCGTTAATATGAATGGTGGACTTGAACTGATTGGAGAGAAAAACGCGCTTGCAGATGGGTCGCAGGTTACTTCCt                                                             <  2:1265529/148‑1 (MQ=255)
                                                                                    tctcGTCAGAGTAATTTATTTGATGGCAGCGCAGAGAAAACGTCAGGTCTGGTTGTGATTGGCGATGGCAATACCGTTAATATGAATGGTGGACTTGAACTGATTGGAGAGAAAAACGCGCTTGCAGATGGGTCGCAGGTTACTTCCt                                                             <  2:3136711/148‑1 (MQ=255)
                                                                                         tCAGAGTAATTTATTTGATGGCAGCGCAGAGAAAACGTCAGGTTTGGTTGTGATTGGCGATGGCAATACCGTTAATATGAATGGTGGACTTGAACTGATTGGAGAGAAAAACGCGCTTGCAGATGGGTCGCAGGTTACTTCCTTgcgc                                                        >  2:49308/1‑148 (MQ=255)
                                                                                         tCAGAGTAATTTATTTGATGGCAGCGCAGAGAAAACGTCAGGTCTGGTTGTGATTGGCGATGGCAATACCGTTAATATGAATGGTGGACTTGAACTGATTGGAGAGAAAAACGCGCTTGCAGATGGGTCGCAGGTTACTTCCTTgcgc                                                        >  1:199019/1‑148 (MQ=255)
                                                                                           agagTAATTTATTTGATGGCAGCGCAGAGAAAACGTCAGGTCTGGTTGTGATTGGCTATGGCAATACCGTTAATATGAATGGTGGACTTGAACTGATTGGAGAGAAAAACGCGCTTGCAGATGGGTCGCAGGTTACTTCCTTGCGcac                                                      <  2:2756934/148‑1 (MQ=255)
                                                                                             agTAATTTATTTGATGGCAGCGCAGAGAAAACGTCAGGTCTGGTTGTGATTGGCGATGGCAATACCGTTAATATGAATGGTGGACTTGAACTGATTGGAGAGAAAAACGCGCTTGCAGATGGGTCGCAGGTTACTTCCTTGCGCACAg                                                    >  2:1977527/1‑148 (MQ=255)
                                                                                             agTAATTTATTTGATGGCAGCGCAGAGAAAACGTCAGGTCTGGTTGTGATTGGCGATGGCAATACCGTTAATATGAATGGTGGACTTGAACTGATTGGAGAGAAAAACGCGCTTGCAGATGGGTCGCAGGTTACTTCCTTGCGCACAg                                                    <  2:2969312/148‑1 (MQ=255)
                                                                                              gTAATTTATTTGATGGCAGCGCAGAGAAAACGTCAGGTCTGGTTGTGATTGGCGATGGCAATACCGTTAATATGAATGGTGGACTTGAACTGATTGGAGAGAAAAACGCGCTTGCAGATGGGTCGCAGGTTACTTCCTTGCGCACAgg                                                   >  2:436260/1‑148 (MQ=255)
                                                                                               tAATTTATTTGATGGCAGCGCAGAGAAAACGTCAGGTCTGGTTGTGATTGGCGATGGCAATACCGTTAATATGAATGGTGGACTTGAACTGATTGGAGAGAAAAACGCGCTTGCAGATGGGTCGCAGGTTACTTCCTTGCGCACAGGa                                                  >  1:3075675/1‑148 (MQ=255)
                                                                                                aaTTTATTTGATGGCAGCGCAGAGAAAACGTCAGGTCTGGTTGTGATTGGCTATGGCAATACCGTTAATATGAATGGTGGACTTGAACTGATTGGAGAGAAAAACGCGCTTGCAGATGGGTCGCAGGTTACTTCCTTGCGCACAGGat                                                 <  2:269573/148‑1 (MQ=255)
                                                                                                 atttatttGATGGCAGCGCAGAGAAAACGTCAGGTCTGGTTGTGATTGGCTATGGCAATACCGTTAATATGAATGGTGGACTTGAACTGATTGGAGAGAAAAACGCGCTTGCAGATGGGTCGCAGGTTACTTCCTTGCGCACAGGata                                                <  2:3260082/148‑1 (MQ=255)
                                                                                                 atttatttGATGGCAGCGCAGAGAAAACGTCAGGTCTGGTTGTGATTGGCTATGGCAATACCGTTAATATGAATGGTGGACTTGAACTGATTGGAGAGAAAAACGCGCTTGCAGATGGGTCGCAGGTTACTTCCTTGCGCACAGGata                                                >  1:1363361/1‑148 (MQ=255)
                                                                                                    tatttGATGGCAGCGCAGAGAAAACGTCAGGTCTGGTTGTGATTGGCTATGGCAATACCGTTAATATGAATGGTGGACTTGAACTGATTGGAGAGAAAAACGCGCTTGCAGATGGGTCGCAGGTTACTTCCTTGCGCACAGGATATAg                                             >  1:3392136/1‑148 (MQ=255)
                                                                                                    tatttGATGGCAGCGCAGAGAAAACGTCAGGTCTGGTTGTGATTGGCTATGGCAATACCGTTAATATGAATGGTGGACTTGAACTGATTGGAGAGAAAAACGCGCTTGCAGATGGGTCGCAGGTTACTTCCTTGCGCACAGGATATAg                                             >  1:3393063/1‑148 (MQ=255)
                                                                                                        tGATGGCAGCGCAGAGAAAACGTCAGGTCTGGTTGTGATTGGCTATGGCAATACCGTTAATATGAATGGTGGACTTGAACTGATTGGAGAGAAAAACGCGCTTGCAGATGGGTCGCAGGTTACTTCCTTGCGCACAGGATATAGTtat                                         >  1:2629951/1‑148 (MQ=255)
                                                                                                           tGGCAGCTCAGAGAAAAAGTCATGTCTGTTTGTGATTTGCTATTGCAATACCGTTTATATGAATGGTGGACTTGAACTGATTTTAGAGAAAAACGCTCTTGCAGATGGGTCTCAGGTTACTTCCTTGCGCACAGGATATAGTTCTAcc                                      <  2:3485657/148‑1 (MQ=255)
                                                                                                                 cgcAGAGAAAACGTCAGGTCTGGTTGTGATTGGCTATGGCAATACCGTTAATATGAATGGTGGACTTGAACTGATTGGAGAGAAAAACGCGCTTGCAGATGGGTCGCAGGTTACTTCCTTGCGCACAGGATATAGTTATACCAGCGtt                                >  2:1062564/1‑148 (MQ=255)
                                                                                                                    agagAAAACGTCAGGTCTGGTTGTGATTGGCGATGGCAATACCGTTAATATGAATGGTGGACTTGAACTGATTGGAGAGAAAAACGCGCTTGCAGATGGGTCGCAGGTTACTTCCTTGCGCACAGGATATAGTTATACCAGCGTTAtt                             >  1:1089932/1‑148 (MQ=255)
                                                                                                                      agaAAACGTCAGGTCTGGTTGTGATTGGCGATGGCAATACCGTTAATATGAATGGTGGACTTGAACTGATTGGAGAGAAAAACGCGCTTGCAGATGGGTCGCAGGTTACTTCCTTGCGCACAGGATATAGTTATACCAGCGTTATTGt                           >  1:1356753/1‑148 (MQ=255)
                                                                                                                          aaCGTCAGGTCTGGTTGTGATTGGCGATGGCAATACCGTTAATATGAATGGTGGACTTGAACTGATTGGAGAGAAAAACGCGCTTGCAGATGGGTCGCAGGTTACTTCCTTGCGCACAGGATATAGTTATACCAGCGTTATTGTCGtt                       >  1:1653567/1‑148 (MQ=255)
                                                                                                                           aCGTCAGGTCTGGTTGTGATTGGCTATGGCAATACCGTTAATATGAATGGTGGACTTGAACTGATTGGAGAGAAAAACGCGCTTGCAGATGGGTCGCAGGTTACTTCCTTGCGCACAGGATATAGTTATACCAGCGTTATTGTCGTTa                      >  1:711467/1‑148 (MQ=255)
                                                                                                                               cAGGTCTGGTTGTGATTGGCTATGGCAATACCGTTAATATGAATGGTGGACTTGAACTGATTGGAGAGAAAAACGCGCTTGCAGATGGGTCGCAGGTTACTTCCTTGCGCACAGGATATAGTTATACCAGCGTTATTGTCGTTAgtgg                  <  1:1427571/148‑1 (MQ=255)
                                                                                                                               cAGGTCTGGTTGTGATTGGCGATGGCAATACCGTTAATATGAATGGTGGACTTGAACTGATTGGAGAGAAAAACGCGCTTGCAGATGGGTCGCAGGTTTGTTCCTTGCGCACAGGATATAGTTATACCAGCGTTATTGTCGTTAgtgg                  <  1:778085/148‑1 (MQ=255)
                                                                                                                               cAGGTCTGGTTGTGATTGGCGATGGCAATACCGTTAATATGAATGGTGGACTTGAACTGATTGGAGAGAAAAACGCGCTTGCAGATGGGTCGCAGGTTACTTCCTTGCGCACAGGATATAGTTATACCAGCGTTATTGTCGTTAgtgg                  <  2:2373905/148‑1 (MQ=255)
                                                                                                                                 ggTCTGGTTGTGATTGGCTATGGCAATACCGTTAATATGAATGGTGGACTTGAACTGATTGGAGAGAAAAACGCGCTTGCAGATGGGTCGCAGGTTACTTCCTTGCGCACAGGATATAGTTATACCAGCGTTATTGTCGTTAgtggtg                <  2:1473615/148‑1 (MQ=255)
                                                                                                                                 ggTCTGGTTGTGATTGGCGATGGCAATACCGTTAATATGAATGGTGGACTTGAACTGATTGGAGAGAAAAACGCGCTTGCAGATGGGTCGCAGGTTACTTCCTTGCGCACAGGATATAGTTATACCAGCGTTATTGTCGTTAgtggtg                <  1:2662066/148‑1 (MQ=255)
                                                                                                                                  gTCTGGTTGTGATTGGCGATGGCAATACCGTTAATATGAATGGTGGACTTGAACTGATTGGAGAGAAAAACGCGCTTGCAGATGGGTCGCAGGTTACTTCCTTGCGCACAGGATATAGTTATACCAGCGTTATTGTCGTTATTGGTGa               >  2:3000344/1‑148 (MQ=255)
                                                                                                                                   tCTGGTTGTGATTGGCGATGGCAATACCGTTAATATGAATGGTGGACTTGAACTGATTGGAGAGAAAAACGCGCTTGCAGATGGGTCGCAGGTTACTTCCTTGCGCACAGGATATAGTTATACCAGCGTTATTGTCGTTAGTGGTGAg              >  2:2905596/1‑148 (MQ=255)
                                                                                                                                       gTTGTGATTGGCGATGGCAATACCGTTAATATGAATGGTGGACTTGAACTGATTGGAGAGAAAAACGCGCTTGCAGATGGGTCGCAGGTTACTTCCTTGCGCACAGGATATAGTTATACCAGCGTTATTGTCGTTAGTGGTGAgtcgt          <  1:3479248/148‑1 (MQ=255)
                                                                                                                                        ttGTGATTGGCTATGGCAATACCGTTAATATGAATGGTGGACTTGAACTGATTGGAGAGAAAAACGCGCTTGCAGATGGGTCGCAGGTTACTTCCTTGCGCACAGGATATAGTTATACCAGCGTTATTGTCGTTAGTGGTGAgtcgtc         >  1:3308787/1‑148 (MQ=255)
                                                                                                                                           tgATTGGCTATGGCAATACCGTTAATATGAATGGTGGACTTGAACTGATTGGAGAGAAAAACGCGCTTGCAGATGGGTCGCAGGTTACTTCCTTGCGCACAGGATATAGTTATACCAGCGTTATTGTCGTTAGTGGTGAGTCGTCGGt      >  1:924770/1‑148 (MQ=255)
                                                                                                                                           tgATTGGCGATGGCAATACCGTTAATATGAATGGTGGACTTGAACTGATTGGAGAGAAAAACGCGCTTGCAGATGGGTCGCAGGTTACTTCCTTGCGCACAGGATATAGTTATACCAGCGTTATTGTCGTTAGTGGTGAGTCGTCGGt      <  1:1877399/148‑1 (MQ=255)
                                                                                                                                            gATTGGCTATGGCAATACCGTTAATATGAATGGTGGACTTGAACTGATTGGAGAGAAAAACGCGCTTGCAGATGGGTCGCAGGTTACTTCCTTGCGCACAGGATATAGTTATACCAGCGTTATTGTCGTTAGTGGTGAGTCGTCGGta     >  1:1184781/1‑148 (MQ=255)
                                                                                                                                            gATTGGCGATGGCAATACCGTTAATATGAATGGTGGACTTGAACTGATTGGAGAGAAAAACGCGCTTGCAGATGGGTCGCAGGTTACTTCCTTGCGCACAGGATATAGTTATACCAGCGTTATTGTCGTTAGTGGTGAGTCGTCGGta     <  1:3489232/148‑1 (MQ=255)
                                                                                                                                             aTTGGCGATGGCAATACCGTTAATATGAATGGTGGACTTGAACTGATTGGAGAGAAAAACGCGCTTGCAGATGGGTCGCAGGTTACTTCCTTGCGCACAGGATATAGTTATACCAGCGTTATTGTCGTTAGTGGTGAGTCGTCGGtat    <  1:2114883/148‑1 (MQ=255)
                                                                                                                                              ttGGCTATGGCAATACCGTTAATATGAATGGTGGACTTGAACTGATTGGAGAGAAAAACGCGCTTGCAGATGGGTCGCAGGTTACTTCCTTGCGCACAGGATATAGTTATACCAGCGTTATTGTCGTTAGTGGTGAGTCGTCGGtata   >  1:655155/1‑148 (MQ=255)
                                                                                                                                              ttGGCGATGGCAATACCGTTAATATGAATGGTGGACTTGAACTGATTGGAGAGAAAAACGCGCTTGCAGATGGGTCGCAGGTTACTTCCTTGCGCACAGGATATAGTTATACCAGCGTTATTGTCGTTAGTGGTGAGTCGTCGGtata   <  1:291820/148‑1 (MQ=255)
                                                                                                                                              ttGGCGATGGCAATACCGTTAATATGAATGGTGGACTTGAACTGATTGGAGAGAAAAACGCGCTTGCAGATGGGTCGCAGGTTACTTCCTTGCGCACAGGATATAGTTATACCAGCGTTATTGTCGTTAGTGGTGAGTCGTCGGtata   <  1:1977628/148‑1 (MQ=255)
                                                                                                                                              ttGGCGATGGCAATACCGTTAATATGAATGGTGGACTTGAACTGATTGGAGAGAAAAACGCGCTTGCAGATGGGTCGCAGGTTACTTCCTTGCGCACAGGATATAGTTATACCAGCGTTATTGTCGTTAGTGGTGAGTCGTCGGtata   <  1:1977402/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                               
TCCGTGGGTATCAACGTTTACGGCAGTGATAATAACGTGACGCTGGATGGAAAGTTAACTGTTGTATCAGACAGTGAGGTTACTTCTCGTCAGAGTAATTTATTTGATGGCAGCGCAGAGAAAACGTCAGGTCTGGTTGTGATTGGCGATGGCAATACCGTTAATATGAATGGTGGACTTGAACTGATTGGAGAGAAAAACGCGCTTGCAGATGGGTCGCAGGTTACTTCCTTGCGCACAGGATATAGTTATACCAGCGTTATTGTCGTTAGTGGTGAGTCGTCGGTATAT  >  CP009273/1461725‑1462015

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: