Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,516,076 T→C 5.4% Q586R (CAA→CGA)  yfeA ← putative diguanylate cyclase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,516,0760TC 5.4% 255.1 / 3.4 74Q586R (CAA→CGA) yfeAputative diguanylate cyclase
Reads supporting (aligned to +/- strand):  ref base T (35/35);  new base C (2/2);  total (37/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 = 6.57e-01

AAACTTATGCAGCTGCTCGATGTTGTACATGCTGCTTTCTGCGTTAGAAAACGCCTGCTCCTCGGTGATCTCAAGAATGACCGCCTGCGGGGAGATGTGATAACGCTTAAACAGACGAATAATCCGCCCGGCAATATTCTTTTGCAGCAGCGTGAGCGGCATTAAATTGACTGAAAAGCGCGGTCCTTTTTTGTCGCAAGGGTGTGTTGCCAGCCACTTCAACAAGGATTCCAGCACTTGCAAATCAAAACGCGCGCTAAGGTTAAACTGAGCAATAAGGG  >  NC_000913/2515934‑2516214
                                                                                                                                              |                                                                                                                                          
aaaCTTATGCAGCTGCTCGATGTTGTACATGCTGCTTTCTGCGTTAGAAAACGCCTGCTCCTCGGTGATCTCAAGAATGACCGCCTGCGGGGAGATGTGATAACGCTTAAACAGACGAATAATCCGCCCGGCAATATTCTTTTgcagc                                                                                                                                       >  3:8334/1‑148 (MQ=255)
       tGCAGCTGCTCGATGTTGTACATGCTGCTTTCTGCGTTAGAAAACGCCTGCTCCTCGGTGATCTCAAGAATGACCGCCTGCGGGGAGATGTGATAACGCTTAAACAGACGAATAATCCGCCCGGCAATATTCTTTTGCAGCAGCGTGa                                                                                                                                >  2:291539/1‑148 (MQ=255)
             tgctCGATGTTGTACATGCTGCTTTCTGCGTTAGAAAACGCCTGCTCCTCGGTGATCTCAAGAATGACCGCCTGCGGGGAGATGTGATAACGCTTAAACAGACGAATAATCCGCCCGGCAATATTCTTTTGCAGCAGCGTGAGCGGCa                                                                                                                          <  4:8334/148‑1 (MQ=255)
                   aTGTTGTACATGCTGCTTTCTGCGTTAGAAAACGCCTGCTCCTCGGTGATCTCAAGAATGACCGCCTGCGGGGAGATGTGATAACGCTTAAACAGACGAATAATCCGCCCGGCAATATTCTTTTGCAGCAGCGTGAGCGGCATTAAAt                                                                                                                    <  3:299898/148‑1 (MQ=255)
                      ttgtACATGCTGCTTTCTGCGTTAGAAAACGCCTGCTCCTCGGTGATCTCAAGAATGACCGCCTGCGGGGAGATGTGATAACGCTTAAACAGACGAATAATCCGCCCGGCAATATTCTTTTGCAGCAGCGTGAGCGGCATTAAATTGa                                                                                                                 >  4:693619/1‑148 (MQ=255)
                         tACATGCTGCTTTCTGCGTTAGAAAACGCCTGCTCCTCGGTGATCTCAAGAATGACCGCCTGCGGGGAGATGTGATAACGCTTAAACAGACGAATAATCCGCCCGGCAATATTCTTTTGCAGCAGCGTGAGCGGCATTAAATTGACTg                                                                                                              >  1:930407/1‑148 (MQ=255)
                            aTGCTGCTTTCTGCGTTAGAAAACGCCTGCTCCTCGGTGATCTCAAGAATGACCGCCTGCGGGGAGATGTGATAACGCTTAAACAGACGAATAATCCGCCCGGCAATATTCTTTTGCAGCAGCGTGAGCGGCATTAAATTGACTGaaa                                                                                                           <  2:930407/148‑1 (MQ=255)
                                 gctTTCTGCGTTAGAAAACGCCTGCTCCTCGGTGATCTCAAGAATGACCGCCTGCGGGGAGATGTGATAACGCTTAAACAGACGAATAATCCGCCCGGCAATATTCTTTTGCAGCAGCGTGAGCGGCATTAAATTGACTGAAAAgcgc                                                                                                      <  1:291539/148‑1 (MQ=255)
                                           ttAGAAAACGCCTGCTCCTCGGTGATCTCAAGAATGACCGCCTGCGGGGAGATGTGATAACGCTTAAACAGACGAATAATCCGCCCGGCAATATTCTTTTGCAGCAGCGTGAGCGGCATTAAATTGACTGAAAAGCGCGGTCCttttt                                                                                            <  3:693619/148‑1 (MQ=255)
                                            tAGAAAACGCCTGCTCCTCGGTGATCTCAAGAATGACCGCCTGCGGGGAGATGTGATAACGCTTAAACAGACGAATAATCCGCCCGGCAATATTCTTTTGCAGCAGCGTGAGCGGCATTaaa                                                                                                                     <  4:671924/122‑1 (MQ=255)
                                            tAGAAAACGCCTGCTCCTCGGTGATCTCAAGAATGACCGCCTGCGGGGAGATGTGATAACGCTTAAACAGACGAATAATCCGCCCGGCAATATTCTTTTGCAGCAGCGTGAGCGGCATTaaa                                                                                                                     >  3:671924/1‑122 (MQ=255)
                                                aaaCGCCTGCTCCTCGGTGATCTCAAGAATGACCGCCTGCGGGGAGATGTGATAACGCTTAAACAGACGAATAATCCGCCCGGCAATATTCTTTTGCAGCAGCGTGAGCGGCATTAAATTGACTGAAAAGCGCGGTCCttttt                                                                                            >  1:336228/1‑143 (MQ=255)
                                                aaaCGCCTGCTCCTCGGTGATCTCAAGAATGACCGCCTGCGGGGAGATGTGATAACGCTTAAACAGACGAATAATCCGCCCGGCAATATTCTTTTGCAGCAGCGTGAGCGGCATTAAATTGACTGAAAAGCGCGGTCCttttt                                                                                            <  2:336228/143‑1 (MQ=255)
                                                           cctcGGTGATCTCAAGAATGACCGCCTGCGGGGAGATGTGATAACGCTTAAACAGACGAATAATCCGCCCGGCAATATTCTTTTg                                                                                                                                           >  1:437313/1‑85 (MQ=255)
                                                           cctcGGTGATCTCAAGAATGACCGCCTGCGGGGAGATGTGATAACGCTTAAACAGACGAATAATCCGCCCGGCAATATTCTTTTg                                                                                                                                           <  2:437313/85‑1 (MQ=255)
                                                           cctcGGTGATCTCAAGAATGACCGCCTGCGGGGAGATGTGATAACGCTTAAACAGACGAATAATCCGCCCGGCAATATTCTTTTGCAGCAGCGTGAGCGGCATTAAATTGa                                                                                                                 <  2:999560/111‑1 (MQ=255)
                                                           cctcGGTGATCTCAAGAATGACCGCCTGCGGGGAGATGTGATAACGCTTAAACAGACGAATAATCCGCCCGGCAATATTCTTTTGCAGCAGCGTGAGCGGCATTAAATTGa                                                                                                                 >  1:999560/1‑111 (MQ=255)
                                                             tcGGTGATCTCAAGAATGACCGCCTGCGGGGAGATGTGATAACGCTTAAACAGACGAATAATCCGCCCGGCAATATTCTTTTGCAGCAGCGTGAGCGGCATTAAATTGa                                                                                                                 <  4:86039/109‑1 (MQ=255)
                                                             tcGGTGATCTCAAGAATGACCGCCTGCGGGGAGATGTGATAACGCTTAAACAGACGAATAATCCGCCCGGCAATATTCTTTTGCAGCAGCGTGAGCGGCATTAAATTGa                                                                                                                 >  3:86039/1‑109 (MQ=255)
                                                              cGGTGATCTCAAGAATGACCGCCTGCGGGGAGATGTGATAACGCTTAAACAGACGAATAATCCGCCCGGCAATATTCTTTTGCAGCAGCGTGAGCGGCATTAAATTGACTGAAAAGCGCGGTCCTTTTTTGTCGCAAGGGTGTGTTgc                                                                         >  3:871718/1‑148 (MQ=255)
                                                                       cAAGAATGACCGCCTGCGGGGAGATGTGATAACGCTTAAACAGACGAATAATCCGCCCGGCAATATTCTTTTGAAGCAGCGTGAGCGGCATTAAATTGa                                                                                                                 >  4:1035516/1‑99 (MQ=255)
                                                                       cAAGAATGACCGCCTGCGGGGAGATGTGATAACGCTTAAACAGACGAATAATCCGCCCGGCAATATTCTTTTGAAGCAGCGTGAGCGGCATTAAATTGa                                                                                                                 <  3:1035516/99‑1 (MQ=255)
                                                                         aGAATGACCGCCTGCGGGGAGATGTGATAACGCTTAAACAGACGAATAATCCGCCCGGCAATATTCTTTTGCAGCAGCGTGAGCGGCATTAAATTGACTGAAAAGCGCGGTCCTTTTTTGTCGCAAGGGTGTGTTGCCAGCCACTTca                                                              >  4:561262/1‑148 (MQ=255)
                                                                         aGAATGACCGCCTGCGGGGAGATGTGATAACGCTTAAACAGACGAATAATCCGCCCGGCAATATTCTTTTGCAGCAGCGTGAGCGGCATTAAATTGACTGAAAAGCGCGGTCCTTTTTTGTCGCAAGGGTGTGTTGCCAGCCACTTca                                                              <  4:871718/148‑1 (MQ=255)
                                                                                  gCCTGCGGGGAGATGTGATAACGCTTAAACAGACGAATAATCCGCCCGGCAATATTCTTTTGCAGCAGCGTGAGCGGCATTAAATTGACTGAAAAGCGCGGTCCTTTTTTGTCGCAAGGGTGTGTTGCCAGCCACTTCAACAAGGAtt                                                     <  3:561262/148‑1 (MQ=255)
                                                                                          ggAGATGTGATAACGCTTAAACAGACGAATAATCCGCCCGGCAATATTCTTTTgc                                                                                                                                          >  4:487729/1‑55 (MQ=255)
                                                                                          ggAGATGTGATAACGCTTAAACAGACGAATAATCCGCCCGGCAATATTCTTTTgc                                                                                                                                          <  3:487729/55‑1 (MQ=255)
                                                                                              aTGTGATAACGCTTAAACAGACGAATAATCCGCCCGGCAATATTCTTTTGCAGCAGCGTGAGCGGCAt                                                                                                                         <  3:841826/68‑1 (MQ=255)
                                                                                              aTGTGATAACGCTTAAACAGACGAATAATCCGCCCGGCAATATTCTTTTGCAGCAGCGTGAGCGGCAt                                                                                                                         >  4:841826/1‑68 (MQ=255)
                                                                                                 tgATAACGCTTAAACAGACGAATAATCCGCCCGGCAATATTCTTTTGCAGCAGCGTGAGCGGCATTAAAt                                                                                                                    >  2:238980/1‑70 (MQ=255)
                                                                                                 tgATAACGCTTAAACAGACGAATAATCCGCCCGGCAATATTCTTTTGCAGCAGCGTGAGCGGCATTAAAt                                                                                                                    <  1:238980/70‑1 (MQ=255)
                                                                                                   aTAACGCTTAAACAGACGAATAATCCGCCCGGCAATATTCTTTTGCAGCAGCGTGAGCGGCATTAAATTGACTGAAAAGCGCGGTCCTTTTTTGTCGCAAGGGTGTGTTGCCAGCCACTTCAACAAGGATTCCAGCACTTGCAAATCa                                    >  4:964854/1‑148 (MQ=255)
                                                                                                     aaCGCTTAAACAGACGAATAATCCGCCCGGCAATATTCTTTTGCAGCAGCGTGAGCGGCATTAAATTGACTGAAAAGCGCGGTCCTTTTTTGTCGCAAGGGTGTGTTGCCAGCCACTTCAACAAGGATTCCAGc                                                >  3:1068674/1‑134 (MQ=255)
                                                                                                     aaCGCTTAAACAGACGAATAATCCGCCCGGCAATATTCTTTTGCAGCAGCGTGAGCGGCATTAAATTGACTGAAAAGCGCGGTCCTTTTTTGTCGCAAGGGTGTGTTGCCAGCCACTTCAACAAGGATTCCAGc                                                <  4:1068674/134‑1 (MQ=255)
                                                                                                      aCGCTTAAACAGACGAATAATCCGCCCGGCAATATTCTTTTGCAGCAGCGTGAGCGGCATTAAATTGACTGAAAAGCGCGGTCCTTTTTTGTCGCAAGGGTGTGTTGCCAGCCACTTCAACAAGGATTCCAGCACTTGCAAATCaaaa                                 >  3:622254/1‑148 (MQ=255)
                                                                                                           tAAACAGACGAATAATCCGCCCGGCAATATTCTTTTGCAGCAGCGTGAGCGGCATTAAATTGACTGAAAAGCGCGGTCCttt                                                                                              >  1:6660/1‑82 (MQ=255)
                                                                                                           tAAACAGACGAATAATCCGCCCGGCAATATTCTTTTGCAGCAGCGTGAGCGGCATTAAATTGACTGAAAAGCGCGGTCCttt                                                                                              <  2:6660/82‑1 (MQ=255)
                                                                                                           tAAACAGACGAATAATCCGCCCGGCAATATTCTTTTGCAGCAGCGTGAGCGGCATTAAATTGACTGAAAAGCGCGGTCCTTTTTTGTCGCAAggg                                                                                 >  4:387271/1‑95 (MQ=255)
                                                                                                           tAAACAGACGAATAATCCGCCCGGCAATATTCTTTTGCAGCAGCGTGAGCGGCATTAAATTGACTGAAAAGCGCGGTCCTTTTTTGTCGCAAggg                                                                                 <  3:387271/95‑1 (MQ=255)
                                                                                                             aaCAGACGAATAATCCGCCCGGCAATATTCTTTTGCAGCAGCGTGAGCGGCATTAAATTGACTGAAAAGCGCGGTCCTTTTTTGTCGCAAGGGTGTGTTGCCAGCCACTTCAACAAGGATTCCAGCACTTGCAAATCAAAACgcgcgc                          >  2:856984/1‑148 (MQ=255)
                                                                                                             aaCAGACGAATAATCCGCCCGGCAATATTCTTTTGCAGCAGCGTGAGCGGCATTAAATTGACTGAAAAGCGCGGTCCTTTTTTGTCGCAAGGGTGTGTTGCCAGCCACTTCAACAAGGATTCCAGCACTTGCAAATCAAAACgcgcgc                          <  3:964854/148‑1 (MQ=255)
                                                                                                               cAGACGAATAATCCGCCCGGCAATATTCTTTTGCAGCAGCGTGAGCGGCATTAAATTGACTGAAAAGCGCGGTCCTTTTTTGTCGCAAGGGTGTGTTGCCAGCCACTTCAACAAGGATTCCAGCACTTGCAAATCAAAACGCGCGCTa                        <  1:856984/148‑1 (MQ=255)
                                                                                                                 gACGAATAATCCGCCCGGCAATATTCTTTTGCAGCAGCGTGAGCGGCATTAAATTGACTGAAAAGCGCGGTCCTTTTTTGTCGCAAGGGTGTGTTGCCAGCCACTTCAACAAGGAt                                                      <  1:1027181/116‑1 (MQ=255)
                                                                                                                 gACGAATAATCCGCCCGGCAATATTCTTTTGCAGCAGCGTGAGCGGCATTAAATTGACTGAAAAGCGCGGTCCTTTTTTGTCGCAAGGGTGTGTTGCCAGCCACTTCAACAAGGAt                                                      >  2:1027181/1‑116 (MQ=255)
                                                                                                                 gACGAATAATCCGCCCGGCAATATTCTTTTGCAGCAGCGAGAGCGGCATTAAATTGACTGAAAAgcgc                                                                                                      >  1:328083/1‑68 (MQ=255)
                                                                                                                 gACGAATAATCCGCCCCGCAATATTCTTTTGCAGCAGCGTGAGCGGCATTAAATTGACTGAAAAgcgc                                                                                                      <  2:328083/68‑1 (MQ=255)
                                                                                                                      ataatCCGCCCGGCAATATTCTTTTGCAGCAGCGTGAGCGGCATTAAATTGACTGAAAAGCGCGGTCCTTTTTTGTCGCAAGGGTGTGTTGCCAGCCACTTCAACAAGGATTCCAGCACTTGCAAATCAAAACGCGCGCTAAGGTTaa                 >  1:988513/1‑148 (MQ=255)
                                                                                                                      ataatCCGCCCGGCAATATTCTTTTGCAGCAGCGTGAGCGGCATTAAATTGACTGAAAAGCGCGGTCCTTTTTTGTCGCAAGGGTGTGTTGCCAGCCACTTCAACAAGGATTCCAGCACTTGCAAATCAAAACGCGCGCTAAGGTTaa                 >  4:194754/1‑148 (MQ=255)
                                                                                                                       taatCCGCCCGGCAATATTCTTTTGCAGCAGCGTGAGCGGCATTAAATTGACTGAAAAGCGCGGTCCTTTTTTGTCGCAAGGGTGTGTTGCCAGCCACTTCAACAAGGATTCCAGCACTTGCAAATCAAAACGCGCGCTAAGGTTaaa                >  4:393188/1‑148 (MQ=255)
                                                                                                                        aatCCGCCCGGCAATATTCTTTTGCAGCAGCGTGAGCGGCATTAAATTGACTGAAAAGCGCGGTCCTTTTTTGTCGCAAGGGTGTGTTGCCAGCCACTTCAACAAGGGTTCCAGCa                                               >  3:548470/1‑116 (MQ=255)
                                                                                                                        aatCCGCCCGGCAATATTCTTTTGCAGCAGCGTGAGCGGCATTAAATTGACTGAAAAGCGCGGTCCTTTTTTGTCGCAAGGGTGTGTTGCCAGCCACTTCAACAAGGGTTCCAGCa                                               <  4:548470/116‑1 (MQ=255)
                                                                                                                          tCCGCCCGGCAATATTCTTTTGCAGCAGCGTGAGCGGCATTAAATTGACTGAAAAGCGCGGTCCTTTTTTGTCGCAAGGGTGTGTTGCCAGCCACTTcaacaa                                                          <  3:274321/103‑1 (MQ=255)
                                                                                                                          tCCGCCCGGCAATATTCTTTTGCAGCAGCGTGAGCGGCATTAAATTGACTGAAAAGCGCGGTCCTTTTTTGTCGCAAGGGTGTGTTGCCAGCCACTTcaacaa                                                          >  4:274321/1‑103 (MQ=255)
                                                                                                                          tCCGCCCGGCAATATTCTTTTGCAGCAGCGTGAGCGGCATTAAATTGACTGAAAAGCGCGGTCCTTTTTTGTCGCAAGGGTGTGTTGCCAGCCACTTCAACAAGGATTCCAGCACTTGCAAATCAAAACGCGCGCTAAGGTTAAACTg             <  3:194754/148‑1 (MQ=255)
                                                                                                                          tCCGCCCGGCAATATTCTTTTGCAGCAGCGTGAGCGGCATTAAATTGACTGAAAAGCGCGGTCCTTTTTTGTCGCAAGGGTGTGTTGCCAGCCACTTCAACAAGGATTCCAGCACTTGCAAATCAAAACGCGCGCTAAGGTCAAACTg             <  2:988513/148‑1 (MQ=255)
                                                                                                                           ccGCCCGGCAATATTCTTTTGCAGCAGCGTGAGCGGCATTAAATTGACTGAAAAGCGCGGTCCTTTTTTGTCGCAAGGGTGTGTTgccagcc                                                                    >  2:384639/1‑92 (MQ=255)
                                                                                                                           ccGCCCGGCAATATTCTTTTGCAGCAGCGTGAGCGGCATTAAATTGACTGAAAAGCGCGGTCCTTTTTTGTCGCAAGGGTGTGTTgccagcc                                                                    <  1:384639/92‑1 (MQ=255)
                                                                                                                           ccGCCCGGCAATATTCTTTTGCAGCAGCGTGAGCGGCATTAAATTGACTGAAAAGCGCGGTCCTTTTTTGTCGCAAGGGTGTGTTGCCAGCCACTTCAACAAGGATTCCAGCACtt                                            <  4:1013484/116‑1 (MQ=255)
                                                                                                                           ccGCCCGGCAATATTCTTTTGCAGCAGCGTGAGCGGCATTAAATTGACTGAAAAGCGCGGTCCTTTTTTGTCGCAAGGGTGTGTTGCCAGCCACTTCAACAAGGATTCCAGCACtt                                            >  3:250896/1‑116 (MQ=255)
                                                                                                                           ccGCCCGGCAATATTCTTTTGCAGCAGCGTGAGCGGCATTAAATTGACTGAAAAGCGCGGTCCTTTTTTGTCGCAAGGGTGTGTTGCCAGCCACTTCAACAAGGATTCCAGCACtt                                            >  3:1013484/1‑116 (MQ=255)
                                                                                                                           ccGCCCGGCAATATTCTTTTGCAGCAGCGTGAGCGGCATTAAATTGACTGAAAAGCGCGGTCCTTTTTTGTCGCAAGGGTGTGTTGCCAGCCACTTCAACAAGGATTCCAGCACtt                                            <  4:250896/116‑1 (MQ=255)
                                                                                                                           ccGCCCGGCAATATTCTTTTGCAGCAGCGTGAGCGGCATTAAATTGACTGAAAAGCGCGGTCCTTTTTTGTCGCAAGGGTGTGTTGCCAGCCACTTCAACAAGGATTCCAGCACTTGCAAATCAAAACGCGCGCTAAGGTTAAACTGa            >  1:614540/1‑148 (MQ=255)
                                                                                                                                     atatTCTTTTGCAGCAGCGTGAGCGGCATTAAATTGACTGAAAAGCGCGGTCCTTTTTTGTCGCAAGGGTGTGTTGCCAGCCACTTCAACAAGGATTCCAGCACTTGCAAATCAAAACGCGCGCTAAGGTTAAACTGAGCAATAAggg  <  3:393188/148‑1 (MQ=255)
                                                                                                                                     atatTCTTTTGCAGCAGCGTGAGCGGCATTAAATTGACTGAAAAGCGCGGTCCTTTTTTGTCGCAAGGGTGTGTTGCCAGCCACTTCAACAAGGATTCCAGCACTTGCAAATCAAAACGCGCGCTAAGGTTAAACTGAGCAATAAggg  <  2:614540/148‑1 (MQ=255)
                                                                                                                                     atatTCTTTCGCAGCAGCGTGAGCGGCATTAAATTGACTGAAAAGCGCGGTCCTTTTTTGTCGCAAGGGTGTGTTGCCAGCCACTTCAACAAGGATTCCAGCACTTGCAAATCAAAACGCGCGCTAAGGTTAAACTGAGCAATAAgg   >  2:307501/1‑147 (MQ=255)
                                                                                                                                     atatTCTTTCGCAGCAGCGTGAGCGGCATTAAATTGACTGAAAAGCGCGGTCCTTTTTTGTCGCAAGGGTGTGTTGCCAGCCACTTCAACAAGGATTCCAGCACTTGCAAATCAAAACGCGCGCTAAGGTTAAACTGAGCAATAAgg   <  1:307501/147‑1 (MQ=255)
                                                                                                                                     atatTCTTTCGCAGCAGCGTGAGCGGCATTAAATTGACTGAAAAGCGCGGTCCTTTTTTGTCGCAAGGGTGTGTTGCCAGCCACTTCAACAAGGATTCCAGCACTTGCAAATCAAAACGCGCGCTAAGGTTAAACTGAGCAATAAgg   >  2:727281/1‑147 (MQ=255)
                                                                                                                                     atatTCTTTCGCAGCAGCGTGAGCGGCATTAAATTGACTGAAAAGCGCGGTCCTTTTTTGTCGCAAGGGTGTGTTGCCAGCCACTTCAACAAGCATTCCAGCACTTGCAAATCAAAACGCGCGCTAAGGTTAAACTGAGCAATAAgg   <  1:727281/147‑1 (MQ=255)
                                                                                                                                       atTCTTTTGCAGCAGCGTGAGCGGCATTAAATTGACTGAAAAGCGCGGTCCTTTTTTGTCGCAAGGGTGTGTTGCCAGCCACTTCAACAAGGATTCCAGCACTTGCAAATCaa                                   <  4:478744/113‑1 (MQ=255)
                                                                                                                                       atTCTTTTGCAGCAGCGTGAGCGGCATTAAATTGACTGAAAAGCGCGGTCCTTTTTTGTCGCAAGGGTGTGTTGCCAGCCACTTCAACAAGGATTCCAGCACTTGCAAATCaa                                   >  3:478744/1‑113 (MQ=255)
                                                                                                                                         tCTTTTGCAGCAGCGTGAGCGGCATTAAATTGACTGAAAAGCGCGGTCCTTTTTTGTCGCAAGGGTGTGTTGCCAGCCACTTcaaca                                                           >  2:345893/1‑87 (MQ=255)
                                                                                                                                         tCTTTTGCAGCAGCGTGAGCGGCATTAAATTGACTGAAAAGCGCGGTCCTTTTTTGTCGCAAGGGTGTGTTGCCAGCCACTTcaaca                                                           <  1:345893/87‑1 (MQ=255)
                                                                                                                                         tCTTTTGCAGCAGCGTGAGCGGCATTAAATTGACTGAAAAGCGCGGTCCTTTTTTGTCGCAAGGGTGTGTTGCCAGCCACTTCAACAAGGATTCCAGCACTTGCAAATCAAAACGCGCGCTAAGGtt                   <  2:373703/127‑1 (MQ=255)
                                                                                                                                         tCTTTTGCAGCAGCGTGAGCGGCATTAAATTGACTGAAAAGCGCGGTCCTTTTTTGTCGCAAGGGTGTGTTGCCAGCCACTTCAACAAGGATTCCAGCACTTGCAAATCAAAACGCGCGCTAAGGtt                   >  1:373703/1‑127 (MQ=255)
                                                                                                                                           ttttGCAGCAGCGTGAGCGGCATTAAATTGACTGAAAAGCGCGGTCCTTTTTTGTCGCAAGGGTGTGTTGCCAGCCACTTCAACAAGGATTCCAGCACTTGCAAATCaaa                                  >  4:260860/1‑110 (MQ=255)
                                                                                                                                           ttttGCAGCAGCGTGAGCGGCATTAAATTGACTGAAAAGCGCGGTCCTTTTTTGTCGCAAGGGTGTGTTGCCAGCCACTTCAACAAGGATTCCAGCACTTGCAAATCaaa                                  <  3:260860/110‑1 (MQ=255)
                                                                                                                                            tttGCAGCAGCGTGAGCGGCATTAAATTGACTGAAAAGCGCGGTCCTTTTTTGTCGCAAGGGTGTGTTGCCAGCCACTTCAACAAGGATTCCAGCACTTGCaaa                                       >  1:166231/1‑104 (MQ=255)
                                                                                                                                            tttGCAGCAGCGTGAGCGGCATTAAATTGACTGAAAAGCGCGGTCCTTTTTTGTCGCAAGGGTGTGTTGCCAGCCACTTCAACAAGGATTCCAGCACTTGCaaa                                       <  2:166231/104‑1 (MQ=255)
                                                                                                                                            tttGCAGCAGCGTGAGCGGCATTAAATTGACTGAAAAGCGCGGTCCTTTTTTGTCGCAAGGGTGTGTTGCCAGCCACTTCAACAAGGATTCCAGCACTTGCAAATCAAAACGCGCGCTAAGGTTAAACTGAGCa         <  1:775976/134‑1 (MQ=255)
                                                                                                                                            tttGCAGCAGCGTGAGCGGCATTAAATTGACTGAAAAGCGCGGTCCTTTTTTGTCGCAAGGGTGTGTTGCCAGCCACTTCAACAAGGATTCCAGCACTTGCAAATCAAAACGCGCGCTAAGGTTAAACTGAGCa         >  2:775976/1‑134 (MQ=255)
                                                                                                                                              |                                                                                                                                          
AAACTTATGCAGCTGCTCGATGTTGTACATGCTGCTTTCTGCGTTAGAAAACGCCTGCTCCTCGGTGATCTCAAGAATGACCGCCTGCGGGGAGATGTGATAACGCTTAAACAGACGAATAATCCGCCCGGCAATATTCTTTTGCAGCAGCGTGAGCGGCATTAAATTGACTGAAAAGCGCGGTCCTTTTTTGTCGCAAGGGTGTGTTGCCAGCCACTTCAACAAGGATTCCAGCACTTGCAAATCAAAACGCGCGCTAAGGTTAAACTGAGCAATAAGGG  >  NC_000913/2515934‑2516214

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