Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,889,129 C→A 90.9% A4S (GCA→TCA)  ygcQ ← putative flavoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,889,1290CA90.9% 395.7 / 32.1 150A4S (GCA→TCA) ygcQputative flavoprotein
Reads supporting (aligned to +/- strand):  ref base C (5/9);  new base A (57/79);  total (62/88)
Fisher's exact test for biased strand distribution p-value = 7.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 154 total aligned reads displayed.

GCCACTGTTCCACCAACGCATCCAGCACCTGTTCCGCCACTACAGGCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGCGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGTCGCTCCGCGTCGTTGCTCTGCTCTTGCCAGTTGCAGACACTGCATCGCTGGCATTTCG  >  CP009273/2888982‑2889273
                                                                                                                                                   |                                                                                                                                                
gCCACTGTTCCACCAACGCATCCAGCACCTGTTCCGCCACTACAGGCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGc                                                                                                                                                  <  1:1739951/148‑1 (MQ=255)
  cACTGTTCCACCAACGCATCCAGCACCTGTTCCGCCACTACAGGCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTgaga                                                                                                                                                >  1:3290831/1‑148 (MQ=255)
  cACTGTTCCACCAACGCATCCAGCACCTGTTCCGCCACTACAGGCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTgaga                                                                                                                                                <  1:598781/148‑1 (MQ=255)
       ttCCACCAACGCATCCAGCACCTGTTCCGCCACTACAGGCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGCGATATTc                                                                                                                                           <  1:1696199/148‑1 (MQ=255)
        tCCACCAACGCATCCAGCACCTGTTCCGCCACTACAGGCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGCGATATTCa                                                                                                                                          >  2:2301739/1‑148 (MQ=255)
           accaACGCATCCAGCACCTGTTCCGCCACTACAGGCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATgg                                                                                                                                       <  1:289476/148‑1 (MQ=255)
           accaACGCATCCAGCACCTGTTCCGCCACTACAGGCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATgg                                                                                                                                       <  1:1891556/148‑1 (MQ=255)
           accaACGCATCCAGCACCTGTTCCGCCACTACAGGCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATgg                                                                                                                                       <  2:2443915/148‑1 (MQ=255)
           accaACGCATCCAGCACCTGTTCCGCCACTACAGGCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATgg                                                                                                                                       <  2:2856155/148‑1 (MQ=255)
           accaACGCATCCAGCACCTGTTCCGCCACTACAGGCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATgg                                                                                                                                       <  2:556410/148‑1 (MQ=255)
           accaACGCATCCAGCACCTGTTCCGCCACTACAGGCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATgg                                                                                                                                       <  2:784196/148‑1 (MQ=255)
            ccaACGCATCCAGCACCTGTTCCGCCACTACAGGCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGCGATATTCATGGc                                                                                                                                      <  2:2347782/148‑1 (MQ=255)
            ccaACGCATCCAGCACCTGTTCCGCCACTACAGGCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGc                                                                                                                                      <  2:1772269/148‑1 (MQ=255)
               aCGCATCCAGCACCTGTTCCGCCACTACAGGCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGc                                                                                                                                   >  2:1188871/1‑148 (MQ=255)
                cGCATCCAGCACCTGTTCCGCCACTACAGGCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCa                                                                                                                                  >  1:921836/1‑148 (MQ=255)
                cGCATCCAGCACCTGTTCCGCCACTACAGGCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCa                                                                                                                                  >  1:874207/1‑148 (MQ=255)
                cGCATCCAGCACCTGTTCCGCCACTACAGGCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCa                                                                                                                                  >  1:2186514/1‑148 (MQ=255)
                    tCCAGCACCTGTTCCGCCACTACAGGCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCt                                                                                                                              <  2:1696063/148‑1 (MQ=255)
                          aCCTGTTCCGCCACTACAGGCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGCGATATTCATGGCTGCATCCTTTGGCg                                                                                                                        <  1:1079257/148‑1 (MQ=255)
                          aCCTGTTCCGCCACTACAGGCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGCGATATTCATGGCTGCATCCTTTGGCg                                                                                                                        <  1:1079482/148‑1 (MQ=255)
                           ccTGTTCCGCCACTACAGGCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGt                                                                                                                       >  2:134195/1‑148 (MQ=255)
                             tGTTCCGCCACTACAGGCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTaa                                                                                                                     <  2:1902632/148‑1 (MQ=255)
                              gTTCCGCCACTACAGGCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTaaa                                                                                                                    >  1:1211060/1‑148 (MQ=255)
                              gTTCCGCCACTACAGGCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTaaa                                                                                                                    >  1:2861074/1‑148 (MQ=255)
                               gtCCGCCACTACAGGCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAAt                                                                                                                   <  2:3277302/147‑1 (MQ=255)
                                tCCGCCACTACAGGCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATa                                                                                                                  <  1:1916255/148‑1 (MQ=255)
                                tCCGCCACTACAGGCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATa                                                                                                                  >  2:3097/1‑148 (MQ=255)
                                  cGCCACTACAGGCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGt                                                                                                                <  1:1185684/148‑1 (MQ=255)
                                   gCCACTACAGGCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTc                                                                                                               <  2:2349085/148‑1 (MQ=255)
                                   gCCACTACAGGCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTc                                                                                                               <  2:379749/148‑1 (MQ=255)
                                    ccACTACAGGCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTcc                                                                                                              <  1:2220297/148‑1 (MQ=255)
                                    ccACTACAGGCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTcc                                                                                                              <  1:2555960/148‑1 (MQ=255)
                                      aCTACAGGCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTg                                                                                                            <  1:1214224/148‑1 (MQ=255)
                                        tACAGGCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGcc                                                                                                          <  1:822395/148‑1 (MQ=255)
                                        tACAGGCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGcc                                                                                                          <  1:821812/148‑1 (MQ=255)
                                            ggCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAg                                                                                                      <  1:231726/148‑1 (MQ=255)
                                            ggCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAg                                                                                                      >  2:1770248/1‑148 (MQ=255)
                                            ggCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAg                                                                                                      <  2:1716324/148‑1 (MQ=255)
                                               tGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGTGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCtt                                                                                                   <  1:38861/148‑1 (MQ=255)
                                               tGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCtt                                                                                                   <  1:2394166/148‑1 (MQ=255)
                                               tGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCtt                                                                                                   <  1:2394724/148‑1 (MQ=255)
                                                 aGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCtttt                                                                                                 <  1:927280/148‑1 (MQ=255)
                                                    ttCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAg                                                                                              >  1:2954481/1‑148 (MQ=255)
                                                     tCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGc                                                                                             <  2:491537/148‑1 (MQ=255)
                                                        aTCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCttt                                                                                          >  1:1272453/1‑148 (MQ=255)
                                                        aTCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCttt                                                                                          <  2:2412362/148‑1 (MQ=255)
                                                        aTCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCttt                                                                                          <  2:1433993/148‑1 (MQ=255)
                                                        aTCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCttt                                                                                          >  1:610547/1‑148 (MQ=255)
                                                         tCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCtttt                                                                                         >  2:939295/1‑148 (MQ=255)
                                                          cTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCttttt                                                                                        <  1:2809444/148‑1 (MQ=255)
                                                           tGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTc                                                                                       <  1:39017/148‑1 (MQ=255)
                                                                 tGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTAc                                                                                 <  2:1354934/148‑1 (MQ=255)
                                                                         ggTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGcc                                                                         <  2:3049817/148‑1 (MQ=255)
                                                                            aCAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCat                                                                      <  1:2680053/148‑1 (MQ=255)
                                                                            aCAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCat                                                                      <  1:2680957/148‑1 (MQ=255)
                                                                            aCAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCat                                                                      <  1:3174162/148‑1 (MQ=255)
                                                                            aCAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCat                                                                      <  1:3174574/148‑1 (MQ=255)
                                                                                ccGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCatcaat                                                                  <  1:374524/148‑1 (MQ=255)
                                                                                 cGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCatcaatc                                                                 >  1:1104488/1‑148 (MQ=255)
                                                                                 cGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCatcaatc                                                                 <  1:2566955/148‑1 (MQ=255)
                                                                                 cGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCatcaatc                                                                 >  1:2803341/1‑148 (MQ=255)
                                                                                     gAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCg                                                             >  1:1180362/1‑148 (MQ=255)
                                                                                      aaaTCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGt                                                            <  1:1912597/148‑1 (MQ=255)
                                                                                      aaaTCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGATACCGTTTGCCCATCAATCAGCGt                                                            >  2:2100161/1‑148 (MQ=255)
                                                                                         tCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCGTAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGTCGc                                                         <  2:653548/148‑1 (MQ=255)
                                                                                         tCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGTCGc                                                         <  2:653442/148‑1 (MQ=255)
                                                                                            tGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGTCGCTcc                                                      >  2:1644813/1‑148 (MQ=255)
                                                                                               gccgccAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGTCGCTCcgcg                                                   <  1:531997/148‑1 (MQ=255)
                                                                                                  gccagccaGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGTCGCTCCGcgtcg                                                >  1:766912/1‑148 (MQ=255)
                                                                                                   ccagccagCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGTCGCTCCGcgtcgt                                               >  2:424089/1‑148 (MQ=255)
                                                                                                       ccagcTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGTCGCTCCGCGTCGTtgct                                           >  2:473199/1‑148 (MQ=255)
                                                                                                        cagcTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTGGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGTCGCTCCGCGTCGTtgctc                                          <  2:1811235/148‑1 (MQ=255)
                                                                                                          gcTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGTCGCTCCGCGTCGTtgctatg                                        >  2:3252236/1‑148 (MQ=255)
                                                                                                           cTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGTCGCTCCGCGTCGTtgctctgc                                       >  1:299157/1‑148 (MQ=255)
                                                                                                             tGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGTCGCTCCGCGTCGTtgctctgctc                                     >  2:1537364/1‑148 (MQ=255)
                                                                                                              gCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGTCGCTCCGCGTCGTTgctctgctct                                    <  2:1656311/148‑1 (MQ=255)
                                                                                                                 aTTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGTCGCTCCGCGTCGTTGCTCTGCTCTTGc                                 >  2:818289/1‑148 (MQ=255)
                                                                                                                   tGCGGCATTCTCCTGATTTATTGTCACAATTGCGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGTCGCTCCGCGTCGTTGCTCTGCTCTTGCCa                               >  2:2435391/1‑148 (MQ=255)
                                                                                                                   tGCGGCATTCTCCTGATTTATTGTCACAATTGCGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGTCGCTCCGCGTCGTTGCTCTGCTCTTGCCa                               >  2:1048262/1‑148 (MQ=255)
                                                                                                                       gCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGTCGCTCCGCGTCGTTGCTCTGCTCTTGCCAGTTg                           <  1:2600963/148‑1 (MQ=255)
                                                                                                                        cATTCTCCTGATTTATTGTAACAATTTAGATATTAATGGCTGCATCCTTTGCCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTATCTTCTACAGTTTGCACATCAATCAGCGTCGCTACGCGTCTTTGCTCTGCTCTTTACAGATGc                          >  2:1781762/1‑148 (MQ=255)
                                                                                                                           tctcCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGTCGCTCCGCGTCGTTGCTCTGCTCTTGCCAGTTGCAGa                       >  1:1095645/1‑148 (MQ=255)
                                                                                                                           tctcCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGTCGCTCCGCGTCGTTGCTCTGCTCTTGCCAGTTGCAGa                       >  2:1450567/1‑148 (MQ=255)
                                                                                                                           tctcCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGTCGCTCCGCGTCGTTGCTCTGCTCTTGCCAGTTGCAGa                       >  1:2747254/1‑148 (MQ=255)
                                                                                                                            ctcCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGTCGCTCCGCGTCGTTGCTCTGCTCTTGCCAGTTGCAGac                      <  1:95513/148‑1 (MQ=255)
                                                                                                                                  aTTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGTCGCTCCGCGTCGTTGCTCTGCTCTTGCCAGTTGCAGACACTGCa                <  2:2861074/148‑1 (MQ=255)
                                                                                                                                  aTTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGTCGCTCCGCGTCGTTGCTCTGCTCTTGCCAGTTGCAGACACTGCa                <  1:3409/148‑1 (MQ=255)
                                                                                                                                   tttATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGTCGCTCCGCGTCGTTGCTCTGCTCTTGCCAGTTGCAGACACTGCAt               >  1:1705345/1‑148 (MQ=255)
                                                                                                                                    ttATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGTCGCTCCGCGTCGTTGCTCTGCTCTTGCCAGTTGCAGACACTGCATc              <  1:2680855/148‑1 (MQ=255)
                                                                                                                                    ttATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGTCGCTCCGCGTCGTTGCTCTGCTCTTGCCAGTTGCAGACACTGCATc              >  1:1612055/1‑148 (MQ=255)
                                                                                                                                        tGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGTCGCTCCGCGTCGTTGCTCTGCTCTTGCCAGTTGCAGACACTGCATCGCTg          >  1:1166916/1‑148 (MQ=255)
                                                                                                                                        tGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGTCGCTCCGCGTCGTTGCTCTGCTCTTGCCAGTTGCAGACACTGCATCGCTg          >  1:1166876/1‑148 (MQ=255)
                                                                                                                                           cacaATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGTCGCTCCGCGTCGTTGCTCTGCTCTTGCCAGTTGCAGACACTGCATCGCTGGCa       <  1:2436015/148‑1 (MQ=255)
                                                                                                                                           cacaATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCATGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGTCGCTCCGCGTCGTTGCTCTGCTCTTGCCAGTTGCAGACACTGCATCGCTGGCa       >  1:1686017/1‑148 (MQ=255)
                                                                                                                                             caATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGTCGCTCCGCGTCGTTGCTCTGCTCTTGCCAGTTGCAGACACTGCATCGCTGGCAtt     >  2:1633715/1‑148 (MQ=255)
                                                                                                                                              aaTTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGTCGCTCCGCGTCGTTGCTCTGCTCTTGCCAGTTGCAGACACTGCATCGCTGGCAttt    <  2:2713981/148‑1 (MQ=255)
                                                                                                                                              aaTTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGTCGCTCCGCGTCGTTGCTCTGCTCTTGCCAGTTGCAGACACTGCATCGCTGGCAttt    <  1:134195/148‑1 (MQ=255)
                                                                                                                                               aTTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGTCGCTCCGCGTCGTTGCTCTGCTCTTGCCAGTTGCAGACACTGCATCGCTGGCATTTc   >  2:2436266/1‑148 (MQ=255)
                                                                                                                                                ttGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGTCGCTCCGCGTCGTTGCTCTGCTCTTGCCAGTTGCAGACACTGCATCGCTGGCATTTCg  <  2:2648775/148‑1 (MQ=255)
                                                                                                                                                ttGAGATATTCATGGCTGCATCCTTGGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGTCGCTCCGCGTCGTTGCTCTGCTCTTGCCAGTTGCAGACACTGCATCGCTGGCATTTCg  <  2:1747560/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
GCCACTGTTCCACCAACGCATCCAGCACCTGTTCCGCCACTACAGGCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGCGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGTCGCTCCGCGTCGTTGCTCTGCTCTTGCCAGTTGCAGACACTGCATCGCTGGCATTTCG  >  CP009273/2888982‑2889273

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: