Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,964,757 A→T 6.7% 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 6.7% 226.3 / 6.1 75Q309L (CAG→CTG) yegEphosphodiesterase domain‑containing diguanylate cyclase
Reads supporting (aligned to +/- strand):  ref base A (40/30);  new base T (2/3);  total (42/33)
Fisher's exact test for biased strand distribution p-value = 6.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.27e-01

CCAGCGCCCTGCAGGAAGAAAAAGAGCGTGCGCAGATTACCCTGCAGGCCATCGGGGATGCGGTGATTACAACCGATGTCGAGGGCAATATCGGCTACATGAACCCGGCCGCCGAGCAGTTGACCCACTGGCAATCCGGCCAGGCCCAGGGCCTGCCGCTGTCGGCACTGTTCAGCCTGGTGGATGAGCACGCCGAGGAGGATGGTCGCAGCTTGGTCGAGCAGGTGCTCAGCGGCAGCCTCAAGGGCGGCTCCGAGCATGCCCGGCTGATCCAGCGCCTGGACGGCAGTACCGTGT  >  NC_002947/1964610‑1964906
                                                                                                                                                   |                                                                                                                                                     
ccAGCGCCCTGCAGGAAGAAAAAGAGCGTGCGCAGATTACCCTGCAGGCCATCGGGGATGCGGTGATTACAACCGATGTCGAGGGCAATATCGGCTACATGAACCCGGCCGCCGCGCAGTTGACCCACTGGCAATCCGGCCAGGCCCAgg                                                                                                                                                     >  2:433480/1‑150 (MQ=255)
ccAGCGCCCTGCAGGAAGAAAAAGAGCGTGCGCAGATTACCCTGCAGGCCATCGGGGATGCGGTGATTACAACCGATGTCGAGGGCAATATCGGCTACATGAACCCGGCCGCCGAGCAGTTGACCCACTGGCAATCCGGCCAGGCCCTgg                                                                                                                                                     <  1:143595/150‑1 (MQ=255)
   ccgcCCTGCAGCAAGAAAAAGAGCGTGCGCAGATTACCCTGCAGGCCATCGGGGATGCGGTTATTATAACCGATGTCGAGGGCAATATCGGCTACATGAACCCGGCCGCCGAGCAGCTGATCCGCAGGCAATCCGGCCAGGCCCAGGGcc                                                                                                                                                  <  1:759101/149‑1 (MQ=255)
     gcCCTGCAGGAAGAAAAAGAGCGTGCGCAGATTACCCTGCAGGCCATCGGGGATGCGGTGATTACAACCGCTGTCTCCTTCAATATCGGCGTGATGAACCCGGCCGCCGAGCAGTTGACCCACTGGCAATCCGGCCAGGCCCAGGGCCTg                                                                                                                                                <  1:411559/150‑1 (MQ=255)
      cccTGCAGGAAGAAAAAGAGCGTGCGCAGATTACCCTGCAGGCCATCGGGGATGCGGTGATTACAACCGATGTCGAGGGCAATATCGGCTACATGAACCCGGCCGCCGAGCAGTTGACCCACTGGCAATCCGGCCAGGCCCAGGGCCTGc                                                                                                                                               <  1:228012/150‑1 (MQ=255)
       ccTGCAGGAAGAAAAAGAGCGTGCGCAGATTACCCTGCAGGCCATCGGGGATGCGGTGATTACAACCGATGTCGAGGGCAATATCGGCTACATGAACCCGGCCGCCGAGCAGTTGACCCACTGGCAATCCGGCCAGGCCCAGGGCCTGcc                                                                                                                                              >  2:232237/1‑150 (MQ=255)
            aGGAAGAAAAAGAGCGTGCGCAGATTACCCTGCAGGCCATCGGGGATGCGGTGATTACAACCGATGTCGAGGGCAATATCGGCTACATGAACCCGGCCGCCGAGCAGTTGACCCACTGGCAATCCGGCCAGGCCCAGGGCCTGCCGCTGt                                                                                                                                         >  2:472033/1‑150 (MQ=255)
                agaaAAAGAGCGTGCGCAGATTACCCTGCAGGCCATCGGGGATGCGGTGATTACAACCGATGTCGAGGGCAATATCGGCTACATGAACCCGGCCGCCGAGCAGTTGACCCACTGGCAATCCGGCCAGGCCCAGGGCCTGCCGCTGTCGGc                                                                                                                                     >  1:676055/1‑150 (MQ=255)
                         gCGTGCGCAGATTACCCTGCAGGCCATCGGGGATGCGGTGATTACAACCGATGTCGAGGGCAATATCGGCTACATGAACCCGGCCGCCGAGCAGTTGACCCACTGGCAATCCGGCCAGGCCCAGGGCCTGCCGCTGTCGGCACTGTTCAg                                                                                                                            >  1:1108386/1‑150 (MQ=255)
                               gcAGATTACCCTGCAGGCCATCGGGGCTGCGGTGATTACAACCGATGTCGAGGGCAATATCGGCTACATGAACCCGGCCGCCGAGCAGTTGACCCACTGGCAATCCGGCCAGGCCCAGGGCCTGCCGCTGTCGGCACTGTTCAGCCtggt                                                                                                                      <  2:988199/150‑1 (MQ=255)
                                   aTTACCCTGCAGGCCATCGGGGATGCGGTGATTACAACCGATGTCGAGGGCAATATCGGCTACATGAACCCGGCCGCCGAGCAGTTGACCCACTGGCAATCCGGCCAGGCCCAGGGCCTGCCGCTGTCGGCACTGTTCAGCCTGGTGGAt                                                                                                                  <  2:370591/150‑1 (MQ=255)
                                     tACCCTGCAGGCCATCGGGGATGCGGTGATTACAACCGATGTCGAGGGCAATATCGGCTACATGAACCCGGCCGCCGAGCAGTTGACCCACTGGCAATCCGGCCAGGCCCAGGGCCTGCCGCTGTCGGCACTGTTCAGCCTGGTGGATGa                                                                                                                <  2:676055/150‑1 (MQ=255)
                                      aCCCTGCAGGCCATCGGGGATGCGGTGATTACAACCGATGTCGAGGGCAATATCGGCTACATGAACCCGGCCGCCGAGCAGTTGACCCACTGGCAATCCGGCCAGGCCCAGGGCCTGCCTCTGTCGGCACTGTTCAGCCTGGTGGATGAg                                                                                                               <  2:1031024/150‑1 (MQ=255)
                                          tGCAGGCCATCGGGGATGCGGTGATTACAACCGATGTCGAGGGCAATATCGGCTACATGAACCCGGCCGCCGAGCAGTTGACCCACTGGCAATCCGGCCAGGCCCAGGGCCTGCCGCTGTCGGCACTGTTCAGCCTGGTGGATGAGCACg                                                                                                           >  1:975579/1‑150 (MQ=255)
                                              ggCCATCGGGGATGCGGTGATTACACCCGATGTCGAGGGCAATATCGGCTACATGAACCCGGCCGCCGAGCAGTTGACCCACTGGCAATCCGGCCAGGCCCAGGGCCTGCCGCTGTCGGCACTGTTCAGCCTGGTGGATGAGCACGCCga                                                                                                       <  2:1107695/150‑1 (MQ=255)
                                                 cATCGGGGATGCGGTGATTACAACCGATGTCGAGGGCAATATCGGCTACATGAACCCGGCCGCCGAGCAGTTGACCCACTGGCAATCCGGCCAGGCCCAGGGCCTGCCGCTGTCGGCACTGTTCAGCCTGGTGGATGAGCACGCCgagga                                                                                                    >  2:572179/1‑150 (MQ=255)
                                                   tCGGGGATGCGGTGATTACAACCGATGTCGAGGGCAATATCGGCTACATGAACCCGGCCGCCGAGCAGTTGACCCACTGGCAATCCGGCCAGGCCCAGGGCCTGCCGCTGTCGGCACTGTTCAGCCTGGTGGATGAGCACGCCGaggagg                                                                                                  >  1:950610/1‑150 (MQ=255)
                                                         aTGCGGTGATTACAACCGATGTCGAGGGCAATATCGGCTACATGAACCCGGCCGCCGAGCAGTTGACCCACTGGCAATCCGGCCAGGCCCAGGGCCTGCCGCTGTCGGCACTGTTCAGCCTGGTGGATGAGCACGCCGAGGAGGATGGTc                                                                                            <  2:1036223/150‑1 (MQ=255)
                                                           gCGGTGATTACAACCGATGTCGAGGGCAATATCGGCTACATGAACCCGGCCGCCGAGCAGTTGACCCACTGGCAATCCGGCCAGGCCCAGGGCCTGCCGCTGTCGGCACTGTTCAGCCTGGTGGAt                                                                                                                  >  2:793558/1‑126 (MQ=255)
                                                           gCGGTGATTACAACCGATGTCGAGGGCAATATCGGCTACATGAACCCGGCCGCCGAGCAGTTGACCCACTGGCAATCCGGCCAGGCCCAGGGCCTGCCGCTGTCGGCACTGTTCAGCCTGGTGGAt                                                                                                                  <  1:793558/126‑1 (MQ=255)
                                                            cGGTGATTACAACCGATGTCGAGGGCAATATCGGCTACATGAACCCGGCCGCCGAGCAGTTGACCCACTGGCAATCCGGCCAGGCCCAGGGCCTGCCGCTGTCGGCACTGTTCAGCCTGGTGGATGAGCACGCCGAGGAGGATGGTCGCa                                                                                         <  2:827901/150‑1 (MQ=255)
                                                                 aTTACAACCGATGTCGAGGGCAATATCGGCTACATGAACCCGGCCGCCGAGCAGTTGACCCACTGGCAATCCGGCCAGGCCCAGGGCCTGCCGCTGTCGGCACTGTTCAGCCTGGTGGATGAGCACGCCGAGGAGGATGGTCGCAGCTTg                                                                                    <  1:232237/150‑1 (MQ=255)
                                                                  ttACAACCGATGTCGAGGGCAATATCGGCTACATTAACCCGGCCGCCGAGCAGTTGACCCACTGGCAATCCGGCCAGGCCCAGGGCCTGCCGCTGTCGGCACTGTTCa                                                                                                                             <  1:625279/108‑1 (MQ=255)
                                                                  ttACAACCGATGTCGAGGGCAATATCGGCTACATGAACCCGGCCGCCGAGCAGTTGACCCACTGGCAATCCGGCCAGGCCCAGGGCCTGCCGCTGTCGGCACTGTTCa                                                                                                                             >  2:625279/1‑108 (MQ=255)
                                                                        ccGATGTCGAGGGCAATATCGGCTACATGAACCCGGCCGCCGAGCAGTTGACCCACTGGCAATCCGGCCAGGCCCAGGGCCTGCCGCTGTCGGCACTGTTCAGCCTGGTGGATGAGCACGCCGAGGAGGATGGTCGCAGCTTGGTCGAGc                                                                             >  2:389949/1‑150 (MQ=255)
                                                                          gATGTCGAGGGCAATATCGGCTACATGAACCCGGCCGCCGAGCAGTTGACCCACTGGCAATCCGGCCAGGCCCAGGGCCTGCCGCTGTCGGCACTGTTCAGCCTGGTGGATGAGCACGCCGAGGAGGATGGTCGCAGCTTGGTCGAGCAg                                                                           >  2:887548/1‑150 (MQ=255)
                                                                           aTGTCGAGGGCAATATCGGCTACATGAACCCGGCCGCCGAGCAGTTGACCCACTGGCAATCCGGCCAGGCCCAGGGCCTGCCGCTGTCGGCACTGTTCAGCCTGGTGGATGAGCACGCCGAGGAGGATGGTCGCAGCTTGGTCGAGCAgg                                                                          <  1:433480/150‑1 (MQ=255)
                                                                              tCGAGGGCAATATCGGCAACATGAACACGGCCGCCGCCCAGTTGACACACAAGCAATCCGACCAGCCCCAGGGCCTGCCCCTGTCGGCACTGTACAGCCTGCTGGATGAGCACGCCGAGGAGGAAGGACGCAGCGTGGTCGCACAGGTGc                                                                       >  2:239593/1‑150 (MQ=255)
                                                                                gAGGGCAATATCGGCTACATGAACCCGGCCGCCGAGCAGTTGACCCACTGGCAATCCGGCCAGGCCCAGGGCCTGCCGCTGTCGGCACTGTTCAGCCTGGTGGATGAGCACGCCGAGGAGGATGGTCGCAGCTTGGTCGAGCAGGTGCTc                                                                     >  2:615005/1‑150 (MQ=255)
                                                                                 aGGGCAATATCGGCTACATGAACCCGGCCGCCGAGCAGTTGACCCACTGGCAATCCGGCCAGGCCCAGGGCCTGCCGCTGTCGGCACTGTTCAGCCTGGTGGATGAGCACGCCGAGGAGGATGGTCGCAGCTTGGTCGAGCAGGTGCTCa                                                                    >  2:1070701/1‑150 (MQ=255)
                                                                                  gggCAATATCGGCTACATGAACCCGGCCGCCGAGCAGTTGACCCACTGGCAATCCGGCCAGGCCCAGGGCCTGCCGCTGTCGGCACTGTTCAGCCTGGTGGATGAGCACGCCGAGGAGGATGGTCGCAGCTTGGTCGAGCAGGTGCTCAg                                                                   >  1:567368/1‑150 (MQ=255)
                                                                                   ggCAATATCGGCTACATGAACCCGGCCGCCGAGCAGTTGACCCACTGGCAATCCGGCCAGGCCCAGGGCCTGCCGCTGTCGGCACTGTTCAGCCTGGTGGATGAGCACGCCGAGGAGGATGGTCGCAGCTTGGTCGAGCAGGTGCTCAGc                                                                  >  1:197111/1‑150 (MQ=255)
                                                                                    gCAATATCGGCTACATGAACCCGGCCGCCGAGCAGTTGACCCACTGGCAATCCGGCCAGGCCCAGGGCCTGCCGCTGTCGGCACTGTTCAGCCTGGTGGATGAGCACGCCGAGGAGGATGGTCGCAGCTTGGTCGAGCAGGTGc                                                                       >  1:887919/1‑144 (MQ=255)
                                                                                    gCAATATCGGCTACATGAACCCGGCCGCCGAGCAGTTGACCCACTGGCAATCCGGCCAGGCCCAGGGCCTGCCGCTGTCGGCACTGTTCAGCCTGGTGGATGAGCACGCCGAGGAGGATGGTCGCAGCTTGGTCGAGCAGGTGc                                                                       <  2:887919/144‑1 (MQ=255)
                                                                                     cAATATCGGCTACATGAACCCGGCCGCCGAGCAGTTGACCCACTGGCAATCCGGCCAGGCCCTGGGCCTGCCGCTGTCGGCACTGTTCAGCCTGGTGGATGAGCACGCCGAGGAGGATGGTCGCAGCTTGGTCGAGCAg                                                                           >  1:1076136/1‑139 (MQ=255)
                                                                                     cAATATCGGCTACATGAACCCGGCCGCCGAGCAGTTGACCCACTGGCAATCCGGCCAGGCCCTGGGCCTGCCGCTGTCGGCACTGTTCAGCCTGGTGGATGAGCACGCCGAGGAGGATGGTCGCAGCTTGGTCGAGCAg                                                                           <  2:1076136/139‑1 (MQ=255)
                                                                                      aaTATCGGCTACATGAACCCGGCCGCCGAGCAGTTGACCCACTGGCAATCCGGCCAGGCCCAGGGCCTGCCGCTGTCGGCACTGTTCAGCCTGGTGGATGAGCACGCCGAGGAGGATGGTCGCAGCTTGGTCGAGCAGGTGCTCAGCGGc                                                               >  2:888219/1‑150 (MQ=255)
                                                                                            ggCTACATGAACCCGGCCGCCGAGCAGTTGACCCACTGGCAATCCGGCCAGGCCCAGGGCCTGCCGCTGTCGGCACTGTTCAGCCTGGTGGATGAGCACGCCGAGGAGGATGGTCGCAGCTTGGt                                                                                  >  1:390259/1‑125 (MQ=255)
                                                                                            ggCTACATGAACCCGGCCGCCGAGCAGTTGACCCACTGGCAATCCGGCCAGGCCCAGGGCCTGCCGCTGTCGGCACTGTTCAGCCTGGTGGATGAGCACGCCGAGGAGGATGGTCGCAGCTTGGt                                                                                  <  2:390259/125‑1 (MQ=255)
                                                                                             gCTACATGAACCCGGCCGCCGAGCAGTTGACCCACTGGCAATCCGGCCAGGCCCAGGGCCTGCCGCTGTCGGCACTGTTCAGCCTGGTGGATGAGCACGCCGAGGAGGATGGTCGCAGCTTGGTCGAGCAGGTGCTCAGCGGCAGCCTCa                                                        <  2:919660/150‑1 (MQ=255)
                                                                                             gCTACATGAACCCGGCCGCCGAGCAGTTGACCCACTGGCAATCCGGCCAGGCCCAGGGCCTGCCGCTGTCGGCACTGTTCAGCCTGGTGGATGAGCACGCCGAGGAGGATGGTCGCAGCTTGGTCGAGCAGGTGCTCAGCGGCAGCCTCa                                                        >  1:1107871/1‑150 (MQ=255)
                                                                                              cTACATGAACCCGGCCGCCGAGCAGTTGACCCACTGGCAATCCGGCCAGGCCCAGGGCCTGCCGCTGTCGGCACTGTTCAGCCTGGTGGATGAGCACGCCGAGGAGGATGGTCGCAGCTTGGTCGAGCAGGTGCTCAGCGGCAGCCTCaa                                                       >  1:975264/1‑150 (MQ=255)
                                                                                              cTACATGAACCCGGCCGCCGAGCAGTTGACCCACTGGCAATCCGGCCAGGCCCAGGGCCTGCCGCTGTCGGCACTGTTCAGCCTGGTGGATGAGCACGCCGAGGAGGATGGTCGCAGCTTGGTCGAGCAGGTGCTCAGCGGCAGCCTCaa                                                       <  2:975579/150‑1 (MQ=255)
                                                                                                       cccGGCCGCCGAGCAGTTGACCCACTGGCAATCCGGCCAGGCCCAGGGCCTGCCGCTGTCGGCACTGTTCAGCCTGGTGGATGAGCACGCCGAGGAGGATGGTCGCAGCTTGGTCGAGCAGGTGCTCAGCGGCAGCCTCAAGGGCGGCTc                                              >  2:940091/1‑150 (MQ=255)
                                                                                                         cGGCCGCCGAGCAGTTGACCCACTGGCAATCCGGCCAGGCCCAGGGCCTGCCGCTGTCGGCACTGTTCAGCCTGGTGGATGAGCACGCCGAGGAGGATGGTCGCAGCTTGGTCGAGCAGGTGCTCAGCGGCAGCCTCAAGGGCGGCTCCg                                            >  1:180201/1‑150 (MQ=255)
                                                                                                         cGGCCGCCGAGCAGTTGACCCACTGGCAATCCGGCCAGGCCCAGGGCCTGCCGCTGTCGGCACTGTTCAGCCTGGTGGATGAGCACGCCGAGGAGGATGGTCGCAGCTTGGTCGAGCAGGTGCTCAGCGGCAGCCTCAAGGGCGGCTCCg                                            >  1:737082/1‑150 (MQ=255)
                                                                                                          ggCCGCCGAGCAGTTGACCCACTGGCAATCCGGCCAGGCCCAGGGCCTGCCGCTGTCGGCACTGTTCAGCCTGGTGGATGAGCACGCCGAGGAGGATGGTCGCAGCTTGGTCGAGCAGGTGCTCAGCGGCAGc                                                            <  1:597462/133‑1 (MQ=255)
                                                                                                          ggCCGCCGAGCAGTTGACCCACTGGCAATCCGGCCAGGCCCAGGGCCTGCCGCTGTCGGCACTGTTCAGCCTGGTGGATGAGCACGCCGAGGAGGATGGTCGCAGCTTGGTCGAGCAGGTGCTCAGCGGCAGCCTCAAGGGCGGCTCCGa                                           <  2:975264/150‑1 (MQ=255)
                                                                                                          ggCCGCCGAGCAGTTGACCCACTGGCAATCCGGCCAGGCCCAGGGCCTGCCGCTGGCGGCACTGTTCAGCCTGGTGGATGAGCACGCCGAGGAGGATGGTCGCAGCTTGGTCGAGCAGGTGCTCAGCGGCAGc                                                            >  2:597462/1‑133 (MQ=255)
                                                                                                              gccgAGCAGTTGACCCACTGGCAATCCGGCCAGGCCCAGGGCCTGCCGCTGTCGGCACTGTTCAGCCTGGTGGATGAGCACGCCGAGGAGGATGGTCGCAGCTTGGTCGAGCAGGTGCTCAGCGGCAGCCTCAAGGGCGGCTCCGAGCAt                                       <  1:389949/150‑1 (MQ=255)
                                                                                                              gccgAGCAGTTGACCCACTGGCAATCCGGCCAGGCCCAGGGCCTGCCGCTGTCGGCACTGTTCAGCCTGGTGGATGAGCACGCCGAGGAGGATGGTCGCAGCTTGGTCGAGCAGGTGCTCAGCGGCAGCCTCAAGGGCGGCTCCGAGCAt                                       >  1:327661/1‑150 (MQ=255)
                                                                                                              cccgAGCAGTTGACCCACTGGCAATCCGGCCAGGCCCAGGGCCTGCCGCCGTCGGCACTGTTCAGCCTGGCGGGTGAGCACGCCCAGCAGGATGGTCGCAGCTTGGTCGAGCAGGTGCTCAGCGGCAGGCTCAAGGGCGGCTCCAAGCAt                                       >  2:368642/2‑150 (MQ=255)
                                                                                                                   gCAGTTGACCCACTGGCAATCCGGCCAGGCCCAGGGCCTGCCGCTGTCGGCACTGTTCAGCCTGGTGGATGAGCACGCCGAGGAGGATGGTCGCAGCTTGGTCGAGCAGGTGCTCAGCGGCAGCCTCAAGGGCGGCTCCGAGCATGCCCg                                  <  2:950610/150‑1 (MQ=255)
                                                                                                                    cAGTTGACCCACTGGCAATCCGGCCAGGCCCAGGGCCTGCCGCTGTCGGCACTTTTCAGCCTGGTGGGTGAGCACGCCGAGGAGGATGGTCGCAGCTTGGTCGAGCAGGTGCTCAGCGGCAGCCTCAAGGGCGGCTCCGAGCATGCCCgg                                 <  2:607679/150‑1 (MQ=255)
                                                                                                                    cAGTTGACCCACTGGCAATCCGGCCAGGCCCAGGGCCTGCCGCTGTCGGCACTGTTCAGCCTGGTGGATGAGCACGCCGAGGAGGATGGTCGCAGCTTGGTCGAGCAGGTGCTCAGCGGCAGCCTCAAGGGCGGCTCCGAGCATGCCCgg                                 >  1:824168/1‑150 (MQ=255)
                                                                                                                            ccACTGGCAATCCGGCCAGGCCCAGGGCCTGCCGCTGTCGGCACTGTTCAGCCTGGTGGATGAGCACGCCGAGGAGGATGGTCGCAGCTTGGTCGAGCAGGTGCTCAGCGGCAGCCTCAAGGGCGGCTCCGAGCATGCCCGGCTGATCCa                         >  1:597821/1‑150 (MQ=255)
                                                                                                                            ccACTGGCAATCCGGCCAGGCCCAGGGCCTGCCGCTGTCGGCACTGTTCAGCCTGGTGGATGAGCACGCCGAGGAGGATGGTCGCAGCTTGGTCGAGCAGGTGCTCAGCGGCAGCCTCAAGGGCGGCTCCGAGCATGCCCGGCTGATCCa                         >  2:356305/1‑150 (MQ=255)
                                                                                                                            ccACTGGCAATCCGGCCAGGCCCAGGGCCTGCCGCTGTCGGCACTGTTCAGCCTGGTGGATGAGCACGCCGAGGAGGATGGTCGCAGCTTGGTCGAGCAGGTGCTCAGCGGCAGCCTCAAGGGCGGCTCCGAGCATGCCCGGCTGATCCa                         <  2:1108386/150‑1 (MQ=255)
                                                                                                                                tGGCAATCCGGCCAGGCCCAGGGCCTTCCGCTGGCGGCACTGTTCAGCCTGGTGGAGGAGCAGATCGGTGCTCATGGTCGCAGCTTGGTCGAGCAGGTGCTCAGCGGCAGCCTCAAGGGCGGCTCCGAGCATGCCCGGCTGATCCAgcgc                     <  1:569642/150‑1 (MQ=255)
                                                                                                                                 ggCAATCCGGCCAGGCCCAGGGCCTGCCGCTGTCGGCACTGTTCAGCCTGGTGGATGAGCACGCCGAGGAGGATGGTCGCAGCTTGGTCGAGCAGGTGCTCAGCGGCAGCCTCAAGGGCGGCTCCGAGCATGCCCGGCTGATCCAGCGcc                    >  1:85366/1‑150 (MQ=255)
                                                                                                                                   cAATCCGGCCAGGCCCAGGGCCTGCCGCTGTCGGCACTGTTCAGCCTGGTGGATGAGCACGCCGAGGAGGATGGTCGCAGCTTGGTCGAGCAGGTGCTCAGCGGCAGCCTCAAGGGCGGCTCCGAGCATGCCCGGCTGATCCAGCGCCTg                  <  2:1107871/150‑1 (MQ=255)
                                                                                                                                    aaTCCGGCCAGGCCCAGGGCCTGCCGCTGTCGGCACTGtt                                                                                                                               <  2:200488/40‑1 (MQ=255)
                                                                                                                                    aaTCCGGCCAGGCCCAGGGCCTGCCGCTGTCGGCACTGtt                                                                                                                               >  1:200488/1‑40 (MQ=255)
                                                                                                                                              ggCCCAGGGCCTGCCGCTGTCGGCACTGTTCAGCCTGGTGGATGAGCACGCCGAGGAGGATGGTCGCAGCTTGGTCGAGCAGGTGCTCAGCGGCAGCCTCAAGGGCGGCTCCGAGCATGCCCGGCTGATCCAGCGCCTGGACGGCAGTAc       <  2:567368/150‑1 (MQ=255)
                                                                                                                                               gCCCTGGGCCTGCCGCTGTCGGCACTGTTCAGCCTGGTGGa                                                                                                                   <  1:151896/41‑1 (MQ=255)
                                                                                                                                               gCCCTGGGCCTGCCGCTGTCGGCACTGTTCAGCCTGGTGGa                                                                                                                   >  2:151896/1‑41 (MQ=255)
                                                                                                                                               gCCCAGGGCCTGCCGCTGTCGGCACTGTTCAGCCTGGTGGATGAGCACGCCGAGGAGGATGGTCGCAGCTTGGTCGAGCAGGTGCTCAGCGGCAGCCTCAAGGGCGGCTCCGAGCATGCCCGGCTGATCCAGCGCCTGGACGGCAGTAcc      <  2:737082/150‑1 (MQ=255)
                                                                                                                                               gCCCAGGGCCTGCCGCTGTCGGCACTGTTCAGCCTGGTGGATGAGCACGCCGAGGAGGATGGTCGCAGCTTGGTCGAGCAGGTGCTCAGCGGCAGCCTCAAGGGCGGCTCCGAGCATGCCCGGCTGATCCAGCGCCTGGACGGCAGTAcc      <  2:180201/150‑1 (MQ=255)
                                                                                                                                               gCCCAGGGCCTGCCGCTGTCGGCACTGTTCAGCCTGGTGGATGAGCACGCCGAGGAGGATGGTCGCAGCTTGGTCGAGCAGGTGCTCAGCGGCAGCCTCAAGGGCGGCTCCGAGCATGCCCGGCTGATCCAGCGCCTGGACGGCAGTAcc      >  1:1074537/1‑150 (MQ=255)
                                                                                                                                                 ccAGGGCCTGCCGCTGTCGGCACTGTTCAGCCTGGTGGATGAGCACGCCGAGGAGGATGGTCGCAGCTTGGTCGAGCAGGTGCTCAGCGGCAGCCTCAAGGGCGGCTCCGAGCATGCCCGGCTGATCCAGCGCCTGGACGGCAGTACCgt    >  2:157617/1‑150 (MQ=255)
                                                                                                                                                 ccAGGGCCTGCCGCTGTCGGCACTGTTCAGCCTGGTGGATGAGCACGCCGAGGAGGATGGTCGCAGCTTGGTCGAGCAGGTGCTCAGCGGCAGCCTCAAGGGCGGCTCCGAGCATGCCCGGCTGATCCAGCGCCTGGACGGCAGTACCgt    >  1:712642/1‑150 (MQ=255)
                                                                                                                                                 ccAGGGCCTGCCGCTGTCGGCACTGTTCAGCCTGGTGGATGAGCACGCCGAGGAGGATGGTCGCAGCTTGGTCGAGCAGGTGCTCAGCCGCAGCCTCAAGGGCGGCTCCGAGCATGACGGAATGATCCAGCTCTGGGACGGCAGTAACgt    >  1:976150/1‑150 (MQ=255)
                                                                                                                                                  cAGGGCCTGCCGCTGTCGGCACTGTTCAGCCTGGTGGATGAGCACGCCGAGGAGGATGGTCGCAGCTTGGTCGAGCAGGTGCTCAGCGGCAGCCTCAAGGGCGGCTCCGAGCATGCCCGGCTGATCCAGCGCCTGGACGGCAGTACCgtg   >  1:1022351/1‑150 (MQ=255)
                                                                                                                                                  cAGGGCCTGCCGCTGTCGGCACTGTTCAGCCTGGTGGATGAGCACGCCGAGGAGGATGGTCGCAGCTTGGTCGAGCAGGGAGTCAGTGGCAGCCTCAAGTGCGGCTACGAGCCTGACAGCCTAACTCAGGTCCTGTTCGTCAGTACCgtg   >  1:20820/1‑150 (MQ=255)
                                                                                                                                                  cAGGGCCTGCCGCTGTCGGCACTGTTCAGCCGGGTGGATGAGCACGCCGAGGAGGATGGTCGCAGCTTGGTCGAGCAGGTGCTCAGCGGCAGCCTCAAGGGCGGCTCCAAGCATGCCCGGCTGATCCAGCGCCTGGACGGCAGTACCgtg   >  2:855913/1‑150 (MQ=255)
                                                                                                                                                   aGGGCCTGCCGCTGTCGGCACTGTTCAGCCTGGTGGATGAGCACGCCGAGGAGGATGGTCGCAGCTTGGTCGAGCAGGTGCTCAGCGGCAGCCTCAAGGGCGGCTCCGAGCATGCCCGGCTGATCCAGCGCCTGGACGGCAGTACCgtgt  <  1:334591/150‑1 (MQ=255)
                                                                                                                                                   aGGGCCTGCCGCTGTCGGCACTGTTCAGCCTGGTGGATGAGCACGCCGAGGAGGATGGTCGCAGCTTGGTCGAGCAGGTGCTCAGCGGCAGCCTCAAGGGCGGCTCCGAGCATGCCCGGCTGATCCAGCGCCTGGACGGCAGTACCgtgt  <  1:615005/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
CCAGCGCCCTGCAGGAAGAAAAAGAGCGTGCGCAGATTACCCTGCAGGCCATCGGGGATGCGGTGATTACAACCGATGTCGAGGGCAATATCGGCTACATGAACCCGGCCGCCGAGCAGTTGACCCACTGGCAATCCGGCCAGGCCCAGGGCCTGCCGCTGTCGGCACTGTTCAGCCTGGTGGATGAGCACGCCGAGGAGGATGGTCGCAGCTTGGTCGAGCAGGTGCTCAGCGGCAGCCTCAAGGGCGGCTCCGAGCATGCCCGGCTGATCCAGCGCCTGGACGGCAGTACCGTGT  >  NC_002947/1964610‑1964906

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: