Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,466,293 C→A 12.7% pseudogene (1520/3497 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,466,2930CA12.7% 833.9 / 63.0 259pseudogene (1520/3497 nt)ydbApseudogene, autotransporter homolog; interrupted by IS2 and IS30
Reads supporting (aligned to +/- strand):  ref base C (96/129);  new base A (15/18);  total (111/148)
Fisher's exact test for biased strand distribution p-value = 8.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.43e-01

Only 100 of 262 total aligned reads displayed.

TAGTAGTGAGGGTAACAACACCGGGATGCTGGAAGTTAATAATAATTCAGCTTTCAATAACCGCGGCGAGTTTATTCTTGATAACGACAAGAATGCTGTGCACATCAACCAGTCCGGTACGCTTTATAATACCGGTCACATGAACATCAGTAATTCTTCCCACAACGGAGCCGTTAATATGTGGGGCGGAAATGGTCGTTTTATCAATGACGGAACGATTGATGTTTCTGCGAAGTCACTGGTAGTCAGCGCTAATAATGCCGGCGATCAGAATGCCTTCTTCTGGAACCAGGAT  >  CP009273/1466146‑1466440
                                                                                                                                                   |                                                                                                                                                   
tagtagtGAGGGTAACAACACCGGGATGCTGGAAGTTAATAATAATTCAGCTTTCAATAACCGCGGCGAGTTTATTCTTGATAACGACAAGAATGCTGTGCACATCAACCAGTCCGGTACGCTTTATAATACCGGTCACATGAACATc                                                                                                                                                     >  1:2505419/1‑148 (MQ=255)
   tagtGAGGGTAACAACACCGGGATGCTGGAAGTTAATAATAATTCAGCTTTCAATAACCGCGGCGAGTTTATTCTTGATAACGACAAGAATGCTGTGCACATCAACCAGTCCTGTACTCTTTATAATACCTGTCACATGAACATCAGt                                                                                                                                                  >  2:3181851/1‑148 (MQ=255)
     gtGAGGGTAACAACACCGGGATGCTGGAAGTTAATAATAATTCAGCTTTCAATAACCGCGGCGAGTTTATTCTTGATAACGACAAGAATGCTGTGCACATCAACCAGTCCGGTACGCTTTATAATACCGGTCACATGAACATCAGTaa                                                                                                                                                <  2:2339393/148‑1 (MQ=255)
      tGAGGGTAACAACACCGGGATGCTGGAAGTTAATAATAATTCAGCTTTCAATAACCGCGGCGAGTTTATTCTTGATAACGACAAGAATGCTGTGCACATCAACCAGTCCGGTACGCTTTATAATACCGGTCACATGAACATCAGTAAt                                                                                                                                               <  1:3095131/148‑1 (MQ=255)
      tGAGGGTAACAACACCGGGATGCTGGAAGTTAATAATAATTCAGCTTTCAATAACCGCGGCGAGTTTATTCTTGATAACGACAAGAATGCTGTGCACATCAACCAGTCCGGTACGCTTTATAATACCGGTCACATGAACATCAGTAAt                                                                                                                                               >  2:1509346/1‑148 (MQ=255)
       gAGGGTAACAACACCGGGATGCTGGAAGTTAATAATAATTCAGCTTTCAATAACCGCGGCGAGTTTATTCTTGATAACGACAAGAATGCTGTGCACATCAACCAGTCCGGTACGCTTTATAATACCGGTCACATGAACATCAGTAAtt                                                                                                                                              >  1:2232792/1‑148 (MQ=255)
          ggTAACAACACCGGGATGCTGGAAGTTAATAATAATTCAGCTTTCAATAACCGCGGCGAGTTTATTCTTGATAACGACAAGAATGCTGTGCACATCAACCAGTCCGGTACGCTTTATAATACCGGTCACATGAACATCAGTAAttctt                                                                                                                                           <  2:2655760/148‑1 (MQ=255)
          ggTAACAACACCGGGATGCTGGAAGTTAATAATAATTCAGCTTTCAATAACCGCGGCGAGTTTATTCTTGATAACGACAAGAATGCTGTGCACATCAACCAGTCCGGTACGCTTTATAATACCGGTCACATGAACATAAGTAAttctt                                                                                                                                           >  1:3785380/1‑148 (MQ=255)
            tAACAACACCGGGATGCTGGAAGTTAATAATAATTCAGCTTTCAATAACCGCGGCGAGTTTATTCTTGATAACGACAAGAATGCTGTGCACATCAACCAGTCCGGTACGCTTTATAATACCGGTCACATGAACATCAGTAATTCTTcc                                                                                                                                         >  1:1541776/1‑148 (MQ=255)
                    ccGTGATGCTGGAAGTTAATAATAATTCAGCTTTCAATAACCGCGGAGTGTTTATTCTTGATAACGACAAGAATGCTGTGCACATCAACCAGTCCGGTAAGCTTTATAATACCGGTCACATGAACATCAGTAATTCTTCCCACAACgg                                                                                                                                 <  2:241581/148‑1 (MQ=255)
                      gggATGCTGGAAGTTAATAATAATTCAGCTTTCAATAACCGCGGCGAGTTTATTCTTGATAACGACAAGAATGCTGTGCACATCAACCAGTCCGGTACGCTTTATAATACCGGTCACATGAACATCAGTAATTCTTCCCACAACGGAg                                                                                                                               <  1:71310/148‑1 (MQ=255)
                          tGCTGGAAGTTAATAATAATTCAGCTTTCAATAACCGCGGCGAGTTTATTCTTGATAACGACAAGAATGCTGTGCACATCAACCAGTCCGGTACGCTTTATAATACCGGTCACATGAACATCAGTAATTCTTCCCACAACGGAGCCGt                                                                                                                           <  2:3703390/148‑1 (MQ=255)
                          tGCTGGAAGTTAATAATAATTCAGCTTTCAATAACCGCGGCGAGTTTATTCTTGATAACGACAAGAATGCTGTGCACATCAACCAGTCCGGTACGCTTTATAATACCGGTCACATGAACATCAGTAATTCTTCCCACAACGGAGCCGt                                                                                                                           <  1:1464624/148‑1 (MQ=255)
                           gCTGGAAGTTAATAATAATTCAGCTTTCAATAACCGCGGCGAGTTTATTCTTGATAACGACAAGAATGCTGTGCACATCAACCAGTCCGGTACGCTTTATAATACCGGTCACATGAACATCAGTAATTCTTCCCACAACGGAGCCGtt                                                                                                                          >  1:3347857/1‑148 (MQ=255)
                            cTGGAAGTTAATAATAATTCAGCTTTCAATAACCGCGGCGAGTTTATTCTTGATAACGACAAGAATGCTGTGCACATCAACCAGTCCGGTACGCTTTATAATACCGGTCACATGAACATCAGTAATTCTTCCCACAACGGAGCCGTTa                                                                                                                         <  2:2824540/148‑1 (MQ=255)
                            cTGGAAGTTAATAATAATTCAGCTTTCAATAACCGCGGCGAGTTTATTCTTGATAACGACAAGAATGCTGTGCACATCAACCAGTCCGGTACGCTTTATAATACCGGTCACATGAACATCAGTAATTCTTCCCACAACGGAGCCGTTa                                                                                                                         <  1:2375165/148‑1 (MQ=255)
                            cTGGAAGTTAATAATAATTCAGCTTTCAATAACCGCGGCGAGTTTATTCTTGATAACGACAAGAATGCTGTGCACATCAACCAGTCCGGTACGCTTTATAATACCGGTCACATGAACATCAGTAATTCTTCCCACAACGGAGCCGTTa                                                                                                                         <  1:367335/148‑1 (MQ=255)
                             tGGAAGTTAATAATAATTCAGCTTTCTATAACCGCGGCGAGTTTATTCTTGATAACGACAAGAATGCTGTGCACATCAACCAGTCCGGTACGCTTTATAATACCGGTCACATGAACATCAGTAATTCTTCCCACAACGGAGCCGTTaa                                                                                                                        <  1:3556038/148‑1 (MQ=255)
                              ggAAGTTAATAATAATTCAGCTTTCAATAACCGCGGCGAGTTTATTCTTGATAACGACAAGAATGCTGTGCACATCAACCAGTCCGGTACGCTTTATAATACCGGTCACATGAACATCAGTAATTCTTCCCACAACGGAGCCGTTaaa                                                                                                                       >  1:998649/1‑147 (MQ=255)
                              ggAAGTTAATAATAATTCAGCTTTCAATAACCGCGGCGAGTTTATTCTTGATAACGACAAGAATGCTGTGCACATCAACCAGTCCGGTACGCTTTATAATACCGGTCACATGAACATAAGTAATTCTTCCCACAACGGAGCCGTTAat                                                                                                                       <  1:2593547/148‑1 (MQ=255)
                               gAAGTTAATAATAATTCAGCTTTCAATAACCGCGGCGAGTTTATTCTTGATAACGACAAGAATGCTGTGCACATCAACCAGTCCGGTACGCTTTATAATACCGGTCACATGAACATCAGTAATTCTTCCCACAACGGAGCCGTTAata                                                                                                                      <  1:1152743/148‑1 (MQ=255)
                               gAAGTTAATAATAATTCAGCTTTCAATAACCGCGGCGAGTTTATTCTTGATAACGACAAGAATGCTGTGCACATCAACCAGTCCGGTACGCTTTATAATACCGGTCACATGAACATAAGTAATTCTTCCCACAACGGAGCCGTTAata                                                                                                                      <  1:3442756/148‑1 (MQ=255)
                               gAAGTGAATAATAATTCAGCTTTCAATAACCGCGGCGAGTTTATTCTTGATAACGACAAGAATGCTGTGCACATCAACCAGTCCGGTACGCTTTATAATACCGGTCACATGAACATCAGTAATTCTTCCCACAACGGAGCCGTTAata                                                                                                                      <  2:2863757/148‑1 (MQ=255)
                                taGTTAATAATAATTCAGCTTTCAATAACCGCGGCGAGTTTATTCTTGATAACGACAAGAATGCTGTGCACATCAACCAGTCCGGTACGCTTTATAATACCGGTCACATGAACATCAGTAATTCTTCCCACAACGGAGCCGTTAatat                                                                                                                     <  1:3716981/147‑1 (MQ=255)
                                aaGTTAATAATAATTCAGCTTTCAATAACCGCGGCGAGTTTATTCTTGATAACGACAAGAATGCTGTGCACATCAACCAGTCCGGTACGCTTTATAATACCGGTCACATGAACATCAGTAATTCTTCCCACAACGGAGCCGTTAatat                                                                                                                     <  1:865969/148‑1 (MQ=255)
                                 aGTTAATAATAATTCAGCTTTCAATAACCGCGGCGAGTTTATTCTTGATAACGACAAGAATGCTGTGCACATCAACCAGTCCGGTACGCTTTATAATACCGGTCACATGAACATCAGTAATTCTTCCCACAACGGAGCCGTTAATAtg                                                                                                                    <  2:1198744/148‑1 (MQ=255)
                                  gTTAATAATAATTCAGCTTTCAATAACCGCGGCGAGTTTATTCTTGATAACGACAAGAATGCTGTGCACATCAACCAGTCCGGTACGCTTTATAATACCGGTCACATGAACATAAGTAATTCTTCCCACAACGGAGCCGTTAATAtgt                                                                                                                   <  2:3674663/148‑1 (MQ=255)
                                    taataataaTTCAGCTTTCAATAACCGCGGCGAGTTTATTCTTGATAACGACAAGAATGCTGTGCACATCAACCAGTCCGGTACGCTTTATAATACCGGTCACATGAACATCAGTAATTCTTCCCACAACGGAGCCGTTAATATGTgg                                                                                                                 >  2:1922334/1‑148 (MQ=255)
                                      ataataatTCAGCTTTCAATAACCGCGGCGAGTTTATTCTTGATAACGACAAGAATGCTGTGCACATCAACCAGTCCGGTACGCTTTATAATACCGGTCACATGAACATCAGTAATTCTTCCCACAACGGAGCCGTTAATATTTgggg                                                                                                               >  1:1260928/1‑148 (MQ=255)
                                         ataatTCAGCTTTCAATAACCGCGGCGAGTTTATTCTTGATAACGACAAGAATGCTGTGCACATCAACCAGTCCGGTACGCTTTATAATACCGGTCACATGAACATCAGTAATTCTTCCCACAACGGAGCCGTTAATATGTGGGGCgg                                                                                                            >  2:1370034/1‑148 (MQ=255)
                                         ataatTCAGCTTTCAATAACCGCGGCGAGTTTATTCTTGATAACGACAAGAATGCTGTGCACATCAACCAGTCCGGTACGCTTTATAATACCGGTCACATGAACATAAGTAATTCTTCCCACAACGGAGCCGTTAATATGTGGGGCgg                                                                                                            >  2:3274398/1‑148 (MQ=255)
                                          taatTCAGCTTTCAATAACCGCGGCGAGTTTATTCTTGATAACGACAAGAATGCTGTGCACATCAACCAGTCCGGTACGCTTTATAATACCGGTCACATGAACATCAGTAATTCTTCCCACAACGGAGCCGTTAATATGTGGGGCGGa                                                                                                           >  2:948608/1‑148 (MQ=255)
                                           aatTCAGCTTTCAATAACCGCGGCGAGTTTATTCTTGATAACGACAAGAATGCTGTGCACATCAACCAGTCCGGTACGCTTTATAATACCGGTCACATGAACATCAGTAATTCTTCCCACAACGGAGCCGTTAATATGTGGGGCGGaa                                                                                                          >  1:2205614/1‑148 (MQ=255)
                                             ttCAGCTTTCAATAACCGCGGCGAGTTTATTCTTGATAACGACAAGAATGCTGTGCACATCAACCAGTCCGGTACGCTTTATAATACCGGTCACATGAACATCAGTAATTCTTCCCACAACGGAGCCGTTAATATGTGGGGCGGAAAt                                                                                                        <  2:3732243/148‑1 (MQ=255)
                                             ttCAGCTTTCAATAACCGCGGCGAGTTTATTCTTAATAACGACAAGAATGCTGTGCACATCAACCAGTCCGGTACGCTTTATAATACCGGTCACATGAACATCAGTAATTCTTCCCACAACGGAGCCGTTAATATGTGGGGCGGAAAt                                                                                                        <  1:707327/148‑1 (MQ=255)
                                              tCAGCTTTCAATAACCGCGGCGAGTTTATTCTTGATAACGACAAGAATGCTGTGCACATCAACCAGTCCGGTACGCTTTATAATACCGGTCACATGAACATCAGTAATTCTTCCCACAACGGAGCCGTTAATATGTGGGGCGGAAATg                                                                                                       <  1:3106002/148‑1 (MQ=255)
                                                aGCTTTCAATAACCGCGGCGAGTTTATTCTTGATAACGACAAGAATGCTGTGCACATCAACCAGTCCGGTACGCTTTATAATACCGGTCACATGAACATCAGTAATTCTTCCCACAACGGAGCCGTTAATATGTGGGGCGGAAATGGt                                                                                                     <  2:888896/148‑1 (MQ=255)
                                                aGCTTTCAATAACCGCGGCGAGTTTATTCTTGATAACGACAAGAATGCTGTGCACATCAACCAGTCCGGTACGCTTTATAATACCGGTCACATGAACATCAGTAATTCTTCCCACAACGGAGCCGTTAATATGTGGGGCGGAAATGGt                                                                                                     <  1:2489897/148‑1 (MQ=255)
                                                aGCTTTCAATAACCGCGGCGAGTTTATTCTTGATAACGACAAGAATGCTGTGCACATCAACCAGTCCGGTACGCTTTATAATACCGGTCACATGAACATCAGTAATTCTTCCCACAACGGAGCCGTTAATATGTGGGGCGGAAATGGt                                                                                                     <  2:3707328/148‑1 (MQ=255)
                                                aGCTTTCAATAACCGCGGCGAGTTTATTCTTGATAACGACAAGAATGCTGTGCACATCAACCAGTCCGGTACGCTTTATAATACCGGTCACATGAACATCAGTAATTCTTCCCACAACGGAGCCGTTAATATGTGGGGCGGAAATGGt                                                                                                     >  2:1087106/1‑148 (MQ=255)
                                                  cTTTCAATAACCGCGGCGAGTTTATTCTTGATAACGACAAGAATGCTGTGCACATCAACCAGTCCGGTACGCTTTATAATACCGGTCACATGAACATAAGTAATTCTTCCCACAACGGAGCCGTTAATATGTGGGGCGGAAATGGTCg                                                                                                   >  2:1820585/1‑148 (MQ=255)
                                                  cTTTCAATAACCGCGGCGAGTTTATTCTTGATAACGACAAGAATGCTGTGCACATCAACCAGTCCGGTACGCTTTATAATACCGGTCACATGAACATAAGTAATTCTTCCCACAACGGAGCCGTTAATATGTGGGGCGGAAATGGTCg                                                                                                   >  2:2344890/1‑148 (MQ=255)
                                                   tttCAATAACCGCGGCGAGTTTATTCTTGATAACGACAAGAATGCTGTGCACATCAACCAGTCCGGTACGCTTTATAATACCGGTCACATGAACATCAGTAATTCTTCCCACAACGGAGCCGTTAATATGTGGGGCGGAAATGGTCGt                                                                                                  >  2:176365/1‑148 (MQ=255)
                                                    ttCAATAACCGCGGCGAGTTTATTCTTGATAACGACAAGAATGCTGTGCACATCAACCAGTCCGGTACGCTTTATAATACCGGTCACATGAACATCAGTAATTCTTCCCACAACGGAGCCGTTAATATGTGGGGCGGAAATGGTCGtt                                                                                                 >  1:3405773/1‑148 (MQ=255)
                                                     tCAATAACCGCGGCGAGTTTATTCTTGATAACGACAAGAATGCTGTGCACATCAACCAGTCCGGTACGCTTTATAATACCGGTCACATGAACATAAGTAATTCTTCCCACAACGGAGCCGTTAATATGTGGGGCGGAAATGGTCGttt                                                                                                <  1:480649/148‑1 (MQ=255)
                                                        aTAACCGCGGCGAGTTTATTCTTGATAACGACAAGAATGCTGTGCACATCAACCAGTCCGGTACGCTTTATAATACCGGTCACATGAACATCAGTAATTCTTCCCACAACGGAGCCGTTAATATGTGGGGCGGAAATGGTCGTTTTAt                                                                                             >  1:412136/1‑148 (MQ=255)
                                                         tAACCGCGGCGAGTTTATTCTTGATAACGACAAGAATGCTGTGCACATCAACCAGTCCGGTACGCTTTATAATACCGGTCACATGAACATCAGTAATTCTTCCCACAACGGAGCCGTTAATATGTGGGGCGGAAATGGTCGTTTTATc                                                                                            >  2:3568825/1‑148 (MQ=255)
                                                            ccGCGGCGAGTTTATTCTTGATAACGACAAGAATGCTGTGCACATCAACCAGTCCGGTACGCTGTATAATACCGGTCACATGAACATAAGTAATTCTTCCCACAACGGAGCCGTGAAAATGTGGGGGGGCAATTGTCGTTTTATCAAt                                                                                         >  1:2881605/1‑148 (MQ=255)
                                                             cgcgGCGAGTTTATTCTTGATAACGACAAGAATGCTGTGCACATCAACCAGTCCGGTACGCTTTATAATACCGGTCACATGAACATCAGTAATTCTTCCCACAACGGAGCCGTTAATATGTGGGGCGGAAATGGTCGTTTTATCAATg                                                                                        <  2:2371126/148‑1 (MQ=255)
                                                               cggcgAGTTTATTCTTGATAACGACAAGAATGCTGTGCACATCAACCAGTCCGGTACGCTTTATAATACCGGTCACATGAACATCAGTAATTCTTCCCACAACGGAGCCGTTAATATGTGGGGCGGAAATGGTCGTTTTATCAATGAc                                                                                      <  2:63075/148‑1 (MQ=255)
                                                                ggcgAGTTTATTCTTGATAACGACAAGAATGCTGTGCACATCAACCAGTCCGGTACGCTTTATAATACCGGTCACATGAACATAAGTAATTCTTCCCACAACGGAGCCGTTAATATGTGGGGCGGAAATGGTCGTTTTATCAATGACg                                                                                     <  1:1209568/148‑1 (MQ=255)
                                                                  cgAGTTTATTCTTTATAACGACAAGAATGCTGTGCACATCCACCAGTCCGGTACGCTTTATAATACCGGTCACATGAACATCAGTAATTCTTCCCACAACGGAGCCGTTAATTTGTGGGGCGGAAATGGTCGTTTTATCAATGACGGa                                                                                   <  1:3114859/148‑1 (MQ=255)
                                                                  cgAGTTTATTCTTGATAACGACAAGAATGCTGTGCACATCAACCAGTCCGGTACGCTTTATAATACCGGTCACATGAACATCAGTAATTCTTCCCACAACGGAGCCGTTAATATGTGGGGCGGAAATGGTCGTTTTATCAATGACGGa                                                                                   >  2:3621925/1‑148 (MQ=255)
                                                                  cgAGTTTATTCTTGATAACGACAAGAATGCTGTGCACATCAACCAGTCCGGTACGCTTTATAATACCGGTCACATGAACATCAGTAATTCTTCCCACAACGGAGCCGTTAATATGTGGGGCGGAAATGGTCGTTTTATCAATGACGGa                                                                                   >  1:1848046/1‑148 (MQ=255)
                                                                     gTTTATTCTTGATAACGACAAGAATGCTGTGCACATCAACCAGTCCGGTACGCTTTATAATACCGGTCACATGAACATAAGTAATTCTTCCCACAACGGAGCCGTTAATATGTGGGGCGGAAATGGTCGTTTTATAAATGACGGAACg                                                                                >  2:3541503/1‑148 (MQ=255)
                                                                       ttATTCTTGATAACGACAAGAATGCTGTGCACATCAACCAGTCCGGTACGCTTTATAATACCGGTCACATGAACATCAGTAATTCTTCCCACAACGGAGCCGTTAATATGTGGGGCGGAACTGGTCGTTTTATCAATGACGGAACGAt                                                                              >  2:1985130/1‑148 (MQ=255)
                                                                       ttATTCTTGATAACGACAAGAATGCTGTGCACATCAACCAGTCCGGTACGCTTTATAATACCGGTCACATGAACATCAGTAATTCTTCCCACAACGGAGCCGTTAATATGTGGGGCGGAAATGGTCGTTTTATCAATGACGGAAAGAt                                                                              >  2:2558674/1‑148 (MQ=255)
                                                                         aTTCTTGATAACGACAAGAATGCTGTGCACATCAACCAGTCCGGTACGCTTTATAATACCGGTCACATGAACATCAGTAATTCTTCCCACAACGGAGCCGTTAATATGTGGGGCGGAAATGGTCGTTTTATCAATGACGGAACGATTg                                                                            <  1:20356/148‑1 (MQ=255)
                                                                         aTTCTTGATAACGACAAGAATGCTGTGCACATCAACCAGTCCGGTACGCTTTATAATACCGGTCACATGAACATCAGTAATTCTTCCCACAACGGAGCCGTTAATATGTGGGGCGGAAATGGTCGTTTTATCAATGACGGAACGATTg                                                                            <  1:3630465/148‑1 (MQ=255)
                                                                           tCTTGATAACGACAAGAATGCTGTGCACATCAACCAGTCCGGTACGCTTTATAATACCGGTCACATGAACATCAGTAATTCTTCCCACAACGGAGCCGTTAATATGTGGGGCGGAAATGGTCGTTTTATCAATGAAGGAACGATTGAt                                                                          >  2:1326280/1‑148 (MQ=255)
                                                                             ttGATAACGACAAGAATGCTGTGCACATCAACCAGTCCGGTACGCTTTATAATACCGGTCACATGAACATCAGTAATTCTTCCCACAACGGAGCCGTTAATATGTGGGGCGGAAATGGTCGTTTTATCAATGACGGAACGATTGATGt                                                                        <  2:1613403/148‑1 (MQ=255)
                                                                              tGATAACGACAAGAATGCTGTGCACATCAACCAGTCCGGTACGCTTTATAATACCGGTCACATGAACATCAGTAATTCTTCCCACAACGGAGCCGTTAATATGTGGGGCGGAAATGGTCGTTTTATCAATGACGGAACGATTGATGtt                                                                       >  2:2044878/1‑148 (MQ=255)
                                                                              tGATAACGACAAGAATGCTGTGCACATCAACCAGTCCGGTACGCTTTATAATACCGGTCACATGAACATCAGTAATTCTTCCCACAACGGAGCCGTTAATATGTGGGGCGGAAATGGTCGTTTTATCAAATACGGAACGAGTGATGtt                                                                       >  1:3198284/1‑148 (MQ=255)
                                                                               gATAACGACAAGAATGCTGTGCACATCAACCAGTCCGGTACGCTTTATAATACCGGTCACATGAACATCAGTAATTCTTCCCACAACGGAGCCGTTAATATGTGGGGCGGAAATGTTCGTTTTATCAATGACGGAACGATTGATGttt                                                                      <  2:3106475/148‑1 (MQ=255)
                                                                               gATAACGACAAGAATGCTGTGCACATCAACCAGTCCGGTACGCTTTATAATACCGGTCACATGAACATCAGTAATTCTTCCCACAACGGAGCCGTTAATATGTGGGGCGGAAATGGTCGTTTTATCAATGACGGAACGATTGATGttt                                                                      <  2:2482179/148‑1 (MQ=255)
                                                                                   aCGACAAGAATGCTGTGCACATCAACCAGTCCGGTACGCTTTATAATACCGGTCACATGAACATCAGTAATTCTTCCCACAACGGAGCCGTGAATATGTGGGGCGGAAATGGTCGTTTTATAAAATACGGAACGATTGATGTTTCTGc                                                                  >  1:2698632/1‑148 (MQ=255)
                                                                                     gACAAGAATGCTGTGCACATCAACCAGTCCGGTACGCTTTATAATACCGGTCACATGAACATCAGTAATTCTTCCCACAACGGAGCCGTTAATATGTGGGGCGGAAATGGTCGTTTTATCAATGACGGAACGATTGATGTTTCTGCGa                                                                >  2:1354875/1‑148 (MQ=255)
                                                                                         aGAATGCTGTGCACATCAACCAGTCCGGTACGCTTTATAATACCGGTCACATGAACATCAGTAATTCTTCCCACAACGGAGCCGTTAATATGTGGGGCGGAAATGGTCGTTTTATCAATGACGGAACGATTGATGTTTCTGCGAAGTc                                                            <  2:3819065/148‑1 (MQ=255)
                                                                                         aGAATGCTGTGCACATCAACCAGTCCGGTACGCTTTATAATACCGGTCACATGAACATCAGTAATTCTTCCCACAACGGAGCCGTTAATATGTGGGGCGGAAATGGTCGTTTTATCAATGACGGAACGATTGATGTTTCTGCGAAGTc                                                            >  2:473957/1‑148 (MQ=255)
                                                                                         aGAATGCTGTGCACATCAACCAGTCCGGTACGCTTTATAATACCGGTCACATGAACATCAGTAATTCTTCCCACAACGGAGCCGTTAATATGTGGGGCGGAAATGGTCGTTTTATCAATGACGGAACGATTGATGTTTCTGCGAAGTc                                                            >  2:1532357/1‑148 (MQ=255)
                                                                                          gAATGCTGTGCACATCAACCAGTCCGGTACGCTTTATAATACCGGTCACATGAACATCAGTAATTCTTCCCACAACGGAGCCGTTAATATGTGGGGCGGAAATGGTCGTTTTAAAAATGACGGCCCGAGTGCTGGTTCTGCGAAGTCa                                                           >  2:2130846/1‑148 (MQ=255)
                                                                                            aTGCTGTGCACATCAACCAGTCCGGTACGCTTTATAATACCGGTCACATGAACATCAGTAATTCTTCCCACAACGGAGCCGTTAATATGTGGGGCGGAAATGGTCGTTTTATCAATGACGGAACGATTGATGTTTCTGCGAAGTCACt                                                         >  2:1658585/1‑148 (MQ=255)
                                                                                              gCTGTGCACATCAACCAGTCCGGTACGCTTTATAATACCGGTCACATGAACATCAGTAATTCTTCCCACAACGGAGCCGTTAATATGTGGGGCGGAAATGGTCGTTTTATCAATGACGGAACGATTGATGTTTCTGCGAAGTCACTgg                                                       <  2:3143676/148‑1 (MQ=255)
                                                                                               cTGTGCACATCAACCAGTCCGGTACGCTTTATAATACCGGTCACATGAACATCAGTAATTCTTCCCACAACGGAGCCGTTAATATGTGGGGCGGAAATGGTCGTTTTATCAATGACGGAACGATTGATGTTTCTGCGAAGTCACTGGt                                                      <  1:949671/148‑1 (MQ=255)
                                                                                                 gtgCACATCAACCAGTCCGGTACGCTTTATAATACCGGTCACATGAACATCAGTAATTCTTCCCACAACGGAGCCGTTAATATGTGGGGCGGAAATGGTCGTTTTATCAATGACGGAACGATTGATGTTTCTGCGAAGTCACTGGTAg                                                    >  2:2905891/1‑148 (MQ=255)
                                                                                                  tgCACATCAACCAGTCCGGTACGCTTTATAATACCGGTCACATGAACATCAGTAATTCTTCCCACAACGGAGCCGTTAATATGTGGGGCGGAAATGGTCGTTTTATCAATGACGGAACGATTGATGTTTCTGCGAAGTCACTGGTAGt                                                   <  2:670359/148‑1 (MQ=255)
                                                                                                    cacaTCAACCAGTCCGGTACGCTTTATAATACCGGTCACATGAACATCAGTAATTCTTCCCACAACGGAGCCGTGAATATGTGGGGCGGAAATGGTCGTTTTATAAATTAAGGAAAGAATGTTTTTTTTTGGAAGTCACAGGTAGTCa                                                 >  1:1617111/1‑148 (MQ=255)
                                                                                                       aTCAACCAGTCCGGTACGCTTTATAATACCGGTCACATGAACATAAGTAATTCTTCCCACAACGGAGCCGTTAATATGTGGGGCGGAAATGGTCGTTTTATCAATGACGGAACGATTGATGTTTCTGCGAAGTCACTGGTAGTCAgcg                                              <  2:1118242/148‑1 (MQ=255)
                                                                                                         cAACCAGTCCGGTACGCTTTATAATACCGGTCACATGAACATCAGTAATTCTTCCCACAACGGAGCCGTTAATATGTGGGGCGGAAATGGTCGTTTTATAAATGACGGAACGATTGATGTTTCTGCGAAGTCACTGGTAGTCAGCGCt                                            >  1:3330870/1‑148 (MQ=255)
                                                                                                           aCCAGTCCGGTACGCTTTATAATACCGGTCACATGAACATCAGTAATTCTTCCCACAACGGAGCCGTTAATATGTGGGGCGGAAATGGTCGTTTTATCAATGACGGAACGATTGATGTTTCTGCGAAGTCACTGGTAGTCAGCGCtaa                                          <  2:3236856/148‑1 (MQ=255)
                                                                                                                      aCGCTTTATAATACCGGTCACATGAACATCAGTAATTCTTCCCACAACGGAGCCGTTAATATGTGGGGCGGAAATGGTCGTTTTATCAATGACGGAACGATTGATGTTTCTGCGAAGTCACTGGTAGTCAGCGCTAATAATGCCGGCg                               >  2:1893522/1‑148 (MQ=255)
                                                                                                                       cGCTTTATAATACCGGTCACATGAACATCAGTAATTCTTCCCACAACGGAGCCGTTAATATGTGGGGCGGAAATGGTCGTTTTATCAATGACGGAACGATTGATGTTTCTGCGAAGTCACTGGTAGTCAGCGCTAATAATGCCGGCGa                              <  1:1738466/148‑1 (MQ=255)
                                                                                                                          tttATAATACCGGTCACATGAACATCAGTAATTCTTCCCACAACGGAGCCGTTAATATGTGGGGCGGAAATGGTCGTTTTATAAATGACGGAACGATTGATGTTTCTGCGAAGTCACTGGTAGTCAGCGCTAATAATGCCGGCGATCa                           >  1:193354/1‑148 (MQ=255)
                                                                                                                            tataATACCGGTCACATGAACATCAGTAATTCTTCCCACAACGGAGCCGTTAATATGTGGGGCGGAAATGGTCGTTTTATCAATGACGGAACGATTGATGTTTCTGCGAAGTCACTGGTAGTCAGCGCTAATAATGCCGGCGATCAGa                         <  2:702150/148‑1 (MQ=255)
                                                                                                                            tataATACCGGTCACATGAACATCAGTAATTCTTCCCACAACGGAGCCGTTAATATGTGGGGCGGAAATGGTCGTTTTATCAATGACGGAACGATTGATGTTTCTGCGAAGTCACTGGTAGTCAGCGCTAATAATGCCGGCGATCAGa                         <  1:3361414/148‑1 (MQ=255)
                                                                                                                            tataAGACCGGTCACATGAACATCAGTAATTCTTCCCACAACGGAGCCGTTAATATGTGGGGCGGAAATGGTCGTTTTATCAATGACGGAACGATTGATGTTTCTGCGAAGTCACTGGTAGTCAGCGCTAATAATGCCGGCGATCAGa                         <  2:2728953/148‑1 (MQ=255)
                                                                                                                             ataataCCGGTCACATGAACATCAGTAATTCTTCCCACAACGGAGCCGTTAATATGTGGGGCGGAAATGGTCGTTTTATCAATGACGGAACGATTGATGTTTCTGCGAAGTCACTGGTAGTCAGCGCTAATAATGCCGGCGATCAGaa                        <  2:3430107/148‑1 (MQ=255)
                                                                                                                               aataCCGGTCACATGAACATCAGTAATTCTTCCCACAACGGAGCCGTTAATATGTGGGGCGGAAATGGTCGTTTTATCAATGACGGAACGAGTGATGTTTCTGCGAAGTCACTGGTAGTCAGCGCTAATAATGCCGGCGATCAGAATg                      >  2:3280767/1‑148 (MQ=255)
                                                                                                                                   ccGGTCACATGAACATCAGTAATTCTTCCCACAACGGAGCCGTTAATATGTGGGGCGGAAATGGTCGTTTTATCAATGACGGAACGATTGATGTTTCTGCGAAGTCACTGGTAGTCAGCGCTAATAATGCCGGCGATCAGAATGCctt                  <  2:2208659/148‑1 (MQ=255)
                                                                                                                                   ccGGTCACATGAACATCAGTAATTCTTCCCACAACGGAGCCGTTAATATGTGGGGCGGAAATGGTCGTTTTATCAATGACGGAACGATTGATGTTTCTGCGAAGTCACTGGTAGTCAGCGCTAATAATGCCGGCGATCAGAATGCctt                  <  2:3419726/148‑1 (MQ=255)
                                                                                                                                   ccGGTCACATGAACATCAGTAATTCTTCACACAACGGAGCCCTTAATATGACGGGCGCACATGGTCGTTTTATAAAAGACGGAACGATATATGTTACTGCGAAGACCCTGGTAGTAAGAGCTAATAATGCCGGCGAACAGaaagaata                  >  2:3176373/1‑142 (MQ=255)
                                                                                                                                   ccGGTCACATGAACATAAGTAATTCTTCCCACAACGGAGCCGTTAATATGTGGGGCGGAAATGGTCGTTTTATCAATGACGGAACGATTGATGTTTCTGCGAAGTCACTGGTAGTCAGCGCTAATAATGCCGGCGATCAGAATGCctt                  <  1:1358986/148‑1 (MQ=255)
                                                                                                                                    cGGTCACATGAACATCAGTAATTCTTCCCACAACGGAGCCGTTAATATGTGGGGCGGAAATGGTCGTTTTATCAATGACGGAACGATTGATGTTTCTGCGAAGTCACTGGTAGTCAGCGCTAATAATGCCGGCGATCAGAATGCcttc                 <  1:352130/148‑1 (MQ=255)
                                                                                                                                       tCACATGAACATCAGTAATTCTTCCCACAACGGAGCCGTTAATATGTGGGGCGGAAATGGTCGTTTTATCAATGACGGAACGATTGATGTTTCTGCGAAGTCACTGGTAGTCAGCGCTAATAATGCCGGCGATCAGAATGCCttcttc              <  2:1848046/148‑1 (MQ=255)
                                                                                                                                        cacaTGAACATAAGTAATTCTTCCCACAACGGAGCCGTTAATATGTGGGGCGGAAATGGTCGTTTTATCAATGACGGAACGATTGATGTTTCTGCGAAGTCACTGGTAGTCAGCGCTAATAATGCCGGCGATCAGAATGCCTtcttct             <  1:327041/148‑1 (MQ=255)
                                                                                                                                          caTGAACATCAGTAATTCTTCCCACAACGGAGCCGTTAATATGTGGGGCGGAAATGGTCGTTTTATCAATGACGGAACGATTGATGTTTCTGCGAAGTCACTGGTAGTCAGCGCTAATAATGCCGGCGATCAGAATGCCTTCTTCTgg           <  1:1509346/148‑1 (MQ=255)
                                                                                                                                           aTGAACATCAGTAATTCTTCCCACCCCGGAGCCGTTAATATGTGGGGCGGAAATGGTCGTTTTATCAATGACGGAACGATTGATGTTTCTGCGAAGTCACTGGTAGTCAGCGCTAATAATGCCGGCGATCAGAATGCCTTCTTCTGGa          <  2:735845/148‑1 (MQ=255)
                                                                                                                                            tGAACATCAGTAATTCTTCCCACAACGGAGCCGTTAATATGTGGGGCGGAAATGGTCGTTTTATCAATGACGGAACGATTGATGTTTCTGCGAAGTCACTGGTAGTCAGCGCTAATAATGCCGGCGATCAGAATGCCTTCTTCTGGaa         <  2:2140870/148‑1 (MQ=255)
                                                                                                                                                 aTCAGTAATTCTTCCCACAACGGAGCCGTTAATATGTGGGGCGGAAATGGTCGTTTTATCAATGACGGAACGATTGATGTTTCTGCGAAGTCACTGGTAGTCCGCGCTAATAATGCCGGCGATCAGAATGCCTTCTTCTGGAACCAgg    >  2:2982770/1‑148 (MQ=255)
                                                                                                                                                   cAGTAATTCTTCCCACAACGGAGCCGTTAATATGTGGGGCGGAAATGGTCGTTTTATCAATGACGGAACGATTGATGTTTCTGCGAAGTCACTGGTAGTCAGCGCTAATAATGCCGGCGATCAGAATGCCTTCTTCTGGAACCAGGAt  <  1:2313881/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
TAGTAGTGAGGGTAACAACACCGGGATGCTGGAAGTTAATAATAATTCAGCTTTCAATAACCGCGGCGAGTTTATTCTTGATAACGACAAGAATGCTGTGCACATCAACCAGTCCGGTACGCTTTATAATACCGGTCACATGAACATCAGTAATTCTTCCCACAACGGAGCCGTTAATATGTGGGGCGGAAATGGTCGTTTTATCAATGACGGAACGATTGATGTTTCTGCGAAGTCACTGGTAGTCAGCGCTAATAATGCCGGCGATCAGAATGCCTTCTTCTGGAACCAGGAT  >  CP009273/1466146‑1466440

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: