Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,230,857 C→G 100% 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,8570CG97.0% 918.3 / 14.0 234A63G (GCC→GGC) ygjRputative NAD(P)‑dependent dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (2/5);  new base G (112/115);  total (114/120)
Fisher's exact test for biased strand distribution p-value = 4.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.74e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 234 total aligned reads displayed.

CAGTTCGTCGAGGCCGCCCATGAGAGCGGTAAATACAAGTTAACCGCCGTATATTCCCGCAGCCTTGAACAGGCCCAGCACTTCGCCAATGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGCCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCATTTTTCCCAGACACAACTTTTCCTTAGCCATAAAATTAATGTGATTTGCGAGAAACCGCTGGCGTCGAATCTGGCGGAAGTGGATGCCGCCATTGCCTGTGCG  >  CP009273/3230712‑3230999
                                                                                                                                                 |                                                                                                                                              
cAGTTCGTCGAGGCCGCCCATGAGAGCGGTAAATACAAGTTAACCGCCGTATATTCCCGCAGCCTTGAACAGGCCCAGCACTTCGCCAATGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCa                                                                                                                                              <  1:442646/148‑1 (MQ=255)
   ttCGTCGAGGCCGCCCATGAGAGCGGTAAATACAAGTTAACCGCCGTATATTCCCGCAGCCTTGAACAGGCCCAGCACTTCGCCAATGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTg                                                                                                                                           >  1:1206624/1‑148 (MQ=255)
   ttCGTCGAGGCCGCCCATGAGAGCGGTAAATACAAGTTAACCGCCGTATATTCCCGCAGCCTTGAACAGGCCCAGCACTTCGCCAATGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTg                                                                                                                                           >  2:688818/1‑148 (MQ=255)
   ttCGTCGAGGCCGCCCATGAGAGCGGTAAATACAAGTTAACCGCCGTATATTCCCGCAGCCTTGAACAGGCCCAGCACTTCGCCAATGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTg                                                                                                                                           <  1:2842516/148‑1 (MQ=255)
   ttCGTCGAGGCCGCCCATGAGAGCGGTAAATACAAGTTAACCGCCGTATATTCCCGCAGCCTTGAACAGGCCCAGCACTTCGCCAATGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTg                                                                                                                                           <  2:2935372/148‑1 (MQ=255)
   ttCGTCGAGGCCGCCCATGAGAGCGGTAAATACAAGTTAACCGCCGTATATTCCCGCAGCCTTGAACAGGCCCAACACTTCGCCAATGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTg                                                                                                                                           >  1:3376360/1‑148 (MQ=255)
       tcgAGGCCGCCCATGAGAGCGGTAAATACAAGTTAACCGCCGTATATTCCCGCAGCCTTGAACAGGCCCAGCACTTCGCCAATGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGAcgc                                                                                                                                       >  2:612133/1‑148 (MQ=255)
        cgAGGCCGCCCATGAGAGCGGTAAATACAAGTTAACCGCCGTATATTCCCGCAGCCTTGAACAGGCCCAGCACTTCGCCAATGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGAcgcg                                                                                                                                      >  1:3156311/1‑148 (MQ=255)
        cgAGGCCGCCCATGAGAGCGGTAAATACAAGTTAACCGCCGTATATTCCCGCAGCCTTGAACAGGCCCAGCACTTCGCCAATGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGAcgcg                                                                                                                                      <  1:3187767/148‑1 (MQ=255)
          aGGCCGCCCATGAGAGCGGTAAATACAAGTTAACCGCCGTATATTCCCGCAGCCTTGAACAGGCCCAGCACTTCGCCAATGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGgt                                                                                                                                    <  2:1472772/148‑1 (MQ=255)
          aGGCCGCCCATGAGAGCGGTAAATACAAGTTAACCGCCGTATATTCCCGCAGCCTTGAACAGGCCCAGCACTTCGCCAATGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGgt                                                                                                                                    <  1:2619984/148‑1 (MQ=255)
           ggCCGCCCATGAGAGCGGTAAATACAAGTTAACCGCCGTATATTCCCGCAGCCTTGAACAGGCCCAGCACTTCGCCAATGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGgtg                                                                                                                                   >  1:1530756/1‑148 (MQ=255)
           ggCCGCCCATGAGAGCGGTAAATACAAGTTAACCGCCGTATATTCCCGCAGCCTTGAACAGGCCCAGCACTTCGCCAATGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGgtg                                                                                                                                   >  2:685186/1‑148 (MQ=255)
           ggCCGCCCATGAGAGCGGTAAATACAAGTTAACCGCCGTATATTCCCGCAGCCTTGAACAGGCCCAGCACTTCGCCAATGATTTTTCTGTCAAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGCCATTGACGCGgtg                                                                                                                                   >  2:3090037/1‑148 (MQ=255)
            gccgccCATGAGAGCGGTAAATACAAGTTAACCGCCGTATATTCCCGCAGCCTTGAACAGGCCCAGCACTTCGCCAATGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGgtgt                                                                                                                                  <  2:795225/148‑1 (MQ=255)
               gccCATGAGAGCGGTAAATACACGTTAACCGCCGTATATTCCCGCAGCCTTGAACAGGCCCAGCACTTCGCCAATGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGtata                                                                                                                               <  1:1708689/148‑1 (MQ=255)
                cccATGAGAGCGGTAAATACAAGTTAACCGCCGTATATTCCCGCAGCCTTGAACAGGCCCAGCACTTCGCCAATGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTatat                                                                                                                              <  2:1516651/148‑1 (MQ=255)
                  cATGAGAGCGGTAAATACAAGTTAACCGCCGTATATTCCCGCAGCCTTGAACAGGCCCAGCACTTCGCCAATGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTg                                                                                                                            >  2:2010548/1‑148 (MQ=255)
                  cATGAGAGCGGTAAATACAAGTTAACCGCCGTATATTCCCGCAGCCTTGAACAGGCCCAGCACTTCGCCAATGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTg                                                                                                                            <  1:1574339/148‑1 (MQ=255)
                    tGAGAGCGGTAAATACAAGTTAACCGCCGTATATTCCCGCAGCCTTGAACAGGCCCAGCACTTCGCCAATGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGcc                                                                                                                          <  1:501065/148‑1 (MQ=255)
                       gagCGGTAAATACAAGTTAACCGCCGTATATTCCCGCAGCCTTGAACAGGCCCAGCACTTCGCCAATGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGGCGCGGTGTATATTGCCAGc                                                                                                                       >  1:3456398/1‑148 (MQ=255)
                       gagCGGTAAATACAAGTTAACCGCCGTATATTCCCGCAGCCTTGAACAGGCCCAGCACTTCGCCAATGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGc                                                                                                                       >  2:185296/1‑148 (MQ=255)
                        agCGGTAAATACAAGTTAACCGCCGTATATTCCCGCAGCCTTGAACAGGCCCAGCACTTCGCCAATGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCTTTGACGCGGTGTATATTGCCAGcc                                                                                                                      <  2:939430/148‑1 (MQ=255)
                            gTAAATACAAGTTAACCGCCGTATATTCCCGCAGCCTTGAACAGGCCCAGCACTTCGCCAATGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGaa                                                                                                                  >  2:445098/1‑148 (MQ=255)
                             tAAATACAAGTTAACCGCCGTATATTCCCGCAGCCTTGAACAGGCCCAGCACTTCGCCAATGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAAt                                                                                                                 >  1:2523933/1‑148 (MQ=255)
                             tAAATACAAGTTAACCGCCGTATATTCCCGCAGCCTTGAACAGGCCCAGCACTTCGCCAATGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAAt                                                                                                                 >  1:350113/1‑148 (MQ=255)
                               aaTACAAGTTAACCGCCGTATATTCCCGCAGCCTTGAACAGGCCCAGCACTTCGCCAATGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTc                                                                                                               <  1:1347875/148‑1 (MQ=255)
                                aTACAAGTTAACCGCCGTATATTCCCGCAGCCTTGAACAGGCCCAGCACTTCGCCAATGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATCTTGCCAGCCCGAATTcc                                                                                                              >  2:2292578/1‑148 (MQ=255)
                                  aCAAGTTAACCGCCGTATATTCCCGCAGCCTTGAACAGGCCCAGCACTTCGCCAATGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCt                                                                                                            >  1:425330/1‑148 (MQ=255)
                                     aGTTAACCGCCGTATATTCCCGCAGCCTTGAACAGGCCCAGCACTTCGCCAATGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCa                                                                                                         <  2:1206624/148‑1 (MQ=255)
                                     aGTTAACCGCCGTATATTCCCGCAGCCTTGAACAGGCCCAGCACTTCGCCAATGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCa                                                                                                         >  1:1049039/1‑148 (MQ=255)
                                     aGTTAACCGCCGTATATTCCCGCAGCCTTGAACAGGCCCAGCACTTCGCCAATGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCa                                                                                                         >  1:2828399/1‑148 (MQ=255)
                                       ttAACCGCCGTATATTCCCGCAGCCTTGAACAGGCCCAGCACTTCGCCAATGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCAtt                                                                                                       >  2:19007/1‑148 (MQ=255)
                                       ttAACCGCCGTATATTCCCGCAGCCTTGAACAGGCCCAGCACTTCGCCAATGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCAtt                                                                                                       >  2:2594363/1‑148 (MQ=255)
                                             gccgTATATTCCCGCAGCCTTGAACAGGCCCAGCACTTCGCCAATGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCATTTTTccc                                                                                                 <  1:1453160/148‑1 (MQ=255)
                                             gccgTATATTCCCGCAGCCTTGAACAGGCCCAGCACTTCGCCAATGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGCCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCATTTTTccc                                                                                                 <  1:2811165/148‑1 (MQ=255)
                                             gccgTATAGTCCCGCAGCCTTGAACAGGCCCAGCACTTCGCCAATGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCATTTTTccc                                                                                                 <  1:799934/148‑1 (MQ=255)
                                                  atatTCCCGCAGCCTTGAACAGGCCCAGCACTTCGCCAATGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCATTTTTCCCAGaca                                                                                            >  1:704765/1‑148 (MQ=255)
                                                    atTCCCGCAGCCTTGAACAGGCCCAGCACTTCGCCAATGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCATTTTTCCCAGAcaca                                                                                          >  1:57178/1‑148 (MQ=255)
                                                      tCCCGCAGCCTTGAACAGGCCCAGCACTTCGCCAATGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCATTTTTCCCAGACACAAc                                                                                        <  1:712198/148‑1 (MQ=255)
                                                       cccGCAGCCTTGAACAGGCCCAGCACTTCGCCAATGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCATTTTTCCCAGACACAACt                                                                                       <  2:344167/148‑1 (MQ=255)
                                                        ccGCAGCCTTGAACAGGCCCAGCACTTCGCCAATGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCATTTTTCCCAGACACAACtt                                                                                      <  2:352369/148‑1 (MQ=255)
                                                        ccGCAGCCTTGAACAGGCCCAGCACTTCGCCAATGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCATTTTTCCCAGACACAACtt                                                                                      <  2:1075467/148‑1 (MQ=255)
                                                        ccGCAGCCTTGAACAGGCCCAGCACTTCGCCAATGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCATTTTTCCCAGACACAACtt                                                                                      <  1:1298866/148‑1 (MQ=255)
                                                         cGCAGCCTTGAACAGGCCCAGCACTTCGCCAATGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCATTTTTCCCAGACACAACttt                                                                                     <  1:855530/148‑1 (MQ=255)
                                                           cAGCCTTGAACAGGCCCAGCACTTCGCCAATGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCATTTTTCCCAGACACAACTTTTc                                                                                   <  1:2507079/148‑1 (MQ=255)
                                                           cAGCCTTGAACAGGCCCAGCACTTCGCCAATGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCATTTTTCCCAGACACAACTTTTc                                                                                   <  1:1994054/148‑1 (MQ=255)
                                                            aGCCTTGAACAGGCCCAGCACTTCGCCAATGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCATTTTTCCCAGACACAACTTTTcc                                                                                  <  1:2718938/148‑1 (MQ=255)
                                                                  gAACAGGCCCAGCACTTCGCCAATGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCATTTTTCCCAGACACAACTTTTCCTTAGcc                                                                            <  2:1111617/148‑1 (MQ=255)
                                                                        gCCCAGCACTTCGCCAATGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCATTTTTCCCAGACACAACTTTTCCTTAGCCATaaaa                                                                      <  2:824148/148‑1 (MQ=255)
                                                                            aGCACTTCGCCAATGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCATTTTTCCCAGACACAACTTTTCCTTAGCCATAAAATTaa                                                                  >  2:2317709/1‑148 (MQ=255)
                                                                             gCACTTCGCCAATGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATTGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCATTTTTCCCAGACACAACTTTTCCTTAGCCATAAAATTAAt                                                                 <  2:2748022/148‑1 (MQ=255)
                                                                              cACTTCGCCAATGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCATTTTTCCCAGACACAACTTTTCCTTAGCCATAAAATTAAtg                                                                >  2:2404302/1‑148 (MQ=255)
                                                                                cTTCGCCAATGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCATTTTTCCCAGACACAACTTTTCCTTAGCCATAAAATTAAtgtg                                                              <  1:2686809/148‑1 (MQ=255)
                                                                                cTTCGCCAATGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCATTTTTCCCAGACACAACTTTTCCTTAGCCATAAAATTAAtgtg                                                              <  2:1065324/148‑1 (MQ=255)
                                                                                cTTCGCCAATGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCATTTTTCCCAGACACAACTTTTCCTTAGCCATAAAATTAAtgtg                                                              <  2:138511/148‑1 (MQ=255)
                                                                                    gCCAATGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCATTTTTCCCAGACACAACTTTTCCTTAGCCATAAAATTAATGTGAttt                                                          >  2:1208252/1‑148 (MQ=255)
                                                                                     caattGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCATTTTTCCCAGACACAACTTTTCCTTATCCATAAAATTAATGTGATTTg                                                         <  1:827307/144‑1 (MQ=255)
                                                                                       aaTGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCATTTTTCCCAGACACAACTTTTCCTTAGCCATAAAATTAATGTGATTTGCg                                                       <  1:1237608/148‑1 (MQ=255)
                                                                                       aaTGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCATTTTTCCCAGACACAACTTTTCCTTAGCCATAAAATTAATGTGATTTGCg                                                       >  2:2486380/1‑148 (MQ=255)
                                                                                       aaTGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCATTTTTCCCAGACACAACTTTTCCTTAGCCATAAAATTAATGTGATTTGCg                                                       >  1:3285967/1‑148 (MQ=255)
                                                                                        aTGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCATTTTTCCCAGACACAACTTTTCCTTAGCCATAAAATTAATGTGATTTGCga                                                      >  2:2513575/1‑148 (MQ=255)
                                                                                         tGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCATTTTTCCCAGACACAACTTTTCCTTAGCCATAAAATTAATGTGATTTGCgag                                                     >  1:265189/1‑148 (MQ=255)
                                                                                         tGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCATTTTTCCCAGACACAACTTTTCCTTAGCCATAAAATTAATGTGATTTGCgag                                                     >  2:1669109/1‑148 (MQ=255)
                                                                                         tGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCATTTTTCCCAGACACAACTTTTCCTTAGCCATAAAATTAATGTGATTTGCgag                                                     >  1:1858149/1‑148 (MQ=255)
                                                                                          gATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCATTTTTCCCAGACACAACTTTTCCTTAGCCATAAAATTAATGTGATTTGCgaga                                                    <  1:1708626/148‑1 (MQ=255)
                                                                                          gATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCATTTTTCCCAGACACAACTTTTCCTTAGCCATAAAATTAATGTGATTTGCgaga                                                    <  2:3005878/148‑1 (MQ=255)
                                                                                          gATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCATTTTTCCCAGACACAACTTTTCCTTAGCCATAAAATTAATGTGATTTGCgaga                                                    <  2:1226187/148‑1 (MQ=255)
                                                                                            tttttCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCATTTTTCCCAGACACAACTTTTCCTTAGCCATAAAATTAATGTGATTTGCGAGaaa                                                  >  1:3363364/1‑148 (MQ=255)
                                                                                            tttttCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCATTTTTCCCAGACACAACTTTTCCTTAGCCATAAAATTAATGTGATTTGCGAGaaa                                                  <  2:1495083/148‑1 (MQ=255)
                                                                                             ttttCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCATTTTTCCCAGACACAACTTTTCCTTAGCCATAAAATTAATGTGATTTGCGAGAAAc                                                 >  1:2383560/1‑148 (MQ=255)
                                                                                               ttCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCATTTTTCCCAGACACAACTTTTCCTTAGCCATAAAATTAATGTGATTTGCGAGAAACCg                                               <  1:2111577/148‑1 (MQ=255)
                                                                                               ttCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCATTTTTCCCAGACACAACTTTTCCTTAGCCATAAAATTAATGTGATTTGCGAGAAACCg                                               <  1:3422660/148‑1 (MQ=255)
                                                                                                tCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCATTTTTCCCAGACACAACTTTTCCTTAGCCATAAAATTAATGTGATTTGCGAGAAACCGc                                              <  2:3398801/148‑1 (MQ=255)
                                                                                                      gAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCATTTTTCCCAGACACAACTTTTCCTTAGCCATAAAATTAATGTGATTTGCGAGAAACCGCTGGCGt                                        <  1:1106357/148‑1 (MQ=255)
                                                                                                       aGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCATTTTTCCCAGACACAACTTTTCCTTAGCCATAAAATTAATGTGATTTGCGAGAAACCGCTGGCGTc                                       >  2:2042288/1‑148 (MQ=255)
                                                                                                         cATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCATTTTTCCCAGACACAACTTTTCCTTAGCCATAAAATTAATGTGATTTGCGAGAAACCGCTGGCGTCGa                                     >  1:3098168/1‑148 (MQ=255)
                                                                                                         cATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCATTTTTCCCAGACACAACTTTTCCTTAGCCATAAAATTAATGTGATTTGCGAGAAACCGCTGGCGTCGa                                     >  1:1323729/1‑148 (MQ=255)
                                                                                                          aTCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCATTTTTCCCAGACACAACTTTTCCTTAGCCATAAAATTAATGTGATTTGCGAGAAACCGCTGGCGTCGaa                                    >  2:1165792/1‑148 (MQ=255)
                                                                                                            cTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCATTTTTCCCAGACACAACTTTTCCTTAGCCATAAAATTAATGTGATTTGCGAGAAACCGCTGGCGTCGAATc                                  >  1:2362779/1‑148 (MQ=255)
                                                                                                              gTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCATTTTTCCCAGACACAACTTTTCCTTAGCCATAAAATTAATGTGATTTGCGAGAAACCGCTGGCGTCGAATCTg                                >  2:2811002/1‑148 (MQ=255)
                                                                                                               tttACCTCGCTGGAAGCGATGGCGGAAAGCGATGCCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCATTTTTCCCAGACACAACTTTTCCTTAGCCATAAAATTAATGTGATTTGCGAGAAACCGCTGGCGTCGAATCTgg                               <  2:327185/148‑1 (MQ=255)
                                                                                                                   ccTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCATTTTTCCCAGACACAACTTTTCCTTAGCCATAAAATTAATGTGATTTGCGAGAAACCGCTGGCGTCGAATCTGGCGGa                           >  2:2846620/1‑148 (MQ=255)
                                                                                                                    cTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCATTTTTCCCAGACACAACTTTTCCTTAGCCATAAAATTAATGTGATTTGCGAGAAACCGCTGGCGTCGAATCTGGCGGaa                          <  1:685536/148‑1 (MQ=255)
                                                                                                                    cTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCATTTTTCCCAGACACAACTTTTCCTTAGCCATAAAATTAATGTGATTTGCGAGAAACCGCTGGCGTCGAATCTGGCGGaa                          <  2:3202331/148‑1 (MQ=255)
                                                                                                                     tCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCATTTTTCCCAGACACAACTTTTCCTTAGCCATAAAATTAATGTGATTTGCGAGAAACCGCTGGCGTCGAATCTGGCGGAAg                         <  2:1068926/148‑1 (MQ=255)
                                                                                                                      cGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCATTTTTCCCAGACACAACTTTTCCTTAGCCATAAAATTAATGTGATTTGCGAGAAACCGCTGGCGTCGAATCTGGCGGAAGt                        >  1:402945/1‑148 (MQ=255)
                                                                                                                      cGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCATTTTTCCCAGACACAACTTTTCCTTAGCCATAAAATTAATGTGATTTGCGAGAAACCGCTGGCGTCGAATCTGGCGGAAGt                        <  1:2764556/148‑1 (MQ=255)
                                                                                                                       gCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCATTTTTCCCAGACACAACTTTTCCTTAGCCATAAAATTAATGTGATTTGCGAGAAACCGCTGGCGTCGAATCTGGCGGAAGTg                       <  2:883007/148‑1 (MQ=255)
                                                                                                                          ggAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCATTTTTCCCAGACACAACTTTTCCTTAGCCATAAAATTAATGTGATTTGCGAGAAACCGCTGGCGTCGAATCTGGCGGAAGTGGAt                    >  1:168951/1‑148 (MQ=255)
                                                                                                                              gCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCATTTTTCCCAGACACAACTTTTCCTTAGCCATAAAATTAATGTGATTTGCGAGAAACCGCTGGCGTCGAATCTGGAGGAAGTGGATgccg                >  1:2248526/1‑148 (MQ=255)
                                                                                                                               cGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCATTTTTCCCAGACACAACTTTTCCTTAGCCATAAAATTAATGTGATTTGCGAGAAACCGCTGGCGTCGAATCTGGCGGAAGTGGATgccgc               >  2:3258572/1‑148 (MQ=255)
                                                                                                                               cGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCATTTTTCCCAGACACAACTTTTCCTTAGCCATAAAATTAATGTGATTTGCGAGAAACCGCTGGCGTCGAATCTGGCGGAAGTGGATgccgc               >  2:2957997/1‑148 (MQ=255)
                                                                                                                                 aTGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCATTTTTCCCAGACACAACTTTTCCTTAGCCATAAAATTAATGTGATTTGCGAGAAACCGCTGGCGTCGAATCTGGCGGAAGTGGATGCCGCCa             <  2:2930562/148‑1 (MQ=255)
                                                                                                                                  tGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCATTTTTCCCAGACACAACTTTTCCTTAGCCATAAAATTAATGTGATTTGCGAGAAACCGCTGGCGTCGAATCTGGCGGAAGTGGATGCCGCCAt            <  2:907919/148‑1 (MQ=255)
                                                                                                                                    gCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCATTTTTCCCAGACACAACTTTTCCTTAGCCATAAAATTAATGTGATTTGCGAGAAACCGCTGGCGTCGAATCTGGCGGAAGTGGATGCCGCCATTg          <  2:704865/148‑1 (MQ=255)
                                                                                                                                     cGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCATTTTTCCCAGACACAACTTTTCCTTAGCCATAAAATTAATGTGATTTGCGAGAAACCGCTGGCGTCGAATCTGGCGGAAGTGGATGCCGCCATTGc         <  2:639592/148‑1 (MQ=255)
                                                                                                                                        aaaGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCATTTTTCCCAGACACAACTTTTCCTTAGCCATAAAATTAATGTGATTTGCGAGAAACCGCTGGCGTCGAATCTGGCGGAAGTGGATGCCGCCATTGCCtg      <  1:644756/148‑1 (MQ=255)
                                                                                                                                        aaaGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCATTTTTCCCAGACACAACTTTTCCTTAGCCATAAAATTAATGTGATTTGCGAGAAACCGCTGGCGTCGAATCTGGCGGAAGTGGATGCCGCCATTGCCtg      >  1:1031014/1‑148 (MQ=255)
                                                                                                                                         aaGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCATTTTTCCCAGACACAACTTTTCCTTAGCCATAAAATTAATGTGATTTGCGAGAAACCGCTGGCGTCGAATCTGGCGGAAGTGGATGCCGCCATTGCCtgt     <  2:617083/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                              
CAGTTCGTCGAGGCCGCCCATGAGAGCGGTAAATACAAGTTAACCGCCGTATATTCCCGCAGCCTTGAACAGGCCCAGCACTTCGCCAATGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGCCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCATTTTTCCCAGACACAACTTTTCCTTAGCCATAAAATTAATGTGATTTGCGAGAAACCGCTGGCGTCGAATCTGGCGGAAGTGGATGCCGCCATTGCCTGTGCG  >  CP009273/3230712‑3230999

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: