Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,461,652 G→T 12.3% pseudogene (2004/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,6520GT12.3% 1169.5 / 105.1 334pseudogene (2004/2513 nt)ydbApseudogene, autotransporter homolog; interrupted by IS2 and IS30
Reads supporting (aligned to +/- strand):  ref base G (147/145);  new base T (19/23);  total (166/168)
Fisher's exact test for biased strand distribution p-value = 6.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 336 total aligned reads displayed.

CGCTGGCAGGCAGTATGCAGGTTGGCGATTTCTCGACCGGGGTAGATCTTAATGGCAACAATAACAGCGTGACGCTGGCGGCAAAAGATCTAAAAGTGGTCGGGCAGAAAGCGACGGGCATAAACGTTTCTGGCGATGCGAATACAGTGAATATCACTGGTAACGTTCTGGTTGATAAGGATAAAACCGCAGACAATGCGGCGGAATATTTCTTCGATCCATCCGTGGGTATCAACGTTTACGGCAGTGATAATAACGTGACGCTGGATGGAAAGTTAACTGTTGTATCAGACAG  >  CP009273/1461504‑1461798
                                                                                                                                                    |                                                                                                                                                  
  cTGGCAGTCAGTATGCAGGTTGGCGATTTCTCGACCGGGGTAGATCTTAATGGCAACAATAACAGCGTGACGCTGGCGGCAAAAGATCTAAAAGTGGTCGGGCAGAAAGCGACGGGCATAAACGTTTCTGGCGATGCGAATACAGTTaa                                                                                                                                                  <  2:2971540/149‑1 (MQ=255)
       aggcagTATGCAGGTTGGCGATTTCTCGACCGGGGTAGATCTTAATGGCAACAATAACAGCGTGACGCTGGCGGCAAAAGATCTAAAAGTGGTCGGGCAGAAAGCGACGGGCATAAACGTTTCTGGCGATGCGAATACAGTGAATATCa                                                                                                                                             >  1:2119256/1‑149 (MQ=255)
        ggcagTATGCAGGTTGGCGATTTCTCGACCGGGGTAGATCTTAATGGCAACAATAACAGCGTGACGCTGGCGGCAAAAGATCTAAAAGTGGTCGGGCAGAAAGCGACGGGCATAAACGTTTCTGGCGATGCGAATACAGTGAATATCAc                                                                                                                                            <  2:3734488/149‑1 (MQ=255)
         gcagTATGCAGGTTGGCGATTTCTCGACCGGGGTAGATCTTAATGGCAACAATAACAGCGTGACGCTGGCGGCAAAAGATCTAAAAGTGGTCGGGCAGAAAGCGACGGGCATAAACGTTTCTGGCGATGCGAATACAGTGAATATCACt                                                                                                                                           <  2:4923526/149‑1 (MQ=255)
            gTATGCAGGTTGGCGATTTCTCGACCGGGGTAGATCTTAATGGCAACAATAACAGCGTGACGCTGGCGGCACAAGATCTAAAAGTGGTCGGGCAGAAAGCGACGGGCATAAACGTTTCTGGCGATGCGAATACAGTTAATATCACTGGt                                                                                                                                        >  2:3637617/1‑149 (MQ=255)
             tATGCAGGTTGGCGATTTCTCGACCGGGGTAGATCTTAATGGCAACAATAACAGCGTGACGCTGGCGGCAAAAGATCTAAAAGTGGTCGGGCAGAAAGCGACGGGCATAAACGTTTCTGGCGATGCGAATACAGTGAATATCACTGGTa                                                                                                                                       <  2:4519728/149‑1 (MQ=255)
             tATGCAGGTTGGCGATTTCTCGACCGGGGTAGATCTTAATGGCAACAATAACAGCGTGACGCGGGCGGCAAAAGATCTAAAAGTGGTCGGGCAGAAAGCGACGGGCATAAACGTTTCTGGCGATGCGAATACAGTTAATATCACTGGTa                                                                                                                                       <  1:929970/149‑1 (MQ=255)
               tGCAGGTTGGCGATTTCTCGACCGGGGTAGATCTTAATGGCAACAATAACAGCGTGACGCTGGCGGCAAAAGATCTAAAAGTGGTCGGGCAGAAAGCGACGGGCATAAACGTTTCTGGCGATGCGAATACAGTGAATATCACTGGTAAc                                                                                                                                     <  2:4194582/149‑1 (MQ=255)
                 cAGGTTGGCGATTTCTCGACCGGGGTATATCTTAATGGCAACAATAACAGCGTGACGCTGGCGGCAAAAGATCTAAAAGTGGTCGGGCAGAAAGCGACGGGCATAAACGTTTCTGGCGATGCGAATACAGTGAATATCACTGGTAACGt                                                                                                                                   <  1:3961144/149‑1 (MQ=255)
                  aGGTTGGCGATTTCTCGACCGGGGTAGATCTTAATGGCAACAATAACAGCGTGACGCTGGCGGCAAAAGATCTAAAAGTGGTCGGGCAGAAAGCGACGGGCATAAACGTTTCTGGCGATGCGAATACAGTGAATATCACTGGTAACGtt                                                                                                                                  <  2:247244/149‑1 (MQ=255)
                  aGGGTGGCGATTTCTCGACCGGGGTAGATCTTAATGGCAACAATAACAGCGTGACGCTGGCGGCAAAAGATCTAAAAGTGGTCGGGCAGAAAGCGACGGGCATAAACGTTTCTGGCGATGCGAATACAGTGAATATCACTGGTAACGtt                                                                                                                                  <  2:3112544/149‑1 (MQ=255)
                   ggTTGGCGATTTCTCGACCGGGGTAGATCTTAATGGCAACAATAACAGCGTGACGCTGGCGGCAAAAGATCTAAAAGTGGTCGGGCAGAAAGCGACGGGCATAAACGTTTCTGGCGATGCGAATACAGTGAATATCACTGGTAACGTTc                                                                                                                                 >  2:1714837/1‑149 (MQ=255)
                    gTTGGCGATTTCTCGACCGGGGTAGATCTTAATGGCAACAATAACAGCGTGACGCTGGCGGCAAAAGATCTAAAAGTGGTCGGGCAGAAAGCGACGGGCATAAACGTTTCTGGCGATGCGAATACAGTGAATATCACTGGTAACGTTCt                                                                                                                                <  1:3593311/149‑1 (MQ=255)
                    gTTGGCGATTTCTCGACCGGGGTAGATCTTAATGGCAACAATAACAGCGTGACGCTGGCGGCAAAAGATCTAAAAGTGGTCGGGCAGAAAGCGACGGGCATAAACGTTTCTGGCGATGCGAATACAGTGAATATCACTGGTAACGTTCt                                                                                                                                <  1:2807286/149‑1 (MQ=255)
                      tGGCGATTTCTCGACCGGGGTAGATCTTAATGGCAACAATAACAGCGTGACGCTGGCGGCAAAAGATCTAAAAGTGGTCGGGCAGAAAGCGACGGGCATAAACGTTTCTGGCGATGCGAATACAGTGAATATCACTGGTAACGTTCTgg                                                                                                                              <  2:2725435/149‑1 (MQ=255)
                           aTTTCTCGACCGGGGTAGATCTTAATGGCAACAATAACAGCGTGACGCTGGCGGCAAAAGATCTAAAAGTGGTCGGGCAGAAAGCGACGGGCATAAACGTTTCTGGCGATGCGAATACAGTGAATATCACTGGTAACGTTCTGGTTGAt                                                                                                                         >  2:4761688/1‑149 (MQ=255)
                           aTTTCTCGACCGGGGTAGATCTTAATGGCAACAATAACAGCGTGACGCTGGCGGCAAAAGATCTAAAAGTGGTCGGGCAGAAAGCGACGGGCATAAACGTTTCTGGCGATGCGAATACAGTGAATATCACTGGTAACGTTCTGGTTGAt                                                                                                                         >  1:1604319/1‑149 (MQ=255)
                           aTTTCTCGACCGGGGTAGATCTTAATGGCAACAATAACAGCGTGACGCTGGCGGCAAAAGATCTAAAAGTGGTCGGGCAGAAAGCGACGGGCATAAACGTTTCTGGCGATGCGAATACAGTGAATATCACTGGTAACGTTCTGGTTGAt                                                                                                                         >  1:3155493/1‑149 (MQ=255)
                            tttCTCGACCGGGGTAGATCTTAATGGCAACAATAACAGCGTGACGCTGGCGGCAAAAGATCTAAAAGTGGTCGGGCAGAAAGCGACGGGCATAAACGTTTCTGGCGATGCGAATACAGTGAATATCACTGGTAACGTTCTGGTTGATa                                                                                                                        <  1:864351/149‑1 (MQ=255)
                                tcGACCGGGGTAGATCTTAATGGCAACAATAACAGCGTGACGCTGGCGGCAAAAGATCTAAAAGTGGTCGGGCAGAAAGCGACGGGCATAAACGTTTCTGGCGATGCGAATACAGTGAATATCACTGGTAACGTTCTGGTTGATAAGGa                                                                                                                    >  2:2867867/1‑149 (MQ=255)
                                tcGACCGGGGTAGATCTTAATGGCAACAATAACAGCGTGACGCTGGCGGCAAAAGATCTAAAAGTGGTCGGGCAGAAAGCGACGGGCATAAACGTTTCTGGCGATGCGAATACAGTGAATATCACTGGTAACGTTCTGGTTGATAAGGa                                                                                                                    >  2:1998091/1‑149 (MQ=255)
                                  gACCGGGGTAGATCTTAATGGCAACAATAACAGCGTGACGCTGGCGGCAAAAGATCTAAAAGTGGTCGGGCAGAAAGCGACGGGCATAAACGTTTCTGGCGATGCGAATACAGTTAATATCACTGGTAACGTTCTGGTTGATAAGGATa                                                                                                                  <  1:1277143/149‑1 (MQ=255)
                                  gACCGGGGTAGATCTTAATGGCAACAATAACAGCGTGACGCTGGCGGCAAAAGATCTAAAAGTGGTCGGGCAGAAAGCGACGGGCATAAACGTTTCTGGCGATGCGAATACAGTGAATATCACTGGTAACGTTCTGGTTGATAAGGATa                                                                                                                  >  2:3675094/1‑149 (MQ=255)
                                         gTAGATCTTAATGGCAACAATAACAGCGTGACGCTGGCGGCAAAAGATCTAAAAGTGGTCGGGCAGAAAGCGACGGGCATAAACGTTTCTGGCGATGCGAATACAGTGAATATCACTGGTAACGTTCTGGTTGATAAGGATAAAACCGc                                                                                                           >  2:1520268/1‑149 (MQ=255)
                                            gATCTTAATGGCAACAATAACAGCGTGACGCTGGCGGCAAAAGATCTAAAAGTGGTCGGGCTGAAAGCGACGGGCATAAACGTTTCTGGCGATGCGAATACAGTGAATATCACTGGTAACGTTCTGGTTGATAAGGATAAAACCGCAGa                                                                                                        >  1:4212687/1‑149 (MQ=255)
                                              tCTTAATGGCAACAATAACAGCGTGACGCTGGCGGCAAAAGATCTAAAAGTGGTCGGGCAGAAAGCGACGGGCATAAACGTTTCTGGCGATGCGAATACAGTGAATATCACTGGTAACGTTCTGGTTGATAAGGATAAAACCGCAGACa                                                                                                      >  1:2327958/1‑149 (MQ=255)
                                              tCTTAATGGCAACAATAACAGCGTGACGCTGGCGGCAAAAGATCTAAAAGTGGTCGGGCAGAAAGCGACGGGCATAAACGTTTCTGGCGATGCGAATACAGTGAATATCACTGGTAACGTTCTGGTTGATAAGGATAAAACCGCAGACa                                                                                                      >  1:3399110/1‑149 (MQ=255)
                                                ttAATGGCAACAATAACAGCGTGACGCTGGCGGCAAAAGATCTAAAAGTGGTCGGGCAGAAAGCGACGGGCATAAACGTTTCTGGCGATGCGAATACAGTGAATATCACTGGTAACGTTCTGGTTGATAAGGATAAAACCGCAGACAAt                                                                                                    >  1:367072/1‑149 (MQ=255)
                                                  aaTGGCAACAATAACAGCGTGACGCTGGCGGCAAAAGATCTAAAAGTGGTCGGGCAGAAAGCGACGGGCATAAACGTTTCTGGCGATGCGAATACAGTTAATATCACTGGTAACGTTCTGGTTGATAAGGATAAAACCGCAGACAATgc                                                                                                  >  2:1133321/1‑149 (MQ=255)
                                                  aaTGGCAACAATAACAGCGTGACGCTGGCGGCAAAAGATCTAAAAGTGGTCGGGCAGAAAGCGACGGGCATAAACGTTTCTGGCGATGCGAATACAGTGAATATCACTGGTAACGTTCTGGTTGATAAGGATAAAACCGCAGACAATgc                                                                                                  >  1:2404774/1‑149 (MQ=255)
                                                  aaTGGCAACAATAACAGCGTGACGCTGGCGGCAAAAGATCTAAAAGTGGTCGGGCAGAAAGCGACGGGCATAAACGTTTCTGGCGATGCGAATACAGTGAATATCACTGGTAACGTTCTGGTTGATAAGGATAAAACCGCAGACAATgc                                                                                                  >  2:2113302/1‑149 (MQ=255)
                                                    tGGCAACAATAACAGCGTGACGCTGGCGGCAAAAGATCTAAAAGTGGTCGGGCAGAAAGCGACGGGCATAAACGTTTCTGGCGATGCGAATACAGTGAATATCACTGGTAACGTTCTGGTTGATAAGGATAAAACCGCAGACAATgcgg                                                                                                >  2:1239533/1‑149 (MQ=255)
                                                          caaTAACAGCGTGACGCTGGCGGCAAAAGATCTAAAAGTGGTCGGGCAGAAAGCGACGGGCATAAACGTTTCTGGCGATGCGAATACAGTTAATATCACTGGTAACGTTCTGGTTGATAAGGATAAAACCGCAGACAATGCGGCGGAat                                                                                          <  2:711814/149‑1 (MQ=255)
                                                          caaTAACAGCGTGACGCTGGCGGCAAAAGATCTAAAAGTGGTCGGGCAGAAAGCGACGGGCATAAACGTTTCTGGCGATGCGAATACAGTGAATATCACTGGTAACGTTCTGGTTGATAAGGATAAAACCGCAGACAATGCGGCGGAat                                                                                          <  2:835826/149‑1 (MQ=255)
                                                                 aGCGTGACGCTGGCGGCAAAAGATCTAAAAGTGGTCGGGCAGAAAGCGACGGGCATAAACGTTTCTGGCGATGCGAATACAGTGAATATCACTGGTAACGTTCTGGTTGATAAGGATAAAACCGCAGACAATGCGGCGGAATATttctt                                                                                   <  2:390197/149‑1 (MQ=255)
                                                                  gCGTGACGCTGGCGGCAAAAGATCTAAATGTGGTCGGGCAGAAAGCGACGGGAATAAACGTTTCTGGCGATGCGAATACAGTGAATATCACTGGTAACGTTCTGGTTGATAAGGATAAAACCGCAGACAATGCGGCGGAATATttcttc                                                                                  <  1:3267120/149‑1 (MQ=255)
                                                                  gCGTGACGCTGGCGGCAAAAGATCTAAAAGTGGTCGGGCAGAAAGCGACGGGCATAAACGTTTCTGGCGATGCGAATACAGTGAATATCACTGGTAACGTTCTGGTTGATAAGGATAAAACCGCAGACAATGCGGCGGAATATttcttc                                                                                  >  1:4251331/1‑149 (MQ=255)
                                                                  gCGTGACGCTGGCGGCAAAAGATCTAAAAGTGGTCGGGCAGAAAGCGACGGGCATAAACGTTTCTGGCGATGCGAATACAGTGAATATCACTGGTAACGTTCTGGTTGATAAGGATAAAACCGCAGACAATGCGGCGGAATATttcttc                                                                                  <  2:2329934/149‑1 (MQ=255)
                                                                  gCGTGACGCTGGCGGCAAAAGATCTAAAAGTGGTCGGGCAGAAAGCGACGGGCATAAACGTTTCTGGCGATGCGAATACAGTGAATATCACTGGTAACGTTCTGGTTGATAAGGATAAAACCGCAGACAATGCGGCGGAATATttcttc                                                                                  >  1:1330556/1‑149 (MQ=255)
                                                                     tGACGCTGGCGGCAAAAGATCTAAAAGTGGTCGGGCAGAAAGCGACGGGCATAAACGTTTCTGGCGATGCGAATACAGTGAATATCACTGGTAACGTTCTGGTTGATAAGGATAAAACCGCAGACAATGCGGCGGAATATTTCTTCGat                                                                               >  2:2157163/1‑149 (MQ=255)
                                                                     tGACGCTGGCGGCAAAAGATCTAAAAGTGGTCGGGCAGAAAGCGACGGGCATAAACGTTTCTGGCGATGCGAATACAGTGAATATCACTGGTAACGTTCTGGTTGATAAGGATAAAACCGCAGACAATGCGGCGGAATATTTCTTCGat                                                                               >  2:1411532/1‑149 (MQ=255)
                                                                      gACGCTGGCGGCAAAAGATCTAAAAGTGGTCGGGCAGAAAGCGACGGGCATAAACGTTTCTGGCGATGCGAATACAGTGAATATCACTGGTAACGTTCTGGTTGATAAGGATAAAACCGCAGACAATGCGGCGGAATATTTCTTCGatc                                                                              <  2:1397779/149‑1 (MQ=255)
                                                                      gACGCTGGCGGCAAAAGATCTAAAAGTGGTCGGGCAGAAAGCGACGGGCATAAACGTTTCTGGCGATGCGAATACAGTGAATATCACTGGTAACGTTCTGGTTGATAAGGATAAAACCGCAGACAATGCGGCGGAATATTTCTTCGatc                                                                              >  1:5145204/1‑149 (MQ=255)
                                                                       aCGCTGGCGGCAAAAGATCTAAAAGTGGTCGGGCAGAAAGCGACGGGCATAAACGTTTCTGGCGATGCGAATACAGTGAATATCACTGGTAACGTTCTGGTTGATAAGGATAAAACCGCAGACAATGCGGCGGAATATTTCTTCGatcc                                                                             <  1:4469215/149‑1 (MQ=255)
                                                                       aCGCTGGCGGCAAAAGATCTAAAAGTGGTCGGGCAGAAAGCGACGGGCATAAACGTTTCTGGCGATGCGAATACAGTGAATATCACTGGTAACGTTCTGGTTGATAAGGATAAAACCGCAGACAATGCGGCGGAATATTTCTTCGatcc                                                                             <  1:1197707/149‑1 (MQ=255)
                                                                           tGGCGGCAAAAGATCTAAAAGTGGTCGGGCAGAAAGCGACGGGCATAAACGTTTCTGGCGATGCGAATACAGTGAATATCACTGGTAACGTTCTGGTTGATAAGGATAAAACCGCAGACAATGCGGCGGAATATTTCTTCGatccatcc                                                                         <  1:111196/149‑1 (MQ=255)
                                                                           tGGCGGCAAAAGATCTAAAAGTGGTCGGGCAGAAAGCGACGGGCATAAACGTTTCTGGCGATGCGAATACAGTGAATATCACTGGTAACGTTCTGGTTGATAAGGATAAAACCGCAGACAATGCGGCGGAATATTTCTTCGatccatcc                                                                         <  1:2858108/149‑1 (MQ=255)
                                                                             gcggcAAAAGATCTAAAAGTGGTCGGGCAGAAAGCGACGGGCATAAACGTTTCTGGCGATGCGAATACAGTTAATATCACTGGTAACGTTCTGGTTGATAAGGATAAAACCGCAGACAATGCGGCGGAATATTTCTTCGATCCATCCGt                                                                       >  2:5115802/1‑149 (MQ=255)
                                                                                 cAAAAGATCTAAAAGTGGTCGGGCAGAAAGCGACGGGCATAAACGTTTCTGGCGATGCGAATACAGTGAATATCACTGGTAACGTTCTGGTTGATAAGGATAAAACCGCAGACAATGCGGCGGAATATTTCTTCGATCCATCCGTGGGt                                                                   <  1:4100210/149‑1 (MQ=255)
                                                                                 cAAAAGATCTAAAAGTGGTCGGGCAGAAAGCGACGGGCATAAACGTTTCTGGCGATGCGAATACAGTGAATATCACTGGTAACGTTCTGGTTGATAAGGATAAAACCGCAGACAATGCGGCGGAATATTTCTTCGATCCATCCGTGGGt                                                                   <  1:2937040/149‑1 (MQ=255)
                                                                                 cAAAAGAACTAAAAGTGGTCGGGCAGAAAGCGACGGGCATAAACGTTTCTGGCGATGCGAATACAGTGAATATCACTGGTAACGTTCTGGTTGATAAGGATAAAACCGCAGACAATGCGGCGGAATATTTCTTCGATCCATCCGTGGGt                                                                   <  2:3235202/149‑1 (MQ=255)
                                                                                    aaGATCTAAAAGTGGTCGGGCAGAAAGCGACGGGCATAAACGTTTCTGGCGATGCGAATACAGTGAATATCACTGGTAACGTTCTGGTTGATAAGGATAAAACCGCAGACAATGCGGCGGAATATTTCTTCGATCCATCCGTGGGTATc                                                                >  1:4679163/1‑149 (MQ=255)
                                                                                      gATCTAAAAGTGGTCGGGCAGAAAGCGACGGGCATAAACGTTTCTGGCGATGCGAATACAGTGAATATCACTGGTAACGTTCTGGTTGATAAGGATAAAACCGCAGACAATGCGGCGGAATATTTCTTCGATCCATCCGTGGGTATCaa                                                              >  1:1045850/1‑149 (MQ=255)
                                                                                      gATCTAAAAGTGGTCGGGCAGAAAGCGACGGGCATAAACGTTTCTGGCGATGCGAATACAGTGAATATCACTGGTAACGTTCTGGTTGATAAGGATAAAACCGCAGACAATGCGGCGGAATATTTCTTCGATCCATCCGTGGGTATCaa                                                              <  2:5001309/149‑1 (MQ=255)
                                                                                         cTAAAAGGGGTCGGGCAGAAAGCGACGGGCATAAACGTTTCTGGCGATGCGAATACAGTGAATATCACTGGTAACGTTCTGGTTGATAAGGATAAAACCGCAGACAATGCGGCGGAATATTTCTTCGATCCATCCGTGGGTATCAACGt                                                           <  2:1623259/149‑1 (MQ=255)
                                                                                          tAAAAGTGGTCGGGCAGAAAGCGACGGGCATAAACGTTTCTGGCGATGCGAATACAGTGAATATCACTGGTAACGTTCTGGTTGATAAGGATAAAACCGCAGACAATGCGGCGGAATATTTCTTCGATCCATCCGTGGGTATCAACGtt                                                          <  1:3007752/149‑1 (MQ=255)
                                                                                             aaGTGGTCGGGCAGAAAGCGACGGGCATAAACGTTTCTGGCGATGCGAATACAGTGAATATCACTGGTAACGTTCTGGTTGATAAGGATAAAACCGCAGACAATGCGGCGGAATATTTCTTCGATCCATCCGTGGGTATCAACGTTTAc                                                       >  2:3495670/1‑149 (MQ=255)
                                                                                             aaGTGGTCGGGCAGAAAGCGACGGGCATAAACGTTTCTGGCGATGCGAATACAGTGAATATCACTGGTAACGTTCTGGTTGATAAGGATAAAACCGCAGACAATGCGGCGGAATATTTCTTCGATCCATCCGTGGGTATCAACGTTTAc                                                       <  1:3027207/149‑1 (MQ=255)
                                                                                             aaGTGGTCGGGCAGAAAGCGACGGGCATAAACGTTTCTGGCGATGCGAATACAGTGAATATCACTGGTAACGTTCTGGTTGATAAGGATAAAACCGCAGACAATGCGGCGGAATATTTCTTCGATCCATCCGTGGGTATCAACGTTTAc                                                       <  2:1851358/149‑1 (MQ=255)
                                                                                              aGTGGTCGGGCAGAAAGCGACGGGCATAAACGTTTCTGGCGATGCGAATACAGTGAATATCACTGGTAACGTTCTGGTTGATAAGGATAAAACCGCAGACAATGCGGCGGAATATTTCTTCGATCCATCCGTGGGTATCAACGTTTACg                                                      >  1:56674/1‑149 (MQ=255)
                                                                                                tGGTCGGGCAGAAAGCGACGGGCATAAACGTTTCTGGCGATGCGAATACAGTGAATATCACTGGTAACGTTCTGGTTGATAAGGATAAAACCGCAGACAATGCGGCGGAATATTTCTTCGATCCATCCGTGGGTATCAACGTTTACGGc                                                    <  1:1475075/149‑1 (MQ=255)
                                                                                                  gTCGGGCAGAAAGCGACGGGCATAAACGTTTCTGGCGATGCGAATACAGTGAATATCACTGGTAACGTTCTGGTTGATAAGGATAAAACCGCAGACAATGCGGCGGAATATTTCTTCGATCCATCCGTGGGTATCAACGTTTACGGCAg                                                  >  1:4935965/1‑149 (MQ=255)
                                                                                                    cGGGCAGAAAGCGACGGGCATAAACGTTTCTGGCGATGCGAATACAGTGAATATCACTGGTAACGTTCTGGTTGATAAGGATAAAACCGCAGACAATGCGGCGGAATATTTCTTCGATCCATCCGTGGGTATCAACGTTTACGGCAGTg                                                >  2:2828027/1‑149 (MQ=255)
                                                                                                     gggCAGAAAGCGACGGGCATAAACGTTTCTGGCGATGCGAATACAGTGAATATCACTGGTAACGTTCTGGTTGATAAGGATAAAACCGCAGACAATGCGGCGGAATATTTCTTCGATCCATCCGTGGGTATCAACGTTTACGGCAGTGa                                               >  1:3348525/1‑149 (MQ=255)
                                                                                                     gggCAGAAAGCGACGGGCATAAACGTTTCTGGCGATGCGAATACAGTGAATATCACTGGTAACGTTCTGGTTGATAAGGATAAAACCGCAGACAATGCGGCGGAATATTTCTTCGATCCATCCGTGGGTATCAACGTTTACGGCAGTGa                                               <  1:4793130/149‑1 (MQ=255)
                                                                                                         aGAAAGCGACGGGCATAAACGTTTCTGGCGATGCGAATACAGTGAATATCACTGGTAACGTTCTGGTTGATAAGGATAAAACCGCAGACAATGCGGCGGAATATTTCTTCGATCCATCCGTGGGTATCAACGTTTACGGCAGTGataat                                           <  1:2239051/149‑1 (MQ=255)
                                                                                                         aGAAAGCGACGGGCATAAACGTTTCTGGCGATGCGAATACAGTGAATATCACTGGTAACGTTCTGGTTGATAAGGATAAAACCGCAGACAATGCGGCGGAATATTTCTTCGATCCATCCGTGGGTATCAACGTTTACGGCAGTGataat                                           >  2:1152059/1‑149 (MQ=255)
                                                                                                           aaaGCGACGGGCATAAACGTTTCTGGCGATGCGAATACAGTGAATATCACTGGTAACGTTCTGGTTGATAAGGATAAAACCGCAGACAATGCGGCGGAATATTTCTTCGATCCATCCGTGGGTATCAACGTTTACGGCAGTGAtaataa                                         <  2:576019/149‑1 (MQ=255)
                                                                                                           aaaGCGACGGGCATAAACGTTTCTGGCGATGCGAATACAGTGAATATCACTGGTAACGTTCTGGTTGATAAGGATAAAACCGCAGACAATGCGGCGGAATATTTCTTCGATCCATCCGTGGGTATCAACGTTTACGGCAGTGAtaataa                                         >  1:777225/1‑149 (MQ=255)
                                                                                                           aaaGCGACGGGCATAAACGTTTCTGGCGATGCGAATACAGTGAATATCACTGGTAACGTTCTGGTTGATAAGGATAAAACCGCAGACAATGCGGCGGAATATTTCTTCGATCCATCCGTGGGTATCAACGTTTACGGCAGTGAtaataa                                         >  2:4578997/1‑149 (MQ=255)
                                                                                                             aGCGACGGGCATAAACGTTTCTGGCGATGCGAATACAGTGAATATCACTGGTAACGTTCTGGTTGATAAGGATAAAACCGCAGACAATGCGGCGGAATATTTCTTCGATCCATCCGTGGGTATCAACGTTTACGGCAGTGATAATAACg                                       >  1:629831/1‑149 (MQ=255)
                                                                                                              gCGACGGGCATAAACGTTTCTGGCGATGCGAATACAGTGAATATCACTGGTAACGTTCTGGTTGATAAGGATAAAACCGCAGACAATGCGGCGGAATATTTCTTCGATCCATCCGTGGGTATCAACGTTTACGGCAGTGATAATAACGt                                      >  1:4129522/1‑149 (MQ=255)
                                                                                                               cGACGGGCATAAACGTTTCTGGCGATGCGAATACAGTGAATATCACTGGTAACGTTCTGGTTGATAAGGATAAAACCGCAGACAATGCGGCGGAATATTTCTTCGATCCATCCGTGGGTATCAACGTTTACGGCAGTGATAATAACGTg                                     <  1:709003/149‑1 (MQ=255)
                                                                                                                gACGGGCATAAACGTTTCTGGCGATGCGAATACAGTGAATATCACTGGTAACGTTCTGGTTGATAAGGATAAAACCGCAGACAATGCGGCGGAATATTTCTTCGATCCATCCGTGGGTATCAACGTTTACGGCAGTGATAATAACGTGa                                    <  2:622954/149‑1 (MQ=255)
                                                                                                                gACGGGCATAAACGTTTCTGGCGATGCGAATACAGTGAATATCACTGGTAACGTTCTGGTTGATAAGGATAAAACCGCAGACAATGCGGCGGAATATTTCTTCGATCCATCCGTGGGTATCAACGTTTACGGCAGTGATAATAACGTGa                                    <  2:5056447/149‑1 (MQ=255)
                                                                                                                 aCGGGCATAAACGTTTCTGGCGATGCGAATACAGTGAATATCACTGGTAACGTTCTGGTTGATAAGGATAAAACCGCAGACAATGCGGCGGAATATTTCTTCGATCCATCCGTGGGTATCAACGTTTACGGCAGTGATAATAACGTGAc                                   <  1:25059/149‑1 (MQ=255)
                                                                                                                   gggCATAAACGTTTCTGGCGATGCGAATACAGTGAATATCACTGGTAACGTTCTGGTTGATAAGGATAAAACCGCAGACAATGCGGCGGAATATTTCTTCGATCCATCCGTGGGTATCAACGTTTACGGCAGTGATAATAACGTGACGc                                 >  1:658115/1‑149 (MQ=255)
                                                                                                                   gggCATAAACGTTTCTGGCGATGCGAATACAGTGAATATCACTGGTAACGTTCTGGTTGATAAGGATAAAACCGCAGACAATGCGGCGGAATATTTCTTCGATCCATCCGTGGGTATCAACGTTTACGGCAGTGATAATAACGTGACGc                                 >  1:1343718/1‑149 (MQ=255)
                                                                                                                    ggCATAAACGTTTCTGGCGATGCGAATACAGTGAATATCACTGGTAACGTTCTGGTTGATAAGGATAAAACCGCAGACAATGCGGCGGAATATTTCTTCGATCCATCCGTGGGTATCAACGTTTACGGCAGTGATAATAACGTGACGCt                                <  2:724723/149‑1 (MQ=255)
                                                                                                                    ggCATAAACGTTTCTGGCGATGCGAATACAGTGAATATCACTGGTAACGTTCTGGTTGATAAGGATAAAACCGCAGACAATGCGGCGGAATATTTCTTCGATCCATCCGTGGGTATCAACGTTTACGGCAGTGATAATAACGTGACGCt                                <  2:2758183/149‑1 (MQ=255)
                                                                                                                     gCATACACGTTTCTGGCGATGCGAATACAGTGAATATCACTGGTAACGTTCTGGTTGATAAGGATAAAACCGCAGACAATGCGGCGGAATATTTCTTCGATCCATCCGTGGGTATCAACGTTTACGGCAGTGATAATAACGTGACGCtg                               <  2:4302901/149‑1 (MQ=255)
                                                                                                                      cATAAACGTTTCTGGCGATGCGAATACAGTGAATATCACTGGTAACGTTCTGGTTGATAAGGATAAAACCGCAGACAATGCGGCGGAATATTTCTTCGATCCATCCGTGGGTATCAACGTTTACGGCAGTGATAATAACGTGACGCtgg                              <  2:3408045/149‑1 (MQ=255)
                                                                                                                        tAAACGTTTCTGGCGATGCGAATACAGTGAATATCACTGGTAACGTTCTGGTTGATAAGGATAAAACCGCAGACAATGCGGCGGAATATTTCTTCGATCCATCCGTGGGTATCAACGTTTACGGCAGTGATAATAACGTGACGCtggat                            >  1:4581045/1‑149 (MQ=255)
                                                                                                                        tAAACGTTTCTGGCGATGCGAATACAGTGAATATCACTGGTAACGTTCTGGTTGATAAGGATAAAACCGCAGACAATGCGGCGGAATATTTCTTCGATCCATCCGTGGGTATCAACGTTTACGGCAGTGATAATAACGTGACGCtggat                            >  1:4596909/1‑149 (MQ=255)
                                                                                                                         aaaCGTTTCTGGCGATGCGAATACAGTGAATATCACTGGTAACGTTCTGGTTGATAAGGATAAAACCGCAGACAATGCGGCGGAATATTTCTTCGATCCATCCGTGGGTATCAACGTTTACGGCAGTGATAATAACGTGACGCtggatg                           >  1:1888149/1‑149 (MQ=255)
                                                                                                                          aaCGTTTCTGGCGATGCGAATACAGTGAATATCACTGGTAACGTTCTGGTTGATAAGGATAAAACCGCAGACAATGCGGCGGAATATTTCTTCGATCCATCCGTGGGTATCAACGTTTACGGCAGTGATAATAACGTGACGCtggatgg                          <  2:3452454/149‑1 (MQ=255)
                                                                                                                             gTTTCTGGCGATGCGAATACAGTGAATATCACTGGTAACGTTCTGGTTGATAAGGATAAAACCGCAGACAATGCGGCGGAATATTTCTTCGATCCATCCGTGGGTATCAACGTTTACGGCAGTGATAATAACGTGACGCTGGATGGaaa                       >  1:3608689/1‑149 (MQ=255)
                                                                                                                              tttCTGGCGATGCGAATACAGTTAATATCACTGGTAACGTTCTGGTTGATAAGGATAACACCGCAGACAATGCGGCGGAATATTTCTTCGATCCATCCGTGGGTATCAACGTTTACGGCAGTGATAATAACGTGACGCTGGATGGAAAg                      <  1:2635490/149‑1 (MQ=255)
                                                                                                                              tttCTGGCGATGCGAATACAGTGAATATCACTGGTAACGTTCTGGTTGATAAGGATAAAACCGCAGACAATGCGGCGGAATATTTCTTCGATCCATCCGTGGGTATCAACGTTTACGGCAGTGATAATAACGTGACGCTGGATGGAAAg                      >  1:207288/1‑149 (MQ=255)
                                                                                                                               ttCTGGCGATGCGAATACAGTGAATATCACTGGTAACGTTCTGGTTGATAAGGATAAAACCGCAGACAATGCGGCGGAATATTTCTTCGATCCATCCGTGGGTATCAACGTTTACGGCAGTGATAATAACGTGACGCTGGATGGAAAGt                     <  1:1506454/149‑1 (MQ=255)
                                                                                                                                  tGGCGATGCGAATACAGTTAATATCACTGGTAACGTTCTGGTTGATAAGGATAAAACCGCAGACAATGCGGCGGAATATTTCTTCGATCCATCCGTGGGTATCAACGTTTACGGCAGTGATAATAACGTGACGCTGGATGGAAAGTTaa                  >  1:719650/1‑149 (MQ=255)
                                                                                                                                     cGATGCGAATACAGTGAATATCACTGGTAACGTTCTGGTTGATAAGGATAAAACCGCAGACAATGCGGCGGAATATTTCTTCGATCCATCCGTGGGTATCAACGTTTACGGCAGTGATAATAACGTGACGCTGGATGGAAAGTTAACtg               <  1:4317404/149‑1 (MQ=255)
                                                                                                                                        tGCGAATACAGTGAATATCACTGGTAACGTTCTGGTTGATAAGGATAAAACCGCAGACAATGCGGCGGAATATTTCTTCGATCCATCCGTGGGTATCAACGTTTACGGCAGTGATAATAACGTGACGCTGGATGGAAAGTTAACtgttg            >  2:2708835/1‑149 (MQ=255)
                                                                                                                                        tGCGAATACAGTGAATATCACTGGTAACGTTCTGGTTGATAAGGATAAAACCGCAGACAATGCGGCGGAATATTTCTTCGATCCATCCGTGGGTATCAACGTTTACGGCAGTGATAATAACGTGACGCTGGATGGAAAGTTAACtgttg            >  1:5096376/1‑149 (MQ=255)
                                                                                                                                         gCGAATACAGTGAATATCACTGGTAACGTTCTGGTTGATAAGGATAAAACCGCAGACAATGCGGCGGAATATTTCTTCGATCCATCCGTGTGTATCAACGTTTACGGCAGTGATAATAACGTGACGCTGGATGGAAAGTTAACtgttgt           <  2:482213/149‑1 (MQ=255)
                                                                                                                                           gAATACAGTGAATATCACTGGTAACGTTCTGGTTGATAAGGATAAAACCGCAGACAATGCGGCGGAATATTTCTTCGATCCATCCGTGGGTATCAACGTTTACGGCAGTGATAATAACGTGACGCTGGATGGAAAGTTAACTGTTGTAt         <  1:612029/149‑1 (MQ=255)
                                                                                                                                             aTACAGTGAATATCACTGGTAACGTTCTGGTTGATAAGGATAAAACCGCAGACAATGCGGCGGAATATTTCTTCGATCCATCCGTGGGTATCAACGTTTACGGCAGTGATAATAACGTGACGCTGGATGGAAAGTTAACTGTTGTATCa       <  1:3790426/149‑1 (MQ=255)
                                                                                                                                                cAGTGAATATCACTGGTAACGTTCTGGTTGATAAGGATAAAACCGCAGACAATGCGGCGGAATATTTCTTCGATCCATCCGTGGGTATCAACGTTTACGGCAGTGATAATAACGTGACGCTGGATGGAAAGTTAACTGTTGTATCAGAc    >  2:2417734/1‑149 (MQ=255)
                                                                                                                                                cAGTGAATATCACTGGTAACGTTCTGGTTGATAAGGATAAAACCGCAGACAATGCGGCGGAATATTTCTTCGATCCATCCGTGGGTATCAACGTTTACGGCAGTGATAATAACGTGACGCTGGATGGAAAGTTAACTGTTGTATCAGAc    >  2:223155/1‑149 (MQ=255)
                                                                                                                                                 aGTGAATATCACTGGTAACGTTCTGGTTGATAAGGATAAAACCGCAGACAATGCGGCGGAATATTTCTTCGATCCATCCGTGGGTATCAACGTTTACGGCAGTGATAATAACGTGACGCTGGATGGAAAGTTAACTGTTGTATCAGACa   <  2:1280173/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
CGCTGGCAGGCAGTATGCAGGTTGGCGATTTCTCGACCGGGGTAGATCTTAATGGCAACAATAACAGCGTGACGCTGGCGGCAAAAGATCTAAAAGTGGTCGGGCAGAAAGCGACGGGCATAAACGTTTCTGGCGATGCGAATACAGTGAATATCACTGGTAACGTTCTGGTTGATAAGGATAAAACCGCAGACAATGCGGCGGAATATTTCTTCGATCCATCCGTGGGTATCAACGTTTACGGCAGTGATAATAACGTGACGCTGGATGGAAAGTTAACTGTTGTATCAGACAG  >  CP009273/1461504‑1461798

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: