Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,964,757 A→T 5.5% Q309L (CAG→CTG)  yegE → phosphodiesterase domain‑containing diguanylate cyclase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,964,7570AT 5.5% 225.6 / 3.9 74Q309L (CAG→CTG) yegEphosphodiesterase domain‑containing diguanylate cyclase
Reads supporting (aligned to +/- strand):  ref base A (37/32);  new base T (2/2);  total (39/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.75e-01

GCAGGAAGAAAAAGAGCGTGCGCAGATTACCCTGCAGGCCATCGGGGATGCGGTGATTACAACCGATGTCGAGGGCAATATCGGCTACATGAACCCGGCCGCCGAGCAGTTGACCCACTGGCAATCCGGCCAGGCCCAGGGCCTGCCGCTGTCGGCACTGTTCAGCCTGGTGGATGAGCACGCCGAGGAGGATGGTCGCAGCTTGGTCGAGCAGGTGCTCAGCGGCAGCCTCAAGGGCGGCTCCGAGCATGCCCGGCTGATCCAGCGCCTGGACGGCAGTACCGT  >  NC_002947/1964620‑1964904
                                                                                                                                         |                                                                                                                                                   
gCAGGAAGAAAAAGAGCGTGCGCAGATTACCCTGCAGGCCATCGGGGATGCGGTGATTACAACCGATGTCGAGGGCAATATCGGCTACATGAACCCGGCCGCCGAGCAGTTGACCCACTGGCAATCCGGCAGGCCCTGGGCCTggccgga                                                                                                                                        <  2:315547/150‑7 (MQ=255)
   ggAAGAAAAAGAGCGTGCGCAGATTACCCTGCAGGCCATCGGGGATGCGGTGATTACAACCGATGTCGAGGGCAATATCGGCTACATGAACCCGGCCGCCGAGCAGTTGACCCACTGGCAATCCGGCCAGGCCCAGGGCCTGCCGCTGTc                                                                                                                                      >  1:1195157/1‑150 (MQ=255)
    gaagaaAAAGAGCGTGCGCAGATTACCCTGCAGGCCATCGGGGATGCGGTGATTACAACCGATGTCGAGGGCAATATCGGCTACATGAACCCGGCCGCCGAGCAGTTGACCCACTGGCAATCCGGCCAGGCCCAGGGCCTGCCGCTGTCg                                                                                                                                     >  2:890703/1‑150 (MQ=255)
      agaaAAAGAGCGTGCGCAGATTACCCTGCAGGCCATCGGGGATGCGGTGATTACAACCGATGTCGAGGGCAATATCGGCTACATGAACCCGGCCGCCGAGCAGTTGACCCACTGGCACCCCGCCCAGGCCCAGGGCCTGCCGCTGTCGGc                                                                                                                                   >  2:654237/1‑150 (MQ=255)
       gaaAAAGAGCGTGCGCAGATTACCCTGCAGGCCATCGGGGATGCGGTGATTACAACCGATGTCGAGGGCAATATCGGCTAGATGAAGGAGGCGCACAAGCATGTGACCCACTGGCAATCCGGCCAGGCCCAGGGCCTGCCGCTGTCGGCa                                                                                                                                  >  1:151898/1‑150 (MQ=255)
          aaaGAGCGTGCGCAGATTACCCTGCAGGCCATCGGGGATGCGGTGATTACAACCGATGTCGAGGGCAATATCGGCTACATGAACCCGGCCGCCGAGCAGTTGACCCACTGGCAATCCGGCCAGGCCCAGGGCCTGCCGCTGTCGGCACTg                                                                                                                               >  1:717466/1‑150 (MQ=255)
          aaaGAGCGTGCGCAGATTACCCTGCAGGCCATCGGGGATGCGGTGATTACAACCGATGTCGAGGGCAATATCGGCTACATGAACCCGGCCGCCGAGCAGTTGACCCACTGGCAATCCGGCCAGGCCCAGGGCCTGCCGCTGTCGGCACTg                                                                                                                               >  1:167691/1‑150 (MQ=255)
           aatagCGTTGGCAGATTACCCTGCAGGCCATCCGGGATGCGGTGATTACAACCGCTGTCGAGGGCAATATCGGCTACAGGACTCCGGCCGCCGAGCAGTTGACCCACTGGCAATCCGGCCCGGCCCAGGGCCTGCCGCTGTCGGCACTGt                                                                                                                              <  2:569650/147‑1 (MQ=255)
             gagCGTGCGCAGATTACCCTGCAGGCCATCGGGGATGCGGTGATTACAACCGATGTCGAGGGCAATATCGGCTACATGAACCCGGCCGCCGAGCAGTTGACCCACTGGCAATCCGGCCAGGCCCAGGGCCTGCCGCTGTCGGCACTGTTc                                                                                                                            <  1:582238/150‑1 (MQ=255)
                     gcAGATTACCCTGCAGGCCATCGGGGANGCGGTGATTACANCCGATGNCGAGGGCAATATCGGCTACATGAACCCGGCCGCCGAGCAGTTGACCCACTGGCAATCCGGCCAGGCCCAGGGCCTGCCGCTGTCGGCACTGTTCAGCCtggt                                                                                                                    >  1:1144882/1‑150 (MQ=255)
                            aCCCTGCAGGCCATCGGGGATGCGGTGATTACAACCGATGTCGAGGGCAATATCGGCTACATGAACCCGGCCGCCGAGCAGTTGACCCACTGGCAATCCGGCCAGGCCCAGGGCCTGCCGCTGTCGTCACTGTTCAGCCTGGTGGATGAg                                                                                                             <  1:772660/150‑1 (MQ=255)
                                tGCAGGCCATCGGGGATGCGGTGATTACAACCGATGTCGAGGGCAATATCGGCTACATGAACCCGGCCGCCGAGCAGTTGACCCACTGGCAATCCGGCCAGGCCCAGGGCCTGCCGCTGTCGGCACTGTTCAGCGTGGTGGATGAGCACg                                                                                                         >  2:895329/1‑150 (MQ=255)
                                   gggCCATCGGGGATGCGGTGATTACAACCGATGTCGAGGGCAATATCGGCTACATGAACCCGGCCGCCGAGCAGTTGACCCACTGGCAATCCGGCCAGGCCCAGGGCCTGCCGCTGTCGGCACTGTTCAGCCTGGTGGATGAGCACGCCg                                                                                                      >  2:1217515/2‑150 (MQ=255)
                                        aTCGGGGATGCGGTGATTACAACCGATGTCGAGGGCAATATCGGCTACATGAACCCGGCCGCCGAGCAGTTGACCCACTGGCAATCCGGCCAGGCCCAGGGCCTGCCGCTGTCGGCACTGTTCAGCCTGGTGGATGAGCACGCCgaggag                                                                                                 <  2:717466/150‑1 (MQ=255)
                                          cGGGGATGCGGTGATTACAACCGATGTCGAGGGCAATATCGGCTACATGAACCCGGCCGCCGAGCAGTTGACCCACTGGCAATCCGGCCAGGCCCAGGGCCTGCCGCTGTCGGCACTGTTCAGCCTGGTGGATGAGCACGCCGAggagga                                                                                               <  1:906660/150‑1 (MQ=255)
                                               aTGCGGTGATTACAACCGATGTCGAGGGCAATATCGGCTACATGAACCCGGCCGCCGAGCAGTTGACCCACTGGCAATCCGGCCAGGCCCAGGGCCTGCCGCTGTCGGCACTGTTCAGCCTGGTGGATGAGCACGCCGAGGAGGATGGTc                                                                                          <  1:890703/150‑1 (MQ=255)
                                                   ggggATTACAACCGATGTCGAGGGCAATATCGGCTACATGAACCCGGCCGCCGAGCAGTTGACCCACTGGCAATCCGGCCAGGCCAAGGGCCTGCCGCTGTCGGCACTGTTCAGCCTGGTGGATGAGCACGCCGAGGAGGATGGTCGCAg                                                                                      >  1:274172/4‑150 (MQ=255)
                                                    gTGATTACAACCGATGTCGAGGGCAATATCGGCTACATGAACCCGGCCGCCGAGCAGTTGACCCACTGGCAATCCGGCCAGGCCCAGGGCCTGCCGCTGTCGGCACTGTTCAGCCTGGTGGATGAGCACGCCGAGGAGGATGGTCGCAGc                                                                                     <  1:170447/150‑1 (MQ=255)
                                                      gatctcAACCGATGTCGAGGGCAATATCGGCTACATGAACCCGGCCGCCGAGCAGTTGACCCACTGGCAATCCGGCCAGGCCCAGGGCCTGCCGCTGTCGGCACTGTTCAGCCTGGTGGATGAGCACGCCGAGGAGGATGGTCGCAGCtt                                                                                   <  2:989552/145‑1 (MQ=255)
                                                          aCAACCGATGTCGAGGGCAATATCGGCTACATGAACCCGGCCGCCGAGCAGTTGACCCACTGGCAATCCGGCCAGGCCCAGGGCCTGCCGCTGTCGGCACTGTTCAGCCTGGTGGATGAGCACGCCGAGGAGGATGGTCGCAGCTTGGTc                                                                               >  2:890941/1‑150 (MQ=255)
                                                            aaCCGATGTCGAGGGCAATATCGGCTACATGAACCCGGCCGCCGAGCAGTTGACCCACTGGCAATCCGGCCCGGCATAGGGCCTGCCGCTGTCGGCCCTGTTAAGCCTGGGAGATGAGCAGGCTTAGGAGGGTGGTCGCAGCTTGGTCGa                                                                             >  1:1248752/1‑150 (MQ=255)
                                                            aaCCGATGTCGAGGGCAATATCGGCTACATGAACCAGGCAGCAGAGCAGCCGACCCACTGGCAATCCGGCCAGGCCCAGGGCCTGCCGCTGTCGGCACTGTTCAGCCTGGTGGATGAGCACGCCGAGGAGGATGGTCGCAGCTTgagatc                                                                             >  1:989552/1‑145 (MQ=255)
                                                                gATGTCTAGTGCAATATCGGCTACATGAACCCGGCCGCCGAGCAGTTGACCCACTGGCAATCCGGCCAGGCCCAGGGCCTGCCGCTGTCGGCACTGTTCAGCCTGGTGGATGAGCACGCCGAGGAGGATGGTCGCAGCTTGGTCGAGCAg                                                                         <  2:753439/150‑1 (MQ=255)
                                                                gATGTCGAGGGCAATATCGGCTACATGAACCCGGCCGCCGAGCAGTTGACCCACTGGCAATCCGGCCAGGCCCAGGGCCTGCCGCTGTCGGCACTGTTCAGCCTGGTGGATGAGCACGCCGAGGAGGATGGTCGCAGCTTGGTCGAGCAg                                                                         <  2:883759/150‑1 (MQ=255)
                                                                gATGTCGAGGGCAATATCGGCTACATGAACCCGGCCGCCGAGCAGTTGACCCACTGGCAATCCGGCCAGGCCCAGGGCCTGCCGCTGTCGGCACTGTTCAGCCTGGTGGATGAGCACGCCGAGGAGGATGGTCGCAGCTTGGTCGAGCAg                                                                         <  1:1251568/150‑1 (MQ=255)
                                                                gATGTCGAGGGCAATATCGGCTACATGAACCCGGCCGCCGAGCAGTTGACCCACTGGCAATCCGGCCAGGCCCAGGGCCTGCCGCTGTCGGCACTGTTCAGCCTGGTGGATGAGCACGCCGAGGAGGATGGTCGCAGCTTGGTCGAGCAg                                                                         <  1:526946/150‑1 (MQ=255)
                                                                gATGTCGAGGGCAATATCGGCTACATGAACCCGGCCGCCGAGCAGTTGACCCACTGGCAATCCGGCCAGGCCAAGGGCCTGCCGCTGTCGGCACTGTTCAGCCTGGTGGATGAGCACGCCGAGGAGGATGGTCGCAGCTTGGTCGAACAg                                                                         <  2:292638/150‑1 (MQ=255)
                                                                            aaTATCGGCTACATGAACCCGGCCGCCGAGCAGTTGACCCACTGGCAATCCGGCCAGGCCCAGGGCCTGCCGCTGTCGGCACTGTTCAGCCTGGTGGATGAGCACGCCGAGGAGGATGGTCGCAGCTTGGTCGAGCAGGTGCTCAGCGGc                                                             <  2:151898/150‑1 (MQ=255)
                                                                                tCGGCTACATGAACCCGGCCGCCGAGCAGTTGACCCACTGGCAATCCGGCCAGGCCCAGGGCCTGCCGCTGTCGGCACTGTTCAGCCTGGTGGATGAGCACGCCGAGGAGGATGGTCGCAGCTTGGTCGAGCAGGTGCTCAGCGGCAGcc                                                         >  2:749623/1‑150 (MQ=255)
                                                                                    cTACATGAACCCGGCCGCCGAGCAGTTGACCCACTGGCAAACCGGCCAGGCCCAGGGCCTGCCGCTGTCGGCACTGTTCAGCCTGGTGGATGAGCACGCCGAGGAGGATGGTCGCAGCTTGGTCGAGCAGGTGCTCAGCGGCAGCCTCaa                                                     <  1:654237/150‑1 (MQ=255)
                                                                                      aCATGAACCCGGCCGCCGAGCAGTTGACCCACTGGCAATCCGGCCAGGCCCAGGGCCTGcc                                                                                                                                            >  2:549435/1‑61 (MQ=255)
                                                                                      aCATGAACCCGGCCGCCGAGCAGTTGACCCACTGGCAATCCGGCCAGGCCCAGGGCCTGcc                                                                                                                                            <  1:549435/61‑1 (MQ=255)
                                                                                      aCATGAACCCGGCCGCCGAGCAGTTGACCCACTGGCAATCCGGCCAGGCCCAGGGCCTGCCGCTGTCGGCACTGTTCAGCCTGGTGGATGAGCACGCCGAGGAGGATGGTCGCAGCTTGGTCGAGCAGGTGCTCAGCGGCAGc                                                          >  1:650726/1‑143 (MQ=255)
                                                                                      aCATGAACCCGGCCGCCGAGCAGTTGACCCACTGGCAATCCGGCCAGGCCCAGGGCCTGCCGCTGTCGGCACTGTTCAGCCTGGTGGATGAGCACGCCGAGGAGGATGGTCGCAGCTTGGTCGAGCAGGTGCTCAGCGGCAGc                                                          <  2:650726/143‑1 (MQ=255)
                                                                                       cATGAACCCGGCCGCCGAGCAGTTGACCCACTGGCAATCCGGCCAGGCCCAGGGCCTGCCGCTGTCGGCACTGTTCAGCCTGGTGGATGAGCACGCCGAGGAGGATGGTCGCAGCTTGGTCGAGCAGGTGCTCAGCGGCAGCCTCAAggg                                                  <  1:787181/150‑1 (MQ=255)
                                                                                        aTGAACCCGGCCGCCGAGCAGTTGACCCACTGGCAATCCGGCCAGGCCCAGGGCCTGCCGCTGTCGGCACTGTTCAGCCTGGTGGATGAGCACGCCGAGGAGGATGGTCGCAGCTTgg                                                                                 >  2:479710/1‑118 (MQ=255)
                                                                                        aTGAACCCGGCCGCCGAGCAGTTGACCCACTGGCAATCCGGCCAGGCCCAGGGCCTGCCGCTGTCGGCACTGTTCAGCCTGGTGGATGAGCACGCCGAGGAGGATGGTCGCAGCTTgg                                                                                 <  1:479710/118‑1 (MQ=255)
                                                                                         tGAACCCGGCCGCCGAGCAGTTGACCCACTGGCAATCCGGCCAGGCCCAGGGCCTGCCGCTGTCGGCACTGTTCAGCCTGGTGGATGAGCACGCCGAGGAGGATGGTCGCAGCTTGGTCGAGCAGGTGCTCAGCGGCAGCCTCAAGggcg                                                >  1:1141122/1‑150 (MQ=255)
                                                                                          gAACCCGCCCGCCGCGGAGTAGCCCCACTGCCATTGTGCCCATGCCCAGGGCCTGCCGCTGTCGGCACTGTTCATCCTGGTGGCTGCGCACGCTCAGGCGGTTGCTCCCAGCTTCCTCCGGCACGTCCCCACCGGCGCCCTCAACtgcgg                                               >  2:4827/1‑150 (MQ=255)
                                                                                            aCCCGGCCGCCGAGCAGTTGACCCACTGGCAATCCGGCCAGGCCCAGGGCCTGCCGCTGTCGGCACTGTTCAGCCTGGTGGATGAGCACGCCGAGGAGGATGGTCGCAGCTTGGTCGAGCAGGTGCTCAGCGGCAGCCTCaa                                                     <  1:1188559/142‑1 (MQ=255)
                                                                                            aCCCGGCCGCCGAGCAGTTGACCCACTGGCAATCCGGCCAGGCCCAGGGCCTGCCGCTGTCGGCACTGTTCAGCCTGGTGGATGAGCACGCCGAGGAGGATGGTCGCAGCTTGGTCGAGCAGGTGCTCAGCGGCAGCCTCaa                                                     >  2:1188559/1‑142 (MQ=255)
                                                                                              ccGGCCGCCGAGCAGTTGACCCACTGGCAATCCGGCCAGGCCCAGGGCCTGCCGCTGTCGGCACTGTTCAGCCTGGTGGATGAGCACGCCGAGGAGGATGGTCGCAGCTTGGTCGAGCAGGTGCTCAGCGGCAGCCTCAAGGGCGGCTcc                                           >  2:746521/1‑150 (MQ=255)
                                                                                                ggCCGCCGAGCAGTTGACCCACTGGCAATCCGGCCAGGCCCAGGGCCTGCCGCTGTCGGCACTGTTCAGCCTGGTGGATGAGCACGCCGAGGAGGATGGt                                                                                           <  1:867537/100‑1 (MQ=255)
                                                                                                ggCCGCCGAGCAGTTGACCCACTGGCAATCCGGCCAGGCCCAGGGCCTGCCGCTGTCGGCACTGTTCAGCCTGGTGGATGAGCACGCCGAGGAGGATGGt                                                                                           >  2:867537/1‑100 (MQ=255)
                                                                                                ggCCGCCGAGCAGTTGACCCACTGGCAATCCGGCCAGGCCCAGGGCCTGCCGCTGTCGGCACTGTTCAGCCTGGTGGATGAGCACGCCGAGGAGGATGGTCGCAGCTTGGTCGAGCAGGTGCTCAGCGGCAGCCTCAAGGGCGGCTCCGa                                         <  2:1195157/150‑1 (MQ=255)
                                                                                                     ccgAGCAGTTGACCCACTGGCAATCCGGCCAGGCCCAGGGCCTGCCGCTGTCGGCACTGTTCAGCCTGGTGGATGAGCACGCCGAGGAGGATGGTCGCAGCTTgg                                                                                 >  2:636187/1‑105 (MQ=255)
                                                                                                     ccgAGCAGTTGACCCACTGGCAATCCGGCCAGGCCCAGGGCCTGCCGCTGTCGGCACTGTTCAGCCTGGTGGATGAGCACGCCGAGGAGGATGGTCGCAGCTTgg                                                                                 <  1:636187/105‑1 (MQ=255)
                                                                                                     ccgAGCAGTTGACCCACTGGCAATCCGGCCAGGCCCAGGGCCTGCCGCTGTCGGCACTGTTCAGCCTGGTGGATGAGCACGCCGAGGAGGATGGTCGCAGCTTGGTCGAGCAGGTGCTCAGCGGCAGCCTCAAGGGCGGCTCCGAGCATg                                    >  1:1208921/1‑150 (MQ=255)
                                                                                                      cgAGCAGTTGACCCACTGGCAATCCGGCCAGGCCCAGGGCCTGCCGCTGTCGGCACTGTTCAGCCTGGTGGATGAGCACGCCGAGGAGGATGGTCGCAGCTTGGTCGAGCAGGTGCTCAGCGGCAGc                                                          <  2:90337/127‑1 (MQ=255)
                                                                                                      cgAGCAGTTGACCCACTGGCAATCCGGCCAGGCCCAGGGCCTGCCGCTGTCGGCACTGTTCAGCCTGGTGGAGGAGCACGCCGAGGAGGATGGTCGCAGCTTGGTCGAGCAGGTGCTCAGCGGCAGc                                                          >  1:90337/1‑127 (MQ=255)
                                                                                                       gAGCAGTTGACCCACTGGCAATCCCGCCAGGCCCAGGGCCTGCCGCTGTCGGCACTGTTCAGCCTGGTGGATGAGCACGCCGAGGAGGATGGTCGCAGCTTGGTCGAGCAGGTGCTCAGCGGCAGCCTCAAGGGCGGCTCCGAGCATGcc                                  <  1:105022/150‑1 (MQ=255)
                                                                                                           agtggccccACTGGCAATCCGGGCAGGCCCGGGGCGCGCGGCTGTCGGCACTGTTCAGCCTGGTGGATGAGCACGCCGAGGAGGATGGTCGCAGCTTGGTCGAGCAGGTGCTCAGCGGCAGCCTCAAGGGCGGCTCCGAGCATGCCCGGc                              <  1:391257/144‑1 (MQ=255)
                                                                                                            gTTGACCCACTGGCAATCCGGCCAGGCCCAGGGCCTGCCGCTGTCGGCACTGTTCAGCCTGGTGGATGAGCACGCCGAggagga                                                                                               >  1:10151/1‑84 (MQ=255)
                                                                                                            gTTGACCCACTGGCAATCCGGCCAGGCCCAGGGCCTGCCGCTGTCGGCACTGTTCAGCCTGGTGGATGAGCACGCCGAggagga                                                                                               <  2:10151/84‑1 (MQ=255)
                                                                                                            gTTGACCCACTGGCAATCCGGCCAGGCCCAGGGCCTGCCGCTGTCGGCACTGTTCAGCCTGGTGGATGAGCACGCCGAGGAGGATGGTCGCAGCTTGGTc                                                                               >  1:868828/1‑100 (MQ=255)
                                                                                                            gTTGACCCACTGGCAATCCGGCCAGGCCCAGGGCCTGCCGCTGTCGGCACTGTTCAGCCTGGTGGATGAGCACGCCGAGGAGGATGGTCGCAGCTTGGTc                                                                               <  2:868828/100‑1 (MQ=255)
                                                                                                             ttGACCCACTGGCAATCCGGCCAGGCCCAGGGCCTGCCGCTGTCGGCACTGTTCAGCCTGGTGGATGAGCACGCCGAGGAGGATGGTCGCAGCTTGGTCGAGCAGGTGCTCAGCGGCAGCCTCAAGGGCGGCTCCGAGCATGCCCGGCTg                            >  1:559379/1‑150 (MQ=255)
                                                                                                                aCCCACTGGCAATCCGGCCAGGCCCAGGGCCTGCCGCTGTCGGCACTGTTCAGCCTGGTGGATGAGCACGCCGAGGAGGATGGTCGCAGCTTGGTCGAGCAGGTGCTCAGCGGCAGCCTCAAGGGCGGCTCCGAGCATg                                    >  2:506457/1‑139 (MQ=255)
                                                                                                                aCCCACTGGCAATCCGGCCAGGCCCAGGGCCTGCCGCTGTCGGCACTGTTCAGCCTGGTGGATGAGCACGCCGAGGAGGATGGTCGCAGCTTGGTCGAGCAGGTGCTCAGCGGCAGCCTCAAGGGCGGCTCCGAGCATg                                    <  1:506457/139‑1 (MQ=255)
                                                                                                                 cccACTGGCAATCCGGCCAGGCCCAGGGCCTGCCGCTGTCGGCACTGTTCAGCCTGGTGGATGAGCACGCCGAGGAGGATGGTCGCAGCTTGGTCTAGCAGGTGCTCAGCGGCAGCCTCAAGGGCGGCTCCGAGCATGCCCGGCTGATcc                        >  2:439385/1‑150 (MQ=255)
                                                                                                                  ccACTGGCAATCCGGCCAGGCCCAGGGCCTGCCGCTGTCGGCACTGTTCAGCCTGGTGGATGAGCACGCCGAGGAGGATGGTCGCAGCTTGGTCGAGCAGGTGCTCAGCGGCAGCCTCAAGGGCGGCTCCGAGCATGCCCGGCTGATCCa                       >  1:785711/1‑150 (MQ=255)
                                                                                                                    aCTGGCAATCCGGCCAGGCCCAGGGCCTGCCGCTGTCGGCACTGTTCAGCCTGGTGGATGAGCACGCCGAGGAGGATGGTCGCAGCTTGGTCGAGCAGGTGCTCAGCGGCAGCCTCAAGGGCGGCTCCGAGCATGCCCGGCTGATCCAgc                     <  1:890941/150‑1 (MQ=255)
                                                                                                                    aCTGGCAATCCGGCCAGGCCCAGGGCCTGCCGCTGTCGGCACTGTTCAGCCTGGTGGATGAGCACGCCGAGGAGGATGGTCGCAGCTTGGTCGAGCAGGTGCTCAGCGGCAGCCTCAAGGGCGGCTCCGAGCATGCCCGGCTGATCCAgc                     >  1:1139414/1‑150 (MQ=255)
                                                                                                                           aTCCGGCCAGGCCCAGGGCCTGCCGCTGTCGGCACTGTTCAGCCTGGTGGATGAGCACGCCGAGGAGGATGGTCGCAGCTTGGTCGAGCAGGTGCTCAGCGGCAGCCTCAAGGGCGGCTCCGAGCATGCCCGGCTGATCCAGCGCCTGGa              >  1:773314/1‑150 (MQ=255)
                                                                                                                              ggtccAGGCCCAGGGCCTGCCGCTGTCGGCACTGTTCAGCCTGGTGGATGAGCACGCCGAGGAGGATGGTCGCAGCTTGGTCGAGCAGGTGCTCAGCGGCAGCCTCAAGGGCGGCTCCGAGCATGCCCGGCTGATCCAGCGCCTGGACgg           >  2:1180533/4‑150 (MQ=255)
                                                                                                                              cGGCCAGGCCCAGGGCCTGCCGCTGTCGGCACTGTTCAGCCTGGTGGATGAGCACGCCGAGGAGGATGGTCGCAGCTTGGTCGAGCAGGTGCTCAGCGGCAGCCTCAAGGGCGGCTCCGAGCATGCCCGGCTGATCCAGCGCCTGGACgg           <  1:746521/150‑1 (MQ=255)
                                                                                                                                ggcAGGCCCTGGGCCTGCCGCTGTCGGCACTGTTCAGCCTGGTGGATGAGCACGCCGAGGAGGATGGTCGCAGCTTGGTCGAGCAGGTGCTCAGCGGCAGCCTCAAGGGCGGCTCCGAGCATGCCCg                                >  1:75999/3‑127 (MQ=255)
                                                                                                                                ggcAGGCCCTGGGCCTGCCGCTGTCGGCACTGTTCAGCCTGGTGGATGAGCACGCCGAGGAGGATGGTCGCAGCTTGGTCGAGCAGGTGCTCAGCGGCAGCCTCAAGGGCGGCTCCGAGCATGCCCg                                <  2:75999/125‑1 (MQ=255)
                                                                                                                                gCCAGGCCCAGGGCCTGCCGCTGTCGGCACTGTTCAGCCTGGTGGATGAGCACGCCGAGGAGGATGGTCGCAGCTTGGTCGAGCAGGTGCTCAGCGGCAGCCTCAAGGGCGGCTCCGAGCATGCCCGGCTGATCCAGCGCCTGGACGGCa         <  1:895329/150‑1 (MQ=255)
                                                                                                                                 gcAGGCCCTGGGCCTGCCGCTGTCGGCACTGTTCAGCCTGGTGGATGAGCACGCCGAGGAGGATGGTCGCAGCTTGGTCGAGCAGGTGCTCAGCGGCAGCCTCAAGGGCGGCTCCGAGCATGCCCGGCTGATCCAGCGCCTGGACGGCAg        >  2:833252/2‑150 (MQ=255)
                                                                                                                                       ccAGGGCCTGCCGCTGTCGGCACTGTTCAGCCTGGTGGATGAGCACGCCGAGGAGGATGGTCGCAGCTTGGTCGAGCAGGTGCTCAGCGGCAGCCTCACGGGCGGCTCCGAGCATGCCCGGCTGATCCAGCGCCTGGACGGCAGTAcccc  >  2:476802/1‑148 (MQ=255)
                                                                                                                                       ccAGGGCCTGCCGCTGTCGGCACTGTTCAGCCTGGTGGATGAGCACGCCGAGGAGGATGGTCGCAGCTTGGTCGAGCAGGTGCTCAGCGGCAGCCTCAAGGGCGGCTCCGAGCATGCCCGGCTGATCCAGCGCCTGGACGGCAGTACCgt  >  1:819773/1‑150 (MQ=255)
                                                                                                                                        cAGGGCCTGCCGCTGTCGGCACTGTTCAGCCTGGTGGATGAGCACGCCGAGGAGGATGGTCGCAGCTTGGTCGAGCAGGTGCTCAGCGGCAGCCTCAAGGGCGGCTCCGAGCATGc                                   >  2:741511/1‑116 (MQ=255)
                                                                                                                                        cAGGGCCTGCCGCTGTCGGCACTGTTCAGCCTGGTGGATGAGCACGCCGAGGAGGATGGTCGCAGCTTGGTCGAGCAGGTGCTCAGCGGCAGCCTCAAGGGCGGCTCCGAGCATGc                                   <  1:741511/116‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                                   
GCAGGAAGAAAAAGAGCGTGCGCAGATTACCCTGCAGGCCATCGGGGATGCGGTGATTACAACCGATGTCGAGGGCAATATCGGCTACATGAACCCGGCCGCCGAGCAGTTGACCCACTGGCAATCCGGCCAGGCCCAGGGCCTGCCGCTGTCGGCACTGTTCAGCCTGGTGGATGAGCACGCCGAGGAGGATGGTCGCAGCTTGGTCGAGCAGGTGCTCAGCGGCAGCCTCAAGGGCGGCTCCGAGCATGCCCGGCTGATCCAGCGCCTGGACGGCAGTACCGT  >  NC_002947/1964620‑1964904

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: