Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,181,504 G→T 88.9% D69Y (GAT→TAT)  pepT → peptidase T

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,181,5040GT88.9% 475.1 / 57.7 189D69Y (GAT→TAT) pepTpeptidase T
Reads supporting (aligned to +/- strand):  ref base G (10/11);  new base T (78/90);  total (88/101)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01

Only 100 of 191 total aligned reads displayed.

GGGGTGAGACAGGTTCCCAGCACGGAAGGCCAATGGAAGTTATTGCATCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCGATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAATTGTTGAAAACTATCGCGGTGGCGATATTGCGCTGGGTATCGGCGATGAAGTTTTATCACCGGTTATG  >  CP009273/1181360‑1181647
                                                                                                                                                |                                                                                                                                               
ggggTGAGACAGGTTCCCAGCACGGAAGGCCAATGGAAGTTATTGCATCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCtata                                                                                                                                              >  1:1874021/1‑148 (MQ=255)
 gggTGAGACAGGTTCCCAGCACGGAAGGCCAATGGAAGTTATTGCATCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTatat                                                                                                                                             >  2:2348569/1‑148 (MQ=255)
 gggTGAGACAGGTTCCCAGCACGGAAGGCCAATGGAAGTTATTGCATCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTatat                                                                                                                                             >  1:2537986/1‑148 (MQ=255)
 gggTGAGACAGGTTCCCAGCACGGAAGGCCAATGGAAGTTATTGCATCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTatat                                                                                                                                             >  2:1907856/1‑148 (MQ=255)
  ggTGAGACAGGTTCCCAGCACGGAAGGCCAATGGAAGTTATTTCATCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATc                                                                                                                                            <  2:1017112/148‑1 (MQ=255)
  ggTGAGACAGGTTCCCAGCACGGAAGGCCAATGGAAGTTATTGCATCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATc                                                                                                                                            <  2:1481096/148‑1 (MQ=255)
  ggTGAGACAGGTTCCCAGCACGGAAGGCCAATGGAAGTTATTGCATCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCGATATc                                                                                                                                            <  2:2020423/148‑1 (MQ=255)
         cAGGTTCCCAGCACGGAAGGCCAATGGAAGTTATTGCATCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGa                                                                                                                                     <  2:1431044/148‑1 (MQ=255)
         cAGGTTCCCAGCACGGAAGGCCAATGGAAGTTATTGCATCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGa                                                                                                                                     <  2:1432690/148‑1 (MQ=255)
          aGGTTCCCAGCACGGAAGGCCAATGGAAGTTATTGCATCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGAt                                                                                                                                    <  1:880256/148‑1 (MQ=255)
          aGGTTCCCAGCACGGAAGGCCAATGGAAGTTATTGCATCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGAt                                                                                                                                    <  2:1191203/148‑1 (MQ=255)
           ggTTCCCAGCACGGAAGGCCAATGGAAGTTATTGCATCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGAtt                                                                                                                                   >  1:3016395/1‑148 (MQ=255)
              tCCCAGCACGGAAGGCCAATGGAAGTTATTGCATCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGc                                                                                                                                <  1:138892/148‑1 (MQ=255)
                ccAGCACGGAAGGCCAATGGAAGTTATTGCATCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCtt                                                                                                                              <  1:555528/148‑1 (MQ=255)
                ccAGCACGGAAGGCCAATGGAAGTTATTGCATCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCtt                                                                                                                              <  1:2478488/148‑1 (MQ=255)
                       ggAAGGCCAATGGAAGTTATTGCATCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTtct                                                                                                                       <  2:154779/148‑1 (MQ=255)
                       ggAAGGCCAATGGAAGTTATTGCATCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTtct                                                                                                                       <  2:1981623/148‑1 (MQ=255)
                        gAAGGCCAATGGAAGTTATTGCATCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTtctc                                                                                                                      >  2:1510661/1‑148 (MQ=255)
                         aaGGCCAATGGAAGTTATTGCATCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCa                                                                                                                     >  1:2083425/1‑148 (MQ=255)
                          aGGCCAATGGAAGTTATTGCATCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCAt                                                                                                                    <  2:893818/148‑1 (MQ=255)
                           ggCCAATGGAAGTTATTGCATCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCAtg                                                                                                                   <  2:2089518/148‑1 (MQ=255)
                            gCCAATGGAAGTTATTGCATCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCAtgt                                                                                                                  <  1:2282183/148‑1 (MQ=255)
                            gCCAATGGAAGTTATTGCATCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCAtgt                                                                                                                  <  1:2390408/148‑1 (MQ=255)
                            gCCAATGGAAGTTATTGCATCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCAtgt                                                                                                                  <  2:504895/148‑1 (MQ=255)
                                aTGGAAGTTATTGCATCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGAt                                                                                                              >  2:2142204/1‑148 (MQ=255)
                                aTGGAAGTTATTGCATCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGAt                                                                                                              >  2:2530509/1‑148 (MQ=255)
                                  ggAAGTTATTGCATCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATAc                                                                                                            <  2:78883/148‑1 (MQ=255)
                                   gAAGTTATTGCATCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATAcc                                                                                                           >  1:2015821/1‑148 (MQ=255)
                                      gTTATTGCATCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCa                                                                                                        >  2:764061/1‑148 (MQ=255)
                                       ttATTGCATCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCAc                                                                                                       <  1:2435819/148‑1 (MQ=255)
                                          ttGCATCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCgg                                                                                                    >  1:2863237/1‑148 (MQ=255)
                                           tGCATCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGa                                                                                                   >  1:870986/1‑148 (MQ=255)
                                           tGCATCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGAGTGGCTTTATTTCTCATGTGGATACCTCACCGGa                                                                                                   >  1:870703/1‑148 (MQ=255)
                                            gCATCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGAt                                                                                                  <  2:2238820/148‑1 (MQ=255)
                                             cATCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGAtt                                                                                                 >  1:1129069/1‑148 (MQ=255)
                                             cATCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCGATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGAtt                                                                                                 >  1:925651/1‑148 (MQ=255)
                                               tCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGc                                                                                               <  2:1001426/148‑1 (MQ=255)
                                                ctgctgAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCa                                                                                              >  2:204166/1‑148 (MQ=255)
                                                    tgAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCgg                                                                                          <  2:2224541/148‑1 (MQ=255)
                                                       aaGAGCAGCTCGAAGAGATGGGGCTTATCACTGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCaa                                                                                       <  2:2799930/148‑1 (MQ=255)
                                                        agagCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCaaa                                                                                      >  1:821093/1‑148 (MQ=255)
                                                             agcTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAAtgt                                                                                 <  2:1486837/148‑1 (MQ=255)
                                                                 cGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAAt                                                                             <  2:875036/148‑1 (MQ=255)
                                                                  gAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATc                                                                            <  2:1133289/148‑1 (MQ=255)
                                                                    agagATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCg                                                                          >  1:1908074/1‑148 (MQ=255)
                                                                    agagATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCg                                                                          >  2:1672798/1‑148 (MQ=255)
                                                                        aTGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCGATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCaaa                                                                      >  2:212266/1‑148 (MQ=255)
                                                                        aTGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCGATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCaaa                                                                      >  2:1882841/1‑148 (MQ=255)
                                                                        aTGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCGATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCaaa                                                                      <  2:2619887/148‑1 (MQ=255)
                                                                         tGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAAt                                                                     <  1:296466/148‑1 (MQ=255)
                                                                          ggggCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCGATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAAtt                                                                    <  1:2420212/148‑1 (MQ=255)
                                                                                 aTCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAATTGTTGaaa                                                             >  1:2763938/1‑148 (MQ=255)
                                                                                   cAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAATTGTTGAAAAc                                                           <  2:1957926/148‑1 (MQ=255)
                                                                                    aaTGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAATTGTTGAAAACt                                                          <  1:1353055/148‑1 (MQ=255)
                                                                                       gtgACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCGATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAATTGTTGATAACTATc                                                       >  2:1793607/1‑148 (MQ=255)
                                                                                        tgACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAATTGTTGAAAACTATcg                                                      <  2:2365916/148‑1 (MQ=255)
                                                                                          aCCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAATTGTTGAAAACTATcgcg                                                    <  1:2559519/148‑1 (MQ=255)
                                                                                          aCCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCGATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAATTGTTGAAAACTATcgcg                                                    <  1:1887350/148‑1 (MQ=255)
                                                                                             ttAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAATTGTTGAAAACTATCGCGGTg                                                 <  1:2142204/148‑1 (MQ=255)
                                                                                             ttAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCGAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAATTGTTGAAAACTATCGCGGTg                                                 <  2:2888054/148‑1 (MQ=255)
                                                                                                  tGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCAATTGGATACCAAACCGGCTTGCAGCGGCAAAAATGTGAATCCGCCCATTGTTTTAAAATATAGCGGTGGCgcg                                            >  1:2810109/1‑146 (MQ=255)
                                                                                                   gagaAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAATTGTTGAAAACTATCGCGGTGGCGata                                           >  2:3118218/1‑148 (MQ=255)
                                                                                                      aaGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAATTGTTGAAAACTATCGCGGTGGCGATATTg                                        <  1:1510660/148‑1 (MQ=255)
                                                                                                        gggCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAATTGTTGAAAACTATCGCGGTGGCGATATTgcg                                      <  2:1570786/148‑1 (MQ=255)
                                                                                                        gggCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAATTGTTGAAAACTATCGCGGTGGCGATATTgcg                                      >  2:386024/1‑148 (MQ=255)
                                                                                                        gggCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAATTGTTGAAAACTATCGCGGTGGCGATATTgcg                                      <  1:589401/148‑1 (MQ=255)
                                                                                                        gggCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAATTGTTGAAAACTATCGCGGTGGCGATATTgcg                                      <  2:1988373/148‑1 (MQ=255)
                                                                                                        gggCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCGTTATTTCTCAATTGGAAACCTCACAGGATTTCACCCGCAAAAATGTGAAACAGCAAATTGATGAAAACAAAAGCGGTGGCGATATTgac                                      >  1:2240463/1‑146 (MQ=255)
                                                                                                          gCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAATTGTTGAAAACTATCGCGGTGGCGATATTGCGCt                                    <  1:2479798/148‑1 (MQ=255)
                                                                                                          gCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAATTGTTGAAAACTATCGCGGTGGCGATATTGCGCt                                    <  1:738931/148‑1 (MQ=255)
                                                                                                           cACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCATCGGCAAAAATGTGAATCCGCAAATTGTTGAAAACTATCGCGGTGGCGATATTGCGCTg                                   >  1:67184/1‑148 (MQ=255)
                                                                                                           cACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAATTGTTGAAAACTATCGCGGTGGCGATATTGCGCTg                                   >  2:905897/1‑148 (MQ=255)
                                                                                                           cACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAATTGTTGAAAACTATCGCGGTGGCGATATTGCGCTg                                   <  1:1639168/148‑1 (MQ=255)
                                                                                                              tttGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAATTGTTGAAAACTATCGCGGTGGCGATATTGCGCTGGGt                                >  2:980534/1‑148 (MQ=255)
                                                                                                              tttGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAATTGTTGAAAACTATAGCGGTGGCGATATTGCGCTGGGt                                <  1:432/148‑1 (MQ=255)
                                                                                                               ttGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGGTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAAGTGTTGAAAACTATCGCGGTGGCGATATTGCGCTGGGTa                               >  1:1788951/1‑148 (MQ=255)
                                                                                                                 gATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAATTGTTGAAAACTATCGCGGTGGCGATATTGCGCTGGGTATc                             <  2:962538/148‑1 (MQ=255)
                                                                                                                  aTGGCGACGTTACCGGCTAACGTCCCTGGCGATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAATTGTTGAAAACTATCGCGGTGGCGATATTGCGCTGGGTATCg                            >  2:2439219/1‑148 (MQ=255)
                                                                                                                      cGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAATTGTTGAAAACTATCGCGGTGGCGATATTGCGCTGGGTATCGGCGa                        >  2:684870/1‑148 (MQ=255)
                                                                                                                        aCGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAATTGTTGAAAACTATCGCGGTGGCGATATTGCGCTGGGTATCGGCGATg                      >  2:2744317/1‑148 (MQ=255)
                                                                                                                            tACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAATTGTTGAAAACTATCGCGGTGGCGATATTGCGCTGGGTATAGGCGATGAAGt                  >  1:684968/1‑148 (MQ=255)
                                                                                                                            tACCGGCTAACGTCCCTGGCGATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAATTGTTGAAAACTATCGCGGTGGCGATATTGCGCTGGGTATCGGCGATGAAGt                  <  1:2519060/148‑1 (MQ=255)
                                                                                                                             aCCGGCTAACGTCCCTGGCGATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAATTGTTGAAAACTATCGCGGTGGCGATATTGCGCTGGGTATCGGCGATGAAGtt                 >  1:2707016/1‑148 (MQ=255)
                                                                                                                              ccGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAATTGTTGAAAACTATCGCGGTGGCGATATTGCGCTGGGTATCGGCGATGAAGttt                <  2:2773171/148‑1 (MQ=255)
                                                                                                                               cGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAATTGTTGAAAACTATCGCGGTGGCGATATTGCGCTGGGTATCGGCGATGAAGtttt               <  2:2680567/148‑1 (MQ=255)
                                                                                                                               cGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAATTGTTGAAAACTATCGCGGTGGCGATATTGCGCTGGGTATCGGCGATGAAGtttt               <  1:261142/148‑1 (MQ=255)
                                                                                                                                ggCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAATTGTTGAAAACTATCGCGGTGGCGATATTGCGCTGGGTATCGGCGATGAAGTTTTa              <  1:2754407/148‑1 (MQ=255)
                                                                                                                                 gCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAATTGTTGAAAACTATCGCGGTGGCGATATTGCGCTGGGTATCGGCGATGAAGTTTTAt             >  2:1119757/1‑148 (MQ=255)
                                                                                                                                  cTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAATTGTTGAAAACTATCGCGGTGGCGATATTGCGCTGGGTATCGGCGATGAAGTTTTATc            <  1:1560950/148‑1 (MQ=255)
                                                                                                                                   tAACGTCCCTGGCGATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAATTGTTGAAAACTATCGCGGTGGCGATATTGCGCTGGGTATCGGCGATGAAGTTTTATCa           >  1:3202809/1‑148 (MQ=255)
                                                                                                                                    aaCGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAATTGTTGAAAACTATCGCGGTGGCGATATTGCGCTGGGTATCGGCGATGAAGTTTTATCAc          >  1:1056427/1‑148 (MQ=255)
                                                                                                                                     aCGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAATTGTTGAAAACTATCGCGGTGGCGATATTGCGCTGGGTATCGGCGATGAAGTTTTATCAcc         <  2:1004841/148‑1 (MQ=255)
                                                                                                                                      cGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAATTGTTGAAAACTATCGCGGTGGCGATATTGCGCTGGGTATCGGCGATGAAGTTTTATCACCg        >  2:37235/1‑148 (MQ=255)
                                                                                                                                       gTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAATTGTTGAAAACTATCGCGGTGGCGATATTGCGCTGGGTATCGGCGATGAAGTTTTATCACCgg       <  1:3280761/148‑1 (MQ=255)
                                                                                                                                         cccTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGACAAAAATGTGAATCCGCAAATTTTTGAAAACTATCGCGGTGGCGATATTGCGCTGGGTATCGGCGCTGTAGTTTTATCACCGGtt     >  1:1991837/1‑148 (MQ=255)
                                                                                                                                         cccTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGACAAAAATGTGAATCCGCAAATTGTTGAAAACTATCGCGGTGGCGATATTGCGCTGGGTATAGGCGATGAAGTTTTATCACCGGtt     >  1:1991832/1‑148 (MQ=255)
                                                                                                                                           cTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAATTGTTGAAAACTATCGCGGTGGCGATATTGCGCTGGGTATCGGCGATGAAGTTTTATCACCGGTTAt   >  2:2326813/1‑148 (MQ=255)
                                                                                                                                            tGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAATTGTTGAAAACTATCGCGGTGGCGATATTGCGCTGGGTATCGGCGATGAAGTTTTATCACCGGTTATg  >  2:503631/1‑148 (MQ=255)
                                                                                                                                            tGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGAGGCAAAAATGTGAATCCGCAAATTGTTGAAAACTATCCCTGTGGCGATATTGCGCTGGGAATAGGAGATGAAGTTTTATCACCGCTTAtt  >  2:2594942/1‑147 (MQ=255)
                                                                                                                                            tGGCTATATCCCGGCGATTGGCTTTATTTATCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATTTGAATCCGCAAATTGTTGAAAAATATCGCGGTGGCGATATTTCGCTGGGTATCGGCGATGAAGTTTTTTAAACGGTTATg  >  1:1990433/1‑148 (MQ=255)
                                                                                                                                                |                                                                                                                                               
GGGGTGAGACAGGTTCCCAGCACGGAAGGCCAATGGAAGTTATTGCATCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCGATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAATTGTTGAAAACTATCGCGGTGGCGATATTGCGCTGGGTATCGGCGATGAAGTTTTATCACCGGTTATG  >  CP009273/1181360‑1181647

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