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)
       gcgGTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAATCGcc                                                                                             <  2:2536883/100‑1 (MQ=255)
           tGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAATCGCctctc                                                                                         >  1:2361684/1‑100 (MQ=255)
              ggCGAAACTTATCTGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAATCGCCTCTCTgg                                                                                      <  2:2038617/100‑1 (MQ=255)
              ggCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAATCGCCTCTCTgg                                                                                      <  1:89152/100‑1 (MQ=255)
               gCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAATCGCCTCTCTGGa                                                                                     >  1:11168/1‑100 (MQ=255)
                cGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAATCGCCTCTCTGGAc                                                                                    <  1:1042806/100‑1 (MQ=255)
                cGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAATCGCCTCTCTGGAc                                                                                    >  2:831429/1‑100 (MQ=255)
                cGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAATCGCCTCTCTGGAc                                                                                    <  2:1537893/100‑1 (MQ=255)
                cGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAATCGCCTCTCTGGAc                                                                                    <  1:1706893/100‑1 (MQ=255)
                 gAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAATCGCCTCTCTGGACg                                                                                   <  1:608753/100‑1 (MQ=255)
                  aaaCTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGCGCTGGCGCTGGCGCAGGCTCGGCATCCAGCCAGGAAAAAGCCTCTCGGGACGa                                                                                  >  2:1705424/1‑100 (MQ=255)
                     cTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAATCGCCTCTCTGGACGAtct                                                                               >  1:322008/1‑100 (MQ=255)
                     cTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAATCGCCTCTCTGGACGAtct                                                                               >  2:1019938/1‑100 (MQ=255)
                      ttATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAATCGCCTCTCTGGACGAtctc                                                                              >  1:1296813/1‑100 (MQ=255)
                       tATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAATCGCCTCTCTGGACGATCTCg                                                                             >  1:1780896/1‑100 (MQ=255)
                       tATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAATCGCCTCTCTGGACGATCTCg                                                                             >  1:213489/1‑100 (MQ=255)
                          ccGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAATCGCCTCTCTGGGCGATCTCGAAg                                                                          >  2:1043558/1‑100 (MQ=255)
                           cGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAATCGCCTCTCTGGACGATCTCGAAGc                                                                         <  1:1808028/100‑1 (MQ=255)
                           cGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGCAAATCGCCTTCCTGGTCGATCTCGaatc                                                                         >  2:775636/1‑98 (MQ=255)
                            gCATTCGTGATAACGGCTGCGGTATCAAAACAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAATCGCCTCTCTGGACGATCTCGAAGcc                                                                        <  2:2061853/100‑1 (MQ=255)
                            gCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAATCGCCTCTCTGGACGATCTCGAAGcc                                                                        <  2:1876182/100‑1 (MQ=255)
                            gCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAATCGCCTCTCTGGACGATCTCGAAGcc                                                                        <  2:2061590/100‑1 (MQ=255)
                            gCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAATCGCCTCTCTGGACGATCTCGAAGcc                                                                        >  1:427790/1‑100 (MQ=255)
                              aTTCGTGATCACGGCTGCGGTATCAAAAAAGATGAGCTGGCGGTGGCGCTGGCTCGTCATGCCACCAGTGAATCGCCTCTCTGGGCGATCTCGAAGCCAt                                                                      <  1:2299748/100‑1 (MQ=255)
                                  gTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAATCGCCTCTCTGGACGATCTCGAAGCCATTATc                                                                  >  2:2469751/1‑100 (MQ=255)
                                  gTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAATCGCCTCTCTGGACGATCTCGAAGCCATTATc                                                                  >  1:1088166/1‑100 (MQ=255)
                                    gATAACGGCTGCGGTCTCAACAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAATCGCCTCTCTGGACGATAGCGCAGCCATTATCAg                                                                >  2:153668/1‑100 (MQ=255)
                                     aTAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGc                                                               <  1:731465/100‑1 (MQ=255)
                                         cGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTgg                                                           <  1:445880/100‑1 (MQ=255)
                                          ggCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAACTCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTggg                                                          <  2:2319949/100‑1 (MQ=255)
                                          ggCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTggg                                                          >  1:1143487/1‑100 (MQ=255)
                                          ggCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTggg                                                          <  1:1018789/100‑1 (MQ=255)
                                           gCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGc                                                         >  1:1043646/1‑100 (MQ=255)
                                           gCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAATCGCCTCTCTAGACGATCTCGAAGCCATTCTCAGCCTGGGc                                                         >  1:1315953/1‑100 (MQ=255)
                                            cTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAATCGCCACTATGGACGATCTCGACGCCATCAACAGCCTCGGca                                                        >  1:2639734/1‑99 (MQ=255)
                                              gCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCttt                                                      <  1:2159743/100‑1 (MQ=255)
                                                 gTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTcgc                                                   <  1:1519283/100‑1 (MQ=255)
                                                 gTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTcgc                                                   <  1:1265837/100‑1 (MQ=255)
                                                 gTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTcgc                                                   <  2:463075/100‑1 (MQ=255)
                                                   aTCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCgg                                                 <  2:1076093/100‑1 (MQ=255)
                                                   aTCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCgg                                                 <  2:322258/100‑1 (MQ=255)
                                                     cAAAAAAGCTGAGGTGGCGATGGCCCTGGCTTGTCATGCCACCAGTAAATCGCCTCTCGGGAGGAGCTCGTAGCCTGTATTAGCCTGTGCTCTCggggag                                               >  2:662128/1‑98 (MQ=255)
                                                     cAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTg                                               >  1:1352838/1‑100 (MQ=255)
                                                     cAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTg                                               >  1:1913874/1‑100 (MQ=255)
                                                      aaaaaaGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGa                                              <  1:926921/100‑1 (MQ=255)
                                                       aaaaaGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAg                                             <  1:1170365/100‑1 (MQ=255)
                                                        aaaaGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAgg                                            <  2:1904588/100‑1 (MQ=255)
                                                              tGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAATCGCCTCTCTGGGCGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTgg                                      >  1:69442/1‑100 (MQ=255)
                                                               gAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGc                                     >  1:1464734/1‑100 (MQ=255)
                                                               gAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGc                                     >  1:872820/1‑100 (MQ=255)
                                                                aGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCg                                    <  2:1710677/100‑1 (MQ=255)
                                                                aGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCg                                    <  1:46488/100‑1 (MQ=255)
                                                                  ctggcgctggcgctggcTCGTCATGCCACCAGTAAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAg                                  <  1:1560593/100‑1 (MQ=255)
                                                                  ctggcgctggcgctggcTCGTCATGCCACCAGTAAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAg                                  <  2:1393219/100‑1 (MQ=255)
                                                                    ggcgctggcgctggcTCGTCATGCCACCAGTAAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGTa                                >  2:1569475/1‑100 (MQ=255)
                                                                     gcggtggcgctggcGCGTCATGCCACCAGTAAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCCGTGAGGCGCTGGCGAGTAt                               <  2:2371250/100‑1 (MQ=255)
                                                                     gcgctggcggtggcTCGTCATGCCACCAGTAAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGTAt                               <  1:986193/100‑1 (MQ=255)
                                                                      cgctggcgctggcTCGTCATGCCACCAGTAAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGTATc                              <  2:1005108/100‑1 (MQ=255)
                                                                         tggcgctggcTCGTCATGCCACCAGTAAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGTATCAGt                           <  2:1071450/100‑1 (MQ=255)
                                                                            cgctggcTCGTCATGCCACCAGTAAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGTATCAGTTCg                        >  1:1733686/1‑100 (MQ=255)
                                                                             gctggcTCGTCATGCCACCAGTAAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGTATCAGTTCgg                       <  2:1310611/100‑1 (MQ=255)
                                                                                ggcTCGTCATGCCACCAGTAAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGTATCAGTTCGGttt                    >  1:776189/1‑100 (MQ=255)
                                                                                 gcTCGTCATGCCACCAGTAAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGTATCAGTTCGGTTTc                   >  1:1244957/1‑100 (MQ=255)
                                                                                 gcTCGTCATGCCACCAGTAAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGTATCAGTTCGGTTTc                   <  1:2128915/100‑1 (MQ=255)
                                                                                  cTCGTCATGCCACCAGTAAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGTATCAGTTCGGTTTcc                  <  2:399031/100‑1 (MQ=255)
                                                                                    cGTCATGCCACCAGTAAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGTATCAGTTCGGTTTCCCg                <  1:2047006/100‑1 (MQ=255)
                                                                                     gTCATGCCACCAGTAAATCGCCTCTCTGGACGATCTCGGAGCCATTATCATCCCGGGCATACCCGGGGGGGCGCCGGCGAATATCAGTTCGGTTTCCCGc               >  2:1604223/1‑100 (MQ=255)
                                                                                     gTCATGCCACCAGTAAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGTATCAGTTCGGTTTCCCGc               >  2:538866/1‑100 (MQ=255)
                                                                                     ctCATGCCACCAGTACATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCCGGGCTTCCGCGGTGAGGCGCTGGCGAGTATCAGTTCGGTTTCCCGc               <  2:1143672/99‑1 (MQ=255)
                                                                                      tCATGCCACCAGTAAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGTATCAGTTCGGTTTCCCGcc              <  2:2060234/100‑1 (MQ=255)
                                                                                      tCATGCCACCAGTAAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGTATCAGTTCGGTTTCCCGcc              <  1:1879017/100‑1 (MQ=255)
                                                                                      tCATGCCACCAGTAAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGTATCAGTTCGGTTTCCCGcc              <  1:716871/100‑1 (MQ=255)
                                                                                          gCCACCAGTAAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGTATCAGTTCGGTTTCCCGCCTGAc          <  1:1297151/100‑1 (MQ=255)
                                                                                          gCCACCAGTAAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGTATCAGTTCGGTTTCCCGCCTGAc          <  1:277407/100‑1 (MQ=255)
                                                                                          gCCACCAGTAAATCGCCTCTCTGGACGATCTCCAAGCCGTTATCAGGCTGGGCTTTCGCGGTGAGGCGCTGGCGAGTATCGGTTCGGTTTCCCGCCTGAc          >  1:2345903/1‑100 (MQ=255)
                                                                                           ccaccaGTAAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGTATCAGTTCGGTTTCCCGCCTGAcn         <  1:2104334/100‑2 (MQ=255)
                                                                                            caccaGTAAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGTATCAGTTCGGTTTCCCGCCTGACGc        <  1:1630530/100‑1 (MQ=255)
                                                                                             accaGTAAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGTATCAGTTCGGTTTCCCGCCTGACGCt       <  1:2353859/100‑1 (MQ=255)
                                                                                              ccaGTAAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGTATCAGTTCGGTTTCCCGCCTGACGCTc      <  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: