Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,368,142 A→T 13.3% M106L (ATG→TTG)  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,1420AT13.3% 225.6 / 24.7 75M106L (ATG→TTG) yhcGDUF1016 family protein in the PD‑(D/E)XK nuclease superfamily
Reads supporting (aligned to +/- strand):  ref base A (32/33);  new base T (5/5);  total (37/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.18e-01

GAGGCCAGGGCTGCGTATGGTGCGCAGCTAATCAAGCGACTATCAAAGGATTTATGTCTAAGGTATAAGCGTGGGTTCTCTGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCT  >  NC_000913/3367995‑3368225
                                                                                                                                                   |                                                                                   
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                                                                               >  2:62766/1‑119 (MQ=255)
                                   gCGACTATCAAAGGATTTATGTCTAAGGTATAAGCGTGGGTTCTCTGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCGCCAGACATTGTCTg                                                                               >  4:570325/1‑119 (MQ=255)
                                   gCGACTATCAAAGGATTTATGTCTAAGGTATAAGCGTGGGTTCTCTGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCGCCAGACATTGTCTg                                                                               <  3:570325/119‑1 (MQ=255)
                                   gCGACTATCAAAGGATTTATGTCTAAGGTATAAGCGTGGGTTCTCTGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCGCCAGACATTGTCTg                                                                               <  1:62766/119‑1 (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)
                                                                 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                                                                          >  2:767550/1‑92 (MQ=255)
                                                                   aGCGTGGGTTCTCTGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACATTGTCTGGTGaa                                                                          <  1:767550/92‑1 (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              >  4:506376/1‑140 (MQ=255)
                                                                               ctGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTgg              <  3:506376/140‑1 (MQ=255)
                                                                                tGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGt    <  2:1128153/149‑1 (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      <  3:1870267/140‑4 (MQ=255)
                                                                                       aaCTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGAcagc      >  4:1870267/1‑137 (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                                    <  1:359897/90‑1 (MQ=255)
                                                                                                          ttttttACCTCTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGa                                    >  2:359897/1‑90 (MQ=255)
                                                                                                                   tctTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTc   <  4:1664960/115‑1 (MQ=255)
                                                                                                                   tctTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTc   >  3:1664960/1‑115 (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)
                                                                                                                                                aCAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTc                                          <  1:705161/47‑1 (MQ=255)
                                                                                                                                                aCAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTc                                          >  2:705161/1‑47 (MQ=255)
                                                                                                                                                   |                                                                                   
GAGGCCAGGGCTGCGTATGGTGCGCAGCTAATCAAGCGACTATCAAAGGATTTATGTCTAAGGTATAAGCGTGGGTTCTCTGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCT  >  NC_000913/3367995‑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: