Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,399,151 Δ1 bp coding (1740/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,399,1510G.99.3% 660.3 / ‑4.1 136coding (1740/1848 nt)mutLmethyl‑directed mismatch repair protein
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base . (65/70);  total (66/70)
Fisher's exact test for biased strand distribution p-value = 4.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 136 total aligned reads displayed.

CAACAAAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTGATTTACATCCGGCGATAAAAGCCCTGAAAGATGAGTGATATCAGTAAGGCGAGCCTGCCTAAGGCGATTTTTTT  >  NC_000913/4399005‑4399295
                                                                                                                                                  |                                                                                                                                                
  acaaAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTCTGGCgg                                                                                                                                           >  4:462230/1‑151 (MQ=255)
    aaaaTTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTCTGGCGGAc                                                                                                                                         <  1:1131376/151‑1 (MQ=255)
    aaaaTTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTCTGGCGGAc                                                                                                                                         <  2:80149/151‑1 (MQ=255)
     aaaTTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTCTGGCGGACg                                                                                                                                        <  4:363632/151‑1 (MQ=255)
        tttACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTCTGGCGGACGTgg                                                                                                                                     >  3:118220/1‑151 (MQ=255)
         ttACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTCTGGCGGACGTgg                                                                                                                                     >  2:1545922/1‑150 (MQ=255)
           aCAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTCTGGCGGACGTGGAAc                                                                                                                                  >  2:1433574/1‑151 (MQ=255)
             aaaTCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTCTGGCGGACGTGGAACgg                                                                                                                                <  2:511163/151‑1 (MQ=255)
               aTCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTCTGGCGGACGTGGAACGGtt                                                                                                                              >  2:357705/1‑151 (MQ=255)
                   tGATTCCTGAACTGATAGGCTGCCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTCTGGCGGACGTGGAACGGTTa                                                                                                                             <  1:672025/148‑1 (MQ=255)
                   tGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTCTGGCGGACGTGGAACGGTTa                                                                                                                             <  1:674999/148‑1 (MQ=255)
                   tGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTCTGGCGGACGTGGAACGGTTATg                                                                                                                           >  1:1108868/1‑150 (MQ=255)
                     aTTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTCTGGCGGACGTGGAACGGTTATGTcc                                                                                                                        >  1:1266522/1‑151 (MQ=255)
                     aTTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTCTGGCGGACGTGGAACGGTTATGTcc                                                                                                                        >  1:26377/1‑151 (MQ=255)
                        ccTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTCTGGCGGACGTGGAACGGTTATGTCCGc                                                                                                                      >  2:569682/1‑150 (MQ=255)
                         cTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTCTGGCGGACGTGGAACGGTTATGTCCGCa                                                                                                                     >  1:1027157/1‑150 (MQ=255)
                         cTGAACTGATAGGCTACATGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTCTGGCGGACGTGGAACGGTTATGTCCGCa                                                                                                                     >  1:1359205/1‑150 (MQ=255)
                          tGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTCTGGCGGACGTGGAACGGTTATGTCCGCa                                                                                                                     >  1:48634/1‑149 (MQ=255)
                          tGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTCTGGCGGACGTGGAACGGTTATGTCCGCAAc                                                                                                                   >  2:145459/1‑151 (MQ=255)
                           gAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTCTGGCGGACGTGGAACGGTTATGTCCGCaa                                                                                                                    <  1:1606745/149‑1 (MQ=255)
                           gAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTCTGGCGGACGTGGAACGGTTATGTCCGCAAc                                                                                                                   <  3:316624/150‑1 (MQ=255)
                           gAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTCTGGCGGACGTGGAACGGTTATGTCCGCAACt                                                                                                                  >  2:1172658/1‑151 (MQ=255)
                              cTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTCTGGCGGACGTGGAACGGTTATGTCCGCAACTtgt                                                                                                               >  1:1521860/1‑151 (MQ=255)
                                gATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTCTGGCGGACGTGGAACGGTTATGTCCGCAACTtgtg                                                                                                              >  1:627119/1‑150 (MQ=255)
                                gATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTCTGGCGGACGTGGAACGGTTATGTCCGCAACTtgtg                                                                                                              <  2:760675/150‑1 (MQ=255)
                                gATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTCTGGCGGACGTGGAACGGTTATGTCCGCAACTtgtg                                                                                                              >  2:156291/1‑150 (MQ=255)
                                gATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGa                                                                                                             <  2:950068/151‑1 (MQ=255)
                                   aGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGaaa                                                                                                           <  1:346897/150‑1 (MQ=255)
                                     gCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAAcg                                                                                                        <  1:146179/151‑1 (MQ=255)
                                        aCCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAAcgcc                                                                                                      >  1:1320166/1‑150 (MQ=255)
                                         ccTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAAcgccgc                                                                                                    >  2:426342/1‑151 (MQ=255)
                                            ggCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGccgccg                                                                                                  <  1:934914/150‑1 (MQ=255)
                                            ggCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCgg                                                                                                 <  2:242577/151‑1 (MQ=255)
                                             gCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCggg                                                                                                <  2:1438590/151‑1 (MQ=255)
                                              cGAAGCAGTCCGTATTCGAACCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCggg                                                                                                >  2:1121097/1‑150 (MQ=255)
                                                aaGCAGTCCGTATTCGAACCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGggtgg                                                                                             >  2:500194/1‑151 (MQ=255)
                                                aaGCAGTCCGTATTCGAACCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGggtgg                                                                                             <  1:666382/151‑1 (MQ=255)
                                                   cAGTCCGTATTCGAACCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCt                                                                                          <  1:1632514/151‑1 (MQ=255)
                                                   cAGTCCGTATTCGAACCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCAAAACCCTCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCt                                                                                          >  1:352626/1‑151 (MQ=255)
                                                     gTCCGTATTCGAACCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGt                                                                                        <  1:426426/151‑1 (MQ=255)
                                                       ccGTATTCGAACCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGACGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGCGAAAACGCCGCCGGGTGGTCTGtt                                                                                       >  4:563935/1‑150 (MQ=255)
                                                              cGAGCCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTATAATCt                                                                                <  1:1242507/150‑1 (MQ=255)
                                                              cGAACCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTg                                                                               >  2:1628470/1‑151 (MQ=255)
                                                              cGAACCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTCTGGCGGACGTGGAACGGTTATGTCAGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCt                                                                                >  1:1665633/1‑150 (MQ=255)
                                                                aaCCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGtt                                                                             >  2:952224/1‑151 (MQ=255)
                                                                  ccTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTg                                                                            >  1:897996/1‑150 (MQ=255)
                                                                  ccTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTg                                                                            >  3:6469/1‑150 (MQ=255)
                                                                  ccTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTg                                                                            <  1:106166/150‑1 (MQ=255)
                                                                   cTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTGAt                                                                          >  2:606600/1‑151 (MQ=255)
                                                                     ggCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTGAttt                                                                        <  3:462247/151‑1 (MQ=255)
                                                                     ggCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTGAttt                                                                        <  3:148506/151‑1 (MQ=255)
                                                                      gCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTGAttt                                                                        >  1:47205/1‑150 (MQ=255)
                                                                       cAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTGATTTAc                                                                      >  4:478443/1‑151 (MQ=255)
                                                                       cAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTGATTTAc                                                                      >  1:1361047/1‑151 (MQ=255)
                                                                        aaTATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTGATTTAc                                                                      >  4:62454/1‑150 (MQ=255)
                                                                        aaTATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTGATTTACa                                                                     <  2:971163/151‑1 (MQ=255)
                                                                        aaTATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTGATTTACa                                                                     <  1:165043/151‑1 (MQ=255)
                                                                         atatTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTGATTTACAt                                                                    <  2:905927/151‑1 (MQ=255)
                                                                         atatTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTGATTTACAt                                                                    >  1:772618/1‑151 (MQ=255)
                                                                           atTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTGATTTACATcc                                                                  >  1:596274/1‑151 (MQ=255)
                                                                           atTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTGATTTACATc                                                                   <  1:452016/150‑1 (MQ=255)
                                                                            ttGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTGATTTACATCCg                                                                 >  4:172654/1‑151 (MQ=255)
                                                                             tGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTGATTTACATCCgg                                                                >  2:1319199/1‑151 (MQ=255)
                                                                              gcgcAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTGATTTACATCCgg                                                                <  2:762149/150‑1 (MQ=255)
                                                                               cgcAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTGATTTACATCCGGCg                                                              >  2:621161/1‑151 (MQ=255)
                                                                                gcAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTGATTTACATCCGGCGa                                                             <  2:669349/151‑1 (MQ=255)
                                                                                 cAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTGATTTACATCCGGCGAt                                                            <  2:682855/151‑1 (MQ=255)
                                                                                 cAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTGATTTACATCCGGCGAt                                                            <  1:794351/151‑1 (MQ=255)
                                                                                  aGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTGATTTACATCCGGCGATa                                                           <  2:1697160/151‑1 (MQ=255)
                                                                                         tGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTGATTTACATCCGGCGATAAAAGcc                                                     <  2:1501587/150‑1 (MQ=255)
                                                                                         tGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTGATTTACATCCGGCGATAAAAGcc                                                     <  2:777169/150‑1 (MQ=255)
                                                                                            aCGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATACCCCTCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTGATTTACATc                                                                   <  2:410110/133‑1 (MQ=255)
                                                                                            aCGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTGATTTACATc                                                                   >  1:410046/1‑133 (MQ=255)
                                                                                            aCGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTGATTTACATCCGGCGATAAAAGCCCTg                                                  >  2:128491/1‑150 (MQ=255)
                                                                                              gAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTGATTTACATCCGGCGATAAAAGCCCTGaa                                                >  2:264174/1‑150 (MQ=255)
                                                                                              gAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTGATTTACATCCGGCGATAAAAGCCCTGaa                                                <  2:1617898/150‑1 (MQ=255)
                                                                                              caaaTCTGATGAGCGAACCTGCGAAGTGGTCAATGGCACTGGCCATAACCCGCTGGCGGACGCGGAACGGGTATGTCCGCAACTTGCGAAAACGCCGCCGGGTGGTCTATTGCAATCTGTTGATTTACATCAGACGACAAAAGCTCTg                                                  >  4:363502/2‑148 (MQ=255)
                                                                                                aaTCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTGATTTACATCCGGCGATAAAAGCCCTGaa                                                >  4:31974/1‑148 (MQ=255)
                                                                                                 aTCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTGATTTACATCCGGCGATAAAAGCCCTGaa                                                >  1:509920/1‑147 (MQ=255)
                                                                                                    tgatgaGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTGATTTACATCCGGCGATAAAAGCCCTGAAAGATg                                           >  1:1018283/1‑149 (MQ=255)
                                                                                                         aGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTGATTTACATCCGGCGATAAAAGCCCTGAAAGATGAGTGata                                    <  2:1249364/151‑1 (MQ=255)
                                                                                                            gAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTGATTTACATCCGGCGATAAAAGCCCTGAAAGATGAGTGATATCa                                 <  1:500115/151‑1 (MQ=255)
                                                                                                              aCATGCGCAGTGGTCAATGGCACAGGCCATAACCCTCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTGATTTACATCCGGCGATAAAAGCCCTGAAAGATGAGTGATATCAg                                >  2:1286311/1‑150 (MQ=255)
                                                                                                                aTGCGCAGTGGTCAATGGCACAGGCCATAACCCTCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTGATTTACATCCGGCGATAAAAGCCCTGAAAGATGAGTGATATCAGTaa                             <  4:8491/151‑1 (MQ=255)
                                                                                                                aTGCGCAGTGGTCAATGGCACAGGCCATAACCCTCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTGATTTACATCCGGCGATAAAAGCCCTGAAAGATGAGTGATATCAGTaa                             <  4:282177/151‑1 (MQ=255)
                                                                                                                   cgcAGTGGTCAATGGCACAGGCCATAACCCTCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTGATTTACATCCGGCGATAAAAGCCCTGAAAGATGAGTGATATCAGTAAGGc                          <  1:145432/151‑1 (MQ=255)
                                                                                                                    gcAGTGGTCAATGGCACAGGCCATAACCCTCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTGATTTACATCCGGCGATAAAAGCCCTGAAAGATGAGTGATATCAGTAAGGCg                         <  2:473447/151‑1 (MQ=255)
                                                                                                                    gcAGTGGTCAATGGCACAGGCCATAACCCTCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTGATTTACATCCGGCGATAAAAGCCCTGAAAGATGAGTGATATCAGAAAGGCg                         >  1:854244/1‑151 (MQ=255)
                                                                                                                        tGGTCAATGGCACAGGCCATAACCCTCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTGATTTACATCCGGCGATAAAAGCCCTGAAAGATGAGTGATATCAGTAAGGCGAgcc                     >  1:586446/1‑151 (MQ=255)
                                                                                                                        tGGTCAATGGCACAGGCCATAACCCTCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTGATTTACATCCGGCGATAAAAGCCCTGAAAGATGAGTGATATCAGTAAGGCGAgcc                     >  4:514833/1‑151 (MQ=255)
                                                                                                                         ggTCAATGGCACAGGCCATAACCCTCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTGATTTACATCCGGCGATAAAAGCCCTGAAAGATGAGTGATATCAGTAAGGCGAgcct                    <  2:1226049/151‑1 (MQ=255)
                                                                                                                         ggTCAATGGCACAGGCCATAACCCTCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTGATTTACATCCGGCGATAAAAGCCCTGAAAGATGAGTGATATCAGTAAGGCGAgcct                    <  1:460280/151‑1 (MQ=255)
                                                                                                                          gTCAATGGCACAGGCCATAACCCTCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTGATTTACATCCGGCGATAAAAGCCCTGAAAGATGAGTGATATCAGTAAGGCGAgcctg                   <  3:337152/151‑1 (MQ=255)
                                                                                                                          gTCAATGGCACAGGCCATAACCCTCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTGATTTACATCCGGCGATAAAAGCCCTGAAAGATGAGTGATATCAGTAAGGCGAgcctg                   <  3:393514/151‑1 (MQ=255)
                                                                                                                             aaTGGCACAGGCCATAACCCTCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTGATTTACATCCGGCGATAAAAGCCCTGAAAGATGAGTGATATCAGTAAGGCGAgcctgcct                <  1:1489914/151‑1 (MQ=255)
                                                                                                                             aaTGGCACAGGCCATAACCCTCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTGATTTACATCCGGCGATAAAAGCCCTGAAAGATGAGTGATATCAGTAAGGCGAgcctgcct                <  2:1323062/151‑1 (MQ=255)
                                                                                                                                  cacaGGCCATAACCCTCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTGATTTACATCCGGCGATAAAAGCCCTGAAAGATGAGTGATATCAGTAAGGCGAGCCTGCCTAAGGc           >  2:716231/1‑151 (MQ=255)
                                                                                                                                      ggCCATAACCCTCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTGATTTACATCCGGCGATAAAAGCCCTGAAAGATGAGTGATATCAGTAAGGCGAGCCTGCCTAAGGCGAtt       <  1:1219223/151‑1 (MQ=255)
                                                                                                                                        ccATAACCCTCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTGATTTACATCCGGCGATAAAAGCCCTGAAAGATGAGTGATATCAGTAAGGCGAGCCTGCCTAAGGCGAtttt     >  1:1129127/1‑151 (MQ=255)
                                                                                                                                           tAACCCTCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTGATTTACATCCGGCGATAAAAGCCCTGAAAGATGAGTGATATCAGTAAGGCGAGCCTGCCTAAGGCGAttttttt  >  3:119770/1‑151 (MQ=255)
                                                                                                                                                  |                                                                                                                                                
CAACAAAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTGATTTACATCCGGCGATAAAAGCCCTGAAAGATGAGTGATATCAGTAAGGCGAGCCTGCCTAAGGCGATTTTTTT  >  NC_000913/4399005‑4399295

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