Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,368,134 A→G 9.5% H103R (CAC→CGC)  yhcG → DUF1016 family protein in the PD‑(D/E)XK nuclease superfamily

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,368,1340AG 9.5% 260.1 / 14.3 84H103R (CAC→CGC) yhcGDUF1016 family protein in the PD‑(D/E)XK nuclease superfamily
Reads supporting (aligned to +/- strand):  ref base A (37/39);  new base G (4/4);  total (41/43)
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.97e-01

CAAGGTGGCGAGGCCAGGGCTGCGTATGGTGCGCAGCTAATCAAGCGACTATCAAAGGATTTATGTCTAAGGTATAAGCGTGGGTTCTCTGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCT  >  NC_000913/3367986‑3368225
                                                                                                                                                    |                                                                                           
cAAGGTGGCGAGGCCAGGGCTGCGTATGGTGCGCAGCTAATCAAGCGACTATCAAAGGATTTATGTCTAAGGTATAAGCGTGGGTTCTCTGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCa                                                                                             <  2:540427/149‑1 (MQ=255)
   ggTGGCGAGGCCAGGGCTGCGTATGGTGCGCAGCTAATCAAGCGACTATCAAAGGATTTATGTCTAAGGTATAAGCGTGGGTTCTCTGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCa                                                                                          <  2:280922/149‑1 (MQ=255)
      ggCGAGGCCAGGGCTGCGTATGGTGCGCAGCTAATCAAGCGACTATCAAAGGATTTATGTCTAAGGTATAAGCGTGGGTTCTCTGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCAc                                                                                            <  2:1842466/144‑1 (MQ=255)
      ggCGAGGCCAGGGCTGCGTATGGTGCGCAGCTAATCAAGCGACTATCAAAGGATTTATGTCTAAGGTATAAGCGTGGGTTCTCTGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCAc                                                                                            >  1:1842466/1‑144 (MQ=255)
         gAGGCCAGGGCTGCGTATGGTGCGCAGCTAATCAAGCGACTATCAAAGGATTTATGTCTAAGGTATAAGCGTGGGTTCTCTGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAAt                                                                                    >  4:689796/1‑149 (MQ=255)
                   cTGCGTATGGTGCGCAGCTAATCAAGCGACTATCAAAGGATTTATGTCTAAGGTATAAGCGTGGGTTCTCTGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGaa                                                                          >  4:1729211/1‑149 (MQ=255)
                        tATGGTGCGCAGCTAATCAAGCGACTATCAAAGGATTTATGTCTAAGGTATAAGCGTGGGTTCTCTGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAac                                                                     <  3:232712/149‑1 (MQ=255)
                              gcgcAGCTAATCAAGCGACTATCAAAGGATTTATGTCTAAGGTATAGGCGTGGGTTCTCTGCAAAAAACTTACGCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTg                                                              <  4:1004062/149‑1 (MQ=255)
                                gcAGCTAATCAAGCGACTATCAAAGGATTTATGTCTAAGGTATAAGCGTGGGTTCTCTGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCAtt                                                               <  4:352327/147‑1 (MQ=255)
                                gcAGCTAATCAAGCGACTATCAAAGGATTTATGTCTAAGGTATAAGCGTGGGTTCTCTGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCAtt                                                               >  3:352327/1‑147 (MQ=255)
                                gcAGCTAATCAAGCGACTATCAAAGGATTTATGTCTAAGGTATAAGCGTGGGTTCTCTGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTgg                                                             <  3:1729211/149‑1 (MQ=255)
                                  aGCTAATCAAGCGACTATCAAAGGATTTATGTCTAAGGTATAAGCGTGGGTTCTCTGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAAtt                                                                        <  4:322673/136‑1 (MQ=255)
                                  aGCTAATCAAGCGACTATCAAAGGATTTATGTCTAAGGTATAAGCGTGGGTTCTCTGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAAtt                                                                        >  3:322673/1‑136 (MQ=255)
                                         cAAGCGACTATCAAAGGATTTATGTCTAAGGTATAAGCGTGGGTTCTCTGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGc                                                    <  3:689796/149‑1 (MQ=255)
                                          aaGCGACTATCAAAGGATTTATGTCTAAGGTATAAGCGTGGGTTCTCTGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCa                                                   >  4:465216/1‑149 (MQ=255)
                                            gCGACTATCAAAGGATTTATGTCTAAGGTATAAGCGTGGGTTCTCTGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCGCCAGACATTGTCTg                                                                               <  1:62766/119‑1 (MQ=255)
                                            gCGACTATCAAAGGATTTATGTCTAAGGTATAAGCGTGGGTTCTCTGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCGCCAGACATTGTCTg                                                                               <  3:570325/119‑1 (MQ=255)
                                            gCGACTATCAAAGGATTTATGTCTAAGGTATAAGCGTGGGTTCTCTGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCGCCAGACATTGTCTg                                                                               >  4:570325/1‑119 (MQ=255)
                                            gCGACTATCAAAGGATTTATGTCTAAGGTATAAGCGTGGGTTCTCTGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCGCCAGACATTGTCTg                                                                               >  2:62766/1‑119 (MQ=255)
                                              gACTATCAAAGGATTTATGTCTAAGGTATAAGCGTGGGTTCTCTGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCGCCAGACATTGTCTGGTGAAtt                                                                        <  4:657784/124‑1 (MQ=255)
                                              gACTATCAAAGGATTTATGTCTAAGGTATAAGCGTGGGTTCTCTGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCGCCAGACATTGTCTGGTGAAtt                                                                        >  3:657784/1‑124 (MQ=255)
                                              gACTATCAAAGGATTTATGTCTAAGGTATAAGCGTGGGTTCTCTGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGaa                                                                          >  1:2885/1‑122 (MQ=255)
                                              gACTATCAAAGGATTTATGTCTAAGGTATAAGCGTGGGTTCTCTGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGaa                                                                          <  2:2885/122‑1 (MQ=255)
                                                  aTCAAAGGATTTATGTCTAAGGTATAAGCGTGGGTTCTCTGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCTGCAGacggtgt                                           >  2:502616/1‑144 (MQ=255)
                                                     aaaGGATTTATGTCTAAGGTATAAGCGTGGGTTCTCTGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGaa                                                          <  4:646098/131‑1 (MQ=255)
                                                     aaaGGATTTATGTCTAAGGTATAAGCGTGGGTTCTCTGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGaa                                                          >  1:296207/1‑131 (MQ=255)
                                                     aaaGGATTTATGTCTAAGGTATAAGCGTGGGTTCTCTGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGaa                                                          >  3:646098/1‑131 (MQ=255)
                                                     aaaGGATTTATGTCTAAGGTATAAGCGTGGGTTCTCTGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGaa                                                          <  2:296207/131‑1 (MQ=255)
                                                     aaaGGATTTATGTCTAAGGTATAAGCGTGGGTTCTCTGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTg                                        >  1:1128153/1‑149 (MQ=255)
                                                      aaGGATTTATGTCTAAGGTATAAGCGTGGGTTCTCTGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGc                                       >  4:1767820/1‑149 (MQ=255)
                                                       aGGATTTATGTCTAAGGTATAAGCGTGGGTTCTCTGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCa                                      <  3:465216/149‑1 (MQ=255)
                                                                 tCTAAGGTATAAGCGTGGGTTCTCTGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACaa                                                                                     >  3:1687619/1‑92 (MQ=255)
                                                                 tCTAAGGTATAAGCGTGGGTTCTCTGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACaa                                                                                     <  4:1687619/92‑1 (MQ=255)
                                                                 tCTAAGGTATAAGCGTGGGTTCTCTGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGaa                                                          <  3:1067520/119‑1 (MQ=255)
                                                                 tCTAAGGTATAAGCGTGGGTTCTCTGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGaa                                                          >  4:1067520/1‑119 (MQ=255)
                                                                  cTAAGGTATAAGCGTGGGTTCTCTGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAAt                                                                         >  3:590818/1‑103 (MQ=255)
                                                                  cTAAGGTATAAGCGTGGGTTCTCTGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAAt                                                                         <  4:590818/103‑1 (MQ=255)
                                                                   tAAGGTATAAGCGTGGGTTCTCTGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCa                                      >  4:303004/1‑137 (MQ=255)
                                                                   tAAGGTATAAGCGTGGGTTCTCTGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCa                                      <  3:303004/137‑1 (MQ=255)
                                                                   tAAGGTATAAGCGTGGGTTCTCTGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGaa                                   >  1:1403398/1‑140 (MQ=255)
                                                                   tAAGGTATAAGCGTGGGTTCTCTGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGaa                                   <  2:1403398/140‑1 (MQ=255)
                                                                        tataAGCGTGGGTTCTCTGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCa                                                                                          <  3:1007248/80‑1 (MQ=255)
                                                                        tataAGCGTGGGTTCTCTGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCa                                                                                          >  4:1007248/1‑80 (MQ=255)
                                                                          taAGCGTGGGTTCTCTGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTg                                                              <  1:1141807/106‑1 (MQ=255)
                                                                          taAGCGTGGGTTCTCTGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTg                                                              >  2:1141807/1‑106 (MQ=255)
                                                                          taAGCGTGGGTTCTCTGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGaa                                   <  2:254224/133‑1 (MQ=255)
                                                                          taAGCGTGGGTTCTCTGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGaa                                   >  1:254224/1‑133 (MQ=255)
                                                                            aGCGTGGGTTCTCTGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACATTGTCTGGTGaa                                                                          <  1:767550/92‑1 (MQ=255)
                                                                            aGCGTGGGTTCTCTGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACATTGTCTGGTGaa                                                                          >  2:767550/1‑92 (MQ=255)
                                                                               gTGGGTTCTCTGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCGAAGATTTgg              <  3:1767820/149‑1 (MQ=255)
                                                                                    ttCTCTGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCTGCAGACGGTGTCTGCGGaa                                   <  4:467134/123‑1 (MQ=255)
                                                                                    ttCTCTGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCTGCAGACGGTGTCTGCGGaa                                   >  3:467134/1‑123 (MQ=255)
                                                                                       tctGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTc                                >  4:1193868/1‑123 (MQ=255)
                                                                                       tctGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTc                                <  3:1193868/123‑1 (MQ=255)
                                                                                        ctGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGGAAAGAt                  <  2:1630950/136‑1 (MQ=255)
                                                                                        ctGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGGAAAGAt                  >  1:1630950/1‑136 (MQ=255)
                                                                                        ctGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTgg              <  3:506376/140‑1 (MQ=255)
                                                                                        ctGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTgg              >  4:506376/1‑140 (MQ=255)
                                                                                         tGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGt    <  2:1128153/149‑1 (MQ=255)
                                                                                            aaaaaaCTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCAc                                                                                            <  4:301156/58‑1 (MQ=255)
                                                                                            aaaaaaCTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCAc                                                                                            >  3:301156/1‑58 (MQ=255)
                                                                                            aaaaaaCTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGt    <  1:125528/146‑1 (MQ=255)
                                                                                            aaaaaaCTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGt    >  2:125528/1‑146 (MQ=255)
                                                                                             aaaaaCTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGacac                                              <  1:1847075/103‑1 (MQ=255)
                                                                                             aaaaaCTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGacac                                              >  2:1847075/1‑103 (MQ=255)
                                                                                                aaCTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTgg              >  4:961580/1‑132 (MQ=255)
                                                                                                aaCTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTgg              <  3:961580/132‑1 (MQ=255)
                                                                                                aaCTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGAcagc      >  4:1870267/1‑137 (MQ=255)
                                                                                                aaCTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGAcagc      <  3:1870267/140‑4 (MQ=255)
                                                                                                        ccAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCg                                                     <  4:1858244/85‑1 (MQ=255)
                                                                                                        ccAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCg                                                     >  3:1858244/1‑85 (MQ=255)
                                                                                                         cAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACaa                                                                                     <  3:468702/52‑1 (MQ=255)
                                                                                                         cAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACaa                                                                                     >  4:468702/1‑52 (MQ=255)
                                                                                                                   ttttttACCTCTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGa                                    >  2:359897/1‑90 (MQ=255)
                                                                                                                   ttttttACCTCTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGa                                    <  1:359897/90‑1 (MQ=255)
                                                                                                                            tctTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTc   >  3:1664960/1‑115 (MQ=255)
                                                                                                                            tctTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTc   <  4:1664960/115‑1 (MQ=255)
                                                                                                                                           ttGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCt  <  2:502193/101‑1 (MQ=255)
                                                                                                                                           ttGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCt  <  4:71529/101‑1 (MQ=255)
                                                                                                                                           ttGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCt  >  1:502193/1‑101 (MQ=255)
                                                                                                                                           ttGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCt  >  3:71529/1‑101 (MQ=255)
                                                                                                                                           tgttaaTTCGCCAGACATTTTCTGGTGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTg     >  2:1259284/5‑98 (MQ=38)
                                                                                                                                           tgttaaTTCGCCAGACATTGTCTGGTGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTg     <  1:1259284/94‑1 (MQ=38)
                                                                                                                                                 ttCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGa          <  2:1150205/87‑1 (MQ=255)
                                                                                                                                                 ttCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGa          >  1:1150205/1‑87 (MQ=255)
                                                                                                                                                    |                                                                                           
CAAGGTGGCGAGGCCAGGGCTGCGTATGGTGCGCAGCTAATCAAGCGACTATCAAAGGATTTATGTCTAAGGTATAAGCGTGGGTTCTCTGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCT  >  NC_000913/3367986‑3368225

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