Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 527,946 A→G 14.4% noncoding (1181/2902 nt) PP_23SC → 23S ribosomal RNA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947527,9460AG14.4% 143.6 / 18.8 63noncoding (1181/2902 nt)PP_23SC23S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base A (23/30);  new base G (3/6);  total (26/37)
Fisher's exact test for biased strand distribution p-value = 7.22e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.16e-01

GGAGGTTGGCTTAGAAGCAGCCACCCTTTAAAGAAAGCGTAATAGCTCACTAGTCGAGTCGGCCTGCGCGGAAGATGTAACGGGGCTCAAACCATGCACCGAAGCTACGGGTATCATCTTA‑‑CGATGATGCGGTAGAGGAGCGTTCTGTAAGCCTGTGAAGGTGAGTTGAGAAGCTTGCTGGAGGTATCAGAAGTGCGAATGCTGACATGAGTAACGACAATGCGAGTGAAAAACTCGCACGCCGAAAG  >  NC_002947/527811‑528058
                                                                                                                                          |                                                                                                                
ggatgTTGGCTTAGAAGCAGCCACCCTTTAAAGAAAGCGTAATAGCTCACTAGTCGAGTCGGCCTGCGCGGAAGATGTAACGGGGCTCAAACCATGCACCGAAGCTACGGGTATCATCTTA‑‑CGATGATGCGGTaga                                                                                                                  >  1:291841/1‑136 (MQ=17)
ggaggTTGGCTTAGAAGCAGCCACCCTTTAAAGAAAGCGTAATAGCTCACTAGTCGAGTCGGCCTGCGCGGAAGATGTAACGGGGCTCAAACCATGCACCGAAGCTACGGGTATCATCTTA‑‑CGATGATGCGGTaga                                                                                                                  >  1:382349/1‑136 (MQ=32)
ggaggTTGGCTTAGAAGCAGCCACCCTTTAAAGAAAGCGTAATAGCTCACTAGTCGAGTCGGCCTGCGCGGAAGATGTAACGGGGCTCAAACCATGCACCGAAGCTACGGGTATCATCTTA‑‑CGATGATGCGGTaga                                                                                                                  >  1:56245/1‑136 (MQ=32)
          ttAGAAGCAGCCACCCTTTAAAGAAAGCGTAATAGCTCACTAGTCGAGTCGGCCTGCGCGGAAGATGTAACGGGGCTCAAACCATGCACCGAAGCTACGGGTATCATCTTA‑‑CGATGATGCGGTAGAGGAGCGTTcg                                                                                                        >  1:122733/1‑135 (MQ=17)
          ttAGAAGCAGCCACCCTTTAAAGAAAGCGTAATAGCTCACTAGTCGAGTCGGCCTGCGCGGAAGATGTAACGGGGCTCAAACCATGCACCGAAGCTACGGGTATCATCTTA‑‑CGATGATGCGGTAGAGGAGCGTTCt                                                                                                        >  2:50036/1‑136 (MQ=32)
          ttAGAAGCAGCCACCCTTTAAAGAAAGCGTAATAGCTCACTAGTCGAGTCGGCCTGCGCGGAAGATGTAACGGGGCTCAAACCATGCACCGAAGCTACGGGTATCATCTTA‑‑CGATGATGCGGTAGAGGAGCGTTCt                                                                                                        >  1:501373/1‑136 (MQ=32)
             gAAGCAGCCATCATTTAAAGAAAGCGTAATAGCTCACTGGTCGAGTCGGCCTGCGCGGAAGATGTAACGGGGCTCAAACCATGCACCGAAGCTACGGGTATCATCTTA‑‑CGATGATGCGGTAGAGGAGCGTTCTGTa                                                                                                     >  1:262167/1‑136 (MQ=12)
             gAAGCAGCCATCATTTAAAGAAAGCGTAATAGCTCACTGGTCGAGTCGGCCTGCGCGGAAGATGTAACGGGGCTCAAACCATGCACCGAAGCTACGGGTATCATCTTA‑‑CGATGATGCGGTAGAGGAGCGTTCTGTa                                                                                                     >  1:13756/1‑136 (MQ=12)
             gAAGCAGCCACCCTTTAAAGAAAGCGTAATAGCTCACTAGTCGAGTCGGCCTGCGCGGAAGATGTAACGGGGCTCAAACCATGCACCGAAGCTACGGGTATCATCTTA‑‑CGATGATGCGGTAGAGGAGCGTTCTGTa                                                                                                     <  2:359802/136‑1 (MQ=32)
               agcagcCACCCTTTAAAGAAAGCGTAATAGCTCACTAGTCGAGTCGGCCTGCGCGGAAGATGTAACGGGGCTCAAACCATGCACCGAAGCTACGGGTATCATCTTA‑‑CGATGATGCGGTAGAGGAGCGTTCTGTAAg                                                                                                   >  1:278815/1‑136 (MQ=32)
                     cACCCTTTAAAGAAAGCGTAATAGCTCACTAGTCGAGTCGGCCTGCGCGGAAGATGTAACGGGGCTCAAACCATGCACCGAAGCTACGGGTATCATCTTA‑‑CGATGATGCGGTAGAGGAGCGTTCTGTAAGCCtgtg                                                                                             <  2:204495/136‑1 (MQ=32)
                      aCCCTTTAAAGAAAGCGTAATAGCTCACTAGTCGAGTCGGCCTGCGCGGAAGATGTAACGGGGCTCAAACCATGCACCGAAGCTACGGGTATCATCTTA‑‑CGATGATGCGGTAGAGGAGCGTTCTGTAAGCCt                                                                                                >  2:301441/1‑132 (MQ=32)
                      aCCCTTTAAAGAAAGCGTAATAGCTCACTAGTCGAGTCGGCCTGCGCGGAAGATGTAACGGGGCTCAAACCATGCACCGAAGCTACGGGTATCATCTTA‑‑CGATGATGCGGTAGAGGAGCGTTCTGTAAGCCt                                                                                                <  1:301441/132‑1 (MQ=32)
                      aCCCTTTAAAGAAAGCGTAATAGCTCACTAGTCGAGTCGGCCTGCGCGGAAGATGTAACGGGGCTCAAACCATGCACCGAAGCTACGGGTATCATCTTA‑‑CGATGATGCGGTAGAGGAGCGTTCTGTAAGCCTGTGa                                                                                            >  1:277712/1‑136 (MQ=32)
                      aCCCTTTAAAGAAAGCGTAATAGCTCACTAGTCGAGTCGGCCTGCGCGGAAGATGTAACGGGGCTCAAACCATGCACCGAAGCTACGGGTATCATCTTA‑‑CGATGATGCGGTAGAGGAGCGTTCTGTAAGCCTGTGa                                                                                            >  1:369344/1‑136 (MQ=32)
                      aCCCTTTAAAGAAAGCGTAATAGCTCACTAGTCGAGTCGGCCTGCGCGGAAGATGTAACGGGGCTAAACCATGCACCGAAGCTGCGGCAGCGACGCTTATGCGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCtgtg                                                                                             >  2:301882/1‑136 (MQ=0)
                           ttAAAGAAAGCGTAATAGCTCACTAGTCGAGTCGGCCTGCGCGGAAGATGTAACGGGGCTCAAACCATGCACCGAAGCTACGGGTATCATCTTA‑‑CGATGATGCGGTAGAGGAGCGTTCTGTAAGCCTGTGAAGGTg                                                                                       >  2:403853/1‑136 (MQ=32)
                           ttAAAGAAAGCGTAATAGCTCACTAGTCGAGTCGGCCTGCGCGGAAGATGTAACGGGGCTCAAACCATGCACCGAAGCTACGGGTATCATCTTA‑‑CGATGATGCGGTAGAGGAGCGTTCTGTAAGCCTGTGAAGGTg                                                                                       >  2:194095/1‑136 (MQ=32)
                           atAAAGAAAGCGTAATAGCTCACTAGTCGAGTCGGCCTGCGCGGAAGATGTAACGGGGCTCAAACCATGCACCGAAGCTACGGGTATCATCTTA‑‑CGATGATGCGGTTGAGGAGCGTTCTGTAAGCCTGTGAAGGTg                                                                                       >  2:114539/2‑136 (MQ=17)
                               agaaagCGTAATAGCTCACTAGTCGAGTCGGCCTGCGCGGAAGATGTAACGGGGCTCAAACCATGCACCGAAGCTACGGGTATCATCTTA‑‑CGATGATGCGGTAGAGGAGCGTTCTGTAAGCCTGTGAAGGTGAGtt                                                                                   <  2:461877/136‑1 (MQ=32)
                                  aagAGTAATAGCTCACTAGTCGAGTCGGCCTGCGCGGAAGATGTAACGGGGCTCAATCCATGCACCGAAGCTACGGGTATCATCTTA‑‑CGATGTTGCGGTAGAGGAGCGTTCTGTAAGCCTGTGAAGGTGAGTTgag                                                                                <  2:11725/136‑1 (MQ=12)
                                     cGTAATAGCTCACTAGTCGAGTCGGCCTGCGCGGAAGATGTAACGGGGCTCAAACCATGCACCGAAGCTACGGGTATCATCTTA‑‑CGATGATGCGGTAGAGGAGCGTTCTGTAAGCCTGTGAAGGTGAGTTGAGAAg                                                                             >  1:483121/1‑136 (MQ=32)
                                     cGTAATAGCTCACTAGTCGAGTCGGCCTGCGCGGAAGATGTAACGGGGCTCAAACCATGCACCGAAGCTACGGGTATCATCTTA‑‑CGATGATGCGGTAGAGGAGCGTTCTGTAAGCCTGTGAAGGTGAGTTGAGAAg                                                                             >  2:102094/1‑136 (MQ=32)
                                     cGTAATAGCTCACTAGTCGAGTCGGCCTGCGCGGAAGATGTAACGGGGCTCAAACCATGCACCGAAGCTACGGGTATCATCTTA‑‑CGATGATGCGGTAGAGGAGCGTTCTGTAAGCCTGTGAAGGTGAGTTGAGAAg                                                                             >  2:441514/1‑136 (MQ=32)
                                           aGCTCACTAGTCGAGTCGGCCTGCGCGGAAGATGTAACGGGGCTCAAACCATGCACCGAAGCTACGGGTATCATCTTA‑‑CGATGATGCGGTAGAGGAGCGTTCTGTAAGCCTGTGAAGGTGAGTTGAGAAGCTTGCt                                                                       <  2:278815/136‑1 (MQ=32)
                                           aGCTCACTAGTCGAGTCGGCCTGCGCGGAAGATGTAACGGGGCTCAAACCATGCACCGAAGCTACGGGTATCATCTTA‑‑CGATGATGCGGTAGAGGAGCGTTCTGTAAGCCTGTGAAGGTGAGTTGAGAAGCTTGCt                                                                       <  2:12824/136‑1 (MQ=32)
                                             cTCACTGGTCGAGTCGGCCTGCGCGGAAGATGTAACGGGGCTCAAACCATGCACCGAAGCTACGGGTATCATCTTA‑‑CGATGATGCGGTAGAGGAGCGTTCTGTAAGCCTGTGAAGGTGAGTTGAGAAGCTTGCTgg                                                                     <  2:396659/136‑1 (MQ=17)
                                                  tAGTCGAGTCGGCCTGCGCGGAAGATGTAACGGGGCTCAAACCATGCACCGAAGCTACGGGTATCATCTTA‑‑CGATGATGCGGTAGAGGAGCGTTCTGTAAGCCTGTGAAGGTGAGTTGAGAAGCTTGCTGGAGGTa                                                                >  1:412347/1‑136 (MQ=32)
                                                   agtcgagtcgGCCTGCGCGGAAGATGTAACGGGGCTCAAACCATGCACCGAAGCTACGGGTATCATCTTA‑‑CGATGATGCGGTAGAGGAGCGTTCTGTAAGCCTGTGAAGGTGAGTTGAGAAGCTTGCTGGAGGTAt                                                               >  1:512455/1‑136 (MQ=32)
                                                         gtcgGCCTGCGCGGAAGATGTAACGGGGCTCAAACCATGCACCGAAGCTACGGGTATCATCTTA‑‑CGATGATGCGGTAGAGGAGCGTTCTGTAAGCCTGTGAAGGTGAGTTGAGAAGCTTGCTGGAGGTATCAGAAg                                                         <  1:323495/136‑1 (MQ=32)
                                                             gCCTGCGCGGAAGATGTAACGGGGCTCAAACCATGCACCGAAGCTACGGGTATCATCTTA‑‑CGATGATGCGGTAGAGGAGCGTTCTGTAAGCCTGTGAAGGTGAGTTGAGAAGCTTGCTGGAGGTATCAGAAGTGCg                                                     >  1:139107/1‑136 (MQ=32)
                                                             gCCTGCGCGGAAGATGTAACGGGGCTAAACCATGCACCGAAGCTGCGGCAGCGACGCTTATGCGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGAGTTGAGAAGCTTGCTGGAGGTATCAGAAGTGc                                                      >  1:347324/1‑136 (MQ=0)
                                                             gCCTGCGCGGAAGATGTAACGGGGCTAAACCATGCACCGAAGCTGCGGCAGCGACGCTTATGCGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGCGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGc                                                      >  1:285840/1‑136 (MQ=0)
                                                               cTGCGCGGAAGATGTAACGGGGCTCAAACCATGCACCGAAGCTACGGGTATCATCTTA‑‑CGATGATGCGGTAGAGGAGCGTTCTGTAAGCCTGTGAAGGTGAGTTGAGAAGCTTGCTGGAGGTATCAGAAGTGCGaa                                                   >  1:296085/1‑136 (MQ=32)
                                                                    cgGAAGATGTAACGGGGCTCAAACCATGCACCGAAGCTACGGGTATCATCTTA‑‑CGATGATGCGGTAGAGGAGCGTTCTGTAAGCCTGTGAAGGTGAGTTGAGAAGCTTGCTGGAGGTATCAGAAGTGCGAATGCTg                                              <  1:320157/136‑1 (MQ=32)
                                                                    cgGAAGATGTAACGGGGCTCAAACCATGCACCGAAGCTACGGGTATCATCTTA‑‑CGATGATGCGGTAGAGGAGCGTTCTGTAAGCCTGTGAAGGTGAGTTGAGAAGCTTGCTGGAGGTATCAGAAGTGCGAATGCTg                                              <  1:71651/136‑1 (MQ=32)
                                                                     ggAAGATGTAACGGGGCTAAACCATGCACCGAAGCTACGGGTATCATCTTA‑‑CGATGATGCGGTAGAGGAGCGTTCTGTAAGCCTGTGAAGGTGAGTTGAGAAGCTTGCTGGAGGTATCAGAAGTGCGAATGCTGAc                                            <  2:136313/136‑1 (MQ=17)
                                                                      gAAGATGTAACGGGGCTCAAACCATGCACCGAAGCTACGGGTATCATCTTA‑‑CGATGATGCGGTAGAGGAGCGTTCTGTAAGCCTGTGAAGGTGAGTTGAGAAGCTTGCTGGAGGTATCAGAAGTGCGAATGCTGAc                                            <  2:501373/136‑1 (MQ=32)
                                                                      gAAGATGTAACGGGGCTCAAACCATGCACCGAAGCTACGGGTATCATCTTA‑‑CGATGATGCGGTAGAGGAGCGTTCTGTAAGCCTGTGAAGGTGAGTTGAGAAGCTTGCTGGAGGTATCAGAAGTGCGAATGCTGAc                                            <  2:510769/136‑1 (MQ=32)
                                                                       aaGATGTAACGGGGCTAAACCATGCACCGAAGCTGCGGCAGCGACGCTTATGCGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTACTGTGAGGTATGCTGGAGGTATCAGAAGTGCGAATGCTGAc                                            <  2:199546/136‑1 (MQ=0)
                                                                       aaGATGTAACGGGGCTAAACCATGCACCGAAGCTGCGGCAGCGACGCTTATGCGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTACTGTGAGGTATGCTGGAGGTATCAGAAGTGCGAATGCTGAc                                            <  2:29172/136‑1 (MQ=0)
                                                                       aaGATGTAACGGGGCTAAACCATGCACCGAAGCTGCGGCAGCGACGCTTATGCGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGCGAAGGTGAGTTGAGAAGCTTGCTGGAGGTATCAGAAGTGCGAATGCTGAc                                            <  2:32206/136‑1 (MQ=0)
                                                                        aGATGTAACGGGGCTCAAACCATGCACCGAAGCTACGGGTATCATCTTA‑‑CGATGATGCGGTAGAGGAGCGTTCTGTAAGCCTGTGAAGGTGAGTTGAGAAGCTTGCTGGAGGTATCAGAAGTGCGAATGCTGACAt                                          <  2:501083/136‑1 (MQ=32)
                                                                        aGATGTAACGGGGCTCAAACCATGCACCGAAGCTACGGGTATCATCTTA‑‑CGATGATGCGGTAGAGGAGCGTTCTGTAAGCCTGTGAAGGTGAGTTGAGAAGCTTGCTGGAGGTATCAGAAGTGCGAATGCTGACAt                                          <  2:475210/136‑1 (MQ=32)
                                                                        aGATGTAACGGGGCTCAAACCATGCACCGAAGCTACGGGTATCATCTTA‑‑CGATGATGCGGTAGAGGAGCGTTCTGTAAGCCTGTGAAGGTGAGTTGAGAAGCTTGCTGGAGGTATCAGAAGTGCGAATGCTGACAt                                          <  2:91363/136‑1 (MQ=32)
                                                                        aGATGTAACGGGGCTCAAACCATGCACCGAAGCTACGGGTATCATCTTA‑‑CGATGATGCGGTAGAGGAGCGTTCTGTAAGCCTGTGAAGGTGAGTTGAGAAGCTTGCTGGAGGTATCAGAAGTGCGAATGCTGACAt                                          <  2:512455/136‑1 (MQ=32)
                                                                        aGATGTAACGGGGCTCAAACCATGCACCGAAGCTACGGGTATCATCTTA‑‑CGATGATGCGGTAGAGGAGCGTTCTGTAAGCCTGTGAAGGTGAGTTGAGAAGCTTGCTGGAGGTATCAGAAGTGCGAATGCTGACAt                                          <  2:388595/136‑1 (MQ=32)
                                                                        aGATGTAACGGGGCTCAAACCATGCACCGAAGCTACGGGTATCATCTTA‑‑CGATGATGCGGTAGAGGAGCGTTCTGTAAGCCTGTGAAGGTGAGTTGAGAAGCTTGCTGGAGGTATCAGAAGTGCGAATGCTGACAt                                          <  2:382349/136‑1 (MQ=32)
                                                                        aGATGTAACGGGGCTCAAACCATGCACCGAAGCTACGGGTATCATCTTA‑‑CGATGATGCGGTAGAGGAGCGTTCTGTAAGCCTGTGAAGGTGAGTTGAGAAGCTTGCTGGAGGTATCAGAAGTGCGAATGCTGACAt                                          <  2:277712/136‑1 (MQ=32)
                                                                        aGATGTAACGGGGCTCAAACCATGCACCGAAGCTACGGGTATCATCTTA‑‑CGATGATGCGGTAGAGGAGCGTTCTGTAAGCCTGTGAAGGTGAGTTGAGAAGCTTGCTGGAGGTATCAGAAGTGCGAATGCTGACAt                                          <  2:291841/136‑1 (MQ=32)
                                                                        aGATGTAACGGGGCTCAAACCATGCACCGAAGCTACGGGTATCATCTTA‑‑CGATGATGCGGTAGAGGAGCGTTCTGTAAGCCTGTGAAGGTGAGTTGAGAAGCTTGCTGGAGGTATCAGAAGTGCGAATGCTGACAt                                          <  2:281105/136‑1 (MQ=32)
                                                                        aGATGTAACGGGGCTCAAACCATGCACCGAAGCTACGGGTATCATCTTA‑‑CGATGATGCGGTAGAGGAGCGTTCTGTAAGCCTGTGAAGGTGAGTTGAGAAGCTTGCTGGAGGTATCAGAAATGCGAATGCTGACAt                                          <  2:369344/136‑1 (MQ=17)
                                                                         gATGTAACGGGGCTAAACCATGCACCGAAGCTGCGGCAGCGACGCTTATGCGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGAGTTGAGAAGCTTGCTGGAGGTATCAGAAGTGCGAATGCTGACAt                                          <  2:347324/136‑1 (MQ=0)
                                                                         gATGTAACGGGGCTAAACCATGCACCGAAGCTGCGGCAGCGACGCTTATGCGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGAGTTGAGAAGCTTGCTGGAGGTATCAGAAGTGCGAATGCTGACAt                                          <  2:473222/136‑1 (MQ=0)
                                                                          ttGTAACGGGGCTAAACCATGCACCGAAGCTGCGGCAGCGACGCTTATGCGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGCGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACAta                                         <  2:109993/135‑2 (MQ=0)
                                                                                      tCAAACCATGCACCGAAGCTACGGGTATCATCTTA‑‑CGATGATGCGGTAGAGGAGCGTTCTGTAAGCCTGTGAAGGTGAGTTGAGAAGCTTGCTGGAGGTATCAGAAGTGCGAATGCTGACATGAGTAACGACAATg                            >  1:78645/1‑136 (MQ=32)
                                                                                      tCAAACCATGCACCGAAGCTACGGGTATCATCTTA‑‑CGATGATGCGGTAGAGGAGCGTTCTGTAAGCCTGTGAAGGTGAGTTGAGAAGCTTGCTGGAGGTATCAGAAGTGCGAATGCTGACATGAGTAACGACAATg                            <  2:400312/136‑1 (MQ=32)
                                                                                        aaaCCATGCACCGAAGCTACGGGTATCATCTTA‑‑CGATGATGCGGTAGAGGAGCGTTCTGTAAGCCTGTGAAGGTGAGTTGAGAAGCTTGCTGGAGGTATCAGAAGTGCGAATGCTGACATGAGTAACGACAATGCg                          >  1:297680/1‑136 (MQ=32)
                                                                                        aaaCCATGCACCGAAGCTACGGGTATCATCTTA‑‑CGATGATGCGGTAGAGGAGCGTTCTGTAAGCCTGTGAAGGTGAGTTGAGAAGCTTGCTGGAGGTATCAGAAGTGCGAATGCTGACATGAGTAACGACAATGCg                          >  1:254043/1‑136 (MQ=32)
                                                                                                 aCCGAAGCTACGGGTATCATCTTA‑‑CGATGATGCGGTAGAGGAGCGTTCTGTAAGCCTGTGAAGGTGAGTTGAGAAGCTTGCTGGAGGTATCAGAAGTGCGAATGc                                                >  1:247808/1‑105 (MQ=32)
                                                                                                 aCCGAAGCTACGGGTATCATCTTA‑‑CGATGATGCGGTAGAGGAGCGTTCTGTAAGCCTGTGAAGGTGAGTTGAGAAGCTTGCTGGAGGTATCAGAAGTGCGAATGc                                                <  2:247808/105‑1 (MQ=32)
                                                                                                              gTATCATCTTA‑‑CGATGATGCGGTAGAGGAGCGTTCTGTAAGCCTGTGAAGGTGAGTTGAGAAGCTTGCTGGAGGTATCAGAAGTGCGAATGTTGACATGAGTAACGACAATGCGAGTGAAAAACTCGCACGCCGaa    <  2:139107/136‑1 (MQ=17)
                                                                                                              gTATCATCTTA‑‑CGATGATGCGGTAGAGGAGCGTTCTGTAAGCCTGTGAAGGTGAGTTGAGAAGCTTGCTGGAGGTATCAGAAGTGCGAATGCTGACATGAGTAACGACAATGCGAGTGAAAAACTCGCACGCCGaa    <  2:48326/136‑1 (MQ=32)
                                                                                                              gTATCATCTTA‑‑CGATGATGCGGTAGAGGAGCGTTCTGTAAGCCTGTGAAGGTGAGTTGAGAAGCTTGCTGGAGGTATCAGAAGTGCGAATGCTGACATGAGTAACGACAATGCGAGTGAAAAACTCGCACGCCGaa    <  2:13756/136‑1 (MQ=32)
                                                                                                              gTATCATCTTA‑‑CGATGATGCGGTAGAGGAGCGTTCTGTAAGCCTGTGAAGGTGAGTTGAGAAGCTTGCTGGAGGTATCAGAAGTGCGAATGCTGACATGAGTAACGACAATGCGAGTGAAAAACTCGCACGCCGaa    <  2:296085/136‑1 (MQ=32)
                                                                                                                atcatcTTA‑‑CGATGATGCGGTAGTGGCGCGTTCTGTAAGCCTGTGATGGTGAGTTGAGAAGCTTGCTGGAGGTATCAGAAGGGCGAATGCTGACATGAGTAACGACAATGCGAGTGAAAAACTCGCACGCCGAAAg  <  2:510045/136‑1 (MQ=255)
                                                                                                                                          |                                                                                                                
GGAGGTTGGCTTAGAAGCAGCCACCCTTTAAAGAAAGCGTAATAGCTCACTAGTCGAGTCGGCCTGCGCGGAAGATGTAACGGGGCTCAAACCATGCACCGAAGCTACGGGTATCATCTTA‑‑CGATGATGCGGTAGAGGAGCGTTCTGTAAGCCTGTGAAGGTGAGTTGAGAAGCTTGCTGGAGGTATCAGAAGTGCGAATGCTGACATGAGTAACGACAATGCGAGTGAAAAACTCGCACGCCGAAAG  >  NC_002947/527811‑528058

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