Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,958,935 A→G 5.5% Q780R (CAA→CGA)  cyaA → adenylate cyclase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,958,9350AG 5.5% 214.1 / 2.0 73Q780R (CAA→CGA) cyaAadenylate cyclase
Reads supporting (aligned to +/- strand):  ref base A (36/33);  new base G (2/2);  total (38/35)
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

CGCTGTACCTGGATGCCAATGCCCTGCAGGATCACGACCTGCGCCTGGTGCTGGAGCAAGGCCGGCCAGCCTGTATTCAGGTGTTCTACCGCCTGCAGGGCGGCTGGGCCGACGTGTATGTGCTGGATGAATACAACGCATTGTGGCAACAACGCCTGCCGCTGCATGACGAAGCGCATCTGCTGTTGCCACTGCAACGCTTTCTGCGCTCGGTGGTGATGCGCCGCGATGCACGCCTGCCCCTCGATAGCCTGCGGCCCGCCTCGCTGGATATCCACTTCGCGCAGCTGTTGCC  >  NC_002947/5958788‑5959082
                                                                                                                                                   |                                                                                                                                                   
cGCTGTACCTGGATGCCAATGCCCTGCAGGATCACGACCTGCGCCTGGTGCTGGAGCAAGGCCGGCCAGCCTGTATTCAGGTGTTCTACCGCCTGCAGGGCGGCTGGGCCGACGTGTATGTGCTGGATGAATACAACGCATTGTGGcaa                                                                                                                                                    <  3:279021/149‑1 (MQ=255)
        cTGGATGCCAATGCCCTGCAGGATCACGACCTGCGCCTGGTGCTGGAGCAAGGCCGGCCAGCCTGTATTCAGGTGTTCTACCGCCTGCAGGGCGGCTGGGCCGACGTGTATGTGCTGGATGAATACAACGCATTGTGGCAACAACGCCt                                                                                                                                            >  6:139613/1‑149 (MQ=255)
        cTGGATGCCAATGCCCTGCAGGATCACGACCTGCGCCTGGTGCTGGAGCAAGGCCGGCCAGCCTGTATTCAGGTGTTCTACCGCCTGCAGGGCGGCTGGGCCGACGTGTATGTGCTGGATGAATACAACGCATTGTGGCAACAACGCCt                                                                                                                                            >  7:248542/1‑149 (MQ=255)
           gATGCCAATGCCCTGCAGGATCACGACCTGCGCCTGGTGCTGGAGCAAGGCCGGCCAGCCTGTATTCAGGTGTTCTACCGCCTGCAGGGCGGCTGGGCCGACGTGTATGTGCTGGATGAATACAACGCATTGTGGCAACAACGCCTGc                                                                                                                                          <  1:11832/148‑1 (MQ=255)
              gCCAATGCCCTGCAGGATCACGACCTGCGCCTGGTGCTGGAGCAAGGCCGGCCAGCCTGTATTCAGGTGTTCTACCGCCTGCAGGGCGGCTGGGCCGACGTGTATGTGCTGGATGAATACAACGCATTGTGGCAACAACGCCTGCCGCt                                                                                                                                      >  5:163987/1‑149 (MQ=255)
                cAATGCCCTGCAGGATCACGACCTGCGCCTGGTGCTGGAGCAAGGCCGGCCAGCCTGTATTCAGGTGTTCTACCGCCTGCAGGGCGGCTGGGCCGACGTGTATGTGCTGGATGAATACAACGCATTGTGGCAACAACGCCTGCCGCTg                                                                                                                                     >  6:171815/1‑148 (MQ=255)
                cAATGCCCTGCAGGATCACGACCTGCGCCTGGTGCTGGAGCAAGGCCGGCCAGCCTGTATTCAGGTGTTCTACCGCCTGCAGGGCGGCTGGGCCGACGTGTATGTGCTGGATGAATACAACGCATTGTGGCAACAACGCCTGCCGCTg                                                                                                                                     >  1:242792/1‑148 (MQ=255)
                cAATGCCCTGCAGGATCACGACCTGCGCCTGGTGCTGGAGCAAGGCCGGCCAGCCTGTATTCAGGTGTTCTACCGCCTGCAGGGCGGCTGGGCCGACGTGTATGTGCTGGATGAATACAACGCATTGTGGCAACAACGCCTGCCGCTGc                                                                                                                                    <  8:136886/149‑1 (MQ=255)
                  aTGCCCTGCAGGATCACGACCTGCGCCTGGTGCTGGAGCAAGGCCGGCCAGCCTGTATTCAGGTGTTCTACCGCCTGCAGGGCGGCTGGGCCGACGTGTATGTGCTGGATGAATACAACGCATTGTGGCAACAACGCCTGCCGCTGCAt                                                                                                                                  <  6:163987/149‑1 (MQ=255)
                       cTGCAGGATCACGACCTGCGCCTGGTGCTGGAGCAAGGCCGGCCAGCCTGTATTCAGGTGTTCTACCGCCTGCAGGGCGGCTGGGCCGACGTGTATGTGCTGGATGAATACAACGCATTGTGGCAACAACGCCTGCCGCTGCAt                                                                                                                                  >  7:267910/1‑144 (MQ=255)
                       cTGCAGGATCACGACCTGCGCCTGGTGCTGGAGCAAGGCCGGCCAGCCTGTATTCAGGTGTTCTACCGCCTGCAGGGCGGCTGGGCCGACGTGTATGTGCTGGATGAATACAACGCATTGTGGCAACAACGCCTGCCGCTGCAt                                                                                                                                  <  8:267910/144‑1 (MQ=255)
                        tGCAGGATCACGACCTGCGCCTGGTGCTGGGGCAAGGCCGGCCAGCCTGTATTCAGGTGTTCTACCGCCTGCAGGGCGGCTGGGCCGACGTGTATGTGCTGGATGAATACAACGCATTGTGGCAACAACGCCTGCCGCTGCATGACg                                                                                                                              >  7:69068/1‑147 (MQ=255)
                        tGCAGGATCACGACCTGCGCCTGGTGCTGGAGCAAGGCCGGCCAGCCTGTATTCAGGTGTTCTACCGCCTGCAGGGCGGCTGGGCGGACGTGTATGTGCTGGATGAATACAACGCATTGTGGCATCAACGCCTGGCGCTGCATGACg                                                                                                                              >  4:58118/1‑147 (MQ=255)
                        tGCAGGATCACGACCTGCGCCTGGTGCTGGAGCAAGGCCGGCCAGCCTGTATTCAGGTGTTCTACCGCCTGCAGGGCGGCTGGGCCGACGTGTATGTGCTGGATGAATACAACGCATTGTGGCAACAACGCCTGCCGCTGCATGACg                                                                                                                              >  5:299963/1‑147 (MQ=255)
                          cAGGATCACGACCTGCGCCTGGTGCTGGAGCAAGGCCGGCCAGCCTGTATTCAGGTGTTCTACCGCCTGCAGGGCGGCTGGGCCGACGTGTATGTGCTGGATGAATACAACGCATTGTGGCAACAACGCCTGCCGCTGCATGACGAAgc                                                                                                                          <  5:155718/149‑1 (MQ=255)
                            ggATCACGACCTGCGCCTGGTGCTGGAGCAAGGCCGGCCAGCCTGTATTCAGGTGTTCTACCGCCTGCAGGGCGGCTGGGCCGACGTGTATGTGCTGGATGAATACAACGCATTGTGGCAACAACGCCTGCCGCTGCATGACGAAgcg                                                                                                                         >  3:261232/1‑148 (MQ=255)
                                cACGACCTGCGCCTGGTGCTGGAGCAAGGCCGGCCAGCCTGTATTCAGGTGTTCTACCGCCTGCAGGGCGGCTGGGCCGACGTGTATGTGCTGGATGAATACAACGCATTGTGGCAACAACGCCTGCCGCTGCATGACGAAGCGCATct                                                                                                                    >  7:249520/1‑149 (MQ=255)
                                  cGACCTGCGCCTGGTGCTGGAGCAAGGCCGGCCAGCCTGTATTCAGGTGTTCTACCGCCTGCAGGGCGGCTGGGCCGACGTGTATGTGCTGGATGAATACAACGCATTGTGGCAACAACGCCTGCCGCTGCATGACGAAGCGCATctgc                                                                                                                  >  2:172639/1‑149 (MQ=255)
                                  cGACCTGCGCCTGGTGCTGGAGCAAGGCCGGCCAGCCTGTATTCAGGTGTTCTACCGCCTGCAGGGCGGCTGGGCCGACGTGTATGTGCTGGATGAATACAACGCATTGTGGCAACAACGCCTGCCGCTGCATGACGAAGCGCATctgc                                                                                                                  >  2:264749/1‑149 (MQ=255)
                                         cgcCTGGTGCTGGAGCAAGGCCGGCCAGCCTGTATTCAGGTGTTCTACCGCCTGCAGGGCGGCTGGGCCGACGTGTATGTGCTGGATGAATACAACGCATTGTGGCAACAACGCCTGCCGCTGCATGACGAAGCGCATCTGCTGTTGc                                                                                                            <  5:171815/148‑1 (MQ=255)
                                           ccTGGTGCTGGAGCAAGGCCGGCCAGCCTGTATTCAGGTGTTCTACCGCCTGCAGGGCGGCTGGGCCGACGTGTATGTGCTGGATGAATACAACGCATTGTGGCAACAACGCCTGCCGCTGCATGACGAAGCGCATCTGCTGTTGCCAc                                                                                                         <  4:4477/149‑1 (MQ=255)
                                           ccTGGTGCTGGAGCAAGGCCGGCCAGCCTGTATTCAGGTGTTCTACCGCCTGCAGGGCGGCTGGGCCGACGTGTATGTGCTGGATGAATACAACGCATTGTGGCAACAACGCCTGCCGCTGCATGACGAAGCGCATCTGCTGTTGCCAc                                                                                                         >  7:233532/1‑149 (MQ=255)
                                            cTGGTGCTGGAGCAAGGCCGGCCAGCCTGTATTCAGGTGTTCTACCGCCTGCAGGGCGGCTGGGCCGACGTGTATGTGCTGGATGAATACAACGCATTGTGGcaac                                                                                                                                                   <  7:115440/106‑1 (MQ=255)
                                            cTGGTGCTGGAGCAAGGCCGGCCAGCCTGTATTCAGGTGTTCTACCGCCTGCAGGGCGGCTGGGCCGACGTGTATGTGCTGGATGAATACAACGCATTGTGGcaac                                                                                                                                                   >  8:115440/1‑106 (MQ=255)
                                               gTGCTGGAGCAAGGCCGGCCAGCCTGTATTCAGGTGTTCTACCGCCTGCAGGGCGGCTGGGCCGACGTGTATGTGCTGGATGAATACAACGCATTGTGGCAACAACGCCTGCCGCTGCATGACGAAGCGCATCTGCTGTTGCCACTg                                                                                                       <  6:299963/147‑1 (MQ=255)
                                               gTGCTGGAGCAAGGCCGGCCAGCCTGTATTCAGGTGTTCTACCGCCTGCAGGGCGGCTGGGCCGACGTGTATGTGCTGGATGAATACAACGCATTGTGGCAACAACGCCTGCCGCTGCATGACGAAGCGCATCTGCTGTTGCCACTGc                                                                                                      <  8:233532/148‑1 (MQ=255)
                                                  cTGGAGCAAGGCCGGCCAGCCTGTATTCAGGTGTTCTACCGCCTGCAGGGCGGCTGGGCCGACGTGTATGTGCTGGATGAATACAACGCATTGTGGCAACAACGCCTGCCGCTGCATGACGAAGCGCATCTGCTGTTGCCACTGCAAc                                                                                                   <  5:139613/148‑1 (MQ=255)
                                                             ccggccAGCCTGTATTCAGGTGTTCTACCGCCTGCAGGGCGGCTGGGCCGACGTGTATGTGCTGGATGAATACAACGCATTGTGGCAACAACGCCTGCCGCTGCATGACGAAGCGCATCTGCTGTTGCCACTGCAACGCTTTCTGCGCt                                                                                       <  3:314378/149‑1 (MQ=255)
                                                                  cAGCCTGTATTCAGGTGTTCTACCGCCTGCAGGGCGGCTGGGCCGACGTGTATGTGCTGGATGAATACAACGCATTGTGGCAACAACGCCTGCCGCTGCATGACGAAGCGCATCTGCTGTTGCCACTGCAACGCTTTCTGCGCTCggtg                                                                                  <  1:264749/149‑1 (MQ=255)
                                                                   aGCCTGTTTTCTGGTGTTGTAGCGCCTGCAGGGCGGCTGGGCCGACGTGTATGTGCTGGATGAATAGAACGCATTGTGGCAACAACGTGTTCGGCTGCATGTCGAAGCGAATCTGCTGTTGCCACTGCAACGGTTTGTGCGCTGggtgg                                                                                 <  2:301515/149‑1 (MQ=255)
                                                                     ccTGTATTCAGGTGTTCTACCGCCTGCAGGGCGGCTGGGCCGACGTGTATGTGCCGGATGAATACAACGCATTGTGGCAACAACGCCTGCCGCTGCATGACGAAGCGCATCTGCTGTTGCCACTGCAACGCTTTCTGCGCTCGgtggtg                                                                               <  8:248542/149‑1 (MQ=255)
                                                                         tATTCAGGTGTTCTACCGCCTGCAGGGCGGCTGGGCCGACGTGTATGTGCTGGATGAATACAACGCATTGTGGCAACAACGCCTGCCGCTGCATGACGAAGCGCATCTGCTGTTGCCACTGCAACGCTTTCTGCGCTCGGTGGTGATg                                                                            <  1:172639/148‑1 (MQ=255)
                                                                            tCAGGTGTTCTACCGCCTGCAGGGCGGCTGGGCCGACGTGTATGTGCTGGATGAATACAACGCATTGTGGCAACAACGCCTGCCGCTGCATGACGAAGCGCATCTGCTGTTGCCACTGCAACGCTTTCt                                                                                            <  6:361219/129‑1 (MQ=255)
                                                                            tCAGGTGTTCTACCGCCTGCAGGGCGGCTGGGCCGACGTGTATGTGCTGGATGAATACAACGCATTGTGGCAACAACGCCTGCCGCTGCATGACGAAGCGCATCTGCTGTTGCCACTGCAACGCTTTCt                                                                                            >  5:361219/1‑129 (MQ=255)
                                                                                gtgtTCTACCGCCTGCAGGGCGGCTGGGCCGACGTGTATGTGCTGGATGAATACAACGCATTGTGGCAACAACGCCTGCCGCTGCATGACGAAGCGCATCTGCTGTTGCCACTGCAACGCTTTCTGCGCTCGGTGGTGATGCGCCGCGa                                                                    >  8:143949/1‑149 (MQ=255)
                                                                                     cTACCGCCTGCAGGGCGGCTGGGCCGACGTGTATGTGCTGGATGAATACAACGCATTGTGGCAACAACGCCTGCCGCTGCATGACGAAGCGCATCTGCTGTTGCCACTGCAACGCTTTCTGCGCTCGGTGGTGATGCGCCGCGATGCAc                                                               >  1:243455/1‑149 (MQ=255)
                                                                                      tACCGCCTGCAGGGCGGCTGGGCCGACGTGTATGTGCTGGATGAATACAACGCATTGTGGCAACATCGCCTGCCGCTGCATGACGAAGCGCATCTGCTGTTGCCACTGCaa                                                                                                    <  1:255677/111‑1 (MQ=255)
                                                                                      tACCGCCTGCAGGGCGGCTGGGCCGACGTGTATGTGCTGGATGAATACAACGCATTGTGGCAACAACGCCTGCCGCTGCATGACGAAGCGCATCTGCTGTTGCCACTGCaa                                                                                                    >  2:255677/1‑111 (MQ=255)
                                                                                       aCCGCCTGCAGGGCGGCTGGGCCGACGTGTATGTGCTGGATGAATACAACGCATTGTGGCAACAACGCCTGCCGCTGCATGACGAAGCGCATc                                                                                                                     <  3:98488/93‑1 (MQ=255)
                                                                                       aCCGCCTGCAGGGCGGCTGGGCCGACGTGTATGTGCTGGATGAATACAACGCATTGTGGCAACAACGCCTGCCGCTGCATGACGAAGCGCATc                                                                                                                     >  4:98488/1‑93 (MQ=255)
                                                                                       aCCGCCTGCAGGGCGGCTGGGCCGACGTGTATGTGCTGGATGAATACAACGCATTGTGGCAACAACGCCTGCCGCTGCATGACGAAGCGCATCTGCTGTTGCCACTGCAACGCTTTCTGCGCTCGGTGGTGATGAGCCGCGATGCACg                                                              >  2:343647/1‑148 (MQ=255)
                                                                                        ccGCCTGCAGGGCGGCTGGGCCGACGTGTATGTGCTGGATGAATACAACGCATTGTGGCAACAACGCCTGCCGCTGCATGACGAAGCGCATCTGCTGTTGCCACTGCAACGCTTTCTGCGCTCGGTGGTGATGCGCCGCGATGCAc                                                               <  2:243455/146‑1 (MQ=255)
                                                                                           ccTGCAGGGCGGCTGGGCCGACGTGTATGTGCTGGATGAATACAACGCATTGTGGCAACAACGCCTGCCGCTGCATGACGAAGCGCATCTGCTGTTGCCACTGCAACGCTTTCt                                                                                            >  5:38772/1‑114 (MQ=255)
                                                                                           ccTGCAGGGCGGCTGGGCCGACGTGTATGTGCTGGATGAATACAACGCATTGTGGCAACAACGCCTGCCGCTGCATGACGAAGCGCATCTGCTGTTGCCACTGCAACGCTTTCt                                                                                            <  6:38772/114‑1 (MQ=255)
                                                                                                 gggCGGCTGGGCCGACGTGTATGTGCTGGATGAATACAACGCATTGTGGCAACAACGCCTGCCGCTGCATGACGAAGCGCATCTGCTGTTGCCACTGCAACGCTTTCTGCGCTCGGTGGTGATGCGCCGCGATGCACGCCTGCCCCTc                                                    >  8:6403/1‑148 (MQ=255)
                                                                                                           gCCGACGTGTATGTGCTGGATGAATACAACGCATTGTGGCAACAACGCCTGCCGCTGCATGACGAAGCGCATCTGCTGTTGCCACTGCAACGCTTTCTGCGCTCGGTGGTGATGCGCCGCGATGCACGCCTGCCCCTCGATAGCCTGCg                                         <  8:249520/149‑1 (MQ=255)
                                                                                                               aCGTGTATGTGCTGGATGAATACAACGCATTGTGGCAACAACGCCTGCCGCTGCATGACGAAGCGCATCTGCTGTTGCCACTGCAACGCTTTCTGCGCTCGGTGGTGATGCGCCGCGATGCACGCCTGCCCCTCGATAGCCTGCGGccc                                     <  4:261232/149‑1 (MQ=255)
                                                                                                                cGTGTATGTGCTGGATGAATACAACGCATTGTGGCAACAACGCCTGCCGCTGCATGACGAAGCGCATCTGCTGTTGCCACTGCAACGCTTTCTGCGCTCGGTGGTGATGCGCCGCGATGCACGCCTGCCCCTCGATAGCCTGCGGCCCg                                    <  2:238875/149‑1 (MQ=255)
                                                                                                                     aTGTGCTGGATGAATACAACGCATTGTGGCAACAACGCCTGCCGCTGCATGACGAAGCGCATCTGCTGTTGCCACTGCAACGCTTTCTGCGCTCGGTGGTGATGCGCCGCGATGCACGc                                                             <  3:150462/119‑1 (MQ=255)
                                                                                                                     aTGTGCTGGATGAATACAACGCATTGTGGCAACAACGCCTGCCGCTGCATGACGAAGCGCATCTGCTGTTGCCACTGCAACGCTTTCTGCGCTCGGTGGTGATGCGCCGCGATGCACGc                                                             >  4:150462/1‑119 (MQ=255)
                                                                                                                     aTGTGCTGGATGAATACAACGCATTGTGGCAACAACGCCTGCCGCTGCATGACGAAGCGCATCTGCTGTTGCCACTGCAACGCTTTCTGCGCTCGGTGGTGATGCGCCGCGATGCACGCCTGCCCCTCGATAGCCTGCGGCCCGCCTCg                               >  3:2710/1‑149 (MQ=255)
                                                                                                                       gtgCTGGATGAATACAACGCATTGTGGCGACAACGCCTGCCGCTGCATGACGAAGCGCATCTGCTGTTGCCACTGCAACGCTTTCTGCGCTCGGTGGTGATGCGCCGCGATGCACGCCTGCCCCTCGATAGCCTGCGGCCCGCCTCGc                              >  8:184540/1‑148 (MQ=255)
                                                                                                                       gtgCTGGATGAATACAACGCATTGTGGCGACAACGCCTGCCGCTGCATGACGAAGCGCATCTGCTGTTGCCACTGCAACGCTTTCTGCGCTCGGTGGTGATGCGCCGCGATGCACGCCTGCCCCTCGATAGCCTGCGGCCCGCCTCGc                              >  8:189763/1‑148 (MQ=255)
                                                                                                                       gtgCTGGATGAATACAACGCATTGTGGCGACAACGCCTGCCGCTGCATGACGAAGCGCATCTGCTGTTGCCACTGCAACGCTTTCTGCGCTCGGTGGTGATGCGCCGCGATGCACGCCTGCCCCTCGATAGCCTGCGGCCCGCCTCGc                              <  7:184540/148‑1 (MQ=255)
                                                                                                                       gtgCTGGATGAATACAACGCATTGTGGCGACAACGCCTGCCGCTGCATGACGAAGCGCATCTGCTGTTGCCACTGCAACGCTTTCTGCGCTCGGTGGTGATGCGCCGCGATGCACGCCTGCCCCTCGATAGCCTGCGGCCCGCCTCGc                              <  7:189763/148‑1 (MQ=255)
                                                                                                                       gtgCTGGATGAATACAACGCATTGTGGCAACAACGCCTGCCGCTGCATGACGAAGCGCATCTGCTGTTGCCACTGCAACGCTTTCTGCGCTCGGTGGTGATGCGCCGCGATGCACGCCTGCCCCTCGATAGCCTGCGGCCCGCCTc                                <  7:337800/146‑1 (MQ=255)
                                                                                                                       gtgCTGGATGAATACAACGCATTGTGGCAACAACGCCTGCCGCTGCATGACGAAGCGCATCTGCTGTTGCCACTGCAACGCTTTCTGCGCTCGGTGGTGATGCGCCGCGATGCACGCCTGCCCCTCGATAGCCTGCGGCCCGCCTc                                >  8:337800/1‑146 (MQ=255)
                                                                                                                       gtgCTGGATGAATACAACGCATTGTGGCAACAACGCCTGCCGCTGCATGACGAAGCGCATCTGCTGTTGCCACTGCAACGCTTTCTGCGCTCGGTGGTGATGCGCCGCGATGCACGCCTGCCCCTCGATAGCCTGCGGCCCGCCTCGCt                             >  7:127737/1‑149 (MQ=255)
                                                                                                                        tgCTGGATGAATACAACGCATTGTGGCAACAACGCCTGCCGCTGCATGACGAAGCGCATCTGCTGTTGCCACTGCAACGCTTTCTGCGCTCGGTGGTGATGCGCCGCGATGCACGCCTGCCCCTCGATAGCCTGCGGCCCGCCTCGCt                             <  8:69068/148‑1 (MQ=255)
                                                                                                                                  aTACAACGCATTGTGGCAACAACGCCTGCCGCTGCATGACGAAGCGCATCTGCTGTTGCCACTGCAACGCTTTCTGCGCTCGGTGGTGATGCGCCGCGATGCACGCCTGCCCCTCGATAGCCTGCGGCCCGCCTCGCTGGATATCCACt                  >  4:350485/1‑149 (MQ=255)
                                                                                                                                  aTACAACGCATTGTGGCAACAACGCCTGCCGCTGCATGACGAAGCGCATCTGCTGTTGCCACTGCAACGCTTTCTGCGCTCGGTGGTGATGCGCCGCGATGCACGCCTGCCCCTCGATAGCCTGCGGCCCGCCTCGCTGGATATCCACt                  >  6:79340/1‑149 (MQ=255)
                                                                                                                                  aTACAACGCATTGTGGCAACAACGCCTGCCGCTGCATGACGAAGCGCATCTGCTGTTGCCACTGCAACGCTTTCTGCGCTCGGTGGTGATGCGCCGCGATGCACGCCTGCCCCTCGATAGCCTGCGGCCCGCCTCGCTGGATATCCACt                  >  8:239726/1‑149 (MQ=255)
                                                                                                                                   tACAACGCATTGTGGCAACAACGCCTGCCGCTGCATGACGAAGCGCATCTGCTGTTGCCACTGCAACGCTTTCTGCGCTCGGTGGTGATGCGCCGCGATGCACGCCTGCCCCTCGATAGCCTGCGGCCCACCTCGCTGGATATCCACtt                 >  2:121386/1‑149 (MQ=255)
                                                                                                                                      aaCGCATTGTGGCAACAACGCCTGCCGCTGCATGACGAAGCGCATCTGCTGTTGCCACTGCAACGCTTTCTGCGCTCGGTGGTGATGcgccg                                                                       <  2:294832/92‑1 (MQ=255)
                                                                                                                                      aaCGCATTGTGGCAACAACGCCTGCCGCTGCATGACGAAGCGCATCTGCTGTTGCCACTGCAACGCTTTCTGCGCTCGGTGGTGATGcgccg                                                                       >  1:294832/1‑92 (MQ=255)
                                                                                                                                      aaCGCATTGTGGCAACAACGCCTGCCGCTGCATGACGAAGCGCATCTGCTGTTGCCACTGCAACGCTTTCTGCGCTCGGTGGTGATGCGCCGCGATGCACGCCTGCCCCTCGATAGCCTGCGGCCCGCCTCGCTGGATATCCACTTcgc              <  7:143949/149‑1 (MQ=255)
                                                                                                                                      aaCGCATTGTGGCAACAACGCCTGCCGCTGCATGACGAAGCGCATCTGCTGTTGCCACTGCAACGCTTTCTGCGCTCGGTGGTGATGCGCCGCGATGCACGCCTGCCCCTCGATAGCCTGCGGCCCGCCTCGCTGGATATCCACTTcgc              >  6:333473/1‑149 (MQ=255)
                                                                                                                                           aTTGTGGCAACAACGCCTGCCGCTGCATGACGAAGCGCATCTGCTGTTGCCACTGCAACGCTTTCTGCGCTCGGTGGTGATGCGCCGCGATGCACGCCTGCCCCTCGATAGCCTGCGGCCCGCCTCGCTGGATATCCACTTCGCGCAGc         >  1:341261/1‑149 (MQ=255)
                                                                                                                                             tgtgGCAACAACGCCTGCCGCTGCATGACGAAGCGCATCTGCTGTTGCCACTGCAACGCTTTCTGCGCTCGGTGGTGATGCGCCGCGATGCACGCCTGCCCCTCGATAGCCTGCGGCCCGCCTCGCTGGATATCCACTTCGCGCAGCt        <  7:6403/148‑1 (MQ=255)
                                                                                                                                               tgGCAACAACGCCTGCCGCTGCATGACGAAGCGCATCTGCTGTTGCCACTGCAACGCTTTCTGCGCTCGGTGGTGATGCGCCGCGATGCACGCCTGCCCCTCGATAGCCTGCGGCCCGCCTCGCTGGATATCCACTTCGCGCAGCTGtt     >  2:242148/1‑149 (MQ=255)
                                                                                                                                                 gCAACAACGCCTGCCGCTTCATGACGAAGCGCATCTGCTGTTGCCACTGCAACTCTTTCTGCGCTCGGTTGTGATGCGCCGCGATGAACGCCTGCCCCTCGATAGCCTGCGGCCCGCCTCGCTGGATATCCACTTcgcg             <  3:300350/139‑1 (MQ=255)
                                                                                                                                                 gCAACAACGCCTGCCGCTGCATGACGAAGCGCATCTGCTGTTGCCACTGCAACGCTTTCTGCGCTCggtg                                                                                  <  3:82292/70‑1 (MQ=255)
                                                                                                                                                 gCAACAACGCCTGCCGCTGCATGACGAAGCGCATCTGCTGTTGCCACTGCAACGCTTTCTGCGCTCggtg                                                                                  >  4:82292/1‑70 (MQ=255)
                                                                                                                                                 gCAACAACGCCTGCCGCTGCATGACGAAGCGCATCTGCTGTTGCCACTGCAACGCTTTCTGCGCTCGGTGGTGATGCGCCGCGATGCACGCCTGCCCCTCGATAGCCTGCGGCCCGCCTCGCTGGATATCCACTTCGCGCAGCTGTTGc   <  5:79340/149‑1 (MQ=255)
                                                                                                                                                 gCAACAACGCCTGCCGCTGCATGACGAAGCGCATCTGCTGTTGCCACTGCAACGCTTTCTGCGCTCGGTGGTGATGCGCCGCGATGCACGCCTGCCCCTCGATAGCCTGCGGCCCGCCTCGCTGGATATCCACTTCGCGCAGCTGTTGc   <  4:2710/149‑1 (MQ=255)
                                                                                                                                                 gCAACAACGCCTGCCGCTGCATCACGAACCGCATCTGCTGTTGCCACTGCAACGCTTTCTGCGCTCGGTGGTGATGCGACGCGATGCACGCCTGCACCTCGATAGCATGCGGCACGCATAGCTGGAAATCCACTACCCGCacatcgca    >  4:300350/1‑141 (MQ=255)
                                                                                                                                                   aacaacGCCTGCCGCTGCATGACGAAGCGCATCTGCTGTTGCCACTGCAACGCTTTCTGCGCTCGGTGGTGATGCGCCGCGATGCACGCCTGCCCCTCGATAGCCTGCGGCCCGCCTCGCTGGATATCCACTTCGCGCAGCTGTTGcc  >  6:212019/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
CGCTGTACCTGGATGCCAATGCCCTGCAGGATCACGACCTGCGCCTGGTGCTGGAGCAAGGCCGGCCAGCCTGTATTCAGGTGTTCTACCGCCTGCAGGGCGGCTGGGCCGACGTGTATGTGCTGGATGAATACAACGCATTGTGGCAACAACGCCTGCCGCTGCATGACGAAGCGCATCTGCTGTTGCCACTGCAACGCTTTCTGCGCTCGGTGGTGATGCGCCGCGATGCACGCCTGCCCCTCGATAGCCTGCGGCCCGCCTCGCTGGATATCCACTTCGCGCAGCTGTTGCC  >  NC_002947/5958788‑5959082

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