Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,230,857 C→G 86.8% A63G (GCC→GGC)  ygjR → putative NAD(P)‑dependent dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,230,8570CG86.8% 737.1 / 90.5 218A63G (GCC→GGC) ygjRputative NAD(P)‑dependent dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (15/14);  new base G (123/66);  total (138/80)
Fisher's exact test for biased strand distribution p-value = 2.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 218 total aligned reads displayed.

CAGTTCGTCGAGGCCGCCCATGAGAGCGGTAAATACAAGTTAACCGCCGTATATTCCCGCAGCCTTGAACAGGCCCAGCACTTCGCCAATGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGCCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCATTTTTCCCAGACACAACTTTTCCTTAGCCATAAAATTAATGTGATTTGCGAGAAACCGCTGGCGTCGAATCTGGCGGAAGTGGATGCCGCCATTGCCTGTGCGCG  >  CP009273/3230712‑3231001
                                                                                                                                                 |                                                                                                                                                
cAGTTCGTCGAGGCCGCCCATGAGAGCGGTAAATCCCCGTACCCCGCCGTATATTCCCGCAGCCTTGAACAGGCCCAGCACTTCGCCAATGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCa                                                                                                                                                <  1:3046282/148‑1 (MQ=255)
cAGTTCGTCGAGGCCGCCCATGAGAGCGGTAAATACAAGTTAACCGCCGTATATTCCCGCAGCCTTGAACAGGCCCAGCACTTCGCCAATGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCa                                                                                                                                                <  1:797051/148‑1 (MQ=255)
    tcgtcgAGGCCGCCCATGAGAGCGGTAAATACAAGTTAACCGCCGTATATTCCCGCAGCCTTGAACAGGCCCAGCACTTCGCCAATGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGa                                                                                                                                            <  1:2692747/148‑1 (MQ=255)
       tcgAGGCCGCCCATGAGAGCGGTAAATACAAGTTAACCGCCGTATATTCCCGCAGCCTTGAACAGGCCCAGCACTTCGCCAATGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGCCATTGAcgc                                                                                                                                         <  2:248112/148‑1 (MQ=255)
            gccgccCATGAGAGCGGTAAATACAAGTTAACCGCCGTATATTCCCGCAGCCTTGAACAGGCCCAGCACTTCGCCAATGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGgtgt                                                                                                                                    >  2:1716102/1‑148 (MQ=255)
              cgccCATGAGAGCGGTAAATACAAGTTAACCGCCGTATATTCCCGCAGCCTTGAACAGGCCCAGCACTTCGCCAATGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGCCATTGACGCGGTGtat                                                                                                                                  <  2:511089/148‑1 (MQ=255)
                  cATGAGAGCGGTAAATACAAGTTAACCGCCGTATATTCCCGCAGCCTTGAACAGGCCCAGCACTTCGCCAATGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTg                                                                                                                              >  2:609987/1‑148 (MQ=255)
                   aTGAGAGCGGTAAATACAAGTTAACCGCCGTATATTCCCGCAGCCTTGAACAGGCCCAGCACTTCGCCAATGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGc                                                                                                                             <  1:2966959/148‑1 (MQ=255)
                    tGAGAGCGGTAAATACAAGTTAACCGCCGTATATTCCCGCAGCCTTGAACAGGCCCAGCACTTCGCCAATGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGcc                                                                                                                            <  2:1559378/148‑1 (MQ=255)
                     gagaGCGGTAAATACAAGTTAACCGCCGTATATTCCTGCAGCCTTGAACAGGCCCAGCACTTCGCCAATGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCa                                                                                                                           <  2:1503817/148‑1 (MQ=255)
                          cGGTAAATACAAGTTAACCGCCGTATATTCCCGCAGCCTTGAACAGGCCCAGCACTTCGCCAATGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCg                                                                                                                      >  1:362302/1‑148 (MQ=255)
                          cGGTAAATACAAGTTAACCGCCGTATATTCCCGCAGCCTTGAACAGGCCCAGCACTTCGCCAATGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCg                                                                                                                      >  2:253458/1‑148 (MQ=255)
                           ggTAAATACAAGTTAACCGCCGTATATTCCCGCAGCCTTGAACAGGCCCAGCACTTCGCCAATGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGa                                                                                                                     <  2:1619562/148‑1 (MQ=255)
                            gTAAATACAAGTTAACCGCCGTATATTCCCGCAGCCTTGAACAGGCCCAGCACTTCGCCAATGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGaa                                                                                                                    >  2:1145669/1‑148 (MQ=255)
                              aaaTACAAGTTAACCGCCGTATATTCCCGCAGCCTTGAACAGGCCCAGCACTTCGCCAATGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGCCATTGACGCGGTGTATATTGCCAGCCCGAAtt                                                                                                                  <  2:1716671/148‑1 (MQ=255)
                                aTACAAGTTAACCGCCGTATATTCCCGCAGCCTTGAACAGGCCCAGCACTTCGCCAATGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTcc                                                                                                                >  2:3065049/1‑148 (MQ=255)
                                aTACAAGTTAACCGCCGTATATTCCCGCAGCCTTGAACAGGCCCAGCACTTCGCCAATGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTcc                                                                                                                >  1:3098694/1‑148 (MQ=255)
                                aTACAAGTTAACCGCCGTATATTCCCGCAGCCTTGAACAGGCCCAGCACTTCGCCAATGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTcc                                                                                                                >  1:1275824/1‑148 (MQ=255)
                                  aCAAGTTAACCGCCGTATATTCCCGCAGCCTTGAACAGGCCCAGCACTTCGCCAATGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCt                                                                                                              >  2:1078771/1‑148 (MQ=255)
                                    aaGTTAACCGCCGTATATTCCCGCAGCCTTGAACAGGCCCAGCACTTCGCCAATGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGc                                                                                                            <  2:2639855/148‑1 (MQ=255)
                                    aaGTTAACCGCCGTATATTCCCGCAGCCTTGAACAGGCCCAGCACTTCGCCAATGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGCCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGc                                                                                                            <  2:3233401/148‑1 (MQ=255)
                                    aaGTTAACCGCCGTATATTCCCGCAGCCTTGAACAGGCCCAGCACTTCGCCAATGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGCCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGc                                                                                                            <  2:1593671/148‑1 (MQ=255)
                                     aGTTAACCGCCGTATATTCCCGCAGCCTTGAACAGGCCCAGCACTTCGCCAATGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCa                                                                                                           >  2:2082317/1‑148 (MQ=255)
                                     aGTTAACCGCCGTATATTCCCGCAGCCTTGAACAGGCCCAGCACTTCGCCAATGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCa                                                                                                           >  2:2080389/1‑148 (MQ=255)
                                      gTTAACCGCCGTATATTCCCGCAGCCTTGAACAGGCCCAGCACTTCGCCAATGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGCCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCAt                                                                                                          <  2:2141591/148‑1 (MQ=255)
                                       ttAACCGCCGTATATTCCCGCAGCCTTGAACAGGCCCAGCACTTCGCCAATGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCAtt                                                                                                         <  1:2278787/148‑1 (MQ=255)
                                       ttAACCGCCGTATATTCCCGCAGCCTTGAACAGGCCCAGCACTTCGCCAATGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCAtt                                                                                                         <  1:826354/148‑1 (MQ=255)
                                        tAACCGCCGTATATTCCCGCAGCCTTGAACAGGCCCAGCACTTCGCCAATGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCAttt                                                                                                        <  2:2792973/148‑1 (MQ=255)
                                        tAACCGCCGTATATTCCCGCAGCCTTGAACAGGCCCAGCACTTCGCCAATGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCAttt                                                                                                        <  2:2425033/148‑1 (MQ=255)
                                          aCCGCCGTATATTCCCGCAGCCTTGAACAGGCCCAGCACTTCGCCAATGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCAttttt                                                                                                      >  1:187372/1‑148 (MQ=255)
                                              ccgTATATTCCCGCAGCCTTGAACAGGCCCAGCACTTCGCCAATGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGCCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCATTTTTCCCa                                                                                                  >  2:1224376/1‑148 (MQ=255)
                                               cgTATATTCCCGCAGCCTTGAACAGGCCCAGCACTTCGCCAATGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCATTTTTCCCAg                                                                                                 >  1:465155/1‑148 (MQ=255)
                                                gTATATTCCCGCAGCCTTGAACAGGCCCAGCACTTCGCCAATGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCATTTTTCCCAGa                                                                                                <  2:1683994/148‑1 (MQ=255)
                                                 tataTTCCCGCAGCCTTGAACAGGCCCAGCACTTCGCCAATGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCATTTTTCCCAGac                                                                                               >  1:1320625/1‑148 (MQ=255)
                                                  atatTCCCGCAGCCTTGAACAGGCCCAGCACTTCGCCAATGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCATTTTTCCCAGaca                                                                                              <  2:2484644/148‑1 (MQ=255)
                                                  atatTCCCGCAGCCTTGAACAGGCCCAGCACTTCGCCAATGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCATTTTTCCCAGaca                                                                                              >  2:758559/1‑148 (MQ=255)
                                                     ttCCCGCAGCCTTGAACAGGCCCAGCACTTCGCCAATGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCATTTTTCCCAGACACaa                                                                                           >  2:1001504/1‑148 (MQ=255)
                                                      tCCCGCAGCCTTGAACAGGCCCAGCACTTCGCCAATGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCATTTTTCCCAGACACAAc                                                                                          >  1:112706/1‑148 (MQ=255)
                                                      tCCCGCAGCCTTGAACAGGCCCAGCACTTCGCCAATGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGCCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCATTTTTCCCAGACACAAc                                                                                          >  2:2035099/1‑148 (MQ=255)
                                                        ccGCAGCCTTGAACAGGCCCAGCACTTCGCCAATGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCATTTTTCCCAGACACAACtt                                                                                        <  2:2839884/148‑1 (MQ=255)
                                                        ccGCAGCCTTGAACAGGCCCAGCACTTCGCCAATGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCATTTTTCCCAGACACAACtt                                                                                        <  2:1354303/148‑1 (MQ=255)
                                                        ccGCAGCCTTGAACAGGCCCAGAACTTCGCCAATGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCATTTTTCCCAGACACAACtt                                                                                        <  2:2975525/148‑1 (MQ=255)
                                                          gCAGCCTTGAACAGGCCCAGCACTTCGCCAATGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCATTTTTCCCAGACACAACtttt                                                                                      <  1:2534246/148‑1 (MQ=255)
                                                            aGCCTTGAACAGGCCCAGCACTTCGCCAATGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCATTTTTCCCAGACACAACTTTTcc                                                                                    >  2:1636540/1‑148 (MQ=255)
                                                              ccTTGAACAGGCCCAGCACTTCGCCAATGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCATTTTTCCCAGACACAACTTTTCCtt                                                                                  >  2:365073/1‑148 (MQ=255)
                                                               cTTGAACAGGCCCAGCACTTCGCCAATGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCATTTTTCCCAGACACAACTTTTCCTTa                                                                                 >  1:3144208/1‑148 (MQ=255)
                                                                       ggCCCAGCACTTCGCCAATGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCATTTTTCCCAGACACAACTTTTCCTTAGCCATaaa                                                                         <  2:2184544/148‑1 (MQ=255)
                                                                       ggCCCAGCACTTCGCCAATGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCATTTTTCCCAGACACAACTTTTCCTTAGCCATaaa                                                                         >  2:179345/1‑148 (MQ=255)
                                                                        gCCCAGCACTTCGCCAATGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCATTTTTCCCAGACACAACTTTTCCTTAGCCATaaaa                                                                        >  2:681157/1‑148 (MQ=255)
                                                                        gCCCAGCACTTCGCCAATGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCATTTTTCCCAGACACAACTTTTCCTTAGCCATaaaa                                                                        >  2:2550960/1‑148 (MQ=255)
                                                                        gCCCAGCACTTCGCCAATGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCATTTTTCCCAGACACAACTTTTCCTTAGCCAGaaaa                                                                        >  2:682398/1‑148 (MQ=255)
                                                                         cccAGCACTTCGCCAATGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCATTTTTCCCAGACACAACTTTTCCTTAGCCATAAAAt                                                                       >  2:1021368/1‑148 (MQ=255)
                                                                         cccAGCACTTCGCCAATGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCATTTTTCCCAGACACAACTTTTCCTTAGCCATAAAAt                                                                       >  2:3043964/1‑148 (MQ=255)
                                                                             gCACTTCGCCAATGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGCCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCATTTTTCCCAGACACAACTTTTCCTTAGCCATAAAATTAAt                                                                   >  1:2472069/1‑148 (MQ=255)
                                                                              cACTTCGCCAATGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCATTTTTCCCAGACACAACTTTTCCTTAGCCATAAAATTAAtg                                                                  <  1:1716109/148‑1 (MQ=255)
                                                                              cACTTCGCCAATGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCATTTTTCCCAGACACAACTTTTCCTTAGCCATAAAATTAAtg                                                                  <  2:2464345/148‑1 (MQ=255)
                                                                              cACTTCGCCAATGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCATTTTTCCCAGACACAACTTTTCCTTAGCCATAAAATTAAtg                                                                  >  2:2286845/1‑148 (MQ=255)
                                                                              cACTTCGCCAATGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCATTTTTCCCAGACACAACTTTTCCTTAGCCATAAAATTAAtg                                                                  >  2:2381533/1‑148 (MQ=255)
                                                                                 ttCGCCAATGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGTCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCATTTTTCCCAGACACAACTTTTCCTTAGCCATAAAATTAATGTGa                                                               >  1:1770054/1‑148 (MQ=255)
                                                                                 ttCGCCAATGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCATTTTTCCCAGACACAACTTTTCCTTAGCCATAAAATTAATGTGa                                                               >  2:2362432/1‑148 (MQ=255)
                                                                                 ttCGCCAATGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCATTTTTCCCAGACACAACTTTTCCTTAGCCATAAAATTAATGTGa                                                               >  1:1665029/1‑148 (MQ=255)
                                                                                   cGCCAATGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCATTTTTCCCAGACACAACTTTTCCTTAGCCATAAAATTAATGTGAtt                                                             <  1:2930927/148‑1 (MQ=255)
                                                                                      cAATGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCATTTTTCCCAGACACAACTTTTCCTTAGCCATAAAATTAATGTGATTTGc                                                          >  1:453533/1‑148 (MQ=255)
                                                                                        aTGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCATTTTTCCCAGACACAACTTTTCCTTAGCCATAAAATTAATGTGATTTGCga                                                        >  2:2607015/1‑148 (MQ=255)
                                                                                        aTGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCATTTTTCCCAGACACAACTTTTCCTTAGCCATAAAATTAATGTGATTTGCga                                                        >  2:1427139/1‑148 (MQ=255)
                                                                                         tGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCATTTTTCCCAGACACAACTTTTCCTTAGCCATAAAATTAATGTGATTTGCgag                                                       >  1:1540069/1‑148 (MQ=255)
                                                                                           aTTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCATTTTTCCCAGACACAACTTTTCCTTAGCCATAAAATTAATGTGATTTGCGAGaa                                                     >  2:3306331/1‑148 (MQ=255)
                                                                                            tttttCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCATTTTTCCCAGACACAACTTTTCCTTAGCCATAAAATTAATGTGATTTGCGAGaaa                                                    >  1:2069680/1‑148 (MQ=255)
                                                                                            tttttCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGCCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCATTTTTCCCAGACACAACTTTTCCTTAGCCATAAAATTAATGTGATTTGCGAGaaa                                                    <  2:1145133/148‑1 (MQ=255)
                                                                                             ttttCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCATTTTTCCCAGACACAACTTTTCCTTAGCCATAAAATTAATGTGATTTGCGAGAAAc                                                   >  1:2198912/1‑148 (MQ=255)
                                                                                                tCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCATTTTTCCCAGACACAACTTTTCCTTAGCCATAAAATTAATGTGATTTGCGAGAAACCGc                                                <  1:3128958/148‑1 (MQ=255)
                                                                                                    tCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCATTTTTCCCAGACACAACTTTTCCTTAGCCATAAAATTAATGTGATTTGCGAGAAACCGCTGGc                                            <  2:1320625/148‑1 (MQ=255)
                                                                                                       aGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCATTTTTCCCAGACACAACTTTTCCTTAGCCATAAAATTAATGTGATTTGCGAGAAACCGCTGGCGTc                                         >  2:1949820/1‑148 (MQ=255)
                                                                                                        gCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCATTTTTCCCAGACACAACTTTTCCTTAGCCATAAAATTAATGTGATTTGCGAGAAACCGCTGGCGTCg                                        >  1:896914/1‑148 (MQ=255)
                                                                                                        gCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCATTTTTCCCAGACACAACTTTTCCTTAGCCATAAAATTAATGTGATTTGCGAGAAACCGCTGGCGTCg                                        >  1:3036566/1‑148 (MQ=255)
                                                                                                         cATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCATTTTTCCCAGACACAACTTTTCCTTAGCCATAAAATTAATGTGATTTGCGAGAAACCGCTGGCGTCGa                                       >  1:1708173/1‑148 (MQ=255)
                                                                                                         cATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCATTTTTCCCAGACACAACTTTTCCTTAGCCATAAAATTAATGTGATTTGCGAGAAACCGCTGGCGTCGa                                       >  1:1974004/1‑148 (MQ=255)
                                                                                                           tCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCATTTTTCCCAGACACAACTTTTCCTTAGCCATAAAATTAATGTGATTTGCGAGAAACCGCTGGCGTCGAAt                                     <  1:1090019/148‑1 (MQ=255)
                                                                                                           tCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCATTTTTCCCAGACACAACTTTTCCTTAGCCATAAAATTAATGTGATTTGCGAGAAACCGCTGGCGTCGAAt                                     >  2:488468/1‑148 (MQ=255)
                                                                                                            cTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCATTTTTCCCAGACACAACTTTTCCTTAGCCATAAAATTAATGTGATTTGCGAGAAACCGCTGGCGTCGAATc                                    <  1:1762371/148‑1 (MQ=255)
                                                                                                               tttACCTCGCTGGAAGCGATGGCGGAAAGCGATGCCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCATTTTTCCCAGACACAACTTTTCCTTAGCCATAAAATTAATGTGATTTGCGAGAAACCGCTGGCGTCGAATCTgg                                 >  1:426626/1‑148 (MQ=255)
                                                                                                                   ccTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCATTTTTCCCAGACACAACTTTTCCTTAGCCATAAAATTAATGTGATTTGCGAGAAACCGCTGGCGTCGAATCTGGCGGa                             >  1:2278209/1‑148 (MQ=255)
                                                                                                                   ccTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCATTTTTCCCAGACACAACTTTTCCTTAGCCATAAAATTAATGTGATTTGCGAGAAACCGCTGGCGTCGAATCTGGCGGa                             <  2:3144208/148‑1 (MQ=255)
                                                                                                                   ccTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCATTTTTCCCAGACACAACTTTTCCTTAGCCATAAAATTAATGTGATTTGCGAGAAACCGCTGGCGTCGAATCTGGCGGa                             >  2:614572/1‑148 (MQ=255)
                                                                                                                   ccTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCATTTTTCCCAGACACAACTTGTCCTTAGCCATAAAATTAATGTGATTTGCGAGACACCGCTGGCGTCGAATCTGGGGGa                             >  2:2031198/1‑148 (MQ=255)
                                                                                                                          ggAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCATTTTTCCCAGACACAACTTTTCCTTAGCCATAAAATTAATGTGATTTGCGAGAAACCGCTGGCGTTGAATCTGGCGGAAGTGGAt                      <  2:721300/148‑1 (MQ=255)
                                                                                                                          ggAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCATTTTTCCCAGACACAACTTTTCCTTAGCCATAAAATTAATGTGATTTGCGAGAAACCGCTGGCGTCGAATCTGGCGGAAGTGGAt                      >  1:2173411/1‑148 (MQ=255)
                                                                                                                            aaGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCATTTTTCCCAGACACAACTTTTCCTTAGCCATAAAATTAATGTGATTTGCGAGAAACCGCTGGCGTCGAATCTGGCGGAAGTGGATgc                    >  1:2534849/1‑148 (MQ=255)
                                                                                                                            aaGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCATTTTTCCCAGACACAACTTTTCCTTAGCCATAAAATTAATGTGATTTGCGAGAAACCGCTGGCGTCGAATCTGGCGGAAGTGGATgc                    <  2:1769423/148‑1 (MQ=255)
                                                                                                                                 aTGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCATTTTTCCCAGACACAACTTTTCCTTAGCCATAAAATTAATGTGATTTGCGAGAAACCGCTGGCGTCGAATCTGGCGGAAGTGGATGCCGCCa               >  1:547399/1‑148 (MQ=255)
                                                                                                                                    gCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCATTTTTCCCAGACACAACTTTTCCTTAGCCATAAAATTAATGTGATTTGCGAGAAACCGCTGGCGTCGAATCTGGCGGAAGTGGATGCCGCCATTg            >  2:2264236/1‑148 (MQ=255)
                                                                                                                                     cGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCATTTTTCCCAGACACAACTTTTCCTTAGCCATAAAATTAATGTGATTTGCGAGAAACCGCTGGCGTCGAATCTGGCGGAAGTGGATGCCGCCATTGc           <  1:2884940/148‑1 (MQ=255)
                                                                                                                                     cGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCATTTTTCCCAGACACAACTTTTCCTTAGCCATAAAATTAATGTGATTTGCGAGAAACCGCTGGCGTCGAATCTGGCGGAAGTGGATGCCGCCATTGc           <  2:1651886/148‑1 (MQ=255)
                                                                                                                                       gAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCATTTTTCCCAGACACAACTTTTCCTTAGCCATAAAATTAATGTGATTTGCGAGAAACCGCTGGCGTCGAATCTGGCGGAAGTGGATGCCGCCATTGCCt         >  1:1008577/1‑148 (MQ=255)
                                                                                                                                        aaaGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCATTTTTCCCAGACACAACTTTTCCTTAGCCATAAAATTAATGTGATTTGCGAGAAACCGCTGGCGTCGAATCTGGCGGAAGTGGATGCCGCCATTGCCtg        >  1:2514976/1‑148 (MQ=255)
                                                                                                                                        aaaGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCATTTTTCCCAGACACAACTTTTCCTTAGCCATAAAATTAATGTGATTTGCGAGAAACCGCTGGCGTCGAATCTGGCGGAAGTGGATGCCGCCATTGCCtg        >  1:3288566/1‑148 (MQ=255)
                                                                                                                                        aaaGCGATGGCATTGACGCGGTGTAGATTGCCAGCCCGAATTCCCTGCATTTTTCCCAGACACAACTTTTCCTTAGCCATAAAATTAATGTGATTTGCGAGAAACCGCTGGCGTCGAATCTGGCGGAAGTGGATGCCGCCATTGCCtg        >  1:2829894/1‑148 (MQ=255)
                                                                                                                                            cGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCATTTTTCCCAGACACAACTTTTCCTTAGCCATAAAATTAATGTGATTTGCGAGAAACCGCTGGCGTCGAATCTGGCGGAAGTGGATGCCGCCATTGCCTGTgcg    >  2:3023780/1‑148 (MQ=255)
                                                                                                                                              aTGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCATTTTTCCCAGACACAACTTTTCCTTAGCCATAAAATTAATGTGATTTGCGAGAAACCGCTGGCGTCGAATCTGGCGGACGTGGATGCCGCCATTGCCTGTGcgcg  >  2:2444965/1‑148 (MQ=255)
                                                                                                                                              aTGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCATTTTTCCCAGACACAACTTTTCCTTAGCCATAAAATTAATGTGATTTGCGAGAAACCGCTGGCGTCGAATCTGGCGGAAGTGGATGCCGCCATTGCCTGTGcgcg  >  1:2117634/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
CAGTTCGTCGAGGCCGCCCATGAGAGCGGTAAATACAAGTTAACCGCCGTATATTCCCGCAGCCTTGAACAGGCCCAGCACTTCGCCAATGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGCCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCATTTTTCCCAGACACAACTTTTCCTTAGCCATAAAATTAATGTGATTTGCGAGAAACCGCTGGCGTCGAATCTGGCGGAAGTGGATGCCGCCATTGCCTGTGCGCG  >  CP009273/3230712‑3231001

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: