Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 527,934 A→T 13.0% noncoding (1169/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,9340AT13.0% 153.9 / 15.3 67noncoding (1169/2902 nt)PP_23SC23S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base A (27/31);  new base T (3/6);  total (30/37)
Fisher's exact test for biased strand distribution p-value = 7.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.91e-01

CTTAGACAGCTAGGAGGTTGGCTTAGAAGCAGCCACCCTTTAAAGAAAGCGTAATAGCTCACTAGTCGAGTCGGCCTGCGCGGAAGATGTAACGGGGCTCAAACCATGCACCGAAGCTACGGGTATCATCTTA‑‑CGATGATGCGGTAGAGGAGCGTTCTGTAAGCCTGTGAAGGTGAGTTGAGAAGCTTGCTGGAGGTATCAGAAGTGCGAATGCTGACATGAGTAACGACAATGCGAGTGAAAAACTCGCACGCCGAAAG  >  NC_002947/527799‑528058
                                                                                                                                          |                                                                                                                            
cTTAGACAGCTAGGAGGTTGGCTTAGAAGCAGCCACCCTTTAAAGAAAGCGTAATAGCTCACTAGTCGAGTCGGCCTGCGCGGAAGATGTAACGGGGCTCAAACCATGCACCGAAGCTACGGGTATCATCTTA‑‑Cga                                                                                                                              <  2:208338/136‑1 (MQ=32)
  tAGACAGCTAGGAGGTTGGCTTAGAAGCAGCCACCCTTTAAAGAAAGCGTAATAGCTCACTAGTCGAGTCGGCCTGCGCGGAAGATGTAACGGGGCTCAAACCATGCACCGAAGCTACGGGTATCATCTTA‑‑Cgatg                                                                                                                            <  1:530552/136‑1 (MQ=32)
   aGACAGCTAGGAGGTTGGCTTAGAAGCAGCCACCCTTTAAAGAAAGCGTAATAGCTCACTAGTCGAGTCGGCCTGCGCGGAAGATGTAACGGGGCTCAAACCATGCACCGAAGCTACGGGTATCATCTTA‑‑Cgatga                                                                                                                           >  1:461877/1‑136 (MQ=32)
            ggatgTTGGCTTAGAAGCAGCCACCCTTTAAAGAAAGCGTAATAGCTCACTAGTCGAGTCGGCCTGCGCGGAAGATGTAACGGGGCTCAAACCATGCACCGAAGCTACGGGTATCATCTTA‑‑Cga                                                                                                                              >  1:128240/1‑124 (MQ=17)
            ggatgTTGGCTTAGAAGCAGCCACCCTTTAAAGAAAGCGTAATAGCTCACTAGTCGAGTCGGCCTGCGCGGAAGATGTAACGGGGCTCAAACCATGCACCGAAGCTACGGGTATCATCTTA‑‑CGATGATGCGGTaga                                                                                                                  >  1:291841/1‑136 (MQ=17)
            ggatgTTGGCTTAGAAGCAGCCACCCTTTAAAGAAAGCGTAATAGCTCACTAGTCGAGCCGGCCTGCGCGGAAGATGTAACGGGGCTCAAACCATGCACCGAAGCTACGGGTATCATCTTA‑‑Cga                                                                                                                              <  2:128240/124‑1 (MQ=17)
            ggatgTTGGCTTAGAAGCAGCCACCCTTTAAAGAAAGCGGAATAGCTCACTAGTCGAGTCGGCCTGCGCGGAAGATGTAACGGGGCTCAACCCATGCACCGAAGCTACGGGTATCATCTTA‑‑CGATGATGCGGTagg                                                                                                                  >  1:278156/1‑135 (MQ=12)
            ggaggTTGGCTTAGAAGCAGCCACCCTTTAAAGAAAGCGTAATAGCTCACTAGTCGAGTCGGCCTGCGCGGAAGATGTAACGGGGCTCAAACCATGCACCGAAGCTACGGGTATCATCTTA‑‑CGATGATGCGGTaga                                                                                                                  >  1:56245/1‑136 (MQ=32)
            ggaggTTGGCTTAGAAGCAGCCACCCTTTAAAGAAAGCGTAATAGCTCACTAGTCGAGTCGGCCTGCGCGGAAGATGTAACGGGGCTCAAACCATGCACCGAAGCTACGGGTATCATCTTA‑‑CGATGATGCGGTaga                                                                                                                  >  1:382349/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:13756/1‑136 (MQ=12)
                         gAAGCAGCCATCATTTAAAGAAAGCGTAATAGCTCACTGGTCGAGTCGGCCTGCGCGGAAGATGTAACGGGGCTCAAACCATGCACCGAAGCTACGGGTATCATCTTA‑‑CGATGATGCGGTAGAGGAGCGTTCTGTa                                                                                                     >  1:262167/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:369344/1‑136 (MQ=32)
                                  aCCCTTTAAAGAAAGCGTAATAGCTCACTAGTCGAGTCGGCCTGCGCGGAAGATGTAACGGGGCTCAAACCATGCACCGAAGCTACGGGTATCATCTTA‑‑CGATGATGCGGTAGAGGAGCGTTCTGTAAGCCTGTGa                                                                                            >  1:277712/1‑136 (MQ=32)
                                  aCCCTTTAAAGAAAGCGTAATAGCTCACTAGTCGAGTCGGCCTGCGCGGAAGATGTAACGGGGCTAAACCATGCACCGAAGCTGCGGCAGCGACGCTTATGCGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCtgtg                                                                                             >  2:301882/1‑136 (MQ=0)
                                       ttAAAGAAAGCGTAATAGCTCACTAGTCGAGTCGGCCTGCGCGGAAGATGTAACGGGGCTCAAACCATGCACCGAAGCTACGGGTATCATCTTA‑‑CGATGATGCGGTAGAGGAGCGTTCTGTAAGCCTGTGAAGGTg                                                                                       >  2:194095/1‑136 (MQ=32)
                                       ttAAAGAAAGCGTAATAGCTCACTAGTCGAGTCGGCCTGCGCGGAAGATGTAACGGGGCTCAAACCATGCACCGAAGCTACGGGTATCATCTTA‑‑CGATGATGCGGTAGAGGAGCGTTCTGTAAGCCTGTGAAGGTg                                                                                       >  2:403853/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:441514/1‑136 (MQ=32)
                                                 cGTAATAGCTCACTAGTCGAGTCGGCCTGCGCGGAAGATGTAACGGGGCTCAAACCATGCACCGAAGCTACGGGTATCATCTTA‑‑CGATGATGCGGTAGAGGAGCGTTCTGTAAGCCTGTGAAGGTGAGTTGAGAAg                                                                             >  2:102094/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:71651/136‑1 (MQ=32)
                                                                                cgGAAGATGTAACGGGGCTCAAACCATGCACCGAAGCTACGGGTATCATCTTA‑‑CGATGATGCGGTAGAGGAGCGTTCTGTAAGCCTGTGAAGGTGAGTTGAGAAGCTTGCTGGAGGTATCAGAAGTGCGAATGCTg                                              <  1:320157/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:91363/136‑1 (MQ=32)
                                                                                    aGATGTAACGGGGCTCAAACCATGCACCGAAGCTACGGGTATCATCTTA‑‑CGATGATGCGGTAGAGGAGCGTTCTGTAAGCCTGTGAAGGTGAGTTGAGAAGCTTGCTGGAGGTATCAGAAGTGCGAATGCTGACAt                                          <  2:512455/136‑1 (MQ=32)
                                                                                    aGATGTAACGGGGCTCAAACCATGCACCGAAGCTACGGGTATCATCTTA‑‑CGATGATGCGGTAGAGGAGCGTTCTGTAAGCCTGTGAAGGTGAGTTGAGAAGCTTGCTGGAGGTATCAGAAGTGCGAATGCTGACAt                                          <  2:382349/136‑1 (MQ=32)
                                                                                    aGATGTAACGGGGCTCAAACCATGCACCGAAGCTACGGGTATCATCTTA‑‑CGATGATGCGGTAGAGGAGCGTTCTGTAAGCCTGTGAAGGTGAGTTGAGAAGCTTGCTGGAGGTATCAGAAGTGCGAATGCTGACAt                                          <  2:388595/136‑1 (MQ=32)
                                                                                    aGATGTAACGGGGCTCAAACCATGCACCGAAGCTACGGGTATCATCTTA‑‑CGATGATGCGGTAGAGGAGCGTTCTGTAAGCCTGTGAAGGTGAGTTGAGAAGCTTGCTGGAGGTATCAGAAGTGCGAATGCTGACAt                                          <  2:501083/136‑1 (MQ=32)
                                                                                    aGATGTAACGGGGCTCAAACCATGCACCGAAGCTACGGGTATCATCTTA‑‑CGATGATGCGGTAGAGGAGCGTTCTGTAAGCCTGTGAAGGTGAGTTGAGAAGCTTGCTGGAGGTATCAGAAGTGCGAATGCTGACAt                                          <  2:475210/136‑1 (MQ=32)
                                                                                    aGATGTAACGGGGCTCAAACCATGCACCGAAGCTACGGGTATCATCTTA‑‑CGATGATGCGGTAGAGGAGCGTTCTGTAAGCCTGTGAAGGTGAGTTGAGAAGCTTGCTGGAGGTATCAGAAGTGCGAATGCTGACAt                                          <  2:291841/136‑1 (MQ=32)
                                                                                    aGATGTAACGGGGCTCAAACCATGCACCGAAGCTACGGGTATCATCTTA‑‑CGATGATGCGGTAGAGGAGCGTTCTGTAAGCCTGTGAAGGTGAGTTGAGAAGCTTGCTGGAGGTATCAGAAGTGCGAATGCTGACAt                                          <  2:281105/136‑1 (MQ=32)
                                                                                    aGATGTAACGGGGCTCAAACCATGCACCGAAGCTACGGGTATCATCTTA‑‑CGATGATGCGGTAGAGGAGCGTTCTGTAAGCCTGTGAAGGTGAGTTGAGAAGCTTGCTGGAGGTATCAGAAGTGCGAATGCTGACAt                                          <  2:277712/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                            <  2:400312/136‑1 (MQ=32)
                                                                                                  tCAAACCATGCACCGAAGCTACGGGTATCATCTTA‑‑CGATGATGCGGTAGAGGAGCGTTCTGTAAGCCTGTGAAGGTGAGTTGAGAAGCTTGCTGGAGGTATCAGAAGTGCGAATGCTGACATGAGTAACGACAATg                            >  1:78645/1‑136 (MQ=32)
                                                                                                    aaaCCATGCACCGAAGCTACGGGTATCATCTTA‑‑CGATGATGCGGTAGAGGAGCGTTCTGTAAGCCTGTGAAGGTGAGTTGAGAAGCTTGCTGGAGGTATCAGAAGTGCGAATGCTGACATGAGTAACGACAATGCg                          >  1:254043/1‑136 (MQ=32)
                                                                                                    aaaCCATGCACCGAAGCTACGGGTATCATCTTA‑‑CGATGATGCGGTAGAGGAGCGTTCTGTAAGCCTGTGAAGGTGAGTTGAGAAGCTTGCTGGAGGTATCAGAAGTGCGAATGCTGACATGAGTAACGACAATGCg                          >  1:297680/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:13756/136‑1 (MQ=32)
                                                                                                                          gTATCATCTTA‑‑CGATGATGCGGTAGAGGAGCGTTCTGTAAGCCTGTGAAGGTGAGTTGAGAAGCTTGCTGGAGGTATCAGAAGTGCGAATGCTGACATGAGTAACGACAATGCGAGTGAAAAACTCGCACGCCGaa    <  2:48326/136‑1 (MQ=32)
                                                                                                                          gTATCATCTTA‑‑CGATGATGCGGTAGAGGAGCGTTCTGTAAGCCTGTGAAGGTGAGTTGAGAAGCTTGCTGGAGGTATCAGAAGTGCGAATGCTGACATGAGTAACGACAATGCGAGTGAAAAACTCGCACGCCGaa    <  2:296085/136‑1 (MQ=32)
                                                                                                                            atcatcTTA‑‑CGATGATGCGGTAGTGGCGCGTTCTGTAAGCCTGTGATGGTGAGTTGAGAAGCTTGCTGGAGGTATCAGAAGGGCGAATGCTGACATGAGTAACGACAATGCGAGTGAAAAACTCGCACGCCGAAAg  <  2:510045/136‑1 (MQ=255)
                                                                                                                                          |                                                                                                                            
CTTAGACAGCTAGGAGGTTGGCTTAGAAGCAGCCACCCTTTAAAGAAAGCGTAATAGCTCACTAGTCGAGTCGGCCTGCGCGGAAGATGTAACGGGGCTCAAACCATGCACCGAAGCTACGGGTATCATCTTA‑‑CGATGATGCGGTAGAGGAGCGTTCTGTAAGCCTGTGAAGGTGAGTTGAGAAGCTTGCTGGAGGTATCAGAAGTGCGAATGCTGACATGAGTAACGACAATGCGAGTGAAAAACTCGCACGCCGAAAG  >  NC_002947/527799‑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: