Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,461,041 A→T 100% pseudogene (1393/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,0410AT96.5% 622.8 / 9.8 198pseudogene (1393/2513 nt)ydbApseudogene, autotransporter homolog; interrupted by IS2 and IS30
Reads supporting (aligned to +/- strand):  ref base A (3/4);  new base T (120/71);  total (123/75)
Fisher's exact test for biased strand distribution p-value = 4.30e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.20e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 198 total aligned reads displayed.

CGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGA  >  CP009273/1460896‑1461185
                                                                                                                                                 |                                                                                                                                                
   tgatgaCGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGc                                                                                                                                             >  2:2515783/1‑148 (MQ=255)
     atgaCGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGa                                                                                                                                           >  1:2244421/1‑148 (MQ=255)
       gaCGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACa                                                                                                                                         <  2:97085/148‑1 (MQ=255)
       gaCGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACa                                                                                                                                         <  2:96651/148‑1 (MQ=255)
        actccACGGCCAACCACAACCGCATCACCCCCGTTGACGGCAAGGAGTCCACCGGTACGGAAATGTCTGGCAATAACGGGAATGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAg                                                                                                                                        <  2:2409506/145‑1 (MQ=14)
          gCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAgcg                                                                                                                                      <  1:2841251/148‑1 (MQ=255)
               ggCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACg                                                                                                                                 >  2:1015387/1‑148 (MQ=255)
               ggCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACg                                                                                                                                 >  2:3093230/1‑148 (MQ=255)
                 cAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGt                                                                                                                               <  1:1719623/148‑1 (MQ=255)
                 cAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGt                                                                                                                               >  1:3093208/1‑148 (MQ=255)
                 cAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGt                                                                                                                               >  1:589740/1‑148 (MQ=255)
                   aacaacaacGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTgg                                                                                                                             >  1:414051/1‑148 (MQ=255)
                   aacaacaacGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTgg                                                                                                                             >  1:1400777/1‑148 (MQ=255)
                   aacaacaacGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTgg                                                                                                                             <  2:169601/148‑1 (MQ=255)
                      aacaacGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATa                                                                                                                          >  1:1302930/1‑148 (MQ=255)
                        caacGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAAc                                                                                                                        >  1:1047000/1‑148 (MQ=255)
                        caacGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAAc                                                                                                                        >  2:3534088/1‑148 (MQ=255)
                        caacGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAAc                                                                                                                        >  1:1136241/1‑148 (MQ=255)
                         aacGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACa                                                                                                                       <  1:3152171/148‑1 (MQ=255)
                         aacGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACa                                                                                                                       <  1:3152227/148‑1 (MQ=255)
                         aacGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACa                                                                                                                       <  1:1526854/148‑1 (MQ=255)
                             gCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAggg                                                                                                                   >  2:2108611/1‑148 (MQ=255)
                             gCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAggg                                                                                                                   >  1:920369/1‑148 (MQ=255)
                             gCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAggg                                                                                                                   >  2:36882/1‑148 (MQ=255)
                                aaaCCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCAc                                                                                                                >  1:552715/1‑148 (MQ=255)
                                 aaCCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCAcc                                                                                                               >  2:3189485/1‑148 (MQ=255)
                                    caccGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTTGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATg                                                                                                            >  2:1698611/1‑148 (MQ=255)
                                    caccGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTTGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATg                                                                                                            >  2:1697614/1‑148 (MQ=255)
                                     accGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGa                                                                                                           <  2:948318/148‑1 (MQ=255)
                                        gTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCg                                                                                                        <  2:1279691/148‑1 (MQ=255)
                                         ttGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGt                                                                                                       >  2:207406/1‑148 (MQ=255)
                                           gACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCa                                                                                                     <  2:522218/148‑1 (MQ=255)
                                            actgctaGGAGTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCAc                                                                                                    <  2:3394077/142‑1 (MQ=255)
                                            aCGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCAc                                                                                                    >  1:3528946/1‑148 (MQ=255)
                                            aCGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCAc                                                                                                    >  1:2605013/1‑148 (MQ=255)
                                                cAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGAt                                                                                                <  1:1779274/148‑1 (MQ=255)
                                                cAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGAt                                                                                                >  2:106421/1‑148 (MQ=255)
                                                 aaGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATc                                                                                               <  1:3534088/148‑1 (MQ=255)
                                                  aGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATcc                                                                                              >  2:1215627/1‑148 (MQ=255)
                                                  aGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATcc                                                                                              >  2:1229034/1‑148 (MQ=255)
                                                  aGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATcc                                                                                              >  2:3347535/1‑148 (MQ=255)
                                                       tCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGt                                                                                         >  2:1987826/1‑148 (MQ=255)
                                                        ccACCGGTACGGAAAGTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTc                                                                                        <  1:3043137/148‑1 (MQ=255)
                                                            cGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATc                                                                                    <  1:466637/148‑1 (MQ=255)
                                                            cGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATc                                                                                    <  1:1961135/148‑1 (MQ=255)
                                                            cGGTACGGAAATGTAGGGCAATAACGGTAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATc                                                                                    <  2:3302207/148‑1 (MQ=11)
                                                              gTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCgg                                                                                  >  1:1186889/1‑148 (MQ=255)
                                                                   gAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATcc                                                                             >  2:754301/1‑148 (MQ=255)
                                                                   gAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATcc                                                                             >  2:1471239/1‑148 (MQ=255)
                                                                      aTTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGa                                                                          >  1:773698/1‑148 (MQ=255)
                                                                         gCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTg                                                                       <  1:1074481/148‑1 (MQ=255)
                                                                              cAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGc                                                                  >  2:2656540/1‑148 (MQ=255)
                                                                              cAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGc                                                                  >  2:876516/1‑148 (MQ=255)
                                                                                 tAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGAc                                                               <  1:2775174/148‑1 (MQ=255)
                                                                                 tAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGt                                                                                                       <  1:3332918/108‑1 (MQ=17)
                                                                                    cGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAg                                                            >  2:2721230/1‑148 (MQ=255)
                                                                                     gggAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAgg                                                           >  2:2401480/1‑148 (MQ=255)
                                                                                       gAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGcc                                                         >  2:3109388/1‑148 (MQ=38)
                                                                                       gAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGcc                                                         >  2:135700/1‑148 (MQ=38)
                                                                                       gAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGcc                                                         >  2:1339976/1‑148 (MQ=38)
                                                                                       gAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGcc                                                         >  2:487232/1‑148 (MQ=38)
                                                                                          ggTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATc                                                      >  1:1272232/1‑148 (MQ=38)
                                                                                           gTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCg                                                     <  1:834297/148‑1 (MQ=38)
                                                                                            tGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGt                                                    >  1:2993903/1‑148 (MQ=38)
                                                                                                tCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAAt                                                >  1:2655484/1‑148 (MQ=255)
                                                                                                 cAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATa                                               <  2:3464276/148‑1 (MQ=255)
                                                                                                   ggACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAAc                                             >  2:2238332/1‑148 (MQ=255)
                                                                                                    gACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACg                                            >  1:3293479/1‑148 (MQ=255)
                                                                                                    gACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACg                                            >  1:96944/1‑148 (MQ=255)
                                                                                                     aCGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGa                                           <  2:2947149/148‑1 (MQ=255)
                                                                                                       ggCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAg                                         <  1:3205144/148‑1 (MQ=255)
                                                                                                        gCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAgg                                        <  2:1905976/148‑1 (MQ=255)
                                                                                                           aTCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCga                                     <  1:376107/148‑1 (MQ=255)
                                                                                                            tCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCgag                                    >  2:1258667/1‑148 (MQ=255)
                                                                                                            tCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGCCAATAACGAAGGCgag                                    >  1:1744072/1‑148 (MQ=255)
                                                                                                               ggATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGc                                 <  2:1128331/148‑1 (MQ=255)
                                                                                                               ggATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGc                                 <  1:2660250/148‑1 (MQ=255)
                                                                                                                gATGTCAGCTGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCa                                <  2:1382877/148‑1 (MQ=255)
                                                                                                                gATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCa                                >  2:2278270/1‑148 (MQ=255)
                                                                                                                gATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCa                                >  1:2729617/1‑148 (MQ=255)
                                                                                                                  tGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACt                              <  2:2773456/148‑1 (MQ=255)
                                                                                                                   gTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTa                             <  2:989327/148‑1 (MQ=255)
                                                                                                                     cAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATc                           <  2:589740/148‑1 (MQ=255)
                                                                                                                     cAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATc                           >  2:1554294/1‑148 (MQ=255)
                                                                                                                        cggcggcggTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCAcc                        <  1:3328914/148‑1 (MQ=255)
                                                                                                                        cggcggcggTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCAcc                        >  2:3009714/1‑148 (MQ=255)
                                                                                                                               ggTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATggcg                 >  1:2443565/1‑148 (MQ=255)
                                                                                                                                 tCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATggcggc               >  2:1613012/1‑148 (MQ=255)
                                                                                                                                   aCGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTTTCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGcggcac             <  2:3093174/148‑1 (MQ=255)
                                                                                                                                   aCGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGcggcac             <  2:2655484/148‑1 (MQ=255)
                                                                                                                                    cGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGcggcacc            <  2:2819507/148‑1 (MQ=255)
                                                                                                                                    cGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGcggcacc            >  2:2393339/1‑148 (MQ=255)
                                                                                                                                         tttatatCTCCGGCGACAGAGGACCGGAGGATAACAAGGGCACCACGTCCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGcggcaccggcac       <  2:2169961/145‑1 (MQ=17)
                                                                                                                                         ttGATATCTCCGGCGACAGCGCAACGGTGGAAAACAAGGGCACCAGGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGcggcaccggcac       <  2:2170582/148‑1 (MQ=255)
                                                                                                                                          tGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACt      >  2:2876714/1‑148 (MQ=255)
                                                                                                                                          tGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACt      >  2:2331479/1‑148 (MQ=255)
                                                                                                                                          tGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACt      >  2:1678982/1‑148 (MQ=255)
                                                                                                                                            atatCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCa    >  2:307401/1‑148 (MQ=255)
                                                                                                                                            atatCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCa    >  2:2100413/1‑148 (MQ=255)
                                                                                                                                              atCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGa  >  2:808689/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
CGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGA  >  CP009273/1460896‑1461185

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: