Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 478,750 C→T 6.6% intergenic (‑133/+31) ylaB ← / ← ylaC putative membrane‑anchored cyclic‑di‑GMP phosphodiesterase/DUF1449 family inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913478,7500CT 6.6% 246.8 / 6.2 76intergenic (‑133/+31)ylaB/ylaCputative membrane‑anchored cyclic‑di‑GMP phosphodiesterase/DUF1449 family inner membrane protein
Reads supporting (aligned to +/- strand):  ref base C (36/35);  new base T (3/2);  total (39/37)
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.83e-01

ATGTCGTGTTCTCACCAGCATATCCTTATTAACACAGTGTAGGGAAAGTATATCTGATTGTGGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCATTTTTTGCAGTTGTTCCTTATGAACATCCGCGACACGCGGAGA  >  NC_000913/478603‑478894
                                                                                                                                                   |                                                                                                                                                
aTGTCGTGTTCTCACCAGCATATCCTTATTAACACAGTGTAGGGAAAGTATATCTGATTGTGGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCa                                                                                                                                                 <  4:14162/149‑1 (MQ=255)
      tgtTCTCACAAGCATATCCTTATTAACACAGTGTAGGGAAAGTATATCTGATTGTGGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACAAGGCGCAGAGCtt                                                                                                                                           <  4:645757/149‑1 (MQ=255)
        ttCTCACCAGCATATCCTTATTAACACAGTGTAGGGAAAGTATATCTGATTGTGGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTa                                                                                                                                         >  3:64924/1‑149 (MQ=255)
           tcACCAGCATATCCTTATTAACACAGTGTAGGGAAAGTATATCTGATTGTGGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGc                                                                                                                                      <  3:63692/149‑1 (MQ=255)
              ccAGCATATCCTTATTAACACAGTGTAGGGAAAGTATATCTGATTGTGGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTg                                                                                                                                   <  4:685886/149‑1 (MQ=255)
               cAGCATATCCTTATTAACACAGTGTAGGGAAAGTATATCTGATTGTGGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTgg                                                                                                                                  >  4:326394/1‑149 (MQ=255)
                 gCATATCCTTATTAACACAGTGTAGGGAAAGTATATCTGATTGTGGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTa                                                                                                                                <  2:350903/149‑1 (MQ=255)
                     atCCTTATTAACACAGTGTAGGGAAAGTATATCTGATTGTGGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTa                                                                                                                            >  2:31889/1‑149 (MQ=255)
                         ttattaACACAGTGTAGGGAAAGTATATCTGATTGTGGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAAccc                                                                                                                        >  3:558486/1‑149 (MQ=255)
                            ttaACACAGTGTAGGGAAAGTATATCTGATTGTGGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTaa                                                                                                                     >  2:122993/1‑149 (MQ=255)
                              tacacAGTGTAGGTAAAGTATATCTGATTGTGGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACt                                                                                                                   >  4:447790/2‑149 (MQ=255)
                               acacAGTGTAGGGAAAGTATATCTGATTGTGGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACtt                                                                                                                  <  4:634882/149‑1 (MQ=255)
                                          ggAAAGTATATCTGATTGTGGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGAt                                                                                                       >  1:553153/1‑149 (MQ=255)
                                           gAAAGTATATCTGATTGTGGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGAtt                                                                                                      >  4:636181/1‑149 (MQ=255)
                                           gAAAGTATATCTGATTGTGGAAAAAGGCTGGCTAGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGAtt                                                                                                      >  3:144803/1‑149 (MQ=255)
                                                    tCTGATTGTGGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGtt                                                                                                          <  3:639446/136‑1 (MQ=255)
                                                    tCTGATTGTGGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGtt                                                                                                          >  4:639446/1‑136 (MQ=255)
                                                    tCTGATTGTGGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGAtt                                                                                                      >  1:174390/1‑140 (MQ=255)
                                                    tCTGATTGTGGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGAtt                                                                                                      <  2:174390/140‑1 (MQ=255)
                                                    tCTGATTGTGGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCgcgag                                                                                             >  4:563761/1‑149 (MQ=255)
                                                    tCTGATTGTGGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCgcgag                                                                                             >  4:164917/1‑149 (MQ=255)
                                                      tGATTGTGGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCgcgagcg                                                                                           >  4:627539/1‑149 (MQ=255)
                                                      tGATTGTGGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCgcgagcg                                                                                           >  2:250175/1‑149 (MQ=255)
                                                            tgGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTa                                                                                     >  4:386240/1‑149 (MQ=255)
                                                            tgGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTa                                                                                     >  3:253188/1‑149 (MQ=255)
                                                               aaaaaGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTaaaa                                                                                  >  2:146943/1‑149 (MQ=255)
                                                               aaaaaGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTaaaa                                                                                  <  1:276727/149‑1 (MQ=255)
                                                                    ggctggctGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATc                                                                             <  3:470928/149‑1 (MQ=255)
                                                                           tggAGATAAGGGGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTaaaaa                                                                      <  3:447790/149‑1 (MQ=255)
                                                                             gagaTAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGa                                                                    <  1:122993/149‑1 (MQ=255)
                                                                                aTAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTa                                                                                                                 >  4:369612/1‑101 (MQ=255)
                                                                                aTAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTa                                                                                                                 <  3:369612/101‑1 (MQ=255)
                                                                                  aaGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGtt                                                               >  2:375821/1‑149 (MQ=255)
                                                                                   aggaggTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTc                                                              >  1:541200/1‑149 (MQ=255)
                                                                                        gTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCtt                                                         >  3:339401/1‑149 (MQ=255)
                                                                                           gagCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTAc                                                      <  2:433662/149‑1 (MQ=255)
                                                                                           gagCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTAc                                                      <  3:326394/149‑1 (MQ=255)
                                                                                              cTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGa                                                   <  3:180046/149‑1 (MQ=255)
                                                                                                ggTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACg                                                 >  4:354067/1‑149 (MQ=255)
                                                                                                   aTTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATc                                              <  2:295555/149‑1 (MQ=255)
                                                                                                   aTTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATc                                              <  3:164917/149‑1 (MQ=255)
                                                                                                        aTGAAGTGATGGGTGGATACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGCACGATCAtttt                                         >  4:578519/1‑149 (MQ=255)
                                                                                                           aaGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCATTTTTTg                                      >  4:156845/1‑149 (MQ=255)
                                                                                                             gTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCATTTTTTGCa                                    <  3:563761/149‑1 (MQ=255)
                                                                                                              tGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACa                                                                  >  2:427550/1‑118 (MQ=255)
                                                                                                              tGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACa                                                                  <  1:427550/118‑1 (MQ=255)
                                                                                                                aTGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCATTTTTTGCAgtt                                 <  1:375821/149‑1 (MQ=255)
                                                                                                                 tGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCATTTTTTGCAgttg                                >  3:120683/1‑149 (MQ=255)
                                                                                                                  gggTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCATTTTTTGCAgttgt                               >  3:219493/1‑149 (MQ=255)
                                                                                                                   ggTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCATTTTTTGCAgttgtt                              <  3:615004/149‑1 (MQ=255)
                                                                                                                     tGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCATTTTTTGCAGTTGTTcc                            <  2:369457/149‑1 (MQ=255)
                                                                                                                       gAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCATTTTTTGCAGGTGTTCCtt                          <  3:636181/149‑1 (MQ=255)
                                                                                                                        aaaCGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCATTTTTTGCAGTTGTTCCTTa                         <  4:3581/149‑1 (MQ=255)
                                                                                                                        aaaCGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCATTTTTTGCAGTTGTTCCTTa                         <  1:146943/149‑1 (MQ=255)
                                                                                                                        aaaCGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCATTTTTTGCAGTTGTTCCTTa                         >  1:425096/1‑149 (MQ=255)
                                                                                                                            gATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCATTTTTTGCAGTTGTTCCTTATGaa                     >  1:533642/1‑149 (MQ=255)
                                                                                                                             aTCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCATTTTTTGCAGTTGTTcc                            <  1:288553/141‑1 (MQ=255)
                                                                                                                             aTCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCATTTTTTGCAGTTGTTcc                            >  2:288553/1‑141 (MQ=255)
                                                                                                                                aTTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCATTTTTTGCAGTTGTTcc                            <  3:41539/138‑1 (MQ=255)
                                                                                                                                aTTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCATTTTTTGCAGTTGTTcc                            >  4:41539/1‑138 (MQ=255)
                                                                                                                                     cATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCATTTTTTGCAGTTGTTCCTTATGAACATCCGCGa            <  4:558486/149‑1 (MQ=255)
                                                                                                                                      aTTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCg                                                                            <  1:177231/84‑1 (MQ=255)
                                                                                                                                      aTTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCg                                                                            >  2:177231/1‑84 (MQ=255)
                                                                                                                                      aTTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCATTTTTTGCAGTTGTTCCTTATGAACATCCGCGac           >  1:503563/1‑149 (MQ=255)
                                                                                                                                       tttGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCATTTTTTGCAGTTGTTCCTTATGAACATCCGCGaca          <  1:250175/149‑1 (MQ=255)
                                                                                                                                         tGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCATTTTTTGCAGTTGTTCCTTATGAACATCCGCGACAcg        <  4:219493/149‑1 (MQ=255)
                                                                                                                                         ataCCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATc                                              >  4:84970/3‑111 (MQ=255)
                                                                                                                                         ataCCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATc                                              <  3:84970/109‑1 (MQ=255)
                                                                                                                                           aCCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCATTTTTTGCAGTTGTTCCTTATGAACATCCGCGACAcgcg      >  4:302377/1‑149 (MQ=255)
                                                                                                                                             cAGGCGTAAAGCTCTGCGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGTTCATTTTTTGCAGTTGTTCCTTATGAACATCCGCGac           <  1:95098/142‑1 (MQ=255)
                                                                                                                                             cAGGCGTAAAGCTCTGCGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCATTTTTTGCAGTTGTTCCTTATGAACATCCGCGac           >  2:95098/1‑142 (MQ=255)
                                                                                                                                             cAGGCGTAAAGCTCTGCGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCATTTTTTGCAGTTGTTCCTTATGAACATCCGCGa            >  1:607402/1‑141 (MQ=255)
                                                                                                                                             cAGGCGTAAAGCTCTGCGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCATTTTTTGCAGTTGTTCCTTATGAACATCCGCGa            <  2:607402/141‑1 (MQ=255)
                                                                                                                                              aGGCGTAAAGCTCTGCGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCATTTTTTGCAGTTGTTCCTTATGAACATCCGCGACACGCGgag   >  2:10939/1‑149 (MQ=255)
                                                                                                                                               ggCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCTTAAAAAGACAGTTCACCTTTACGGACGATCATTTTTTGCAGTTGTTCCTTATGAACATCCGCGACACGCGgaga  <  2:503563/149‑1 (MQ=255)
                                                                                                                                               ggCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCATTTTTTGCAGTTGTTCCTTATGAACATCCGCGACACGCGgaga  <  1:31889/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
ATGTCGTGTTCTCACCAGCATATCCTTATTAACACAGTGTAGGGAAAGTATATCTGATTGTGGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCATTTTTTGCAGTTGTTCCTTATGAACATCCGCGACACGCGGAGA  >  NC_000913/478603‑478894

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