Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,397,649 Δ1 bp coding (238/1848 nt) mutL → methyl‑directed mismatch repair protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,397,6460A.98.7% 359.5 / ‑0.2 79coding (235/1848 nt)mutLmethyl‑directed mismatch repair protein
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base . (30/47);  total (32/47)
Fisher's exact test for biased strand distribution p-value = 3.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATCGAACGCGGTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGTATCAGTTCGGTTTCCCGCCTGACGCTCACTT  >  NC_000913/4397547‑4397745
                                                                                                   |                                                                                                   
atCGAACGCGGTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTa                                                                                                     >  2:1365834/1‑100 (MQ=255)
atCGAACGCGGTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTa                                                                                                     >  1:460337/1‑100 (MQ=255)
atCGAACGCGGTGGGGCGAAACTTATCCGCATTCGTAATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGGCCTGCCACCAGTa                                                                                                     >  2:1262303/1‑100 (MQ=255)
atCGAACGCGGCGGGGCGAAACTTATCCGCATTCGTGATATCGGCTGCGGTATCAAAAAAGATGAGCTGGCGGTGGCGCTGGCTCGGCATGCCACCAGTa                                                                                                     <  1:2329576/100‑1 (MQ=255)
 tCGAACGCGGTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTaa                                                                                                    >  2:788245/1‑100 (MQ=255)
 tCGAACGCGGTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTaa                                                                                                    >  2:1148377/1‑100 (MQ=255)
      cgcgGTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTaaatcgc                                                                                               <  2:1051891/100‑5 (MQ=255)
       gcgGTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGT‑AAATCGcc                                                                                             <  2:2536883/100‑1 (MQ=255)
           tGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGT‑AAATCGCctctc                                                                                         >  1:2361684/1‑100 (MQ=255)
              ggCGAAACTTATCTGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGT‑AAATCGCCTCTCTgg                                                                                      <  2:2038617/100‑1 (MQ=255)
              ggCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGT‑AAATCGCCTCTCTgg                                                                                      <  1:89152/100‑1 (MQ=255)
               gCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGT‑AAATCGCCTCTCTGGa                                                                                     >  1:11168/1‑100 (MQ=255)
                cGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGT‑AAATCGCCTCTCTGGAc                                                                                    <  1:1042806/100‑1 (MQ=255)
                cGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGT‑AAATCGCCTCTCTGGAc                                                                                    >  2:831429/1‑100 (MQ=255)
                cGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGT‑AAATCGCCTCTCTGGAc                                                                                    <  2:1537893/100‑1 (MQ=255)
                cGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGT‑AAATCGCCTCTCTGGAc                                                                                    <  1:1706893/100‑1 (MQ=255)
                 gAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGT‑AAATCGCCTCTCTGGACg                                                                                   <  1:608753/100‑1 (MQ=255)
                  aaaCTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGCGCTGGCGCTGGCGCAGGCTCGGCATCCAGCCAGGAAAA‑AGCCTCTCGGGACGa                                                                                  >  2:1705424/1‑100 (MQ=255)
                     cTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGT‑AAATCGCCTCTCTGGACGAtct                                                                               >  1:322008/1‑100 (MQ=255)
                     cTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGT‑AAATCGCCTCTCTGGACGAtct                                                                               >  2:1019938/1‑100 (MQ=255)
                      ttATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGT‑AAATCGCCTCTCTGGACGAtctc                                                                              >  1:1296813/1‑100 (MQ=255)
                       tATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGT‑AAATCGCCTCTCTGGACGATCTCg                                                                             >  1:1780896/1‑100 (MQ=255)
                       tATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGT‑AAATCGCCTCTCTGGACGATCTCg                                                                             >  1:213489/1‑100 (MQ=255)
                          ccGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGT‑AAATCGCCTCTCTGGGCGATCTCGAAg                                                                          >  2:1043558/1‑100 (MQ=255)
                           cGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGT‑AAATCGCCTCTCTGGACGATCTCGAAGc                                                                         <  1:1808028/100‑1 (MQ=255)
                           cGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAG‑CAAATCGCCTTCCTGGTCGATCTCGaatc                                                                         >  2:775636/1‑98 (MQ=255)
                            gCATTCGTGATAACGGCTGCGGTATCAAAACAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGT‑AAATCGCCTCTCTGGACGATCTCGAAGcc                                                                        <  2:2061853/100‑1 (MQ=255)
                            gCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGT‑AAATCGCCTCTCTGGACGATCTCGAAGcc                                                                        <  2:1876182/100‑1 (MQ=255)
                            gCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGT‑AAATCGCCTCTCTGGACGATCTCGAAGcc                                                                        <  2:2061590/100‑1 (MQ=255)
                            gCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGT‑AAATCGCCTCTCTGGACGATCTCGAAGcc                                                                        >  1:427790/1‑100 (MQ=255)
                              aTTCGTGATCACGGCTGCGGTATCAAAAAAGATGAGCTGGCGGTGGCGCTGGCTCGTCATGCCACCAGT‑GAATCGCCTCTCTGGGCGATCTCGAAGCCAt                                                                      <  1:2299748/100‑1 (MQ=255)
                                  gTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGT‑AAATCGCCTCTCTGGACGATCTCGAAGCCATTATc                                                                  >  2:2469751/1‑100 (MQ=255)
                                  gTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGT‑AAATCGCCTCTCTGGACGATCTCGAAGCCATTATc                                                                  >  1:1088166/1‑100 (MQ=255)
                                    gATAACGGCTGCGGTCTCAACAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGT‑AAATCGCCTCTCTGGACGATAGCGCAGCCATTATCAg                                                                >  2:153668/1‑100 (MQ=255)
                                     aTAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGT‑AAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGc                                                               <  1:731465/100‑1 (MQ=255)
                                         cGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGT‑AAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTgg                                                           <  1:445880/100‑1 (MQ=255)
                                          ggCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGT‑AACTCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTggg                                                          <  2:2319949/100‑1 (MQ=255)
                                          ggCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGT‑AAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTggg                                                          >  1:1143487/1‑100 (MQ=255)
                                          ggCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGT‑AAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTggg                                                          <  1:1018789/100‑1 (MQ=255)
                                           gCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGT‑AAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGc                                                         >  1:1043646/1‑100 (MQ=255)
                                           gCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGT‑AAATCGCCTCTCTAGACGATCTCGAAGCCATTCTCAGCCTGGGc                                                         >  1:1315953/1‑100 (MQ=255)
                                            cTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGT‑AAATCGCCACTATGGACGATCTCGACGCCATCAACAGCCTCGGca                                                        >  1:2639734/1‑99 (MQ=255)
                                              gCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGT‑AAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCttt                                                      <  1:2159743/100‑1 (MQ=255)
                                                 gTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGT‑AAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTcgc                                                   <  1:1519283/100‑1 (MQ=255)
                                                 gTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGT‑AAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTcgc                                                   <  1:1265837/100‑1 (MQ=255)
                                                 gTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGT‑AAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTcgc                                                   <  2:463075/100‑1 (MQ=255)
                                                   aTCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGT‑AAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCgg                                                 <  2:1076093/100‑1 (MQ=255)
                                                   aTCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGT‑AAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCgg                                                 <  2:322258/100‑1 (MQ=255)
                                                     cAAAAAAGCTGAGGTGGCGATGGCCCTGGCTTGTCATGCCACCAGT‑AAATCGCCTCTCGGGAGGAGCTCGTAGCCTGTATTAGCCTGTGCTCTCggggag                                               >  2:662128/1‑98 (MQ=255)
                                                     cAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGT‑AAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTg                                               >  1:1352838/1‑100 (MQ=255)
                                                     cAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGT‑AAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTg                                               >  1:1913874/1‑100 (MQ=255)
                                                      aaaaaaGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGT‑AAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGa                                              <  1:926921/100‑1 (MQ=255)
                                                       aaaaaGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGT‑AAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAg                                             <  1:1170365/100‑1 (MQ=255)
                                                        aaaaGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGT‑AAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAgg                                            <  2:1904588/100‑1 (MQ=255)
                                                              tGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGT‑AAATCGCCTCTCTGGGCGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTgg                                      >  1:69442/1‑100 (MQ=255)
                                                               gAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGT‑AAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGc                                     >  1:1464734/1‑100 (MQ=255)
                                                               gAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGT‑AAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGc                                     >  1:872820/1‑100 (MQ=255)
                                                                aGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGT‑AAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCg                                    <  2:1710677/100‑1 (MQ=255)
                                                                aGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGT‑AAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCg                                    <  1:46488/100‑1 (MQ=255)
                                                                  ctggcgctggcgctggcTCGTCATGCCACCAGT‑AAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAg                                  <  1:1560593/100‑1 (MQ=255)
                                                                  ctggcgctggcgctggcTCGTCATGCCACCAGT‑AAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAg                                  <  2:1393219/100‑1 (MQ=255)
                                                                    ggcgctggcgctggcTCGTCATGCCACCAGT‑AAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGTa                                >  2:1569475/1‑100 (MQ=255)
                                                                     gcggtggcgctggcGCGTCATGCCACCAGT‑AAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCCGTGAGGCGCTGGCGAGTAt                               <  2:2371250/100‑1 (MQ=255)
                                                                     gcgctggcggtggcTCGTCATGCCACCAGT‑AAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGTAt                               <  1:986193/100‑1 (MQ=255)
                                                                      cgctggcgctggcTCGTCATGCCACCAGT‑AAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGTATc                              <  2:1005108/100‑1 (MQ=255)
                                                                         tggcgctggcTCGTCATGCCACCAGT‑AAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGTATCAGt                           <  2:1071450/100‑1 (MQ=255)
                                                                            cgctggcTCGTCATGCCACCAGT‑AAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGTATCAGTTCg                        >  1:1733686/1‑100 (MQ=255)
                                                                             gctggcTCGTCATGCCACCAGT‑AAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGTATCAGTTCgg                       <  2:1310611/100‑1 (MQ=255)
                                                                                ggcTCGTCATGCCACCAGT‑AAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGTATCAGTTCGGttt                    >  1:776189/1‑100 (MQ=255)
                                                                                 gcTCGTCATGCCACCAGT‑AAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGTATCAGTTCGGTTTc                   >  1:1244957/1‑100 (MQ=255)
                                                                                 gcTCGTCATGCCACCAGT‑AAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGTATCAGTTCGGTTTc                   <  1:2128915/100‑1 (MQ=255)
                                                                                  cTCGTCATGCCACCAGT‑AAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGTATCAGTTCGGTTTcc                  <  2:399031/100‑1 (MQ=255)
                                                                                    cGTCATGCCACCAGT‑AAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGTATCAGTTCGGTTTCCCg                <  1:2047006/100‑1 (MQ=255)
                                                                                     gTCATGCCACCAGT‑AAATCGCCTCTCTGGACGATCTCGGAGCCATTATCATCCCGGGCATACCCGGGGGGGCGCCGGCGAATATCAGTTCGGTTTCCCGc               >  2:1604223/1‑100 (MQ=255)
                                                                                     gTCATGCCACCAGT‑AAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGTATCAGTTCGGTTTCCCGc               >  2:538866/1‑100 (MQ=255)
                                                                                     ctCATGCCACCAGT‑ACATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCCGGGCTTCCGCGGTGAGGCGCTGGCGAGTATCAGTTCGGTTTCCCGc               <  2:1143672/99‑1 (MQ=255)
                                                                                      tCATGCCACCAGT‑AAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGTATCAGTTCGGTTTCCCGcc              <  2:2060234/100‑1 (MQ=255)
                                                                                      tCATGCCACCAGT‑AAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGTATCAGTTCGGTTTCCCGcc              <  1:1879017/100‑1 (MQ=255)
                                                                                      tCATGCCACCAGT‑AAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGTATCAGTTCGGTTTCCCGcc              <  1:716871/100‑1 (MQ=255)
                                                                                          gCCACCAGT‑AAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGTATCAGTTCGGTTTCCCGCCTGAc          <  1:1297151/100‑1 (MQ=255)
                                                                                          gCCACCAGT‑AAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGTATCAGTTCGGTTTCCCGCCTGAc          <  1:277407/100‑1 (MQ=255)
                                                                                          gCCACCAGT‑AAATCGCCTCTCTGGACGATCTCCAAGCCGTTATCAGGCTGGGCTTTCGCGGTGAGGCGCTGGCGAGTATCGGTTCGGTTTCCCGCCTGAc          >  1:2345903/1‑100 (MQ=255)
                                                                                           ccaccaGT‑AAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGTATCAGTTCGGTTTCCCGCCTGAcn         <  1:2104334/100‑2 (MQ=255)
                                                                                            caccaGT‑AAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGTATCAGTTCGGTTTCCCGCCTGACGc        <  1:1630530/100‑1 (MQ=255)
                                                                                             accaGT‑AAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGTATCAGTTCGGTTTCCCGCCTGACGCt       <  1:2353859/100‑1 (MQ=255)
                                                                                              ccaGT‑AAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGTATCAGTTCGGTTTCCCGCCTGACGCTc      <  2:669423/100‑1 (MQ=255)
                                                                                                   aaaaTCGCCTATCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGTAACAGTTCGGTTTCCCGCCTGACGCTCACtt  >  2:2388660/1‑100 (MQ=255)
                                                                                                   |                                                                                                   
ATCGAACGCGGTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGTATCAGTTCGGTTTCCCGCCTGACGCTCACTT  >  NC_000913/4397547‑4397745

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑