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,8570CG99.2% 1514.1 / 0.3 388A63G (GCC→GGC) ygjRputative NAD(P)‑dependent dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (3/0);  new base G (223/162);  total (226/162)
Fisher's exact test for biased strand distribution p-value = 2.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 388 total aligned reads displayed.

GCCAGTTCGTCGAGGCCGCCCATGAGAGCGGTAAATACAAGTTAACCGCCGTATATTCCCGCAGCCTTGAACAGGCCCAGCACTTCGCCAATGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGCCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCATTTTTCCCAGACACAACTTTTCCTTAGCCATAAAATTAATGTGATTTGCGAGAAACCGCTGGCGTCGAATCTGGCGGAAGTGGATGCCGCCATTGCCTGTGCGCGGG  >  CP009273/3230710‑3231003
                                                                                                                                                   |                                                                                                                                                  
gCCAGTTCGTCGAGGCCGCCCATGAGAGCGGTAAATACAAGTTAACCGCCGTATATTCCCGCAGCCTTGAACAGGCCCAGCACTTCGCCAATGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCa                                                                                                                                                  <  1:974344/150‑1 (MQ=255)
gCCAGTTCGTCGAGGCCGCCCATGAGAGCGGTAAATACAAGTTAACCGCCGTATATTCCCGCAGCCTTGAACAGGCCCAGCACTTCGCCAATGAATTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCa                                                                                                                                                  >  1:1893781/1‑150 (MQ=255)
   aGTTCGTCGAGGCCGCCCATGAGAGCGGTAAATACAAGTTAACCGCCGTATATTCCCGCAGCCTTGAACAGGCCCAGCACTTCGCCAATGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTg                                                                                                                                               >  1:1112497/1‑150 (MQ=255)
   aGTTCGTCGAGGCCGCCCATGAGAGCGGTAAATACAAGTTAACCGCCGTATATTCCCGCAGCCTTGAACAGGCCCAGCACTTCGCCAATGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTg                                                                                                                                               >  1:486845/1‑150 (MQ=255)
        gtcgAGGCCGCCCATGAGAGCGGTAAATACAAGTTAACCGCCGTATATTCCCGCAGCCTTGAACAGGCCCAGCACTTCGCCAATGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGAcgcg                                                                                                                                          >  1:407111/1‑150 (MQ=255)
        gtcgAGGCCGCCCATGAGAGCGGTAAATACAAGTTAACCGCCGTATATTCCCGCAGCCTTGAACAGGCCCAGCACTTCGCCAATGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGAcgcg                                                                                                                                          <  1:2813145/150‑1 (MQ=255)
          cgAGGCCGCCCATGAGAGCGGTAAATACAAGTTAACCGCCGTATATTCCCGCAGCCTTGAACAGGCCCAGCACTTCGCCAATGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGgt                                                                                                                                        >  1:1274906/1‑150 (MQ=255)
          cgAGGCCGCCCATGAGAGCGGTAAATACAAGTTAACCGCCGTATATTCCCGCAGCCTTGAACAGGCCCAGCACTTCGCCAATGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGgt                                                                                                                                        >  1:3886129/1‑150 (MQ=255)
           gAGGCCGCCCATGAGAGCGGTAAATACAAGTTAACCGCCGTATATTCCCGCAGCCTTGAACAGGCCCAGCACTTCGCCAATGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGgtg                                                                                                                                       >  2:503967/1‑150 (MQ=255)
                cgccCATGAGAGCGGTAAATACAAGTTAACCGCCGTATATTCCCGCAGCCTTGAACAGGCCCAGCACTTCGCCAATGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTatat                                                                                                                                  >  2:918536/1‑150 (MQ=255)
                cgccCATGAGAGCGGTAAATACAAGTTAACCGCCGTATATTCCCGCAGCCTTGAACAGGCCCAGCACTTCGCCAATGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTatat                                                                                                                                  >  1:205950/1‑150 (MQ=255)
                 gccCATGAGAGCGGTAAATACAAGTTAACCGCCGTATATTCCCGCAGCCTTGAACAGGCCCAGCACTTCGCCAATGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATAtt                                                                                                                                 <  1:1074696/150‑1 (MQ=255)
                  cccATGAGAGCGGTAAATACAAGTTAACCGCCGTATATTCCCGCAGCCTTGAACAGGCCCAGCACTTCGCCAATGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTg                                                                                                                                >  1:1706222/1‑150 (MQ=255)
                  cccATGAGAGCGGTAAATACAAGTTAACCGCCGTATATTCCCGCAGCCTTGAACAGGCCCAGCACTTCGCCAATGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGGatatat                                                                                                                                >  2:3034987/1‑148 (MQ=255)
                    cATGAGAGCGGTAAATACAAGTTAACCGCCGTATATTCCCGCAGCCTTGAACAGGCCCAGCACTTCGCCAATGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGcc                                                                                                                              <  2:2763340/150‑1 (MQ=255)
                       gagaGCGGTAAATACAAGTTAACCGCCGTATATTCCCGCAGCCTTGAACAGGCCCAGCACTTCGCCAATGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGc                                                                                                                           <  1:3851971/150‑1 (MQ=255)
                            cGGTAAATACAAGTTAACCGCCGTATATTCCCGCAGCCTTGAACAGGCCCAGCACTTCGCCAATGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGaa                                                                                                                      <  1:2765356/150‑1 (MQ=255)
                                aaaTACAAGTTAACCGCCGTATATTCCCGCAGCCTTGAACAGGCCCAGCACTTCGCCAATGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTcc                                                                                                                  >  2:1514449/1‑150 (MQ=255)
                                  aTACAAGTTAACCGCCGTATATTCCCGCAGCCTTGAACAGGCCCAGCACTTCGCCAATGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCt                                                                                                                <  2:3860603/150‑1 (MQ=255)
                                    aCAAGTTAACCGCCGTATATTCCCGCAGCCTTGAACAGGCCCAGCACTTCGCCAATGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGc                                                                                                              >  1:3553797/1‑150 (MQ=255)
                                      aaGTTAACCGCCGTATATTCCCGCAGCCTTGAACAGGCCCAGCACTTCGCCAATGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCAt                                                                                                            <  2:1943851/150‑1 (MQ=255)
                                       aGTTAACCGCCGTATATTCCCGCAGCCTTGAACAGGCCCAGCACTTCGCCAATGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCAtt                                                                                                           >  2:2554474/1‑150 (MQ=255)
                                        gTTAACCGCCGTATATTCCCGCATCCTTGAACATGCCCAGCACTTCGCCAATGATTTTTCTGTCTAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGATAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCAttt                                                                                                          <  1:1028678/150‑1 (MQ=255)
                                         ttAACCGCCGTATATTCCCGCAGCCTTGAACAGGCCCAGCACTTCGCCAATGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCAtttt                                                                                                         <  2:3352527/150‑1 (MQ=255)
                                         ttAACCGCCGTATATTCCCGCAGCCTTGAACAGGCCCAGCACTTCGCCAATGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCAtttt                                                                                                         <  1:3741690/150‑1 (MQ=255)
                                         ttAACCGCCGTATAGTCCCGCAGCCTTGAACAGGCCCAGCACTTCGCCAATGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCAtttt                                                                                                         <  1:2621156/150‑1 (MQ=255)
                                          tAACCGCCGTATATTCCCGCAGCCTTGAACAGGCCCAGCACTTCGCCAATGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCAttttt                                                                                                        >  1:2150849/1‑150 (MQ=255)
                                          tAACCGCCGTATATTCCCGCAGCCTTGAACAGGCCCAGCACTTCGCCAATGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCAttttt                                                                                                        <  2:1962227/150‑1 (MQ=255)
                                          tAACCGCCGTATATTCCCGCAGCCTTGAACAGGCCCAGCACTTCGCCAATGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCAttttt                                                                                                        >  2:2327307/1‑150 (MQ=255)
                                          tAACCGCCGTATATTCCCGCAGCCTTGAACAGGCCCAGCACTTCGCCAATGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCAttttt                                                                                                        >  2:4688676/1‑150 (MQ=255)
                                             ccgccgTATATTCCCGCAGCCTTGAACAGGCCCAGCACTTCGCCAATGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCATTTTTccc                                                                                                     >  1:30001/1‑150 (MQ=255)
                                             ccgccgTATATTCCCGCAGCCTTGAACAGGCCCAGCACTTCGCCAATGATTTTCCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCATTTTTccc                                                                                                     <  1:3402753/150‑1 (MQ=255)
                                              cgccgTATATTCCCGCAGCCTTGAACAGGCCCAGCACTTCGCCAATGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCATTTTTCCCa                                                                                                    <  1:4706105/150‑1 (MQ=255)
                                              cgccgTATATTCCCGCAGCCTTGAACAGGCCCAGCACTTCGCCAATGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCATTTTTCCCa                                                                                                    >  1:3316097/1‑150 (MQ=255)
                                                 cgTATATTCCCGCAGCCTTGAACAGGCCCAGCACTTCGCCAATGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCATTTTTCCCAGac                                                                                                 >  1:778931/1‑150 (MQ=255)
                                                 cgTATATTCCCGCAGCCTTGAACAGGCCCAGCACTTCGCCAATGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCATTTTTCCCAGac                                                                                                 <  1:458077/150‑1 (MQ=255)
                                                     tatTCCCGCAGCCTTGAACAGGCCCAGCACTTCGCCAATGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCATTTTTCCCAGACACaa                                                                                             <  2:2526121/150‑1 (MQ=255)
                                                       ttCCCGCAGCCTTGAACAGGCCCAGCACTTCGCCAATGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCATTTTTCCCAGACACAACt                                                                                           <  2:278662/150‑1 (MQ=255)
                                                       ttCCCGCAGCCTTGAACAGGCCCAGCACTTCGCCAATGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCATTTTTCCCAGACACAACt                                                                                           >  2:2146438/1‑150 (MQ=255)
                                                           cGCAGCCTTGAACAGGCCCAGCACTTCGCCAATGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCATTTTTCCCAGACACAACTTTTc                                                                                       >  2:1730591/1‑150 (MQ=255)
                                                            gCAGCCTTGAACAGGCCCAGCACTTCGCCAATGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCATTTTTCCCAGACACAACTTTTcc                                                                                      >  1:2287542/1‑150 (MQ=255)
                                                            gCAGCCTTGAACAGGCCCAGCACTTCGCCAATGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCATTTTTCCCAGACACAACTTTTcc                                                                                      >  2:1466265/1‑150 (MQ=255)
                                                            gCAGCCTTGAACAGGCCCAGCACTTCGCCAATGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCATTTTTCCCAGACACAACTTTTcc                                                                                      >  2:1054031/1‑150 (MQ=255)
                                                              aGCCTTGATCAGGCCCAGCACTTCGCCAATGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCATTTTTCCCAGACACAACTTTTCCtt                                                                                    <  2:133661/150‑1 (MQ=255)
                                                                ccTTGAACAGGCCCAGCACTTCGCCAATGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCATTTTTCCCAGACACAACTTTTCCTTAg                                                                                  <  2:2593846/150‑1 (MQ=255)
                                                                ccTTGAACAGGCCCAGCACTTCGCCAATGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCATTTTTCCCAGACACAACTTTTCCTTAg                                                                                  <  1:2394895/150‑1 (MQ=255)
                                                                ccTTGAACAGGCCCAGCACTTCGCCAATGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCATTTTTCCCAGACACAACTTTTCCTTAg                                                                                  >  1:3171639/1‑150 (MQ=255)
                                                                  ttGATCAGGCCCAGCACTTCGCCAATGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCATTTTTCCCAGACACAACTTTTCCTTAGcc                                                                                <  2:408863/150‑1 (MQ=255)
                                                                  ttGAACAGGCCCAGCACTTCGCCAATGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCATTTTTCCCAGACACAACTTTTCCTTAGcc                                                                                <  2:2371854/150‑1 (MQ=255)
                                                                   tGAACAGGCCCAGCACTTCGCCAATGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCATTTTTCCCAGACACAACTTTTCCTTAGCCa                                                                               <  2:4852351/150‑1 (MQ=255)
                                                                      aCAGGCCCAGCACTTCGCCAATGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCATTTTTCCCAGACACAACTTTTCCTTAGCCATaa                                                                            >  2:4909413/1‑150 (MQ=255)
                                                                       cAGGCCCAGCACTTCGCCAATGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCATTTTTCCCAGACACAACTTTTCCTTAGCCATaaa                                                                           >  2:321354/1‑150 (MQ=255)
                                                                       cAGGCCCAGCACTTCGCCAATGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCATTTTTCCCAGACACAACTTTTCCTTAGCCATaaa                                                                           >  1:3652733/1‑150 (MQ=255)
                                                                         ggCCCAGCACTTCGCCAATGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCATTTTTCCCAGACACAACTTTTCCTTAGCCATAAAAt                                                                         >  2:3644457/1‑150 (MQ=255)
                                                                              aGCACTTCGCCAATGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCATTTTTCCCAGACACAACTTTTCCTTAGCCATAAAATTAAtg                                                                    >  2:2472451/1‑150 (MQ=255)
                                                                                cACTTCGCCAATGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCATTTTTCCCAGACACAACTTTTCCTTAGCCATAAAATTAAtgtg                                                                  >  2:4056196/1‑150 (MQ=255)
                                                                                cACTTCGCCAATGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCATTTTTCCCAGACACAACTTTTCCTTAGCCATAAAATTAAtgtg                                                                  >  2:2669809/1‑150 (MQ=255)
                                                                                cACTTCGCCAATGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCATTTTTCCCAGACACAACTTTTCCTTAGCCATAAAATTAAtgtg                                                                  >  1:4195953/1‑150 (MQ=255)
                                                                                 aCTTCGCCAATGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCATTTTTCCCAGACACAACTTTTCCTTAGCCATAAAATTAATGTGa                                                                 >  1:448890/1‑150 (MQ=255)
                                                                                     cGCCAATGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCATTTTTCCCAGACACAACTTTTCCTTAGCCATAAAATTAATGTGATTTg                                                             <  1:1188514/150‑1 (MQ=255)
                                                                                       ccAATGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCATTTTTCCCAGACACAACTTTTCCTTAGCCATAAAATTAATGTGATTTGCg                                                           <  2:4617668/150‑1 (MQ=255)
                                                                                        cAATGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCATTTTTCCCAGACACAACTTTTCCTTAGCCATAAAATTAATGTGATTTGCga                                                          >  2:737941/1‑150 (MQ=255)
                                                                                        cAATGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCATTTTTCCCAGACACAACTTTTCCTTAGCCATAAAATTAATGTGATTTGCga                                                          >  2:1830027/1‑150 (MQ=255)
                                                                                         aaTGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCATTTTTCCCAGACACAACTTTTCCTTAGCCATAAAATTAATGTGATTTGCgag                                                         <  2:4335493/150‑1 (MQ=255)
                                                                                         aaTGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCATTTTTCCCAGACACAACTTTTCCTTAGCCATAAAATTAATGTGATTTGCgag                                                         >  2:1264838/1‑150 (MQ=255)
                                                                                          aTGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCATTTTTCCCAGACACAACTTTTCCTTAGCCATAAAATTAATGTGATTTGCgaga                                                        <  1:946864/150‑1 (MQ=255)
                                                                                           tGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCATTTTTCCCAGACACAACTTTTCCTTAGCCATAAAATTAATGTGATTTGCGAGaa                                                       >  1:3826122/1‑150 (MQ=255)
                                                                                           tGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCATTTTTCCCAGACACAACTTTTCCTTAGCCATAAAATTAATGTGATTTGCGAGaa                                                       >  1:1331339/1‑150 (MQ=255)
                                                                                              tttttCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCATTTTTCCCAGACACAACTTTTCCTTAGCCATAAAATTAATGTGATTTGCGAGAAAcc                                                    >  2:2215896/1‑150 (MQ=255)
                                                                                              tttttCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCATTTTTCCCAGACACAACTTTTCCTTAGCCATAAAATTAATGTGATTTGCGAGAAAcc                                                    <  1:4252341/150‑1 (MQ=255)
                                                                                                tttCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCATTTTTCCCAGACACAACTTTTCCTTAGCCATAAAATTAATGTGATTTGCGAGAAACCGc                                                  <  2:1589279/150‑1 (MQ=255)
                                                                                                   cTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCATTTTTCCCAGACACAACTTTTCCTTAGCCATAAAATTAATGTGATTTGCGAGAAACCGCTgg                                               <  1:1947018/150‑1 (MQ=255)
                                                                                                   cTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCATTTTTCCCAGACACAACTTTTCCTTAGCCATAAAATTAATGTGATTTGCGAGAAACCGCTgg                                               >  1:1503383/1‑150 (MQ=255)
                                                                                                        gAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCATTTTTCCCAGACACAACTTTTCCTTAGCCATAAAATTAATGTGATTTGCGAGAAACCGCTGGCGTCg                                          >  2:4229827/1‑150 (MQ=255)
                                                                                                         aGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCATTTTTCCCAGACACAACTTTTCCTTAGCCATAAAATTAATGTGATTTGCGAGAAACCGCTGGCGTCGa                                         >  1:4017072/1‑150 (MQ=255)
                                                                                                            aTCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCATTTTTCCCAGACACAACTTTTCCTTAGCCATAAAATTAATGTGATTTGCGAGAAACCGCTGGCGTCGAATc                                      <  2:2872505/150‑1 (MQ=255)
                                                                                                              cTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCATTTTTCCCAGACACAACTTTTCCTTAGCCATAAAATTAATGTGATTTGCGAGAAACCGCTGGCGTCGAATCTg                                    <  1:5092517/150‑1 (MQ=255)
                                                                                                               tGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCAAATTCCCTGCATTTTTCCCAGACACAACTTTTCCTTGGCCATAAAATTAATGTGATTTGCGAGAAACCGCTGGCGTCGAATCTgg                                   >  1:1361758/1‑150 (MQ=255)
                                                                                                                gTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCATTTTTCCCAGACACAACTTTTCCTTAGCCATAAAATTAATGTGATTTGCGAGAAACCGCTGGCGTCGAATCTGGc                                  <  1:1463923/150‑1 (MQ=255)
                                                                                                                  ttACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCATTTTTCCCAGACACAACTTTTCCTTAGCCATAAAATTAATGTGATTTGCGAGAAACCGCTGGCGTCGAATCTGGCgg                                >  2:1115720/1‑150 (MQ=255)
                                                                                                                   tACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCATTTTTCCCAGACACAACTTTTCCTTAGCCATAAAATTAATGTGATTTGCGAGAAACCGCTGGCGTCGAATCTGGCGGa                               >  2:3722228/1‑150 (MQ=255)
                                                                                                                   tACCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCATTTTTCCCAGACACAACTTTTCCTTAGCCATAAAATTAATGTGATTTGCGAGAAACCGCTGGCGTCGAATCTGGCGGa                               <  1:4829959/150‑1 (MQ=255)
                                                                                                                    aCCTCGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCATTTTTCCCAGACACAACTTTTCCTTAGCCATAAAATTAATGTGATTTGCGAGAAACCGCTGGCGTCGAATCTGGCGGaa                              >  1:4332376/1‑150 (MQ=255)
                                                                                                                        cGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCATTTTTCCCAGACACAACTTTTCCTTAGCCATAAAATTAATGTGATTTGCGAGAAACCGCTGGCGTCGAATCTGGCGGAAGTgg                          <  2:3924954/150‑1 (MQ=255)
                                                                                                                        cGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCATTTTTCCCAGACACAACTTTTCCTTAGCCATAAAATTAATGTGATTTGCGAGAAACCGCTGGCGTCGAATCTGGCGGAAGTgg                          >  2:536137/1‑150 (MQ=255)
                                                                                                                        cGCTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCATTTTTCCCAGACACAACTTTTCCTTAGCCATAAAATTAATGTGATTTGCGAGAAACCGCTGGCGTCGAATCTGGCGGAAGTgg                          >  2:814789/1‑150 (MQ=255)
                                                                                                                          cTGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCATTTTTCCCAGACACAACTTTTCCTTAGCCATAAAATTAATGTGATTTGCGAGAAACCGCTGGCGTCGAATCTGGCGGAAGTGGAt                        >  1:3054112/1‑150 (MQ=255)
                                                                                                                           tGGAAGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCATTTTTCCCAGACACAACTTTTCCTTAGCCATAAAATTAATGTGATTTGCGAGAAACCGCTGGCGTCGAATCTGGCGGAAGTGGATg                       >  2:3998439/1‑150 (MQ=255)
                                                                                                                              aaGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCATTTTTCCCAGACACAACTTTTCCTTAGCCATAAAATTAATGTGATTTGCGAGAAACCGCTGGCGTCGAATCTGGCGGAAGTGGATgccg                    >  2:4825889/1‑150 (MQ=255)
                                                                                                                              aaGCGATGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCATTTTTCCCAGACACAACTTTTCCTTAGCCATAAAATTAATGTGATTTGCGAGAAACCGCTGGCGTCGAATCTGGCGGAAGTGGATgccg                    >  2:2020264/1‑150 (MQ=255)
                                                                                                                                    tGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCATTTTTCCCAGACACAACTTTTCCTTAGCCATAAAATTAATGTGATTTGCGAGAAACCGCTGGCGTCGAATCTGGCGGAAGTGGATGCCGCCATTg              >  2:3543656/1‑150 (MQ=255)
                                                                                                                                    tGGCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCATTTTTCCCAGACACAACTTTTCCTTAGCCATAAAATTAATGTGATTTGCGAGAAACCGCTGGCGTCGAATCTGGCGGAAGTGGATGCCGCCATTg              >  1:611770/1‑150 (MQ=255)
                                                                                                                                      gCGGAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCATTTTTCCCAGACACAACTTTTCCTTAGCCATAAAATTAATGTGATTTGCGAGAAACCGCTGGCGTCGAATCTGGCGGAAGTGGATGCCGCCATTGcc            <  1:2528427/150‑1 (MQ=255)
                                                                                                                                         gAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCATTTTTCCCAGACACAACTTTTCCTTAGCCATAAAATTAATGTGATTTGCGAGAAACCGCTGGCGTCGAATCTGGCGGAAGTGGATGCCGCCATTGCCtgt         <  1:5017869/150‑1 (MQ=255)
                                                                                                                                         gAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCATTTTTCCCAGACACAACTTTTCCTTAGCCATAAAATTAATGTGATTTGCGAGAAACCGCTGGCGTCGAATCTGGCGGAAGTGGATGCCGCCATTGCCtgt         <  2:2990544/150‑1 (MQ=255)
                                                                                                                                         gAAAGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCATTTTTCCCAGACACAACTTTTCCTTAGCCATAAAATTAATGTGATTTGCGAGAAACCGCTGGCGTCGAATCTGGCGGAAGTGGATGCCGCCATTGCCtgt         >  1:2109076/1‑150 (MQ=255)
                                                                                                                                          aaaGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCATTTTTCCCAGACACAACTTTTCCTTAGCCATAAAATTAATGTGATTTGCGAGAAACCGCTGGCGTCGAATCTGGCGGAAGTGGATGCCGCCATTGCCtgtg        >  2:5005013/1‑150 (MQ=255)
                                                                                                                                           aaGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCATTTTTCCCAGACACAACTGTTCCTTAGCCATCAAATTAATGTGATTTGCGAGAAACCGCTGGCGTCGAATCTGGCGGAAGTGGATGCCGCCATAGCCGGTgc       >  1:2946287/1‑150 (MQ=255)
                                                                                                                                            aGCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCATTTTTCCCAGACACAACTTTTCCTTAGCCATAAAATTAATGTGATTTGCGAGAAACCGCTGGCGTCGAATCTGGCGGAAGTGGATGCCGCCATTGCCTGTgcg      <  2:1158352/150‑1 (MQ=255)
                                                                                                                                             gCGATGGCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCATTTTTCCCAGACACAACTTTTCCTTAGCCATAAAATTAATGTGATTTGCGAGAAACCGCTGGCGTCGAATCTGGCGGAAGTGGATGCCGCCATTGCCTGTgcgc     <  2:2276531/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
GCCAGTTCGTCGAGGCCGCCCATGAGAGCGGTAAATACAAGTTAACCGCCGTATATTCCCGCAGCCTTGAACAGGCCCAGCACTTCGCCAATGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGCCATTGACGCGGTGTATATTGCCAGCCCGAATTCCCTGCATTTTTCCCAGACACAACTTTTCCTTAGCCATAAAATTAATGTGATTTGCGAGAAACCGCTGGCGTCGAATCTGGCGGAAGTGGATGCCGCCATTGCCTGTGCGCGGG  >  CP009273/3230710‑3231003

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: