Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,143,783 T→C L173P (CTG→CCG)  yegE → putative diguanylate cyclase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,143,7830TC97.3% 274.5 / ‑2.6 75L173P (CTG→CCG) yegEputative diguanylate cyclase
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base C (43/30);  total (44/31)
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.30e-01

GCACTCGGCAGCGCCATTGTTCCGCCTCTGTTGGGGGGTGTTCTGGTTGTCCTGCTGACGCCCGGAGACGATCCTCTCAGGGCATTTTTGATATGGGTACTGTCAGAATCCATCGGCGCTCTGGCACTGGTGCCGCTGGGATTGTTATTTAAACCACACTATCTGCTGCGCCATCGCAACCCACGGTTGCTTTTTGAGTCGCTGCTCACGTTAGCCATCACACTGACGTTAAGCTGGCTTTCGATGCTGTATCTGCCGTGGCCTTTTACTTT  >  NC_000913/2143647‑2143918
                                                                                                                                        |                                                                                                                                       
gCACTCGGCAGCGCCATTGTTCCGCCTCTGTTGGGGGGTGTTCTGGTTGTCCTGCTGACGCCCGGAGACGATCCTCTCAGGGCATTTTTGATATGGGTACTGTCAGAATCCATCGGCGCTCTGGCACTGGTGCCGCCgg                                                                                                                                       >  2:949032/1‑139 (MQ=255)
gCACTCGGCAGCGCCATTGTTCCGCCTCTGTTGGGGGGTGTTCTGGTTGTCCTGCTGACGCCCGGAGACGATCCTCTCAGGGCATTTTTGATATGGGTACTGTCAGAATCCATCGGCGCTCTGGCACTGGTGCCGCCgg                                                                                                                                       >  1:1495511/1‑139 (MQ=255)
      ggCAGCGCCATTGTTCCGCCTCTGTTGGGGGGTGTTCTGGTTGTCCTGCTGACGCCCGGAGACGATCCTCTCAGGGCATTTTTGATATGGGTACTGTCAGAATCCATCGGCGCTCTGGCACTGGTGCCGCCGGGATTg                                                                                                                                  >  2:1611423/1‑138 (MQ=255)
      ggCAGCGCCATTGTTCCGCCTCTGTTGGGGGGTGTTCTGGTTGTCCTGCTGACGCCCGGAGACGATCCTCTCAGGGCATTTTTGATATGGGTACTGTCAGAATCCATCGGCGCTCTGGCACTGGTGCCGCCGGGATTg                                                                                                                                  <  1:1611423/138‑1 (MQ=255)
                ttGTTCCGCCTCTGTTGGGGGGTGTTCTGGTTGTCCTGCTGACGCCCGGAGACGATCCTCTCAGGGCATTTTTGATATGGGTACTGTCAGAATCCATCGGCGCTCTGGCACTGGTGCCGCCGGGATTg                                                                                                                                  <  1:794351/128‑1 (MQ=255)
                ttGTTCCGCCTCTGTTGGGGGGTGTTCTGGTTGTCCTGCTGACGCCCGGAGACGATCCTCTCAGGGCATTTTTGATATGGGTACTGTCAGAATCCATCGGCGCTCTGGCACTGGTGCCGCCGGGATTg                                                                                                                                  >  2:794351/1‑128 (MQ=255)
                  gTTCCGCCTCTGTTGGGGGGTGTTCTGGTTGTCCTGCTGACGCCCGGAGACGATCCTCTCAGGGCATTTTTGATATGGGTACTGTCAGAATCCATCGGCGCTCTGGCACTGGTGCCGCCGGGATTGTTATTTAAACcac                                                                                                                     <  1:717667/139‑1 (MQ=255)
                             gTTGGGGGGTGTTCTGGTTGTCCTGCTGACGCCCGGAGACGATCCTCTCAGGGCATTTTTGATATGGGTACTGTCAGAATCCATCGGCGCTCTGGCACTGGTGCCGCCGGGATTGTTATTTa                                                                                                                           <  2:996536/122‑1 (MQ=255)
                             gTTGGGGGGTGTTCTGGTTGTCCTGCTGACGCCCGGAGACGATCCTCTCAGGGCATTTTTGATATGGGTACTGTCAGAATCCATCGGCGCTCTGGCACTGGTGCCGCCGGGATTGTTATTTa                                                                                                                           >  1:996536/1‑122 (MQ=255)
                                      tgtTCTGGTTGTCCTGCTGACGCCCGGAGACGATCCTCTCAGGGCATTTTTGATATGGGTACTGTCAGAATCCATCGGCGCTCTGGCACTGGTGCCGCCGGGTTTGTTATTTAAACCACACTATCTGCTGCGCCATCg                                                                                                  >  1:268843/1‑138 (MQ=255)
                                      tgtTCTGGTTGTCCTGCTGACGCCCGGAGACGATCCTCTCAGGGCATTTTTGATATGGGTACTGTCAGAATCCATCGGCGCTCTGGCACTGGTGCCGCCGGGTTTGTTATTTAAACCACACTATCTGCTGCGCCATCg                                                                                                  <  2:268843/138‑1 (MQ=255)
                                        ttCTGGTTGTCCTGCTGACGCCCGGAGACGATCCTCTCAGGGCATTTTTGATATGGGTACTGTCAGAATCCATCGGCGCTCTGGCACTGGTGCCGCCGGGATTGTTATTTAAACCACACTATCTGCTGCGCCATCGCaa                                                                                               >  1:608835/1‑139 (MQ=255)
                                            ggTTGTCCTGCTGACGCCCGGAGACGATCCTCTCAGGGCATTTTTGATATGGGTACTGTCAGAATCCATCGGCGCTCTGGCACTGGTGCCGCCGGGATTGTTATTTAAACCACACTATCTGCTGCGCCATCGCAACCCa                                                                                           <  2:362853/139‑1 (MQ=255)
                                                gTCCTGCTGACGCCCGGAGACGATCCTCTCAGGGCATTTTTGATATGGGTACTGTCAGAATCCATCGGCGCTCTGGCACTGGTGCCGCCGGGATTGTTATTTAAACCACACTATCTGCTGCGCCATCGCAACCCACGGt                                                                                       >  2:1560623/1‑139 (MQ=255)
                                                 tCCTGCTGACGCCCGGAGACGATCCTCTCAGGGCATTTTTGATATGGGTACTGTCAGAATCCATCGGCGCTCTGGCACTGGTGCCGCCGGGATTGTTATTTAAACCACACTATCTGCTGCGCCATCGCAACCCACGGtt                                                                                      >  2:1439974/1‑139 (MQ=255)
                                                   ctgctgACGCCCGGAGACGATCCTCTCAGGGCATTTTTGATATGGGTACTGTCAGAATCCATCGGCGCTCTGGCACTGGTGCCGCCGGGATTGTTATTTAAACCACACTATCTGCTGCGCCATCGCAACCCACGGTTGc                                                                                    >  2:451800/1‑139 (MQ=255)
                                                    tgctgACGCCCGGAGACGATCCTCTCAGGGCATTTTTGATATGGGTACTGTCCGAATCCATCGGCGCTCTGGCACTGGTGCCGCCGGGATTGTTATTTAAACCACACTACCTGCTGCGCCATCGCAACCCACGGTTGCt                                                                                   >  2:1115204/1‑139 (MQ=255)
                                                           gCCCGGAGACGATCCTCTCAGGGCATTTTTGATATGGGTACTGTCAGAATCCATCGGCGCTCTGGCACTGGTGCCGCCGGGATTGTTATTTAAACCACACTATCTGCTGCGCCATCGCAACCCACGGTTGCTTTTTGAg                                                                            >  2:1576766/1‑139 (MQ=255)
                                                             ccGGAGACGATCCTCTCAGGGCATTTTTGATATGGGTACTGTCAGAATCCATCGGCGCTCTGGCACTGGTGCCGCCGGGATTGTTATTTAAACCACACTATCTGCTGCGCCATCGCAACCCACGGTTGCTTTTTGAGTc                                                                          <  1:494621/139‑1 (MQ=255)
                                                              cGGAGACGATCCTCTCAGGGCATTTTTGATATGGGTACTGTCAGAATCCATCGGCGCTCTGGCACTGGTGCCGCCGGGATTGTTATTTAAACCACACTATCTGCTGCGCCATCGCAACCCACGGTTGCTTTTTGAGTCg                                                                         <  1:1576766/139‑1 (MQ=255)
                                                                gagaCGATCCTCTCAGGGCATTTTTGATATGGGTACTGTCAGAATCCATCGGCGCTCTGGCACTGGTGCCGCCGGGATTGTTATTTAAACCACACTATCTGCTGCGCCATCGCAACCCACGGTTGCTTTTTGAGTCgct                                                                       >  1:333526/1‑139 (MQ=255)
                                                                 agaCGATCCTCTCAGGGCATTTTTGATATGGGTACTGTCAGAATCCATCGGCGCTCTGGCACTGGTGCCGCCGGGATTGTTATTTAAACCACACTATCTGCTGCGCCATCGCAACCCACGGTTGCTTTTTGAGTCgctg                                                                      <  1:1022475/139‑1 (MQ=255)
                                                                    cGATCCTCTCAGGGCATTTTTGATATGGGTACTGTCAGAATCCATCGGCGCTCTGGCACTGGTGCCGCCGGGATTGTTATTTAAACCACACTATCTGCTGCGCCATCGCAACCCACGGTTGCTTTTTGAGTCGCTGCTc                                                                   <  2:608835/139‑1 (MQ=255)
                                                                           ctcAGGGCATTTTTGATATGGGTACTGTCAGAATCCATCGGCGCTCTGGCACTGGTGCCGCCGGGATTGTTAt                                                                                                                              >  2:783605/1‑73 (MQ=255)
                                                                           ctcAGGGCATTTTTGATATGGGTACTGTCAGAATCCATCGGCGCTCTGGCACTGGTGCCGCCGGGATTGTTAt                                                                                                                              <  1:783605/73‑1 (MQ=255)
                                                                            tcAGGGCATTTTTGATATGGGTACTGTCAGAATCCATCGGCGCTCTGGCACTGGTGCCGCCGGGATTGTTATTTAAACCACACTATCTGCTGCGCCATCGCAACCCACGGTTGCTTTTTGAGTCGCTGCTCACGTTAGc                                                           <  2:1495511/139‑1 (MQ=255)
                                                                               gggCATTTTTGATATGGGTACTGTCAGAATCCATCGGCGCTCTGGCACTGGTGCCGCCGGGATTGTTATTTAAACCACACTATCTGCTGCGCCATCGCAACCCACGGTTGCTTTTTGAGTCGCTGCTCACGTTAGCCAt                                                        >  2:197678/1‑139 (MQ=255)
                                                                                 gCATTTTTGATATGGGTACTGTCAGAATCCATCGGCGCTCTGGCACTGGTGCCGCCGGGATTGTTATTTAAACCACACTATCTGCTGCGCCATCGCAACCCACGGTTGCTTTTTGAGTCGCTGCTCACGTTAGCCATca                                                      <  1:949032/139‑1 (MQ=255)
                                                                                     ttttGATATGGGTACTGTCAGAATCCATCGGCGCTCTGGCACTGGTGCCGCCGGGATTGTTATTTAAACCACACTATCTGCTGCGCCATCGCAACCCACGGTTGCTTTTTGAGTCGCTGCTCACGTTAGCCATCACACt                                                  >  1:53169/1‑139 (MQ=255)
                                                                                     ttttGATATGGGTACTGTCAGAATCCATCGGCGCTCTGGCACTGGTGCCGCCGGGATTGTTATTTAAACCACACTATCTGCTGCGCCATCGCAACCCACGGTTGCTTTTTGAGTCGCTGCTCACGTTAGCCATCACACt                                                  >  1:980294/1‑139 (MQ=255)
                                                                                              gggTACTGTCAGAATCCATCGGCGCTCTGGCACTGGTGCCGCCGGGATTGTTATTTAAACCACACTATCTGCTGCGCCATCGCAACCCACGGTTGCTTTTTGAGTCGCTGCTCACGTTAGCCATCACACTGACGTTAAg                                         >  2:7818/1‑139 (MQ=255)
                                                                                               ggTACTGTTAGAATCCATCGGCGCTCTGGCACTGGTGCCGCCGGGATTGTTATTTAAACCACACTATCTGCTGCGCGATCGCAACCCACGGTTGCTTTTTGAGTCGCTGCTCACGTTAGCCATCACACTGACGTTAAGc                                        <  2:980294/139‑1 (MQ=255)
                                                                                               ggTACTGTCAGAATCCATCGGCGCTCTGGCACTGGTGCCGCCGGGATTGTTATTTAAACCACACTATCTGCTGCGCCATCGCAACCCACGGTTGCTTTTTGAGTCGCTGCTCACGTTAGCCATCACACTGACGTTAAGc                                        >  2:1376117/1‑139 (MQ=255)
                                                                                               ggTACTGTCAGAATCCATCGGCGCTCTGGCACTGGTGCCGCCGGGATTGTTATTTAAACCACACTATCTGCTGCGCCATCGCAACCCACGGTTGCTTTTTGAGTCGCTGCTCACGTTAGCCATCACACTGACGTTAAGc                                        >  2:330020/1‑139 (MQ=255)
                                                                                                    tGTCAGAATCCATCGGCGCTCTGGCACTGGTGCCGCCGGGATTGTTATTTAAACCACACTATCTGCTGCGCCATCGCAACCCACGGTTGCTTTTTGAGTCGCTGCTCACGtt                                                              <  1:456306/112‑1 (MQ=255)
                                                                                                    tGTCAGAATCCATCGGCGCTCTGGCACTGGTGCCGCCGGGATTGTTATTTAAACCACACTATCTGCTGCGCCATCGCAACCCACGGTTGCTTTTTGAGTCGCTGCTCACGtt                                                              >  2:456306/1‑112 (MQ=255)
                                                                                                     gTCAGAATCCATCGGCGCTCTGGCACTGGTGCCGCCGGGATTGTTCTTTAAAcc                                                                                                                       <  2:1070593/54‑1 (MQ=255)
                                                                                                     gTCAGAATCCATCGGCGCTCTGGCACTGGTGCCGCCGGGATTGTTATTTAAAcc                                                                                                                       >  1:1070593/1‑54 (MQ=255)
                                                                                                      tCAGAATCCATCGGCGCTCTGGCACTGGTGCCGCCGGGATTGTTATTTAAACCACACTATCTGCTGCGCCATCGCAACCCACGGTTGCTTTTTGAGt                                                                           <  2:19630/97‑1 (MQ=255)
                                                                                                      tCAGAATCCATCGGCGCTCTGGCACTGGTGCCGCCGGGATTGTTATTTAAACCACACTATCTGCTGCGCCATCGCAACCCACGGTTGCTTTTTGAGt                                                                           >  1:19630/1‑97 (MQ=255)
                                                                                                         gAATCCATCGGCGCTCTGGCACTGGTGCCGCCGGGATTGTTATTTAAACCACACTATCTGCTGCGCCATCGCAACCCACGGTTGCTTTTTGAGTCGCTGCTCACGTTAGCCATCACACTGAc                                               <  2:228993/122‑1 (MQ=255)
                                                                                                         gAATCCATCGGCGCTCTGGCACTGGTGCCGCCGGGATTGTTATTTAAACCACACTATCTGCTGCGCCATCGCAACCCACGGTTGCTTTTTGAGTCGCTGCTCACGTTAGCCATCACACTGAc                                               >  1:228993/1‑122 (MQ=255)
                                                                                                         gAATCCATCGGCGCTCTGGCACTGGTGCCGCCGGGATTGTTATTTAAACCACACTATCTGCTGCGCCATCGCAACCCACGGTTGCTTTTTGAGTCGCTGCTCACGTTAGCCATCACACTGACGTTAAGCTGGCTTTCg                               >  1:487228/1‑138 (MQ=255)
                                                                                                         gAATCCATCGGCGCTCTGGCACTGGTGCCGCCGGGATTGTTATTTAAACCACACTATCTGCTGCGCCATCGCAACCCACGGTTGCTTTTTGAGTCGCTGCTCACGTTAGCCATCACACTGACGTTAAGCTGGCTTTCg                               <  2:487228/138‑1 (MQ=255)
                                                                                                         gAATCCATCGGCGCTCTGGCACTGGTGCCGCCGGGATTGTTATTTAAACCACACTATCTGCTGCGCCATCGCAACCCACGGTTGCTTTTTGAGTCGCTGCTCACGTTAGCCATCACACTGACGTTAAGCTGGCTTTCGa                              >  1:545716/1‑139 (MQ=255)
                                                                                                          aaTCCATCGGCGCTCTGGCACTGGTGCCGCCGGGATTGTTATTTAAACCACACTATCTGCTGCGCCATCGCAACCCACGGTTGCTTTTTGAGTCgctgc                                                                     >  1:304329/1‑99 (MQ=255)
                                                                                                          aaTCCATCGGCGCTCTGGCACTGGTGCCGCCGGGATTGTTATTTAAACCACACTATCTGCTGCGCCATCGCAACCCACGGTTGCTTTTTGAGTCgctgc                                                                     <  2:304329/99‑1 (MQ=255)
                                                                                                          aaTCCATCGGCGCTCTGGCACTGGTGCCGCCGGGATTGTTATTTAAACCACACTATCTGCTGCGCCATCGCAACCCACGGTTGCTTTTTGAGTCGCTGCTCACGTTAGCCATCACACTGACGTTAAGCTGGCTTTCGAt                             >  1:399340/1‑139 (MQ=255)
                                                                                                          aaTCCATCGGCGCTCTGGCACTGGTGCCGCCGGGATTGTTATTTAAACCACACTATCTGCTGCGCCATCGCAACCCACGGTTGCTTTTTGAGTCGCTGCTCACGTTAGCCATCACACTGACGTTAAGCTGGCTTTCGAt                             >  2:1213148/1‑139 (MQ=255)
                                                                                                           aTCCATCGGCGCTCTGGCACTGGTGCCGCCGGGATTGTTATTTAAACCACACTATCTGCTGCGCCATCGCAACCCACGGTTGCTTTTTGAGTCGCTGCTCACGTTAGCCATCACACTGACGTTAAGCTGGCTTTCGATg                            >  1:1345067/1‑139 (MQ=255)
                                                                                                           aTCCATCGGCGCTCTGGCACTGGTGCCGCCGGGATTGTTATTTAAACCACACTATCTGCTGCGCCATCGCAACCCACGGTTGCTTTTTGAGTCGCTGCTCACGTTAGCCATCACACTGACGTTAAGCTGGCTTTCGATg                            >  2:757824/1‑139 (MQ=255)
                                                                                                            tCCATCGGCGCTCTGGCACTGGTGCCGCCGGGATTGTTATTTAAACCACACTATCTGCTGCGCCATCGCAACCCACGGTTGCTTTTTGAGTCGCTGCTCACGTTAGCCATCACACTGACGTTAAGCTGGCTTTCGATGc                           <  1:197678/139‑1 (MQ=255)
                                                                                                             ccATCGGCGCTCTGGCACTGGTGCCGCCGGGATTGTTATTTAAACCACACTATCTGCTGCGCCATCGCAACCCACGGtt                                                                                      >  2:1450592/1‑79 (MQ=255)
                                                                                                             ccATCGGCGCTCTGGCACTGGTGCCGCCGGGATTGTTATTTAAACCACACTATCTGCTGCGCCATCGCAACCCACGGtt                                                                                      <  1:1450592/79‑1 (MQ=255)
                                                                                                              cATCGGCGCTCTGGCACTGGTGCCGCCGGGATTGTTATTTAAACCACACTATCTGCTGCGCCATCGCAACCCACGGTTGCTTTTTGAGTCGCTGCTCACGTTAGCCATCACACTGACGTTAAGCTGGCTTTCGATGCTg                         >  1:1367133/1‑139 (MQ=255)
                                                                                                               aTCGGCGCTCTGGCACTGGTGCCGCCGGGATTGTTATTTAAACCACACTATCTGCTGCGCCATCGCAACCCACGGTTGCTTTTTGAGTCGCTGCTCACGTTAGCCATCACACTGACGTTAAGCTGGCTTTCGATGCTGt                        >  2:260384/1‑139 (MQ=255)
                                                                                                                tCGGCGCTCTGGCACTGGTGCCGCTGGGATTGTTATTTAAACCACACTATCTGCTGCGCCATCGCAACCCACGGTTGCTTTTTGAGTCGCTGCTCACGTTAGCCATCACACTGACg                                              <  1:761875/116‑1 (MQ=255)
                                                                                                                tCGGCGCTCTGGCACTGGTGCCGCTGGGATTGTTATTTAAACCACACTATCTGCTGCGCCATCGCAACCCACGGTTGCTTTTTGAGTCGCTGCTCACGTTAGCCATCACACTGACg                                              >  2:761875/1‑116 (MQ=255)
                                                                                                                tCGGCGCTCTGGCACTGGTGCCGCCGGGATTGTTATTTAAAc                                                                                                                        <  1:570306/42‑1 (MQ=255)
                                                                                                                tCGGCGCTCTGGCACTGGTGCCGCCGGGATTGTTATTTAAAc                                                                                                                        >  2:570306/1‑42 (MQ=255)
                                                                                                                tCGGCGCTCTGGCACTGGTGCCGCCGGGATTGTTATTTAAACCACACTATCTGCTGCGCCATCGCAACCCACGGTTGCTTTTTGAGTCGCTGCTCACGTTAGCCATCACACTGACGTTAAGCTGGCttt                                 <  1:680940/129‑1 (MQ=255)
                                                                                                                tCGGCGCTCTGGCACTGGTGCCGCCGGGATTGTTATTTAAACCACACTATCTGCTGCGCCATCGCAACCCACGGTTGCTTTTTGAGTCGCTGCTCACGTTAGCCATCACACTGACGTTAAGCTGGCttt                                 >  2:680940/1‑129 (MQ=255)
                                                                                                                tCGGCGCTCTGGCACTGGTGCCGCCGGGATTGTTATTTAAACCACACTATCTGCTGCGCCATCGCAACCCACGGTTGCTTTTTGAGTCGCTGCTCACGTTAGCCATCACACTGACGTTAAGCTGGCTTTCGATGCTGTa                       >  1:1297880/1‑139 (MQ=255)
                                                                                                                 cGGCGCTCTGGCACTGGTGCCGCCGGGATTGTTATTTAAACCACACTATCTGCTGCGCCATCGCAACCCACGGTTGCTTTTTGAGTCGCTGCTCACGTTAGCCATCACACTGACGTTAAGCTGGCTTTCGATGCTGTAt                      <  2:1345067/139‑1 (MQ=255)
                                                                                                                      ctctGGCACTGGTGCCGCCGGGATTGTTATTTAAACCACACTATCTGCTGCGCCATCGCAACCCACGGTTGCTTTTTGAGTCGCTGCTCACGTTAGCCATCACACTGACGTTAAGCTg                                      <  1:116233/118‑1 (MQ=255)
                                                                                                                      ctctGGCACTGGTGCCGCCGGGATTGTTATTTAAACCACACTATCTGCTGCGCCATCGCAACCCACGGTTGCTTTTTGAGTCGCTGCTCACGTTAGCCATCACACTGACGTTAAGCTg                                      >  2:116233/1‑118 (MQ=255)
                                                                                                                      ctctGGCACTGGTGCCGCCGGGATTGTTATTTAAACCACACTATCTGCTGCGCCATCGCAACCCACGGTTGCTTTTTGAGTCGCTGCTCACGTTAGCCATCACACTGACGTTAAGCTGGCTTTCGATGCTGTATCTGcc                 >  1:779270/1‑139 (MQ=255)
                                                                                                                        ctGGCACTGGTGCCGCCGGGATTGTTATTTAAACCACACTATCTGCTGCGCCATCGCAACCCACGGTTGCTTTTTGAGTCGCTGCTCACGTTAGCCATCACACTGACGTTAAGCTGGCTTTCGATGCTGTATCTGCCGt               >  2:1512895/1‑139 (MQ=255)
                                                                                                                         tGGCACTGGTGCCGCCGGGATTGTTATTTAAACCACACTATCTGCTGCGCCATCGCAACCCACGGTTGCTTTTTGAGTCGCTGCTCACGTTAGCCATCACACTGACGTTAAGCTGGCTTTCGATGCTGTATCTGCCGTg              <  2:779270/139‑1 (MQ=255)
                                                                                                                          ggCACTGGTGCCGCCGGGATTGTTATTTAAACCACACTATCTGCTGCGCCATCGCAACCCACGGTTGCTTTTTGAGTCGCTGCTCACGTTAGCCATCACACTGACGTTAAGCTGGCTTTCGATGCTGTATCTGCCGTgg             <  1:260384/139‑1 (MQ=255)
                                                                                                                                 gTGCCGCCGGGATTGTTATTTAAACCACACTATCTGCTGCGCCATCGCAACCCACGGTTGCTTTTTGAGTCGCTGCTCACGTTAGCCATCACACTGACg                                              >  2:825663/1‑99 (MQ=255)
                                                                                                                                 gTGCCGCCGGGATTGTTATTTAAACCACACTATCTGCTGCGCCATCGCAACCCACGGTTGCTTTTTGAGTCGCTGCTCACGTTAGCCATCACACTGACg                                              <  1:825663/99‑1 (MQ=255)
                                                                                                                                 gTGCCGCCGGGATTGTTATTTAAACCACACTATCTGCTGCGCCATCGCAACCCACGGTTGCTTTTTGAGTCGCTGCTCACGTTAGCCATCACACTGACGTTAAGCTGGCTTTCGATGCTGTATCTGCCGTGGCCTTTTa      >  2:38111/1‑139 (MQ=255)
                                                                                                                                 gTGCCGCCGGGATTGTTATTTAAACCACACTATCTGCTGCGCCATCGCAACCCACGGTTGCTTTTTGAGTCGCTGCTCACGTTAGCCATCACACTGACGTTAAGCTGGCTTTCGATGCTGTATCTGCCGTGGCCTTTTa      >  1:943721/1‑139 (MQ=255)
                                                                                                                                     cGCCGGGATTGTTATTTAAACCACACTATCTGCTGCGCCATCGCAACCCACGGTTGCTTTTTGAGTCGCTGCTCACGTTAGCCATCACACTGACGTTAAGCTGGCTTTCGATGCTGTATCTGCCGTGGCCTTTTACttt  <  1:1560623/139‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                       
GCACTCGGCAGCGCCATTGTTCCGCCTCTGTTGGGGGGTGTTCTGGTTGTCCTGCTGACGCCCGGAGACGATCCTCTCAGGGCATTTTTGATATGGGTACTGTCAGAATCCATCGGCGCTCTGGCACTGGTGCCGCTGGGATTGTTATTTAAACCACACTATCTGCTGCGCCATCGCAACCCACGGTTGCTTTTTGAGTCGCTGCTCACGTTAGCCATCACACTGACGTTAAGCTGGCTTTCGATGCTGTATCTGCCGTGGCCTTTTACTTT  >  NC_000913/2143647‑2143918

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: