Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 868,897 G→T 17.1% D155Y (GAC→TAC)  yliE → putative membrane‑anchored cyclic‑di‑GMP phosphodiesterase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273868,8970GT17.1% 519.8 / 66.5 174D155Y (GAC→TAC) yliEputative membrane‑anchored cyclic‑di‑GMP phosphodiesterase
Reads supporting (aligned to +/- strand):  ref base G (70/74);  new base T (13/17);  total (83/91)
Fisher's exact test for biased strand distribution p-value = 6.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 177 total aligned reads displayed.

ATAAATATTATCTTTTTTCGCGGAGGGCGAATGATAAAATCGTTGTTTATCATCAGGAACAAGTAAATGGACCGCTTGTCGACGAGTCAGGGCGGGTTATTTTTGCCGATTTTAACCCATCGAAACGACCGTGGTCGGTGGCTTCAGATGACTCTAACAACAGCTGGAATCCGGCATACAATTGCTTTGATCGTCCGGGTAAAAAATGTATCTCTTTTACGCTACACATCAACGGCAAAGATCACGATTTGTTAGCGGTGGATAAAATTCATGTCGATTTAAACTGGCGATATCTGA  >  CP009273/868748‑869044
                                                                                                                                                      |                                                                                                                                                   
ataaataTTATCTTTTTTCGCGGAGGGCGAATGATAAAATCGTTGTTTATCATCAGGAACAAGTAAATGGACCGCTTGTCGACGAGTCAGGGCGGGTTATTTTTGCCGATTTTAACCCATCGAAACGACCGTGGTCGGTGGCTTCAGATg                                                                                                                                                     >  1:4346518/1‑150 (MQ=255)
 taaataTTATCTTTTTTCGCGGAGGGCGAATGATAAAATCGTTGTTTATCATCAGGAACAAGTAAATGGACCGCTTGTCGACGAGTCAGGGCGGGTTATTTTTGCCGATTTTAACCCATCGAAACGACCGTGGTCGGTGGCTTCAGATGa                                                                                                                                                    >  2:354951/1‑150 (MQ=255)
 taaataTTATCTTTTTTCGCGGAGGGCGAATGATAAAATCGTTGTTTATCATCAGGAACAAGTAAATGGACCGCTTGTCGACGAGTCAGGGCGGGTTATTTTTGCCGATTTTAACCCATCGAAACGACCGTGGTCGGTGGCTTCAGATGa                                                                                                                                                    >  2:354871/1‑150 (MQ=255)
   aataTTATCTTTTTTCGCGGAGGGCGAATGATAAAATCGTTGTTTATCATCAGGAACAAGTAAATGGACCGCTTGTCGACGAGTCAGGGCGGGTTATTTTTGCCGATTTTAACCCATCGAAACGACCGTGGTCGGTGGCTTCAGATGAct                                                                                                                                                  >  2:2406164/1‑150 (MQ=255)
   aataTTATCTTTTTTCGCGGAGGGCGAATGATAAAATCGTTGTTTATCATCAGGAACAAGTAAATGGACCGCTTGTCGACGAGTCAGGGCGGGTTATTTTTGCCGATTTTAACCCATCGAAACGACCGTGGTCGGTGGCTTCAGATGAct                                                                                                                                                  >  2:4099706/1‑150 (MQ=255)
   aataTTATCTTTTTTCGCGGAGGGCCAATGATAAAATCGTTGTTTATCATCAGGAACATGTAAATGGACCGATTGACGACGCGTCAGGGCGGGTTATTTTTGCCGATTTTAACCAATCGAAACGACCGTGGTCGGTGGCTTCAGATTAct                                                                                                                                                  <  2:223781/150‑1 (MQ=255)
   aataTTATCTTTTGTCGCGGTGGGCTAATGATAAAATGGTTTTTTATACTCCTCAACAAGTAACTGGACCTCTGTTCTACGAGTCTGGGCGGGTTATTTTTGCCGATTTTAACCCAACTAAACGACCGTGGTCGGTGGCTTCAGATGAct                                                                                                                                                  <  2:1477115/150‑1 (MQ=255)
     tattatCTTTTTTCGCGGAGGGCGAATGATAAAATCGTTGTTTATCATCAGGAACAAGTAAATGGACCGCTTGTCGACGAGTCAGGGCGGGTTATTTTTGCCGATTTTAACCCATCGAAACGACCGTGGTCGGTGGCTTCAGATTActct                                                                                                                                                >  1:1605512/1‑150 (MQ=255)
     tattatATTTTTTCGCGGCGGGCGAATGATAAAATCGTTGTTTATCTTCTGGAACAATTAAATGGACCGCTTGTCGACGAGTCTGGGGGGTTTATTTTTGCCGATTTTAACCCATCTAAACGACCGTGGTCGGTGGCTTCAGATGActct                                                                                                                                                <  2:3627198/150‑1 (MQ=255)
      attatCTTTTTGCGCGGTGGGCTAATTATAAAATCGTTGTTTATCAACAGTAACAATTAAAGGGCCCGCTTGTCGACGCGTCTGGGCGGGTTATTTTGGCCGTTTTTACCCCATCGAAACGACCGTGGTCGGTGGCTTCAGATGACTCTa                                                                                                                                               <  1:2156555/150‑1 (MQ=255)
       ttttattttttGGCGGTGGGCGATTGATAAAATGGTTGTTTATCATCTGTAACAATAAATGGGACCGCTTGTCGACGTGTCGGGGGGGGTTTTTTTTGCCGATTTTAACCCATCGAAACGACCGTGGTCGGTGGCTTCAGATGACTCTaa                                                                                                                                              <  2:1225407/145‑1 (MQ=255)
       ttatATTTTTTCGGGGTGGGCTATTTATAAATTGGTTTTTTATAATCTGGAAAAATTAAATGGCCCGTTTGTCGACGTGTCTGGGGGGTTTTTTTTTGCCGATTTTAACCAATCGAAACGACCGTGGTCGGTGGCTTCAGATGACTCTaa                                                                                                                                              <  2:3296419/150‑1 (MQ=255)
        tttCTTTTTTCGCGGAGGGCTAATGATAAAATCTTTTTTTAAAATCTGTAACAAGTAAATGGCCCGCTTGTCGACGAGTCAGGGCGGGTTTTTTTTGCCTATTTTAACCCATCGAAACGACCGTGGTCGGTGGCTTCAGATGACTCTaac                                                                                                                                             <  2:1334954/148‑1 (MQ=255)
        tatCTTTTTTCGCGGAGGGCGAATGATAAAATCGTTGTTTATCATCAGGAACAAGTAAATGGACCGCTTGTCGACGAGTCAGGGCGGGTTATTTTTGCCGATTTTAACCCATCGAAACGACCGTGGTCGGTGGCTTCAGATGACTCTaac                                                                                                                                             >  1:1398723/1‑150 (MQ=255)
        tatCTTTTTTCGCGGAGGGCGAATGATAAAATCGTTGTTTATCATCAGGAACAAGTAAATGGACCGCTTGTCGACGAGTCAGGGCGGGTTATTTTTGCCGATTTTAACCCATCGAAACGACCGTGGTCGGTGGCTTCAGATGACTCTaac                                                                                                                                             >  2:2016643/1‑150 (MQ=255)
           cTTTTTTCGCGGAGGGCGAATGATAAAATCGTTGTTTATCATCAGGAACAAGTAAATGGACCGCTTGTCGACGAGTCAGGGCGGGTTATTTTTGCCGATTTTAACCCATCGAAACGACCGTGGTCGGTGGCTTCAGATGACTCTaacaac                                                                                                                                          >  2:1581583/1‑150 (MQ=255)
           cTTTTTTCGCGGAGGGCGAATGATAAAATCGTTGTTTATCATCAGGAACAAGTAAATGGACCGCTTGTCGACGAGTCAGGGCGGGTTATTTTTGCCGATTTTAACCCATCGAAACGACCGTGGTCGGTGGCTTCAGATGACTCTaacaac                                                                                                                                          >  2:3769195/1‑150 (MQ=255)
             tttttGGGGGAGGGCTATTTTTAAATTGGTTGTTTATCATCTGGAACATGTAAATGGACGGGTTGTCGACGTGTCTGGGGGTGTTTTTTTTGCCGTTTTTAACACTTCTAAACGACCGTGGTCGGTGGCTTCTGATTTCTCTAACAACAg                                                                                                                                        <  2:1649410/150‑1 (MQ=255)
                 tCGCGGTGGGCGATTTATAAATTGGTTTTTTTTAATCTGTAACAAGTAATTGGACCGCTGTTCGACGTGTCAGGGGGGGTTTTTTTTGCCGTTTTTAACCCATCGAAACGACCGTGGTCGGTGGCTTCAGATGACTCTAACAACAGCTgg                                                                                                                                    <  2:4234630/150‑1 (MQ=255)
                  cgcgGAGGGCGAATGATAAAATCGTTGTTTATCATCAGGAACAAGTAAATGGACCGCTTGTCGACGAGTCAGGGCGGGTTATTTTTGCCGATTTTAACCCATCGAAACGACCGTGGTCGGTGGCTTCAGATGACTCTAACAACAGCTGGa                                                                                                                                   >  1:919116/1‑150 (MQ=255)
                     ggAGGGCGAATGATAAAATCGTTGTTTATCATCAGGAACAAGTAAATGGACCGCTTGTCGACGAGTCAGGGCGGGTTATTTTTGCCGATTTTAACCCATCGAAACGACCGTGGTCGGTGGCTTCAGATGACTCTAACAACAGCTGGAATc                                                                                                                                >  2:2372249/1‑150 (MQ=255)
                     ggAGGGCGAATGATAAAATCGTTGTTTATCATCAGGAACAAGTAAATGGACCGCTTGTCGACGAGTCAGGGCGGGTTATTTTTGCCGATTTTAACCCATCGAAACGACCGTGGTCGGTGGCTTCAGATGACTCTAACAACAGCTGGAATc                                                                                                                                >  2:1662713/1‑150 (MQ=255)
                      gAGGGCGAATGATAAAATCGTTGTTTATCATCAGGAACAAGTAAATGGACCGCTTGTCGACGAGTCAGGGCGGGTTATTTTTGCCGATTTTAACCCATCGAAACGACCGTGGTCGGTGGCTTCAGATGACTCTAACAACAGCTGGAATcc                                                                                                                               >  1:569101/1‑150 (MQ=255)
                       tgggataaTTATAAATTGTTTTTTTATAATCAGTAAAATTTAATTGGACGGTTTTTCTACGTTTCTGGGGGGTTTTTTTTTGCCTATTTAAACCAATCTAAACGACCGTGGTCGGTGGCTTCAGATTACTCTAACAACAGCTGGAATCCg                                                                                                                              <  1:3019724/144‑1 (MQ=255)
                       aGGGCGAATGATAAAATCGTTGTTTATCATCAGGAACAAGTAAATGGACCGCTTGTCGACGAGTCAGGGCGGGTTATTTTTGCCGATTTTAACCCATCGAAACGACCGTGGTCGGTGGCTTCAGATTACTCTAACAACAGCTGGAATCCg                                                                                                                              >  2:1510217/1‑150 (MQ=255)
                       aGGGCGAATGATAAAATCGTTGTTTATCATCAGGAACAAGTAAATGGACCGCTTGTCGACGAGTCAGGGCGGGTTATTTTTGCCGATTTTAACCCATCGAAACGACCGTGGTCGGTGGCTTCAGATGACTCTAACAACAGCTGGAATCCg                                                                                                                              >  1:3668690/1‑150 (MQ=255)
                        gggCGAATGATAAAATGGTTGTTTATCATCAGGAACAAGTAAATGGACCGCTTGTCGACGAGTCAGGGCGGGTTATTTTTGCCGATTTTAACCCATCGAAACGACCGTGGTCGGTGGCTTCAGATTACTCTAACAACAGCTGGAATCCgg                                                                                                                             <  2:2702623/150‑1 (MQ=255)
                          gctaaTTATAAAATCTTTGTTTAAAATCAGTAACAAGTAAATGGACCGCTTGTCGACGAGTCTGGGGGGGTTATTTTTGCCGATTTTAACCCAACGAAACGACCGTGGTCGGTGGCTTCAGATGACTCTAACAACAGCTGGAATCCGGCa                                                                                                                           <  2:1416047/147‑1 (MQ=255)
                          gatattGAAAAAATCTTTGTTTATAATCAGGAAAATTTAAATGGCCGGATTGTCGACGTGTCTGGGGGGTTTTTTTTTGCCTTTTTTAACCCATGTAAACGACCGTGGTCGGTGGCTTCAGATGACTCTAACAACAGCTGGAATCCGGCa                                                                                                                           <  1:635386/145‑1 (MQ=255)
                          gCGAATGATAAAATCGTTTTTTTTAATCTGGAACAAGTAAATGGACCGCTTGTCGACGAGTCAGGGCGGTTTATTTTTGCCGATTTTAACCCAACGAAACGACCGTGGTCGGTGGCTTCAGATGACTCTAACAACAGCTGGAATCCGGCa                                                                                                                           <  2:1416417/150‑1 (MQ=255)
                             aaTGATAAAATCGTTGTTTATCATCAGGAACAAGTAAATGGACCGCTTGTCGACGAGTCAGGGCGGGTTATTTTTGCCGATTTTAACCCATCGAAACGACCGTGGTCGGTGGCTTCAGATGACTCTAACAACAGCTGGAATCCGGCATAc                                                                                                                        >  1:2365704/1‑150 (MQ=255)
                                  tAAAATCGTTGTTTATCATCAGGAACAAGTAAATGGACCGCTTGTCGACGAGTCAGTGAGGGTTATTTTTGCAGATTTTAACCCATCGAAACGACAGTGGTCGGTGGCTTCAGATTACTCTAACAACAGCTGGAATACGTCATACAAttt                                                                                                                   >  2:1947862/1‑149 (MQ=255)
                                  tAAAATCGTTGTTTATCATCAGGAACAAGTAAATGGACCGCTTGTCGACGAGTCAGGGCGGGTTATTTTTGCCGATTTTAACCCATCGAAACGACCGTGGTCGGTGGCTTCAGATTACTCTAACAACAGCTGGAATCCGGCATACAATTg                                                                                                                   >  2:131713/1‑150 (MQ=255)
                                  tAAAATCGTTGTTTATCATCAGGAACAAGTAAATGGACCGCTTGTCGACGAGTCAGGGCGGGTTATTTTTGCCGATTTTAACCCATCGAAACGACCGTGGTCGGTGGCTTCAGATGACTCTAACAACAGCTGGAATCCGGCATACAATTg                                                                                                                   >  2:4178802/1‑150 (MQ=255)
                                      atttttgttTATAATCATGAACTATTAATTGGACGTGTTGTCGACGAGTGTGGGGGGTTTATTTTTTCCTATTTTAACCCACCGAAACGACCGTGGTCGGTGGCTTCAGATGACTCTAACAACAGCTGGAATCCGGCATACAATTGCttt                                                                                                               <  1:3763638/146‑1 (MQ=255)
                                         gttgttTATCATCAGGAACAAGTAAATGGACCGCTTGTCGACGAGTCAGGGCGGTTTATTTTTGCCGATTTTAACCCATCGAAACGACCGTGGTCGGTGGCTTCAGATGACTCTAACAACAGCTGGAATCCGGCATACAATTGCTTTGAt                                                                                                            <  2:696564/150‑1 (MQ=255)
                                          ttttttATAATCAGGAACAATTAAATGGACCGCTTGTCGACGCGTCAGGGCGGGTTATTTTTGCCTTTTTTAACCCATCGACACGACCGTGGTCGGTGGCTTCAGATTACTCTAACAACAGCTGGAATCCGGCATACAATTGCTTTGATc                                                                                                           <  1:3538199/147‑1 (MQ=255)
                                          ttgttTAACATCAGCAACAAGAAAATGGACCGCTTGTCGACGAGTCAGGGCGGTTTTTTTTTGCCTATTTTAACCCATCGAAACGACCGTGGTCGGTGGCTTCAGATGACTCTAACAACAGCTGGAATCCGGCATACAATTGCTTTGATc                                                                                                           <  2:449061/150‑1 (MQ=255)
                                           tgttTATCATCAGTAACAAGTAAATGGCCCGCTTGTCGACGAGTCTGGGCGGGTTATTTTTGCCTATTTTAACCCATCGAAACGACCGTGGTCGGTGGCTTCAGATTACTCTAACAACAGCTGGAATCCGGCATACAATTGCTTTGATCg                                                                                                          <  2:1670435/150‑1 (MQ=255)
                                           tgttTATCATCAGGAACAAGTAAATGGACCGCTTGTCGCCGAGTCAGGGCGGTTTATTTTTGCCGATTTTAACCCATCGAAACGACCGTGGTCGGTGGCTTCAGATTACTCTAACAACAGCTGGAATCCGGCATACAATTGCTTTGATCg                                                                                                          <  2:3859971/150‑1 (MQ=255)
                                           tgttTATCATCAGGAACAAGTAAATGGACCGCTTGTCGACGAGTCAGGGCGGGTTATTTTTGCCGATTTTAACCCATCGAAACGACCGTGGTCGGTGGCTTCAGATTACTCTAACAACAGCTGGAATCCGGCATACAATTGCTTTGATCg                                                                                                          >  1:2633569/1‑150 (MQ=255)
                                           tgttTATCATCAGGAACAAGTAAATGGACCGCTTGTCGACGAGTCAGGGCGGGTTATTTTTGCCGATTTTAACCCATCGAAACGACCGTGGTCGGTGGCTTCAGATGACTCTAACAACAGCTGGAATCCGGCATACAATTGCTTTGATCg                                                                                                          >  2:3662042/1‑150 (MQ=255)
                                             tttATCATCAGGAACAAGTAAATGGACCGCTTGTCGACGAGTCAGGGCGGGTTATTTTTGCCGATTTTAACCCATCGAAACGACCGTGGTCGGTGGCTTCAGATGACTCTAACAACAGCTGGAATCCGGCATACAATTGCTTTGATCGTc                                                                                                        <  2:2825937/150‑1 (MQ=255)
                                               tATCATCAGGAACAAGTAAATGGACCGCTTGTCGACGAGTCAGGGCGGGTTATTTTTGCCGATTTTAACCCATCGAAACGACCGTGGTCGGTGGCTTCAGATGACTCTAACAACAGCTGGAATCCGGCATACAATTGCTTTGATCGTCCg                                                                                                      >  1:4200669/1‑150 (MQ=255)
                                                 tcatcaGGAACAAGTAAATGGACCGCTTGTCGACGAGTCAGGGCGGGTTATTTTTGCCGATTTTAACCCATCGAAACGACCGTGGTCGGTGGCTTCAGATGACTCTAACAACAGCTGGAATCCGGCATACAATTGCTTTGATCGTCCggg                                                                                                    >  2:469529/1‑150 (MQ=255)
                                                   atcaGGAACAAGAAAATGGACCGCTTGTCGACGAGTCTGGGCGGGTTATTTTTGCCTATTTTAACCCATCGAAACGACCGTGGTCGGTGGCTTCAGATGACTCTAACAACAGCTGGAATCCGGCATACAATTGCTTTGATCGTCCGGGTa                                                                                                  <  2:2342782/150‑1 (MQ=255)
                                                    tcaGGAACAAGTAAATGGACCGCTTGTCGACGAGTCAGGGCGGGTTATTTTTGCCGATTTTAACCCATCGAAACGACCGTGGTCGGTGGCTTCAGATGACTCTAACAACAGCTGGAATCCGGCATACAATTGCTTTGATCGTCCGGGTaa                                                                                                 >  1:3595749/1‑150 (MQ=255)
                                                         aaCAAGTAAATGGACCGCTTGTCGACGAGTCAGGGCGGGTTATTTTTGCCGATTTTAACCCATCGAAACGACCGTGGTCGGTGGCTTCAGATTACTCTAACAACAGCTGGAATCCGGCATACAATTGCTTTGATCGTCCGGGTAAAAAAt                                                                                            >  1:2280240/1‑150 (MQ=255)
                                                          aCAAGTAAATGGACCGCTTGTCGACGAGTCAGGGCGGGTTATTTTTGCCGATTTTAACCCATCGAAACGACCGTGGTCGGTGGCTTCAGATGACTCTAACAACAGCTGGAATCCGGCATACAATTGCTTTGATCGTCCGGGTAAAAAATg                                                                                           >  2:2529722/1‑150 (MQ=255)
                                                             attAAATGCCCCGATTGTCGACGAGTCAGGGGGGTTTATTTTTGCCGATTTAACCCCATCGACACGACCGTGGTCGGTGGCTTCAGATGACTCTAACAACAGCTGGAATCCGGCATACAATTGCTTTGATCGTCCGGGTAAAAAATGTAt                                                                                        <  1:101616/148‑1 (MQ=255)
                                                             aGTAAAGGGACCGCTTGTCGACGAGTCAGGGCGGGTTATTTTTGCCGATTTTAACCCATCGAAACGACCGTGGTCGGTGGCTTCAGATGACTCTAACAACAGCTGGAATCCGGCATACAATTGCTTTGATCGTCCGGGTAAAAAATGTAt                                                                                        <  2:2850536/150‑1 (MQ=255)
                                                              ttAAATGGCCCGCTTGTGGACGTGTCTGGGCGGTTTATTTTTGCCTATTTTAACCCATCTAAACGACCGTGGTCGGTGGCTTCAGATGACTCTAACAACAGCTGGAATCCGGCATACAATTGCTTTGATCGTCCGGGTAAAAAATGTAtt                                                                                       <  2:3120966/149‑2 (MQ=255)
                                                               tAATTGGACCGCTTGTCGACGAGTCAGGGCGGTTTATTTTTGCCGATTTTAACCCATCTAAACGACCGTGGTCGGTGGCTTCAGATGACTCTAACAACAGCTGGAATCCGGCATACAATTGCTTTGATCGTCCGGGTAAAAAATGTAtct                                                                                      <  1:2504827/150‑1 (MQ=255)
                                                               tAAATGGACCGCTTGTCGAGGTGTCTGGGGGGTTTTTTTTTGCGTATTTTAACCCATCGAAACGACCGTGGTCGGTGGCTTCAGATGACTCTAACAACAGCTGGAATCCGGCATACAATTGCTTTGATCGTCCGGGTAAAAAATGTAtct                                                                                      <  1:2504729/150‑1 (MQ=255)
                                                               tAAATGGACCGCTTGTCGACGAGTCAGGGCGGGTTATTTTTGCCGATTTTAACCCATCGAAACGACCGTGGTCGGTGGCTTCAGATGACTCTAACAACAGCTGGAATCCGGCATACAATTGCTTTGATCGTCCGGGTAAAAAATGTAtct                                                                                      >  2:4000449/1‑150 (MQ=255)
                                                                aattGGACGGCTTGTCGCCGAGTCTGGGCGGGTTTTTTTGGCCTATTTAAACCCATCGAAACGACCGTGGTCGGTGGCTTCAGATGACTCTAACAACAGCTGGAATCCGGCATACAATTGCTTTGATCGTCCGGGTAAAAAATGTAtctc                                                                                     <  2:793542/147‑1 (MQ=255)
                                                                aaaTGGACCGCTTGTCGACGAGTCAGGGCGGGTTATTTTTGCCGATTTTAACCCATCGAAACGACCGTGGTCGGTGGCTTCAGATGACTCTAACAACAGCTGGAATCCGGCATACAATTGCTTTGATCGTCCGGGTAAAAAATGTAtctc                                                                                     <  2:2862545/150‑1 (MQ=255)
                                                                aaaTGGACCGCTTGTCGACGAGTCAGGGCGGGTTATTTTTGCCGATTTTAACCCATCGAAACGACCGTGGTCGGTGGCTTCAGATGACTCTAACAACAGCTGGAATCCGGCATACAATTGCTTTGATCGTCCGGGTAAAAAATGTAtctc                                                                                     >  2:3147808/1‑150 (MQ=255)
                                                                  ttGGCCCGCTTGGCGACGTGTCCGGGGGGGTTATTTTTGCCGATTTTACCCCATCGACACGACCGTGGTCGGTGGCTTCAGATGACTCTAACAACAGCTGGAATCCGGCATACAATTGCTTTGATCGTCCGGGTAAAAAATGTATCTCtt                                                                                   <  2:1494312/149‑1 (MQ=255)
                                                                   tGGACCGCTTGTCGACGAGTCTGGGCGGGTTATTTTTGCCTATTTTAACCCATCGAAACGACCGTGGTCGGTGGCTTCAGATGACTCTAACAACAGCTGGAATCCGGCATACAATTGCTTTGATCGTCCGGGTAAAAAATGTATCTCttt                                                                                  <  2:1059601/150‑1 (MQ=255)
                                                                     ggcccGCTTGTCGACGAGTCTGGGGGGTTTATTTTTGCCTATTTTAACCCATCGAAACGACCGTGGTCGGTGGCTTCAGATGACTCTAACAACAGCTGGAATCCGGCATACAATTGCTTTGATCGTCCGGGTAAAAAATGTATCTCtttt                                                                                 <  1:722954/147‑1 (MQ=255)
                                                                     ggACCGTTTGTCGACGTGTCTGGGGGGGTTATTTTTGCCGTTTTTAACCAATGGAAACGACCGTGGTCGGTGGCTTCAGATGACTCTAACAACAGCTGGAATCCGGCATACAATTGCTTTGATCGTCCGGGTAAAAAATGTATCTCtttt                                                                                 <  1:1386422/150‑1 (MQ=255)
                                                                              gTCGACGAGTCAGGGCGGGTTATTTTTGCCGATTTTAACCCATCGAAACGACCGTGGTCGGTGGCTTCAGATGACTCTAACAACAGCTGGAATCCGGCATACAATTGCTTTGATCGTCCGGGTAAAAAATGTATCTCTTTTACGCTacac                                                                        >  1:502911/1‑150 (MQ=255)
                                                                                cgacgaGTCAGGGCGGGTTATTTTTGCCGATTTTAACCCATCGAAACGACCGTGGTCGGTGGCTTCAGATGACTCTAACAACAGCTGGAATCCGGCATACAATTGCTTTGATCGTCCGGGTAAAAAATGTATCTCTTTTACGCTACACAt                                                                      >  1:106779/1‑150 (MQ=255)
                                                                                cgacgaGTCAGGGCGGGTTATTTTTGCCGATTTTAACCCATCGAAACGACCGTGGTCGGTGGCTTCAGATGACTCTAACAACAGCTGGAATCCGGCATACAATTGCTTTGATCGTCCGGGTAAAAAATGTATCTCTTTTAAGCTAAACAt                                                                      >  1:4009062/1‑150 (MQ=255)
                                                                                   cttgTCTGGGGGGGTTATTTTTGCCTTTTTTACCCCATCGACACGACCGTGGTCGGTGGCTTCAGATTACTCTAACAACAGCTGGAATCCGGCATACAATTGCTTTGATCGTCCGGGTAAAAAATGTATCTCTTTTACGCTACACATCaa                                                                   <  1:3743562/147‑1 (MQ=255)
                                                                                        cAGGGCGGGTTATTTTTGCCGATTTTAACCCATCGAAACGACCGTGGTCGGTGGCTTCAGATGACTCTAACAACAGCAGGAATCCGGCATACAATTGCTTTGATCGTCCGGGTAAAAAATGTATCTCTTTTACGCTACACATCAACGGCa                                                              >  1:153990/1‑150 (MQ=255)
                                                                                          ggggggTTTTTTTTTGCCTATTTTACCCTTTGGAAACGACCGTGGTCGGTGGCTTCAGATGTCTCTAACAACAGGTGGAATCCGGCATAGTTTTGCTTTGATCGTCCGGGTAAAAAATGTATCTCTTTTACGCTAGAGTTCATCGGCatt                                                            <  2:2904270/150‑3 (MQ=255)
                                                                                              gggTTATTTTTGCCTATTTTAACCCATCGAAACGACCGTGGTCGGTGGCTTCAGATGACTCTAACAACAGCTGGAATCCGGCATACAATTGCTTTGATCGTCCGGGTAAAAAATGTATCTCTTTTACGCTACACATCAACGGCAAAGATc                                                        <  1:3192893/150‑1 (MQ=255)
                                                                                               gtttATTTTTGCCGATTTAAACCCAACGAAACGACCGTGGTCGGTGGCTTCAGATGACTCTAACAACAGCTGGAATCCGGCATACAATTGCTTTGATCGTCCGGGTAAAAAATGTATCTCTTTTACGCTACACATCAACGGCAAAGATCa                                                       <  1:3590168/148‑1 (MQ=255)
                                                                                               ggTTATTTTTGCCGATTTTAACCCATCGAAACGACCGTGGTCGGTGGCTTCAGATGACTCTAACAACAGCTGGAATCCGGCATACAATTGCTTTGATCGTCCGGGTAAAAAATGTATCTCTTTTACGCTACACATCAACGGCAAAGATCa                                                       >  1:532043/1‑150 (MQ=255)
                                                                                               ggTTATTTTTGCCGATTTTAACCCATCGAAACGACCGTGGTCGGTGGCTACAGATGACTCTAACAACAGCTGGAATCCGGCATACAATTGCTTTGATCGTCCGGGTAAAAAATGTATCTCTTTTACGCTACACATCAACGGCAAAGATCa                                                       <  2:303936/150‑1 (MQ=255)
                                                                                                tttATTTTTGCCGATTTTAACCCATCGAAACGACCGTGGTCGGTGGCTTCAGATGACTCTAACAACAGCTGGAATCCGGCATACAATTGCTTTGATCGTCCGGGTAAAAAATGTATCTCTTTTACGCTACACATCAACGGCAAAGATCAc                                                      <  2:2176992/149‑1 (MQ=255)
                                                                                                     ttttGCCGATTTTAACCCATCGAAACGACCGTGGTCGGTGGCTTCAGATGACTCTAACAACAGCTGGAATCCGGCATACAATTGCTTTGATCGTCCGGGTAAAAAATGTATCTCTTTTACGCTACACATCAACGGCAAAGATCACGAttt                                                 <  2:2975084/150‑1 (MQ=255)
                                                                                                           cGATTTTAACCCATCGAAACGACCGTGGTCGGTGGCTTCAGATGACTCTAACAACAGCTGGAATCCGGCATACAATTGCTTTGATCGTCCGGGTAAAAAATGTATCTCTTTTACGCTACACATCAACGGCAAAGATCACGATTTGTTAGc                                           <  2:2506939/150‑1 (MQ=255)
                                                                                                             aTTTTAACCCATCGAAACGACCGTGGTCGGTGGCTTCAGATGACTCTAACAACAGCTGGAATCCGGCATACAATTGCTTTGATCGTCCGGGTAAAAAATGTATCTCTTTTACGCTACACATCAACGGCAAAGATCACGATTTGTTAGCgg                                         <  1:3769195/150‑1 (MQ=255)
                                                                                                                ttAACCCATCGAAACGACCGTGGTCGGTGGCTTCAGATGACTCTAACAACAGCTGGAATCCGGCATACAATTGCTTTGATCGTCCGGGTAAAAAATGTATCTCTTTTACGCTACACATCAACGGCAAAGATCACGATTTGTTAGCGGTgg                                      <  2:907669/150‑1 (MQ=255)
                                                                                                                 taccccATCTAAACGACCGTGGTCGGTGGCTTCAGATGACTCTAACAACAGCTGGAATCCGGCATACAATTGCTTTGATCGTCCGGGTAAAAAATGTATCTCTTTTACGCTACACATCAACGGCAAAGATCACGATTTGTTAGCGGTGGa                                     <  1:2356184/147‑1 (MQ=255)
                                                                                                                 tAACCCATCGACACGACCGTGGTCGGTGGCTTCAGATGACTCTAACAACAGCTGGAATCCGGCATACAATTGCTTTGATCGTCCGGGTAAAAAATGTATCTCTTTTACGCTACACATCAACGGCAAAGATCACGATTTGTTAGCGGTGGa                                     <  1:182284/150‑1 (MQ=255)
                                                                                                                      cATCGAAACGACCGTGGTCGGTGGCTTCAGATGACTCTAACAACAGCTGGAATCCGGCATACAATTGCTTTGATCGTCCGGGTAAAAAATGTATCTCTTTTACGCTACACATCAACGGCAAAGATCACGATTTGTTAGCGGTGGATaaaa                                >  1:1376400/1‑150 (MQ=255)
                                                                                                                       atctacaCGACCGTGGTCGGTGGCTTCAGATGACTCTAACAACAGCTGGAATCCGGCATACAATTGCTTTGATCGTCCGGGTAAAAAATGTATCTCTTTTACGCTACACATCAACGGCAAAGATCACGATTTGTTAGCGGTGGATAAAAt                               <  1:2651035/144‑1 (MQ=255)
                                                                                                                       aTCGAAACGACCGTGGTCGGTGGCTTCAGATGACTCTAACAACAGCTGGAATCCGGCATACAATTGCTTTGATCGTCCGGGTAAAAAATGTATCTCTTTTACGCTACACATCAACGGCAAAGATCACGATTTGTTAGCGGTGGATAAAAt                               >  2:2852724/1‑150 (MQ=255)
                                                                                                                         cGAAACGACCGTGGTCGGTGGCTTCAGATGACTCTAACAACAGCTGGAATCCGGCATACAATTGCTTTGATCGTCCGGGTAAAAAATGTATCTCTTTTACGCTACACATCAACGGCAAAGATCACGATTTGTTAGCGGTGGATAAAATTc                             >  1:2119050/1‑150 (MQ=255)
                                                                                                                          gAAACGACCGTGGTCGGTGGCTTCAGATTACTCTAACAACAGCTGGAATCCGGCATACAATTGCTTTGATCGTCCGGGTAAAAAATGTATCTCTTTTACGCTACACATCAACGGCAAAGATCACGATTTGTTAGCGGTGGATAAAATTCa                            >  1:16653/1‑150 (MQ=255)
                                                                                                                          gAAACGACCGTGGTCGGTGGCTTCAGATGACTCTAACAACAGCTGGAATCCGGCATACAATTGCTTTGATCGTCCGGGTAAAAAATGTATCTCTTTTACGCTACACATCAACGGCAAAGATCACGATTTGTTAGCGGTGGATAAAATTCa                            >  2:210181/1‑150 (MQ=255)
                                                                                                                            aaCGACCGTGGTCGGTGGCTTCAGATGACTCTAACAACAGCTGGAATCCGGCATACAATTGCTTTGATCGTCCGGGTAAAAAATGTATCTCTTTTACGCTACACATCAACGGCAAAGATCACGATTTGTTAGCGGTGGATAAAATTCATg                          <  1:3934698/150‑1 (MQ=255)
                                                                                                                               gACCGTGGTCGGTGGCTTCAGATGACTCTAACAACAGCTGGAATCCGGCATACAATTGCTTTGATCGTCCGGGTAAAAAATGTATCTCTTTTACGCTACACATCAACGGCAAAGATCACGATTTGTTAGCGGTGGATAAAATTCATGTCg                       >  2:556884/1‑150 (MQ=255)
                                                                                                                               gACCGTGGTCGGTGGCTTCAGATGACTCTAACAACAGCTGGAATCCGGCATACAATTGCTTTGATCGTCCGGGTAAAAAATGTATCTCTTTTACGCTACACATCAACGGCAAAGATCACGATTTGTTAGCGGTGGATAAAATTCATGTCg                       >  2:2371036/1‑150 (MQ=255)
                                                                                                                               gACCGTGGTCGGTGGCTTCAGATGACTCTAACAACAGCTGGAATCCGGCATACAATTGCTTTGATCGTCCGGGTAAAAAATGTATCTCTTTTACGCTACACATCAACGGCAAAGATCACGATTTGTTAGCGGTGGATAAAATTCATGTCg                       >  1:2243888/1‑150 (MQ=255)
                                                                                                                                   gTGGTCGGTGGCTTCAGATTACTCTAACAACAGCTGGAATCCGGCATACAATTGCTTTGATCGTCCGGGTAAAAAATGTATCTCTTTTACGCTACACATCAACGGCAAAGATCACGATTTGTTAGCGGTGGATAAAATTCATGTCGAttt                   <  1:1315956/150‑1 (MQ=255)
                                                                                                                                        cGGTGGCTTCAGATGACTCTAACAACAGCTGGAATCCGGCATACAATTGCTTTGATCGTCCGGGTAAAAAATGTATCTCTTTTACGCTACACATCAACGGCAAAGATCACGATTTGTTAGCGGTGGATAAAATTCATGTCGATTTAAACt              >  2:3836976/1‑150 (MQ=255)
                                                                                                                                        cGGTGGCTTCAGATGACTCTAACAACAGCTGGAATCCGGCATACAATTGCTTTGATCGTCCGGGTAAAAAATGTATCTCTTTTACGCTACACATCAACGGCAAAGATCACGATTTGTTAGCGGTGGATAAAATTCATGTCGATTTAAACt              <  2:2116691/150‑1 (MQ=255)
                                                                                                                                         ggTGGCTTCAGATGACTCTAACAACAGCTGGAATCCGGCATACAATTGCTTTGATCGTCCGGGTAAAAAATGTATCTCTTTTACGCTACACATCAACGGCAAAGATCACGATTTGTTAGCGGTGGATAAAATTCATGTCGATTTAAACTg             <  1:1657206/150‑1 (MQ=255)
                                                                                                                                           tGGCTTCAGATGACTCTAACAACAGCTGGAATCCGGCATACAATTGCTTTGATCGTCCGGGTAAAAAATGTATCTCTTTTACGCTACACATCAACGGCAAAGATCACGATTTGTTAGCGGTGGATAAAATTCATGTCGATTTAAACTGGc           <  2:3320722/150‑1 (MQ=255)
                                                                                                                                            ggCTTCAGATGACTCTAACAACAGCTGGAATCCGGCATACAATTGCTTTGATCGTCCGGGTAAAAAATGTATCTCTTTTACGCTACACATCAACGGCAAAGATCACGATTTGTTAGCGGTGGATAAAATTCATGTCGATTTAAACTGGCg          >  1:4301106/1‑150 (MQ=255)
                                                                                                                                              cTTCAGATGACTCTAACAACAGCTGGAATCCGGCATACAATTGCTTTGATCGTCCGGGTAAAAAATGTATCTCTTTTACGCTACACATCAACGGCAAAGATCACGATTTGTTAGCGGTGGATAAAATTCATGTCGATTTAAACTGGCGat        <  1:4053230/150‑1 (MQ=255)
                                                                                                                                              cTTCAGATGACTCTAACAACAGCTGGAATCCGGCATACAATTGCTTTGATCGTCCGGGTAAAAAATGTATCTCTTTTACGCTACACATCAACGGCAAAGATCACGATTTGTTAGCGGTGGATAAAATTCATGTCGATTTAAACTGGCGat        <  1:990543/150‑1 (MQ=255)
                                                                                                                                               ttCAGATGACTCTAACAACAGCTGGAATCCGGCATACAATTGCTTTGATCGTCCGGGTAAAAAATGTATCTCTTTTACGCTACACATCAACGGCAAAGATCACGATTTGTTAGCGGTGGATAAAATTCATGTCGATTTAAACTGGCGata       >  1:3032985/1‑150 (MQ=255)
                                                                                                                                               ttCAGATGACTCTAACAACAGCTGGAATCCGGCATACAATTGCTTTGATCGTCCGGGTAAAAAATGTATCTCTTTTACGCTACACATCAACGGCAAAGATCACGATTTGTTAGCGGTGGATAAAATTCATGTCGATTTAAACTGGCGata       <  2:917694/150‑1 (MQ=255)
                                                                                                                                                    aTGACTCTAACAACAGCTGGAATCCGGCATACAATTGCTTTGATCGTCCGGGTAAAAAATGTATCTCTTTTACGCTACACATCAACGGCAAAGATCACGATTTGTTAGCGGTGGATAAAATTCATGTCGATTTAAACTGGCGATATCTGa  >  1:3105138/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                   
ATAAATATTATCTTTTTTCGCGGAGGGCGAATGATAAAATCGTTGTTTATCATCAGGAACAAGTAAATGGACCGCTTGTCGACGAGTCAGGGCGGGTTATTTTTGCCGATTTTAACCCATCGAAACGACCGTGGTCGGTGGCTTCAGATGACTCTAACAACAGCTGGAATCCGGCATACAATTGCTTTGATCGTCCGGGTAAAAAATGTATCTCTTTTACGCTACACATCAACGGCAAAGATCACGATTTGTTAGCGGTGGATAAAATTCATGTCGATTTAAACTGGCGATATCTGA  >  CP009273/868748‑869044

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