Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,510,074 C→A 15.9% R413L (CGA→CTA)  yfeA ← putative diguanylate cyclase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,510,0740CA15.9% 861.5 / 113.5 283R413L (CGA→CTA) yfeAputative diguanylate cyclase
Reads supporting (aligned to +/- strand):  ref base C (129/109);  new base A (22/23);  total (151/132)
Fisher's exact test for biased strand distribution p-value = 5.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.79e-01

Only 100 of 292 total aligned reads displayed.

CCGCCAGCCAGCTTAACTGCCCCAACAAGGGTTGCAGGGTTTCCTGATTTCCATCAAAACGCCCCCAGGCAGCACCATAGCCCATATCCAGCCCGGTATTGTTCCAGTGAATTTGCCGACTATTCAGGATGTTAACCATATGCTGGAGTCGCCCTTCCGTTTCCGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGTACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACT  >  CP009273/2509925‑2510222
                                                                                                                                                     |                                                                                                                                                    
ccGCCAGCCAGCTTAACTGCCCCAACAAGGGTTGCAGGGTTTCCTGATTTCCATCAAAACGCCCCCAGGCAGCACCATAGCCCATATCCAGCCCGGTATTGTTCCAGTGAATTTGCCGACTATTCAGGATGTTAACCATATGCTGGAGTc                                                                                                                                                      <  1:977889/150‑1 (MQ=255)
ccGCCAGCCAGCTTAACTGCCCCAACAAGGGTTGCAGGGTTTCCTGATTTCCATCAAAACGCCCCCAGGCAGCACCATAGCCCATATCCAGCCCGGTATTGTTCCAGTGAATTTGCCGACTATTCAGGATGTTAACCATATGCTGGAGTc                                                                                                                                                      <  1:176861/150‑1 (MQ=255)
ccGCCAGCCAGCTTAACTGCCCCAACAAGGGTTGCAGGGTTTCCTGATTTCCATCAAAACGCCCCCAGGCAGCACCATAGCCCATATCCAGCCCGGTATTGTTCCAGTGAATTTGCCGACTATTCAGGATGTTAACCATATGCTGGAGTc                                                                                                                                                      <  1:1548139/150‑1 (MQ=255)
  gccagccaGCTTAACTGCCCCAACATGGGTTGCAGGGTTTCCTGATTTCCATCAAAACGCCCCCAGGCAGCACCATAGCCCATATCCAGCCCGGTATTGTTCCAGTGAATTTGCCGACTATTCAGGATGTTAACCATATGCTGGAGTCGc                                                                                                                                                    <  1:408697/150‑1 (MQ=255)
   ccagccagCTTAACTGCCCCAACAAGGGTTGCAGGGTTTCCTGATTTCCATCAAAACGCCCCCAGGCAGCACCATAGCCCATATCCAGCCCGGTATTGTTCCAGTGAATTTGCCGACTATTCAGGATGTTAACCATATGCTGGAGTCGcc                                                                                                                                                   <  1:2341642/150‑1 (MQ=255)
       ccagcTTAACTGCCCCAACAAGGGTTGCAGGGTTTCCTGATTTCCATCAAAACGCCCCCAGGCAGCACCATAGCCCATATCCAGCCCGGTATTGTTCCAGTGAATTTGCCGACTATTCAGGATGTTAACCATATGCTGGAGTCGCCCTTc                                                                                                                                               >  2:479226/1‑150 (MQ=255)
       ccagcTTAACTGCCCCAACAAGGGTTGCAGGGTTTCCTGATTTCCATCAAAACGCCCCCAGGCAGCACCATAGCCCATATCCAGCCCGGTATTGTTCCAGTGAATTTGCCGACTATTCAGGATGTTAACCATATGCTGGAGTAGCCCTTc                                                                                                                                               >  1:1515705/1‑150 (MQ=255)
        cagcTTAACTGCCCCAACAAGGGTTGCAGGGTTTCCTGATTTCCATCAAAACGCCCCCAGGCAGCACCATAGCCCATATCCAGCCCGGTATTGTTCCAGTGAATTTGCCGACTATTCAGGATGTTAACCATATGCTGGAGTCGCCCTTcc                                                                                                                                              <  2:4114344/150‑1 (MQ=255)
        cagcTTAACTGCCCCAACAAGGGTTGCAGGGTTTCCTGATTTCCATCAAAACGCCCCCAGGCAGCACCATAGCCCATATCCAGCCCGGTATTGTTCCAGTGAATTTGCCGACTATTCAGGATGTTAACCATATGCTGGAGTCGCCCTTcc                                                                                                                                              <  2:1841111/150‑1 (MQ=255)
         agcTTAACTGCCCCAACAAGGGTTGCAGGGTTTCCTGATTTCCATCAAAACGCCCCCAGGCAGCACCATAGCCCATATCCAGCCCGGTATTGTTCCAGTGAATTTGCCGACTATTCAGGATGTTAACCATATGCTGGAGTCGCCCTTCCg                                                                                                                                             >  2:272119/1‑150 (MQ=255)
           cTTAACTGCCCCAACAAGGGTTGCAGGGTTTCCTGATTTCCAACAAAACGCCCCCAGGCAGCACCATAGCCCATATCCAGCCCGGTATTGTTCCAGTGAATTTGCCGACTATTCAGGATGTTAACCATATGCTGGAGTCGCCCTTCCGtt                                                                                                                                           >  2:2761730/1‑150 (MQ=255)
            ttAACTGCCCCAACAAGGGTTGCAGGGTTTCCTGATTTCCATCAAAACGCCCCCAGGCAGCACCATAGCCCATATCCAGCCCGGTATTGTTCCAGTGAATTTGCCGACTATTCAGGATGTTAACCATATGCTGGAGTCGCCCTTCCGttt                                                                                                                                          >  2:4501166/1‑150 (MQ=255)
             tAACTGCCCCAACAAGGGTTGCAGGGTTTCCTGATTTCCATCAAAACGCCCCCAGGCAGCACCATAGCCCATATCCAGCCCGGTATTGTTCCAGTGAATTTGCCGACTATTCGGGATGTTAACCATATGCTGGAGTCGCCCTTCCGTTTc                                                                                                                                         <  2:167803/150‑1 (MQ=255)
              aaCTGCCCCAACAAGGGTTGCAGGGTTTCCTGATTTCCATCAAAACGCCCCCAGGCAGCACCATAGCCCATATCCAGCCCGGTATTGTTCCAGTGAATTTGCCGACTATTCAGGATGTTAACCATATGCTGGAGTCGCCCTTCCGTTTcc                                                                                                                                        >  1:2394468/1‑150 (MQ=255)
              aaCTGCCCCAACAAGGGTTGCAGGGTTTCCTGATTTCCATCAAAACGCCCCCAGGCAGCACCATAGCCCATATCCAGCCCGGTATTGTTCCAGTGAATTTGCCGACTATTCAGGATGTTAACCATATGCTGGAGTCGCCCTTCCGTTTcc                                                                                                                                        <  2:1239234/150‑1 (MQ=255)
               aCTGCCCCAACAAGGGTTGCAGGGTTTCCTGATTTCCATCAAAACGCCCCCAGGCAGCACCATAGCCCATATCCAGCCCGGTATTGTTCCAGTGAATTTGCCGACTATTCAGGATGTTAACCATATGCTGGAGTCGCCCTTCCGTTTCCg                                                                                                                                       <  2:312589/150‑1 (MQ=255)
                          aaGGGTTGCAGGGTTTCCTGATTTCCATCAAAACGCCCCCAGGCAGCACCATAGCCCATATCCAGCCCGGTATTGTTCCAGTGAATTTGCCGACTATTCAGGATGTTAACCATATGCTGGAGTCGCCCTTCCGTTTCCGGCCCGCTCAGc                                                                                                                            >  1:598411/1‑150 (MQ=255)
                          aaGGGTTGCAGGGTTTCCTGATTTCCATCAAAACGCCCCCAGGCAGCACCATAGCCCATATCCAGCCCGGTATTGTTCCAGTGAATTTGCCGACTATTCAGGATGTTAACCATATGCTGGAGTCGCCCTTCCGTTTCCGGCCCGCTCAGc                                                                                                                            >  2:518072/1‑150 (MQ=255)
                          aaGGGTTGCAGGGTTTCCTGATTTCCATCAAAACGCCCCCAGGCAGCACCATAGCCCATATCCAGCCCGGTATTGTTCCAGTGAATTTGCCGACTATTCAGGATGTTAACCATATGCTGGAGTAGCCCTTCCGTTTCCGGCCCGCTCAGc                                                                                                                            >  1:2197081/1‑150 (MQ=255)
                           aGGGTTGCAGGGTTTCCTGATTTCCATCAAAACGCCCCCAGGCAGCACCATAGCCCATATCCAGCCCGGTATTGTTCCAGTGAATTTGCCGACTATTCAGGATGTTAACCATATGCTGGAGTCGCCCTTCCGTTTCCGGCCCGCTCAGCa                                                                                                                           <  2:2674249/150‑1 (MQ=255)
                             ggTTGCAGGGTTTCCTGATTTCCATCAAAACGCCCCCAGGCAGCACCATAGCCCATATCCAGCCCGGTATTGTTCCAGTGAATTTGCCGACTATTCAGGATGTTAACCATATGCTGGAGTCGCCCTTCCGTTTCCGGCCCGCTCAGCACt                                                                                                                         >  1:308574/1‑150 (MQ=255)
                                tGCAGGGTTTCCTGATTTCCATCAAAACGCCCCCAGGCAGCACCATAGCCCATATCCAGCCCGGTATTGTTCCAGTGAATTTGCCGACTATTCAGGATGTTAACCATATGCTGGAGTCGCCCTTCCGTTTCCGGCCCGCTCAGCACTAAc                                                                                                                      <  2:1414307/150‑1 (MQ=255)
                                     ggTTTCCTGATTTCCATCAAAACGCCCCCAGGCAGCACCATAGCCCATATCCAGCCCGGTATTGTTCCAGTGAATTTGCCGACTATTCAGGATGTTAACCATATGCTGGAGTCGCCCTTCCGTTTCCGGCCCGCTCAGCACTAAcagcag                                                                                                                 <  2:3966985/150‑1 (MQ=255)
                                      gTTTCCTGATTTCCATCAAAACGCCCCCAGGCAGCACCATAGCCCATATCCAGCCCGGTATTGTTCCAGTGAATTTGCCGACTATTCAGGATGTTAACCATATGCTGGAGTCGCCCTTCCGTTTCCGGCCCGCTCAGCACTAACAGCAGt                                                                                                                <  2:1506208/150‑1 (MQ=255)
                                      gTTTCCTGATTTCCATCAAAACGCCCCCAGGCAGCACCATAGCCCATATCCAGCCCGGTATTGTTCCAGTGAATTTGCCGACTATTCAGGATGTTAACCATATGCTGGAGTCGCCCTTCCGTTTCCGGCCCGCTCAGCACTAACAGCAGt                                                                                                                >  2:3452611/1‑150 (MQ=255)
                                        ttCCTGATTTCCATCAAAACGCCCCCAGGCAGCACCATAGCCCATATCCAGCCCGGTATTGTTCCAGTGAATTTGCCGACTATTCAGGATGTTAACCATATGCTGGAGTCGCCCTTCCGTTTCCGGCCCGCTCAGCACTAACAGCAGTTc                                                                                                              <  2:1714713/150‑1 (MQ=255)
                                           cTGATTTCCATCAAAACGCCCCCAGGCAGCACCATAGCCCATATCCAGCCCGGTATTGTTCCAGTGAATTTGCCGACTATTCAGGATGTTAACCATATGCTGGAGTCGCCCTTCCGTTTCCGGCCCGCTCAGCACTAACAGCAGTTCACt                                                                                                           >  1:1465314/1‑150 (MQ=255)
                                           cTGATTTCCATCAAAACGCCCCCAGGCAGCACCATAGCCCATATCCAGCCCGGTATTGTTCCAGTGAATTTGCCGACTATTCAGGATGTTAACCATATGCTGGAGTCGCCCTTCCGTTTCCGGCCCGCTCAGCACTAACAGCAGTTCACt                                                                                                           >  1:2773895/1‑150 (MQ=255)
                                            tGATTTCCATCAAAACGCCCCCAGGCAGCACCATAGCCCATATCCAGCCCGGTATTGTTCCAGTGAATTTGCCGACTATTCAGGATGTTAACCATATGCTGGAGTCGCCCTTCCGTTTCCGGCCCGCTCAGCACTAACAGCAGTTCACTa                                                                                                          >  1:3186847/1‑150 (MQ=255)
                                             gATTTCCATCAAAACGCCCCCAGGCAGCACCATAGCCCATATCCAGCCCGGTATTGTTCCAGTGAATTTGCCGACTATTCAGGATGTTAACCATATGCTGGAGTCGCCCTTCCGTTTCCGGCCCGCTCAGCACTAACAGCAGTTCACTAc                                                                                                         <  1:1345470/150‑1 (MQ=255)
                                                ttCCATCAAAACGCCCCCAGGCAGCACCATAGCCCATATCCAGCCCGGTATTGTTCCAGTGAATTTGCCGACTATTCAGGATGTTAACCATATGCTGGAGTCGCCCTTCCGTTTCCGGCCCGCTCAGCACTAACAGCAGTTCACTACCCg                                                                                                      >  2:1537616/1‑150 (MQ=255)
                                                 tCCATCAAAACGCCCCCAGGCAGCACCATAGCCCATATCCAGCCCGGTATTGTTCCAGTGAATTTGCCGACTATTCAGGATGTTAACCATATGCTGGAGTCGCCCTTCCGTTTCCGGCCCGCTCAGCACTAACAGCAGTTCACTACCCgg                                                                                                     <  1:3940285/150‑1 (MQ=255)
                                                   cATCAAAACGCCCCCAGGCAGCACCATAGCCCATATCCAGCCCGGTATTGTTCCAGTGAATTTGCCGACTATTCAGGATGTTAACCATATGCTGGAGTCGCCCTTCCGTTTCCGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCa                                                                                                   <  2:3764065/150‑1 (MQ=255)
                                                    aTCAAAACGCCCCCAGGCAGCACCATAGCCCATATCCAGCCCGGTATTGTTCCAGTGAATTTGCCGACTATTCAGGATGTTAACCATATGCTGGAGTCGCCCTTCCGTTTCCGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCaa                                                                                                  <  1:1630439/150‑1 (MQ=255)
                                                     tCAAAACGCCCCCAGGCAGCACCATAGCCCATATCCAGCCCGGTATTGTTCCAGTGAATTTGCCGACTATTCAGGATGTTAACCATATGCTGGAGTCGCCCTTCCGTTTCCGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAAt                                                                                                 <  2:3931461/150‑1 (MQ=255)
                                                     tCAAAACGCCCCCAGGCAGCACCATAGCCCATATCCAGCCCGGTATTGTTCCAGTGAATTTGCCGACTATTCAGGATGTTAACCATATGCTGGAGTCGCCCTTCCGTTTCCGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAAt                                                                                                 <  1:4168221/150‑1 (MQ=255)
                                                     tCAAAACGCCCCCAGGCAGCACCATAGCCCATATCCAGCCCGGTATTGTTCCAGTGAATTTGCCGACTATTCAGGATGTTAACCATATGCTGGAGTCGCCCTTCCGTTTCCGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAAt                                                                                                 >  1:4322769/1‑150 (MQ=255)
                                                      cAAAACGCCCCCAGGCAGCACCATAGCCCATATCCAGCCCGGTATTGTTCCAGTGAATTTGCCGACTATTCAGGATGTTAACCATATGCTGGAGTCGCCCTTCCGTTTCCGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAAtt                                                                                                <  1:1040881/150‑1 (MQ=255)
                                                      cAAAACGCCCCCAGGCAGCACCATAGCCCATATCCAGCCCGGTATTGTTCCAGTGAATTTGCCGACTATTCAGGATGTTAACCATATGCTGGAGTAGCCCTTCCGTTTCCGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAAtt                                                                                                >  1:3416331/1‑150 (MQ=255)
                                                      cAAAACGCCCCCAGGCAGCACCATAGCCCATATCCAGCCCGGTATTGTTCCAGTGAATTTGCCGACTATTCAGGATGTTAACCATATGCTGGAGTAGCCCTTCCGTTTCCGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAAtt                                                                                                <  1:1604122/150‑1 (MQ=255)
                                                         aaCGCCCCCAGGCAGCACCATAGCCCATATCCAGCCCGGTATTGTTCCAGTGAATTTGCCGACTATTCAGGATGTTAACCATATGCTGGAGTCGCCCTTCCGTTTCCGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGAt                                                                                             >  1:2754002/1‑150 (MQ=255)
                                                         aaCGCCCCCAGGCAGCACCATAGCCCATATCCAGCCCGGTATTGTTCCAGTGAATTTGCCGACTATTCAGGATGTTAACCATATGCTGGAGTCGCCCTTCCGTTTCCGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGAt                                                                                             >  2:2178180/1‑150 (MQ=255)
                                                           cGCCCCCAGGCAGCACCATAGCCCATATCCAGCCCGGTATTGTTCCAGTGAATTTGCCGACTATTCAGGATGTTAACCATATGCTGGAGTCGCCCTTCCGTTTCCGGCCCGCt                                                                                                                                <  1:2013578/113‑1 (MQ=255)
                                                             cccccAGGCAGCACCATAGCCCATATCCAGCCCGGTATTGTTCCAGTGAATTTGCCGACTATTCAGGATGTTAACCATATGCTGGAGTCGCCCTTCCGTTTCCGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACaa                                                                                         <  1:1814848/150‑1 (MQ=255)
                                                                 cAGGCAGCACCATAGCCCATATCCAGCCCGGTATTGTTCCAGTGAATTTGCCGACTATTCAGGATGTTAACCATATGCTGGAGTCGCCCTTCCGTTTCCGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACttt                                                                                     >  2:1155672/1‑150 (MQ=255)
                                                                  aGGCAGCACCATAGCCCATATCCAGCCCGGTATTGTTCCAGTGACTTTGCCGACTATTCAGGATGTTAACCATATGCTGGAGTCGCCCTTCCGTTTCCGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACtttt                                                                                    >  2:3445914/1‑150 (MQ=255)
                                                                  aGGCAGCACCATAGCCCATATCCAGCCCGGTATTGTTCCAGTGAATTTGCCGACTATTCAGGATGTTAACCATATGCTGGAGTCGCCCTTCCGTTTCCGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACtttt                                                                                    >  1:1069296/1‑150 (MQ=255)
                                                                  aGGCAGCACCATAGCCCATATCCAGCCCGGTATTGTTCCAGTGAATTTGCCGACTATTCAGGATGTTAACCATATGCTGGAGTCGCCCTTCCGTTTCCGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACtttt                                                                                    <  1:2970961/150‑1 (MQ=255)
                                                                    gcagcaCCATAGCCCATATCCAGCCCGGTATTGTTCCAGTGAATTTGCCGACTATTCAGGATGTTAACCATATGCTGGAGTAGCCCTTCCGTTTCCGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCg                                                                                  <  2:2706277/150‑1 (MQ=255)
                                                                         aCCATAGCCCATATCCAGCCCGGTATTGTTCCAGTGAATTTGCCGACTATTCAGGATGTTAACCATATGCTGGAGTCGCCCTTCCGTTTCCGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTc                                                                             >  2:4081565/1‑150 (MQ=255)
                                                                          ccATAGCCCATATCCAGCCCGGTATTGTTCCAGTGAATTTGCCGACTATTCAGGATGTTAACCATATGCTGGAGTCGCCCTTCCGTTTCCGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTcc                                                                            <  1:4510718/150‑1 (MQ=255)
                                                                           cATAGCCCATATCCAGCCCGGTATTGTTCCAGTGAATTTGCCGACTATTCAGGATGTTAACCATATGCTGGAGTCGCCCTTCCGTTTCCGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCt                                                                           >  1:4393253/1‑150 (MQ=255)
                                                                              aGCCCATATCCAGCCCGGTATTGTTCCAGTGAATTTGCCGACTATTCAGGATGTTAACCATATGCTGGAGTCGCCCTTCCGTTTCCGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCa                                                                        >  1:2085996/1‑150 (MQ=255)
                                                                                cccATATCCAGCCCGGTATTGTTCCAGTGAATTTGCCGACTATTCAGGATGTTAACCATATGCTGGAGTCGCCCTTCCGTTTCCGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCACTTGATACAACTTTTCGTTTTCCTGCCTc                                                                      >  2:3536391/1‑150 (MQ=255)
                                                                                   atatCCAGCCCGGTATTGTTCCAGTGAATTTGCCGACTATTCAGGATGTTAACCATATGCTGGAGTCGCCCTTCCGTTTCCGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAAc                                                                   <  2:2038706/150‑1 (MQ=255)
                                                                                   atatCCAGCCCGGTATTGTTCCAGTGAATTTGCCGACTATTCAGGATGTTAACCATATGCTGGAGTCGCCCTTCCGTTTCCGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAAc                                                                   <  2:232269/150‑1 (MQ=255)
                                                                                   atatCCAGCCCGGTATTGTTCCAGTGAATTTGCCGACTATTCAGGATGTTAACCATATGCTGGAGTCGCCCTTCCGTTTCCGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAAc                                                                   >  2:302464/1‑150 (MQ=255)
                                                                                   atatACTGCCCTGTATTGTTCCAGTGAATTTGCCTAATTTTCAGTTTGTTAAACATATTATGGAGTCTCCATTGCTTTTCCTGCACTATAAGCACTAACATCATTTCACTACCCTTAATTTGATACAACTTTTAGTTTTCCTTCATCaaa                                                                   <  2:3062205/150‑2 (MQ=255)
                                                                                    tatCCAGCCCGGTATTGTTCCAGTGAATTTGCCGACTATTCAGGATGTTAACCATATGCTGGAGTCGCCCTTCCGTTTCCGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACg                                                                  <  1:836897/150‑1 (MQ=255)
                                                                                    tatCCAGCCCGGTATTGTTCCAGTGAATTTGCCGACTATTCAGGATGTTAACCATATGCTGGAGTCGCCCTTCCGTTTCCGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACg                                                                  >  2:3843009/1‑150 (MQ=255)
                                                                                       ccAGCCCGGTATTGTTCCAGTGAATTTGCCGACTATTCAGGATGTTAACCATATGCTGGAGTCGCCCTTCCGTTTCCGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGgca                                                               <  1:688026/150‑1 (MQ=255)
                                                                                       ccAGCCCGGTATTGTTCCAGTGAATTTGCCGACTATTCAGGATGTTAACCATATGCTGGAGTCGCCCTTCCGTTTCCGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGgca                                                               <  1:688471/150‑1 (MQ=255)
                                                                                       ccAGCCCGGTATTGTTCCAGTGAATTTGCCGACTATTCAGGATGTTAACCATATGCTGGAGTAGCCCTTCCGTTTCCGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGgca                                                               >  1:3712499/1‑150 (MQ=255)
                                                                                        cAGCCCGGTATTGTTCCAGTGAATTTGCCGACTATTCAGGATGTTAACCATATGCTGGAGTAGCCCTTCCGTTTCCGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGgcag                                                              <  2:2806098/150‑1 (MQ=255)
                                                                                         aGCCCGGTATTGTTCCAGTGAATTTGCCGACTATTCAGGATGTTAACCATATGCTGGAGTCGCCCTTCCGTTTCCGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGgcagc                                                             >  2:67542/1‑150 (MQ=255)
                                                                                          gCCCGGTATTGTTCCAGTGAATTTGCCGACTATTCAGGATGTTAACCATATGCTGGAGTCGCCCTTCCGTTTCCGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGgcagca                                                            >  1:4263964/1‑150 (MQ=255)
                                                                                            ccGGTATTGTTCCAGTGAATTTGCCGACTATTCAGGATGTTAACCATATGCTGGAGTCGCCCTTCCGTTTCCGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCagcagc                                                          >  1:877545/1‑150 (MQ=255)
                                                                                              ggTATTGTTCCAGTGAATTTGCCGACTATTCAGGATGTTAACCATATGCTGGAGTAGCCCTTCCGTTTCCGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGt                                                        >  1:445934/1‑150 (MQ=255)
                                                                                               gTATTGTTCCAGTGAATTTGCCGACTATTCAGGATGTTAACCATATGCTGGAGTAGCCCTTCCGTTTCCGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGTa                                                       <  1:1179941/150‑1 (MQ=255)
                                                                                                    gTTCCAGTGAATTTGCCGACTATTCAGGATGTTAACCATATGCTGGAGTCGCCCTTCCGTTTCCGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGTACGGCa                                                  >  1:2427482/1‑150 (MQ=255)
                                                                                                    gTTCCAGTGAATTTGCCGACTATTCAGGATGTTAACCATATGCTGGAGTCGCCCTTCCGTTTCCGGCCCGCTAAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACTGCAGCAGCGTACGGCa                                                  <  1:4186897/150‑1 (MQ=255)
                                                                                                      tCCAGTGAATTTGCCGACTATTCAGGATGTTAACCATATGCTGGAGTCGCCCTTCCGTTTCCGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGTACGGCaaa                                                >  2:423843/1‑150 (MQ=255)
                                                                                                       ccAGTGAATTTGCCGACTATTCAGGATGTTAACCATATGCTGGAGTCGCCCTTCCGTTTCCGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGTACGGCAAAt                                               <  1:1925247/150‑1 (MQ=255)
                                                                                                          gTGAATTTGCCGACTATTCAGGATGTTAACCATATGCTGGAGTCACCCTTCCGTTTCCGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGTACGGCAAATTGa                                            >  2:1145696/1‑150 (MQ=255)
                                                                                                            gAATTTGCCGACTATTCAGGATGTTAACCATATGCTGGAGTCGCCCTTCCGTTTCCGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGTACGGCAAATTGAGc                                          <  2:949424/150‑1 (MQ=255)
                                                                                                             aaTTTGCCGACTATTCAGGATGTTAACCATATGCTGGAGTCGCCCTTCCGTTTCCGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGTACGGCAAATTGAGCg                                         >  2:1447952/1‑150 (MQ=255)
                                                                                                             aaTTTGCCGACTATTCAGGATGTTAACCATATGCTGGAGTCGCCCTTCCGTTTCCGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGTACGGCAAATTGAGCg                                         >  1:44850/1‑150 (MQ=255)
                                                                                                                ttGCCGACTATTCAGGATGTTAACCATATGCTGGAGTAGCCCTTCCGTTTCCGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGTACGGCAAATTGAGCGGAt                                      <  1:3482600/150‑1 (MQ=255)
                                                                                                                 tGCCGACTATTCAGGATGTTAACCATATGCTGGAGTCGCCCTTCCGTTTCCGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGTACGGCAAATTGAGCGGATa                                     <  2:3299417/150‑1 (MQ=255)
                                                                                                                   ccGACTATTCAGGATGTTAACCATATGCTGGAGTCGCCCTTCCGTTTCCGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGTACGGCAAATTGAGCGGATACa                                   <  1:3129089/150‑1 (MQ=255)
                                                                                                                     gACTATTCAGGATGTTAACCATATGCTGGAGTCGCCCTTCCGTTTCCGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGTACGGCAAATTGAGCGGATACAGt                                 <  2:1587767/150‑1 (MQ=255)
                                                                                                                      aCTATTCAGGATGTTAACCATATGCTGGAGTCGCCCTTCCGTTTCCGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGTACGGCAAATTGAGCGGATACAGTg                                <  1:2770619/150‑1 (MQ=255)
                                                                                                                       cTATTCAGGATGTTAACCATATGCTGGAGTCGCCCTTCCGTTTCCGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGTACGGCAAATTGAGCGGATACAGTGa                               <  2:2899109/150‑1 (MQ=255)
                                                                                                                         aTTCAGGATGTTAACCATATGCTGGAGTAGCCCTTCCGTTTCCGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGTACGGCAAATTGAGCGGATACAGTGAAc                             >  2:464304/1‑150 (MQ=255)
                                                                                                                         aTTCAGGATGTTAACCATATGCTGGAGTAGCCCTTCCGTTTCCGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGTACGGCAAATTGAGCGGATACAGTGAAc                             <  1:3874056/150‑1 (MQ=255)
                                                                                                                                 tGTTAACCATATGCTGGAGTCGCCCTTCCGTTTCCGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGTACGGCAAATTGAGCGGATACAGTGAACGCGcatca                     >  2:2923279/1‑150 (MQ=255)
                                                                                                                                   ttAACCATATGCTGGAGTAGCCCTTCCGTTTCCGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGTACGGCAAATTGAGCGGATACAGTGAACGCGCATCAtt                   <  2:2423275/150‑1 (MQ=255)
                                                                                                                                    tAACCATATGCTGGAGTAGCCCTTCCGTTTCCGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGTACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTa                  >  2:729594/1‑150 (MQ=255)
                                                                                                                                    tAACCATATGCTGGAGTAGCCCTTCCGTTTCCGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGTACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTa                  <  2:3728576/150‑1 (MQ=255)
                                                                                                                                      aCCATATGCTGGAGTCGCCCTTCCGTTTCCGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGTACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAg                >  1:986394/1‑150 (MQ=255)
                                                                                                                                       ccATATGCTGGAGTCGCCCTTCCGTTTCCGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGTACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGc               >  1:2581607/1‑150 (MQ=255)
                                                                                                                                       ccATATGCTGGAGTAGCCCTTCCGTTTCCGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGTACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGc               >  2:109485/1‑150 (MQ=255)
                                                                                                                                         atatGCTGGAGTCGCCCTTCCGTTTCCGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGTACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCg             <  1:3547100/150‑1 (MQ=255)
                                                                                                                                         atatGCTGGAGTCGCCCTTCCGTTTCCGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGTACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCg             <  2:1314339/150‑1 (MQ=255)
                                                                                                                                         atatGCTGGAGTCGCCCTTCCGTTTCCGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGTACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCg             <  2:2019284/150‑1 (MQ=255)
                                                                                                                                          tatGCTGGAGTCGCCCTTCCGTTTCCGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGTACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGt            <  2:1264687/150‑1 (MQ=255)
                                                                                                                                          tatGCTGGAGTCGCCCTTCCGTTTCCGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGTACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGt            >  1:2708831/1‑150 (MQ=255)
                                                                                                                                          tatGCTGGAGTCGCCCTTCCGTTTCCGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGTACGGCAAATTGAGCGGATACAGAGAACGCGCATCATTAAGCCGt            >  1:3153950/1‑150 (MQ=255)
                                                                                                                                            tGCTGGAGTCGCCCTTCCGTTTCCGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGTACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTaa          <  2:4488186/150‑1 (MQ=255)
                                                                                                                                             gCTGGAGTCGCCCTTCCGTTTCCGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGTACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAAt         >  2:1437489/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CCGCCAGCCAGCTTAACTGCCCCAACAAGGGTTGCAGGGTTTCCTGATTTCCATCAAAACGCCCCCAGGCAGCACCATAGCCCATATCCAGCCCGGTATTGTTCCAGTGAATTTGCCGACTATTCAGGATGTTAACCATATGCTGGAGTCGCCCTTCCGTTTCCGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGTACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACT  >  CP009273/2509925‑2510222

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: