Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 545,522 C→A 27.4% A147E (GCG→GAG)  ylbF → putative anaerobic allantoin catabolic oxamate carbamoyltransferase; DUF2877 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273545,5220CA27.4% 807.1 / 297.3 372A147E (GCG→GAG) ylbFputative anaerobic allantoin catabolic oxamate carbamoyltransferase; DUF2877 family protein
Reads supporting (aligned to +/- strand):  ref base C (127/143);  new base A (52/50);  total (179/193)
Fisher's exact test for biased strand distribution p-value = 5.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

Only 100 of 378 total aligned reads displayed.

GCGTATTACGCCGCCTGCGCATCCTCAACCACTTGCAGCTGCATGGATGCAACGCGCGGAGGAAACCGGGCTTTTCGGGCCACTGGCGTTGGCGGCAAGCGATCCGCTGCCTGCTGAGTTACGCCAGTTTCGTCACTGTTTTCAGGCCGCGCTCAATGGCGTTAAGACCGACTGGCGGCACTGGCTGGGTAAAGGCCCCGGATTAACGCCGAGTCATGATGACACGCTGAGCGGAATGCTGCTGGCGGCCTGGTATTATGGCGCTTTAGATGCGCGCTCCGGTCGTCCGTTTTTTGCCT  >  CP009273/545373‑545671
                                                                                                                                                     |                                                                                                                                                     
   tATTACGCCGCCTGCGCATCCTCAACCACTTGCAGCTGCATGGATGCAACGCGCGGAGGAAACCGGGCTTTTCGGGCCACTGGCGTTGGCGGCAAGCGATCCGCTGCCTGCTGAGTTACGCCAGTTTCGTCACTGTTTTCAGGCCGCGCt                                                                                                                                                    <  1:1114454/150‑1 (MQ=255)
   tATTACGCCGCCTGCGCATCCTCAACCACTTGCAGCTGCATGGATGCAACGCGCGGAGGAAACCGGGCTTTTCGGGCCACTGGCGTTGGCGGCAAGCGATCCGCTGCCTGCTGAGTTACGCCAGTTTCGTCACTGTTTTCAGGCCGCGCt                                                                                                                                                    <  1:3366455/150‑1 (MQ=255)
     ttACGCCGCCTGCGCATCCTCAACCACTTGCAGCTGCATGGATGCAACGCGCGGAGGAAACCGGGCTTTTCGGGCCACTGGCGTTGGCGGCAAGCGATCCGCTGCCTGCTGAGTTACGCCAGTTTCGTCACTGTTTTCAGGCCGAGCTCa                                                                                                                                                  <  2:399633/150‑1 (MQ=255)
       aCGCCGCCTGCGCATCCTCAACCACTTGCAGCTGCATGGATGCAACGCGCGGAGGAAACCGGGCTTTTCGGGCCACTGGCGTTGGCGGCAAGCGATCCGCTGCCTGCTGAGTTACGCCAGTTTCGTCACTGTTTTCAGGCCGCGCTCAAt                                                                                                                                                >  2:4496707/1‑150 (MQ=255)
       aCGCCGCCTGCGCATCCTCAACCACTTGCAGCTGCATGGATGCAACGCGCGGAGGAAACCGGGCTTTTCGGGCCACTGGCGTTGGCGGCAAGCGATCCGCTGCCTGCTGAGTTACGCCAGTTTCGTCACTGTTTTCAGGCCGCGCTCAAt                                                                                                                                                >  1:1682638/1‑150 (MQ=255)
       aCGCCGCCTGCGCATCCTCAACCACTTGCAGCTGCATGGATGCAACGCGCGGAGGAAACCGGGCTTTTCGGGCCACTGGCGTTGGCGGCAAGCGATCCGCTGCCTGCTGAGTTACGCCAGTTTCGTCACTGTTTTCAGGCCGAGCTCAAt                                                                                                                                                <  1:4791059/150‑1 (MQ=255)
        cgccgcCTGCGCATCCTCAACCACTTGCAGCTGCATGGATGCAACGCGCGGAGGAAACCGGGCTTTTCGGGCCACTGGCGTTGGCGGCAAGCGATCCGCTGCCTGCTGAGTTACGCCAGTTTCGTCACTGTTTTCAGGCCGCGCTCAATg                                                                                                                                               <  1:2874263/150‑1 (MQ=255)
          ccgccTGCGCATCCTCAACCACTTGCAGCTGCATGGATGCAACGCGCGGAGGAAACCGGGCTTTTCGGGCCACTGGCGTTGGCGGCAAGCGATCCGCTGCCTGCTGAGTTACGCCAGTTTCGTCACTGTTTTCAGGCCGAGCTCAATGGc                                                                                                                                             >  1:925902/1‑150 (MQ=255)
           cgccTGCGCATCCTCAACCACTTGCAGCTGCATGGATGCAACGCGCGGAGGAAACCGGGCTTTTCGGGCCACTGGCGTTGGCGGCAAGCGATCCGCTGCCTGCTGAGTTACGCCAGTTTCGTCACTGTTTTCAGGCCGCGCTCAATGGCg                                                                                                                                            >  2:4023704/1‑150 (MQ=255)
            gccTGCGCATCCTCAACCACTTGCAGCTGCATGGATGCAACGCGCGGAGGAAACCGGGCTTTTCGGGCCACTGGCGTTGGCGGCAAGCGATCCGCTGCCTGCTGAGTTACGCCAGTTTCGTCACTGTTTTCAGGCCGCGCTCAATGGCGt                                                                                                                                           >  2:868658/1‑150 (MQ=255)
            gccTGCGCATCCTCAACCACTTGCAGCTGCATGGATGCAACGCGCGGAGGAAACCGGGCTTTTCGGGCCACTGGCGTTGGCGGCAAGCGATCCGCTGCCTGCTGAGTTACGCCAGTTTCGTCACTGTTTTCAGGCCGCGCTCAATGGCGt                                                                                                                                           <  1:2590373/150‑1 (MQ=255)
             ccTGCGCATCCTCAACCACTTGCAGCTGCATGGATGCAACGCGCGGAGGAAACCGGGCTTTTCGGGCCACTGGCGTTGGCGGCAAGCGATCCGCTGCCTGCTGAGTTACGCCAGTTTCGTCACTGTTTTCAGGCCGCGCTCAATGGCGtt                                                                                                                                          <  1:4163658/150‑1 (MQ=255)
              cTGCGCATCCTCAACCACTTGCAGCTGCATGGATGCAACGCGCGGAGGAAACCGGGCTTTTCGGGCCACTGGCGTTGGCGGCAAGCGATCCGCTGCCTGCTGAGTTACGCCAGTTTCGTCACTGTTTTCAGGCCGCGCTCAATGGCGTTa                                                                                                                                         >  2:4361981/1‑150 (MQ=255)
              cTGCGCATCCTCAACCACTTGCAGCTGCATGGATGCAACGCGCGGAGGAAACCGGGCTTTTCGGGCCACTGGCGTTGGCGGCAAGCGATCCGCTGCCTGCTGAGTTACGCCAGTTTCGTCACTGTTTTCAGGCCGCGCTCAATGGCGTTa                                                                                                                                         <  1:241073/150‑1 (MQ=255)
               tGCGCATCCTCAACCACTTGCAGCTGCATGGATGCAACGCGCGGAGGAAACCGGGCTTTTCGGGCCACTGGCGTTGGCGGCAAGCGATCCGCTGCCTGCTGAGTTACGCCAGTTTCGTCACTGTTTTCAGGCCGCGCTCAATGGCGTTaa                                                                                                                                        >  1:3647790/1‑150 (MQ=255)
                    aTCCTCAACCACTTGCAGCTGCATGGATGCAACGCGCGGAGGAAACCGGGCTTTTCGGGCCACTGGCGTTGGCGGCAAGCGATCCGCTGCCTGCTGAGTTACGCCAGTTTCGTCACTGTTTTCAGGCCGAGCTCAATGGCGTTAAGACCg                                                                                                                                   >  2:4327815/1‑150 (MQ=255)
                      ccTCAACCACTTGCAGCTGCATGGATGCAATGCGCGGAGGAAACCGGGCTTTTCGGGCCACTGGCGTTGGCGGCAAGCGATCCGCTGCCTGCTGAGTTACGCCAGTTTCGTCACTGTTTTCAGGCCGCGCTCAATGGCGTTAAGACCGAc                                                                                                                                 >  1:740331/1‑150 (MQ=255)
                       cTCAACCACTTGCAGCTGCATGGATGCAACGCGCGGAGGAAACCGGGCTTTTCGGGCCACTGGCGTTGGCGGCAAGCGATCCGCTGCCTGCTGAGTTACGCCAGTTTCGTCACTGTTTTCAGGCCGCGCTCAATGGCGTTAAGACCGACt                                                                                                                                <  2:4724870/150‑1 (MQ=255)
                       cTCAACCACTTGCAGCTGCATGGATGCAACGCGCGGAGGAAACCGGGCTTTTCGGGCCACTGGCGTTGGCGGCAAGCGATCCGCTGCCTGCTGAGTTACGCCAGTTTCGTCACTGTTTTCAGGCCGCGCTCAATGGCGTTAAGACCGACt                                                                                                                                >  2:1652464/1‑150 (MQ=255)
                        tCAACCACTTGCAGCTGCATGGATGCAACGCGCGGAGGAAACCGGGCTTTTCGGGCCACTGGCGTTGGCGGCAAGCGATCCGCTGCCTGCTGAGTTACGCCAGTTTCGTCACTGTTTTCAGGCCGCGCTCAATGGCGTTAAGACCGACTg                                                                                                                               <  1:211370/150‑1 (MQ=255)
                               cTTGCAGCTGCATGGATGCAACGCGCGGAGGAAACCGGGCTTTTCGGGCCACTGGCGTTGGCGGCAAGCGATCCGCTGCCTGCTGAGTTACGCCAGTTTCGTCACTGTTTTCAGGCCGCGCTCAATGGCGTTAAGACCGACTGGCGGCAc                                                                                                                        >  2:3264042/1‑150 (MQ=255)
                               cTTGCAGCTGCATGGATGCAACGCGCGGAGGAAACCGGGCTTTTCGGGCCACTGGCGTTGGCGGCAAGCGATCCGCTGCCTGCTGAGTTACGCCAGTTTCGTCACTGTTTTCAGGCCGCGCTCAATGGCGTTAAGACCGACTGGCGGCAc                                                                                                                        >  1:4750467/1‑150 (MQ=255)
                                ttGCTGCTGGATGGTTGCAACGCGCGGTTGAATGCGGGCTTTTCGGGCCACTGGCGTTGGCGGCAAGAGATCCGCTGCCTGCTGAGTTACGCCAGTTTCGTCACTTTTTTCAGGCCGCGCTAAATGGCGTTAATACCGACTGGCGGCAct                                                                                                                       <  2:381925/150‑1 (MQ=255)
                                  gCAGCTGCATGGATGCAACGCGCGGAGGAAACCGGGCTTTTCGGGCCACTGGCGTTGGCGGCAAGCGATCCGCTGCCTGCTGAGTTACGCCAGTTTCGTCACTGTTTTCAGGCCGCGCTCAATGGCGTTAAGACCGACTGGCGGCActgg                                                                                                                     <  1:868658/150‑1 (MQ=255)
                                   cAGCTGCATGGATGCAACGCGCGGAGGAAACCGGGCTTTTCGGGCCACTGGCGTTGGCGGCAAGCGATCCGCTGCCTGCTGAGTTACGCCAGTTTCGTCACTGTTTTCAGGCCGCGCTCAATGGCGTTAAGACCGACTGGCGGCActggc                                                                                                                    >  2:1219919/1‑150 (MQ=255)
                                    aGCTGCATGGATGCAACGCGCGGAGGAAACCGGGCTTTTCGGGCCACTGGCGTTGGCGGCAAGCGATCCGCTGCCTGCTGAGTTACGCCAGTTTCGTCACTGTTTTCAGGCCGAGCTCAATGGCGTTAAGACCGACTGGCGGCActggct                                                                                                                   <  1:2707071/150‑1 (MQ=255)
                                        gCATGGATGCAACGCGCGGAGGAAACCGGGCTTTTCGGGCCACTGGCGTTGGCGGCAAGCGATCCGCTGCCTGCTGAGTTACGCCAGTTTCGTCACTGTTTTCAGGCCGCGCTCAATGGCGTTAAGACCGACTGGCGGCACTGGCTGGGt                                                                                                               >  1:3025682/1‑150 (MQ=255)
                                        gCATGGATGCAACGCGCGGAGGAAACCGGGCTTTTCGGGCCACTGGCGTTGGCGGCAAGCGATCCGCTGCCTGCTGAGTTACGCCAGTTTCGTCACTGTTTTCAGGCCGCGCTCAATGGCGTTAAGACCGACTGGCGGCACTGGCTGGGt                                                                                                               >  2:1887459/1‑150 (MQ=255)
                                          aTGGATGCAACGCGCGGAGGAAACCGGGCTTTTCGGGCCACTGGCGTTGGCGGCAAGCGATCCGCTGCCTGCTGAGTTACGCCAGTTTCGTCACTGTTTTCAGGCCGCGCTCAATGGCGTTAAGACCGACTGGCGGCACTGGCTGGGTaa                                                                                                             <  2:3567550/150‑1 (MQ=255)
                                            ggATGCAACGCGCGGAGGAAACCGGGCTTTTCGGGCCACTGGCGTTGGCGGCAAGCGATCCGCTGCCTGCTGAGTTACGCCAGTTTCGTCACTGTTTTCAGGCCGCGCTCAATGGCGTTAAGACCGACTGGCGGCACTGGCTGGGTAAAg                                                                                                           >  1:3161190/1‑150 (MQ=255)
                                                 cAACGCGCGGAGGAAACCGGGCTTTTCGGGCCACTGGCGTTGGCGGCAAGCGATCCGCTGCCTGCTGAGTTACGCCAGTTTCGTCACTGTTTTCAGGCCGCGCTCAATGGCGTTAAGACCGACTGGCGGCACTGGCTGGGTAAAGGcccc                                                                                                      >  2:2280582/1‑150 (MQ=255)
                                                  aaCGCGCGGAGGAAACCGGGCTTTTCGGGCCACTGGCGTTGGCGGCAAGCGATCCGCTGCCTGCTGAGTTACGCCAGTTTCGTCACTGTTTTCAGGCCGCGCTCAATGGCGTTAAGACCGACTGGCGGCACTGGCTGGGTAAAGGCCCCg                                                                                                     <  2:555283/150‑1 (MQ=255)
                                                  aaCGCGCGGAGGAAACCGGGCTTTTCGGGCCACTGGCGTTGGCGGCAAGCGATCCGCTGCCTGCTGAGTTACGCCAGTTTCGTCACTGTTTTCAGGCCGCGCTCAATGGCGTTAAGACCGACTGGCGGCACTGGCTGGGTAAAGGCCCCg                                                                                                     <  2:3757084/150‑1 (MQ=255)
                                                    cgcgcgGAGGAAACCGGGCTTTTCGGGCCACTGGCGTTGGCGGCAAGCGATCCGCTGCCTGCTGAGTTACGCCAGTTTCGTCACTGTTTTCAGGCCGCGCTCAATGGCGTTAAGACCGACTGGCGGCACTGGCTGGGTAAAGGCCCCGGa                                                                                                   <  1:963438/150‑1 (MQ=255)
                                                      cgcgGAGGAAACCGGGCTTTTCGGGCCACTGGCGTTGGCGGCAAGCGATCCGCTGCCTGCTGAGTTACGCCAGTTTCGTCACTGTTTTCAGGCCGCGCTCAATGGCGTTTAGACCGACTGGCGGCACTGGCTGGGTAAAGGCCCCGGAtt                                                                                                 <  1:1622096/150‑1 (MQ=255)
                                                         ggaggaAACCGGGCTTTTCGGGCCACTGGCGTTGGCGGCAAGCGATCCGCTGCCTGCTGAGTTACGCCAGTTTCGTCACTGTTTTCAGGCCGCGCTCAATGGCGTTAAGACCGACTGGCGGCACTGGCTGGGTAAAGGCCCCGGATTAAc                                                                                              <  1:3938771/150‑1 (MQ=255)
                                                         ggaggaAACCGGGCTTTTCGGGCCACTGGCGTTGGCGGCAAGCGATCCGCTGCCTGCTGAGTTACGCCAGTTTCGTCACTGTTTTCAGGCCGCGCTCAATGGCGTTAAGACCGACTGGCGGCACTGGCCGGGTAAAGGCCCCGGATTAAc                                                                                              >  2:297926/1‑150 (MQ=255)
                                                           aggaAACCGGGCTTTTCGGGCCACTGGCGTTGGCGGCAAGCGATCCGCTGCCTGCTGAGTTACGCCAGTTTCGTCACTGTTTTCAGGCCGCGCTCAATGGCGTTAAGACCGACTGGCGGCACTGGCTGGGTAAAGGCCCCGGATTAACGc                                                                                            >  1:1038097/1‑150 (MQ=255)
                                                            ggaAACCGGGCTTTTCGGGCCACTGGCGTTGGCGGCAAGCGATCCGCTGCCTGCTGAGTTACGCCAGTTTCGTCACTGTTTTCAGGCCGCGCTCAATGGCGTTAAGACCGACTGGCGGCACTGGCTGGGTAAAGGCCCCGGATTAACGcc                                                                                           >  2:618933/1‑150 (MQ=255)
                                                             gaAACCGGGCTTTTCGGGCCACTGGCGTTGGCGGCAAGCGATCCGCTGCCTGCTGAGTTACGCCAGTTTCGTCACTGTTTTCAGGCCGCGCTCAATGGCGTTAAGACCGACTGGCGGCACTGGCTGGGTAAAGGCCCCGGATTAACGCCg                                                                                          <  1:596637/150‑1 (MQ=255)
                                                             gaAACCGGGCTTTTCGGGCCACTGGCGTTGGCGGCAAGCGATCCGCTGCCTGCTGAGTTACGCCAGTTTCGTCACTGTTTTCAGGCCGAGCTCAATGGCGTTAAGACCGACTGGCGGCACTGGCTGGGTAAAGGCCCCGGATTAACGCCg                                                                                          <  2:1063253/150‑1 (MQ=255)
                                                             gaAACCGGGCTTTTCGGGCCACTGGCGTTGGCGGCAAGCGATCCGCTGCCTGCTGAGTTACGCCAGTTTCGTCACTGTTTTCAGGCCGAGCTCAATGGCGTTAAGACCGACTGGCGGCACTGGCTGGGTAAAGGCCCCGGATTAACGCCg                                                                                          >  1:997105/1‑150 (MQ=255)
                                                             gaAACCGGGCTTTTCGGGCCACTGGCGTTGGCGGCAAGCGATCCGCTGCCTGCTGAGTTACGCCAGTTTCGTCACTGTTTTCAGGCCGAGCTCAATGGCGTTAAGACCGACTGGCGGCACTGGCTGGGTAAAGGCCCCGGATTAACGCCg                                                                                          <  2:3396682/150‑1 (MQ=255)
                                                              aaaCCGGGCTTTTCGGGCCACTGGCGTTGGCGGCAAGCGATCCGCTGCCTGCTGAGTTACGCCAGTTTCGTCACTGTTTTCAGGCCGCGCTCAATGGCGTTAAGACCGACTGGCGGCACTGGCTGGGTAAAGGCCCCGGATTAACGCCGa                                                                                         >  2:2719399/1‑150 (MQ=255)
                                                              aaaCCGGGCTTTTCGGGCCACTGGCGTTGGCGGCAAGCGATCCGCTGCCTGCTGAGTTACGCCAGTTTCGTCACTGTTTTCAGGCCGAGCTCAATGGCGTTAAGACCGACTGGCGGCACTGGCTGGGTAAAGGCCCCGGATTAACGCCGa                                                                                         >  2:3918512/1‑150 (MQ=255)
                                                              aaaCCGGGCTTTTCGGGCCACTGGCGTTGGCGGCAAGCGATCCGCTGCCTGCTGAGTTACGCCAGTTTCGTCACTGTTTTCAGGCCGAGCTCAATGGCGTTAAGACCGACTGGCGGCACTGGCTGGGTAAAGGCCCCGGATTAACGCCGa                                                                                         >  2:175161/1‑150 (MQ=255)
                                                              aaaCCGGGCTTTTCGGGCCACTGGCGTTGGCGGCAAGCGATCCGCTGCCTGCTGAGTTACGCCAGTTTCGTCACTGTTTTCAGGCCGAGCTCAATGGCGTTAAGACCGACTGGCGGCACTGGCTGGGTAAAGGCCCCGGATTAACGCCGa                                                                                         >  2:4563191/1‑150 (MQ=255)
                                                               aaCCGGGCTTTTCGGGCCACTGGCGTTGGCGGCAAGCGATCCGCTGCCTGCTGAGTTACGCCAGTTTCGTCACTGTTTTCAGGCCGCGCTCAATGGCGTTAAGACCGACTGGCGGCACTGGCTGGGTAAAGGCCCCGGATTAACGCCGAg                                                                                        <  2:3208481/150‑1 (MQ=255)
                                                                    ggCTTTTCGGGCCACTGGCGTTGGCGGCAAGCGATCCGCTGCCTGCTGAGTTACGCCAGTTTCGTCACTGTTTTCAGGCCGCGCTCAATGGCGTTAAGACCGACTGGCGGCACTGGCTGGGTAAAGGCCCCGGATTAACGCCGAGTCatg                                                                                   <  1:2051570/150‑1 (MQ=255)
                                                                    ggCTTTTCGGGCCACTGGCGTTGGCGGCAAGCGATCCGCTGCCTGCTGAGTTACGCCAGTTTCGTCACTGTTTTCAGGCCGCGCTCAATGGCGTTAAGACCGACTGGCGGCACTGGCTGGGTAAAGGCCCCGGATTAACGCCGAGTCatg                                                                                   >  1:3419540/1‑150 (MQ=255)
                                                                     ccTTTTCGGGACACTGGCGTTGGCGGCAAGCGATCAGCAGCCTGCTGAGTTACGCCAGTTTCGTCACTGTTTTCAGGCCGCGCTCAATGGCGATAAGACCGACTGGCGGCACTGGCTGGGTAAATGACCCGGATAAACTCCGAGTCatga                                                                                  >  2:2592125/2‑150 (MQ=255)
                                                                      cTTTTCGGGCCACTGGCGTTGGCGGCAAGCGATCCGCTGCCTGCTGAGTTACGCCAGTTTCGTCACTGTTTTCAGGCCGCGCTCAATGGCGTTAAGACCGACTGGCGGCACTGGCTGGGTAAAGGCCCCGGATTAACGCCGAGTCatgat                                                                                 >  2:1337932/1‑150 (MQ=255)
                                                                       ttttCGGGCCACTGGCGTTGGCGGCAAGCGATCCGCTGCCTGCTGAGTTACGCCAGTTTCGTCACGGTTTTCAGGCCGCGCTCAATGGCGTTAAGACCGACTGGCGGCACTGGCTGGGTAAAGGCCCCGGATTAACGCCGAGTCatgatg                                                                                >  2:3673946/1‑150 (MQ=255)
                                                                         ttCGGGCCACTGGCGTTGGCGGCAAGCGATCCGCTGCCTGCTGAGTTACGCCAGTTTCGTCACTGTTTTCAGGCCGAGCTCAATGGCGTTAAGACCGACTGGCGGCACTGGCTGGGTAAAGGCCCCGGATTAACGCCGAGTCATGATGac                                                                              <  1:1711714/150‑1 (MQ=255)
                                                                            gggCCACTGGCGTTGGCGGCAAGCGATCCGCTGCCTGCTGAGTTACGCCAGTTTCGTCACTGTTTTCAGGCCGCGCTCAATGGCGTTAAGACCGACTGGCGGCACTGGCTGGGTAAAGGCCCCGGATTAACGCCGAGTCATGATGACACg                                                                           >  1:3786409/1‑150 (MQ=255)
                                                                               ccACTGGCGTTGGCGGCAAGCGATCCGCTGCCTGCTGAGTTACGCCAGTTTCGTCACTGTTTTCAGGCCGCGCTCAATGGCGTTAAGACCGACTGGCGGCACTGGCTGGGTAAAGGCCCCGGATTAACGCCGAGTCATGATGACACGCTg                                                                        >  1:3181952/1‑150 (MQ=255)
                                                                                cACTGGCGTTGGCGGCAAGCGATCCGCTGCCTGCTGAGTTACGCCAGTTTCGTCACTGTTTTAAGGCCGCGCTCAATGGCGTTAAGACCGACTGGCGGCACTGGCTGGGTAAAGGCCCCGGATTAACGCCGAGTCATGATGACACGCTGa                                                                       >  1:4491634/1‑150 (MQ=255)
                                                                                  cTGGCGTTGGCGGCAAGCGATCCGCTGCCTGCTGAGTTACGCCAGTTTCGTCACTGTTTTCAGGCCGAGCTCAATGGCGTTAAGACCGACTGGCGGCACTGGCTGGGTAAAGGCCCCGGATTAACGCCGAGTCATGATGACACGCTGAGc                                                                     <  1:1295965/150‑1 (MQ=255)
                                                                                       gTTGGCGGCAAGCGATCCGCTGCCTGCTGAGTTACGCCAGTTTCGTCACTGTTTTCAGGCCGCGCTCAATGGCGTTAAGACCGACTGGCGGCACTGGCTGGGTAAAGGCCCCGGATTAACGCCGAGTCATGATGACACGCTGAGCGGAAt                                                                <  2:1682638/150‑1 (MQ=255)
                                                                                       gTTGGCGGCAAGCGATCCGCTGCCTGCTGAGTTACGCCAGTTTCGTCACTGTTTTCAGGCCGAGCTCAATGGCGTTAAGACCGACTGGCGGCACTGGCTGGGTAAAGGCCCCGGATTAACGCCGAGTCATGATGACACGCTGAGCGGAAt                                                                <  1:2559365/150‑1 (MQ=255)
                                                                                         tGGCGGCAAGCGATCCGCTGCCTGCTGAGTTACGCCAGTTTCGTCACTGTTTTCAGGCCGCGCTCAATGGCGTTAAGACCGACTGGCGGCACTGGCTGGGTAAAGGCCCCGGATTAACGCCGAGTCATGATGACACGCTGAGCGGAAtgc                                                              >  1:2318771/1‑150 (MQ=255)
                                                                                           gcggcAAGCGATCCGCTGCCTGCTGAGTTACGCCAGTTTCGTCACTGTTTTCAGGCCGCGCTCAATGGCGTTAAGACCGACTGGCGGCACTGGCTGGGTAAAGGCCCCGGATTAACGCCGAGTCATGATGACACGCTGAGCGGAAtgctg                                                            <  2:2886164/150‑1 (MQ=255)
                                                                                             ggcAAGCGATCCGCTGCCTGCTGAGTTACGCCAGTTTCGTCACTGTTTTCAGGCCGCGCTCAATGGCGTTAAGACCGACTGGCGGCACTGGCTGGGTAAAGGCCCCGGATTAACGCCGAGTCATGATGACACGCTGAGCGGAATgctgct                                                          >  1:4668744/1‑150 (MQ=255)
                                                                                             ggcAAGCGATCCGCTGCCTGCTGAGTTACGCCAGTTTCGTCACTGTTTTCAGGCCGAGCTCAATGGCGTTAAGACCGACTGGCGGCACTGGCTGGGTAAAGGCCCCGGATTAACGCCGAGTCATGATGACACGCTGAGCGGAATgctgct                                                          <  1:617492/150‑1 (MQ=255)
                                                                                             ggcAAGCGATCCGCTGCCTGCTGAGTTACGCCAGTTTCGTCACTGTTTTCAGGCCGAGCTCAATGGCGTTAAGACCGACTGGCGGCACTGGCTGGGTAAAGGCCCCGGATTAACGCCGAGTCATGATGACACGCTGAGCGGAATgctgct                                                          <  1:4729998/150‑1 (MQ=255)
                                                                                                   cGATCCGCTGCCTGCTGAGTTACGCCAGTTTCGTCACTGTTTTCAGGCCGCGCTCAATGGCGTTAAGACCGACTGGCGGCACTGGCTGGGTAAAGGCCCCGGATTAACGCCGAGTCATGATGACACGCTGAGCGGACTGCTGCTggcggc                                                    >  1:3308326/1‑150 (MQ=255)
                                                                                                    gATCCGCTGCCTGCTGAGTTACGCCAGTTTCGTCACTGTTTTCAGGCCGAGCTCAATGGCGTTAAGACCGACTGGCGGCACTGGCTGGGTAAAGGCCCCGGATTAACGCCGAGTCATGATGACACGCTGAGCGGAATGCTGCTGGCGGcc                                                   >  1:4541448/1‑150 (MQ=255)
                                                                                                    gATCCGCTGCCTGCTGAGTTACGCCAGTTTCGTCACTGTTTTCAGGCCGAGCTCAATGGCGTTAAGACCGACTGGCGGCACTGGCTGGGTAAAGGCCCCGGATTAACGCCGAGTCATGATGACACGCTGAGCGGAATGCTGCTGGCGGcc                                                   >  1:2096999/1‑150 (MQ=255)
                                                                                                     aTCCGCTGCCTGCTGAGTTACGCCAGTTTCGTCACTGTTTTCAGGCCGCGCTCAATGGCGTTAAGACCGACTGGCGGCACTGGCTGGGTAAAGGCCCCGGATTAACGCCGAGTCATGATGACACGCTGAGCGGAATGCTGCTGGCGGCCt                                                  <  1:224959/150‑1 (MQ=255)
                                                                                                      tCCGCTGCCTGCTGAGTTACGCCAGTTTCGTCACTGTTTTCAGGCCGCGCTCAATGGCGTTAAGACCGACTGGCGGCACTGGCTGGGTAAAGGCCCCGGATTAACGCCGAGTCATGATGACACGCTGAGCGGAATGCTGCTGGCGGCCTg                                                 <  2:804446/150‑1 (MQ=255)
                                                                                                        cGCTGCCTGCTGAGTTACGCCAGTTTCGTCACTGTTTTCAGGCCGCGCTCAATGGCGTTAAGACCGACTGGCGGCACTGGCTGGGTAAAGGCCCCGGATTAACGCCGAGTCATGATGACACGCTGAGCGGAATGCTGCTGGCGGCCTGGt                                               >  1:1239123/1‑150 (MQ=255)
                                                                                                        cGCTGCCTGCTGAGTTACGCCAGTTTCGTCACTGTTTTCAGGCCGAGCTCAATGGCGTTAAGACCGACTGGCGGCACTGGCTGGGTAAAGGCCCCGGATTAACGCCGAGTCATGATGACACGCTGAGCGGAATGCTGCTGGCGGCCTGGt                                               >  2:316300/1‑150 (MQ=255)
                                                                                                              ctgaagAGTTACGCCAGTTTCGACACTGTTTTCAGGCCGAGCTCAATGGCGTTAAGACCGACTGGCGGCACTGGCTGGGTAAAGCACCCGGATAAACGCCGAGTCATGATGACACGCTGAGCGGAATGCTGCTGGCGGCCGGGAATTATg                                         >  2:4865958/6‑150 (MQ=255)
                                                                                                                gctgAGTTACGCCAGTTTCGTCACTGTTTTCAGGCCGCGCTCAATGGCGTTAAGACCGACTGGCGGCACTGGCTGGGTAAAGGCCCCGGATTAACGCCGAGTCATGATGACACGCTGAGCGGAATGCTGCTGGCGGCCTGGTATTATGgc                                       >  2:708093/1‑150 (MQ=255)
                                                                                                                   gAGTTACGCCAGTTTCGTCACTGTTTTCAGGCCGCGCTCAATGGCGTTAAGACCGACTGGCGGCACTGGCTGGGTAAAGGCCCCGGATTAACGCCGAGTCATGATGACACGCTGAGCGGAATGCTGCTGGCGGCCTGGTATTATGGCGCt                                    <  1:4502291/150‑1 (MQ=255)
                                                                                                                   gAGTTACGCCAGTTTCGTCACTGTTTTCAGGCCGCGCTCAATGGCGTTAAGACCGACTGGCGGCACTGGCTGGGTAAAGGCCCCGGATTAACGCCGAGTCATGATGACACGCTGAGCGGAATGCTGCTGGCGGCCTGGTATTATGGCGCt                                    <  2:3348834/150‑1 (MQ=255)
                                                                                                                    aGTTACGCCAGTTTCGTCACTGTTTTCAGGCCGCGCTCAATGGCGTTAAGACCGACTGGCGGCACTGGCTGGGTAAAGGCCCCGGATTAACGCCGAGTCATGATGACACGCTGAGCGGAATGCTGCTGGCGGCCTGGTATTATGGCGCtt                                   <  2:2091720/150‑1 (MQ=255)
                                                                                                                     gTTACGCCAGTTTCGTCACTGTTTTCAGGCCGCGCTCAATGGCGTTAAGACCGACTGGCGGCACTGGCTGGGTAAAGGCCCCGGATTAACGCCGAGTCATGATGACACGCTGAGCGGAATGCTGCTGGCGGCCTGGTATTATGGCGCttt                                  <  1:825245/150‑1 (MQ=255)
                                                                                                                     gTTACGCCAGTTTCGTCACTGTTTTCAGGCCGAGCTCAATGGCGTTAAGACCGACTGGCGGCACTGGCTGGGTAAAGGCCCCGGATTAACGCCGAGTCATGATGACACGCTGAGCGGAATGCTGCTGGCGGCCTGGTATTATGGCGCttt                                  <  2:1156962/150‑1 (MQ=255)
                                                                                                                      ttACGCCAGTTTCGTCACTGTTTTCAGGCCGCGCTCAATGGCGTTAAGACCGACTGGCGGCACTGGCTGGGTAAAGGCCCCGGATTAACGCCGAGTCATGATGACACGCTGAGCGGAATGCTGCTGGCGGCCTGGTATTATGGCGCTTTa                                 <  2:2484575/150‑1 (MQ=255)
                                                                                                                         cGCCAGTTTCGTCACTGTTTTCAGGCCGCGCTCAATGGCGTTAAGACCGACTGGCGGCACTGGCTGGGTAAAGGCCCCGGATTAACGCCGAGTCATGATGACACGCTGAGCGGAATGCTGCTGGCGGCCTGGTATTATGGCGCTTTAGAt                              <  1:2369488/150‑1 (MQ=255)
                                                                                                                         cGCCAGTTTCGTCACTGTTTTCAGGCCGCGCTCAATGGCGTTAAGACCGACTGGCGGCACTGGCTGGGTAAAGGCCCCGGATTAACGCCGAGTCATGATGACACGCTGAGCGGAATGCTGCTGGCGGCCTGGTATTATGGCGCTTTAGAt                              >  2:1461355/1‑150 (MQ=255)
                                                                                                                         cGCCAGTTTCGTCACTGTTTTCAGGCCGCGCTCAATGGCGTTAAGACCGACTGGCGGCACTGGCTGGGTAAAGGCCCCGGATTAACGCCGAGTCATGATGACACGCTGAGCGGAATGCTGCTGGCGGCCTGGTATTATGGCGCTTTAGAt                              >  2:1261004/1‑150 (MQ=255)
                                                                                                                          gCCAGTTTCGTCACTGTTTTCAGGCCGCGCTCAATGGCGTTAAGACCGACTGGCGGCACTGGCTGGGTAAAGGCCCCGGATTAACGCCGAGTCATGATGACACGCTGAGCGGAATGCTGCTGGCGGCCTGGTATTATGGCGCTTTAGATg                             <  1:3494376/150‑1 (MQ=255)
                                                                                                                          gCCAGTTTCGTCACTGTTTTCAGGCCGAGCTCAATGGCGTTAAGACCGACTGGCGGCACTGGCTGGGTAAAGGCCCCGGATTAACGCCGAGTCATGATGACACGCTGAGCGGAATGCTGCTGGCGGCCTGGTATTATGGCGCTTTAGATg                             <  2:2373487/150‑1 (MQ=255)
                                                                                                                             aGTTTCGTCACTGTTTTCAGGCCGCGCTCAATGGCGTTAAGACCGACTGGCGGCACTGGCTGGGTAAAGGCCCCGGATTAACGCCGAGTCATGATGACACGCTGAGCGGAATGCTGCTGGCGGCCTGGTATTATGGCGCTTTAGATgcgc                          >  2:2453220/1‑150 (MQ=255)
                                                                                                                                ttCGTCACTGTTTTCAGGCCGCGCTCAATGGCGTTAAGACCGACTGGCGGCACTGGCTGGGTAAAGGCCCCGGATTAACGCCGAGTCATGATGACACGCTGAGCGGAATGCTGCTGGCGGCCTGGTATTATGGCGCTTTAGATGCGCGCt                       >  1:4857498/1‑150 (MQ=255)
                                                                                                                                ttCGTCACTGTTTTCAGGCCGAGCTCAATGGCGTTAAGACCGACTGGCGGCACTGGCTGGGTAAAGGCCCCGGATTAACGCCGAGTCATGATGACACGCTGAGCGGAATGCTGCTGGCGGCCTGGTATTATGGCGCTTTAGATGCGCGCt                       <  1:1951712/150‑1 (MQ=255)
                                                                                                                                 tCGTCACTGTTTTCAGGCCGCGCTCAATGGCGTTAAGACCGACTGGCGGCACTGGCTGGGTAAAGGCCCCGGATTAACGCCGAGTCATGATGACACGCTGAGCGGAATGCTGCTGGCGGCCTGGTATTATGGCGCTTTAGATGCGCGCTc                      <  2:3492718/150‑1 (MQ=255)
                                                                                                                                   gTCACTGTTTTCAGGCCGCGCTCAATGGCGTTAAGACCGACTGGCGGCACTGGCTGGGTAAAGGCCCCGGATTAACGCCGAGTCATGATGACACGCTGAGCGGAATGCTGCTGGCGGCCTGGTATTATGGCGCTTTAGATGCGCGCTCCg                    <  1:3264042/150‑1 (MQ=255)
                                                                                                                                   gTCACTGTTTTCAGGCCGAGCTCAATGGCGTTAAGACCGACTGGCGGCACTGGCTGGGTAAAGGCCCCGGATTAACGCCGAGTCATGATGACACGCTGAGCGGAATGCTGCTGGCGGCCTGGTATTATGGCGCTTTAGATGCGCGCTCCg                    <  1:4327815/150‑1 (MQ=255)
                                                                                                                                        tGTTTTCAGGCCGCGCTCAATGGCGTTAAGACCGACTGGCGGCACTGGCTGGGTAAAGGCCCCGGATTAACGCCGAGTCATGATGACACGCTGAGCGGAATGCTGCTGGCGGCCTGGTATTATGGCGCTTTAGATGCGCGCTCCGgtcgt               >  2:4606318/1‑150 (MQ=255)
                                                                                                                                        tGTTTTCAGGCCGCGCTCAATGGCGTTAAGACCGACTGGCGGCACTGGCTGGGTAAAGGCCCCGGATTAACGCCGAGTCATGATGACACGCTGAGCGGAATGCTGCTGGCGGCCTGGTATTATGGCGCTTTAGATGCGCGCTCCGgtcgt               >  2:2155662/1‑150 (MQ=255)
                                                                                                                                         gTTTTCAGGCCGCGCTCAATGGCGTTAAGACCGACTGGCGGCACTGGCTGGGTAAAGGCCCCGGATTAACGCCGAGTCATGATGACACGCTGAGCGGAATGCTGCTGGCGGCCTGGTATTATGGCGCTTTAGATGCGCGCTCCGgtcgtc              <  1:511321/150‑1 (MQ=255)
                                                                                                                                         gTTTTCAGGCCGAGCTCAATGGCGTTAAGACCGACTGGCGGCACTGGCTGGGTAAAGGCCCCGGATTAACGCCGAGTCATGATGACACGCTGAGCGGAATGCTGCTGGCGGCCTGGTATTATGGCGCTTTAGATGCGCGCTCCGgtcgtc              >  1:1544554/1‑150 (MQ=255)
                                                                                                                                           tttCAGGCCGCGCTCAATGGCGTTAAGACCGACTGGCGGCACTGGCTGGGTAAAGGCCCCGGATTAACGCCGAGTCATGATGACACGCTGAGCGGAATGCTGCTGGCGGCCTGGTATTATGGCGCTTTAGATGCGCGCTCCGGTCGTCCg            <  1:1652593/150‑1 (MQ=255)
                                                                                                                                               aGGCCGCGCTCAATGGCGTTAAGACCGACTGGCGGCACTGGCTGGGTAAAGGCCCCGGATTAACGCCGAGTCATGATGACACGCTGAGCGGAATGCTGCTGGCGGCCTGGTATTATGGCGCTTTAGATGCGCGCTCCGGTCGTCCGtttt        >  2:423023/1‑150 (MQ=255)
                                                                                                                                                 gCCGCGCTCAATGGCGTTAAGACCGACTGGCGGCACTGGCTGGGTAAAGGCCCCGGATTAACGCCGAGTCATGATGACACGCTGAGCGGAATGCTGCTGGCGGCCTGGTATTATGGCGCTTTAGATGCGCGCTCCGGTCGTCCGtttttt      <  1:4574127/150‑1 (MQ=255)
                                                                                                                                                 gCCGCGCTCAATGGCGTTAAGACCGACTGGCGGCACTGGCTGGGTAAAGGCCCCGGATTAACGCCGAGTCATGATGACACGCTGAGCGGAATGCTGCTGGCGGCCTGGTATTATGGCGCTTTAGATGCGCGCTCCGGTCGTCCGtttttt      <  1:4023704/150‑1 (MQ=255)
                                                                                                                                                     cgcTCAATGGCGTTAAGACCGACTGGCGGCACTGGCTGGGTAAAGGCCCCGGATTAACGCCGAGTCATGATGACACGCTGAGCGGAATGCTGCTGGCGGCCTGGTATTATGGCGCTTTAGATGCGCGCTCCGGTCGTCCGTTTTTTGCCt  <  2:2032690/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GCGTATTACGCCGCCTGCGCATCCTCAACCACTTGCAGCTGCATGGATGCAACGCGCGGAGGAAACCGGGCTTTTCGGGCCACTGGCGTTGGCGGCAAGCGATCCGCTGCCTGCTGAGTTACGCCAGTTTCGTCACTGTTTTCAGGCCGCGCTCAATGGCGTTAAGACCGACTGGCGGCACTGGCTGGGTAAAGGCCCCGGATTAACGCCGAGTCATGATGACACGCTGAGCGGAATGCTGCTGGCGGCCTGGTATTATGGCGCTTTAGATGCGCGCTCCGGTCGTCCGTTTTTTGCCT  >  CP009273/545373‑545671

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: