Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,104,065 G→A *351* (TAG→TAA mutY → adenine DNA glycosylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,104,0650GA100.0% 282.8 / NA 86*351* (TAG→TAAmutYadenine DNA glycosylase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (40/46);  total (40/46)

TCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGTTACAGCAGTTACGCACTGGCGCGCCGGTTTAGCGCGTGAGTCGATAAAGAGGATGATTTATGAGCAGAACGATTTTTTGTACTTTCCTGCAACGTGAAGCAGAAGGTCAGGATTTTCAGCTGTACCCCGGCGAGCTGGGAAAACGCATCTATAACGAG  >  NC_000913/3103930‑3104191
                                                                                                                                       |                                                                                                                              
tCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGTTACAGCAGTTACGCACTGGCGCGCCGGTTTAacgc                                                                                                                             >  2:344740/1‑139 (MQ=255)
   cGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGTTACAGCAGTTACGCACTGGCGCGCCGGTTTAACGCGTg                                                                                                                          <  1:106370/139‑1 (MQ=255)
      gtcgtcATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGTTACAGCAGTTACGCACTGGCGCGCCGGTTTAACGCGTGAGt                                                                                                                       <  2:828956/139‑1 (MQ=255)
       tcgtcATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGTTACAGCAGTTACGCACTGGCGCGCCGGTTTAACGCGTGAGTc                                                                                                                      <  1:1117358/139‑1 (MQ=255)
        cgtcATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGTTACAGCAGTTACGCACTGGCGCGCCGGTTTAACGCGTGAGTCg                                                                                                                     >  2:1654382/1‑139 (MQ=255)
            aTTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGTTACAGCAGTTACGCACTGGCGCGCCGGTTTAACGCGTGAg                                                                                                                        >  1:1038759/1‑132 (MQ=255)
            aTTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGTTACAGCAGTTACGCACTGGCGCGCCGGTTTAACGCGTGAg                                                                                                                        <  2:1038759/132‑1 (MQ=255)
            aTTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGTTACAGCAGTTACGCACTGGCGCGCCGGTTTAACGCGTGAGTCGATaa                                                                                                                 >  2:304044/1‑139 (MQ=255)
               cACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGTTACAGCAGTTACGCACTGGCGCGCCGGTTTAACGCGTGAGTCGAt                                                                                                                   <  1:394434/134‑1 (MQ=255)
               cACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGTTACAGCAGTTACGCACTGGCGCGCCGGTTTAACGCGTGAGTCGAt                                                                                                                   >  2:394434/1‑134 (MQ=255)
               cACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGACCTAGCGGCTCCCGTGGAGCGTTTGTTACAGCAGTTACGCACTGGCGCGCCGGTTTAACGCGTg                                                                                                                          <  1:466077/127‑1 (MQ=255)
               cACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGACCTAGCGGCTCCCGTGGAGCGTTTGTTACAGCAGTTACGCACTGGCGCGCCGGTTTAACGCGTg                                                                                                                          >  2:466077/1‑127 (MQ=255)
                   ggCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGTTACAGCAGTTACGCACTGGCGCGCCGGTTTAACGCGTGAGTCGATAAAGAGgat                                                                                                          <  2:1594896/139‑1 (MQ=255)
                   ggCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGTTACAGCAGTTACGCACTGGCGCGCCGGTTTAACGCGTGAGTCGATAAAGAGgat                                                                                                          <  1:1654382/139‑1 (MQ=255)
                    gCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGTTACAGCAGTTACGCACTGGCGCGCCGGTTTAACGCGTGAGTCGATAAAGAGgatg                                                                                                         >  2:1089248/1‑139 (MQ=255)
                     cTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGTTACAGCAGTTACGCACTGGCGCGCCGGTTTAACGCGTGAg                                                                                                                        >  2:1260589/1‑123 (MQ=255)
                     cTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGTTACAGCAGTTACGCACTGGCGCGCCGGTTTAACGCGTGAg                                                                                                                        <  1:1260589/123‑1 (MQ=255)
                     cTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGTTACAGCAGTTACGCACTGGCGCGCCGGTTTAACGCGTGAGTCGATAAAGAGgatga                                                                                                        >  2:1608006/1‑139 (MQ=255)
                        cATGGATGAAGGCAATGCGCTCTCGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGTTACAGCAGTTACGCACTGGCGCGCCGGTTTAACGCGTGAGTCGATAAAGAGGATGAttt                                                                                                     <  2:1348888/139‑1 (MQ=255)
                                    cAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGTTACAGCAGTTACGCACTGGCGCGCCGGTTTAACGCGTGAGTCGATAAAGAGGATGATTTATGAGCAGAACg                                                                                         <  1:637792/139‑1 (MQ=255)
                                       tGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGTTACAGCAGTTACGCACTGGCGCGCCGGTTTAACGCGTGAGTCGATAAAGAGGATGCTTTATGAGCAGAACGAtt                                                                                      >  2:1123834/1‑139 (MQ=255)
                                         cgcTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGTTACAGCAGTTACGCACTGGCGCGCCGGTTTAACGCGTGAGTCGATAAAGAGGATGATTTATGAGCAGAACGAtttt                                                                                    <  1:268029/139‑1 (MQ=255)
                                                 tataACTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGTTACAGCAGTTACGCACTGGCGCGCCGGTTTAAcgcg                                                                                                                            >  2:1137114/1‑91 (MQ=255)
                                                 tataACTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGTTACAGCAGTTACGCACTGGCGCGCCGGTTTAAcgcg                                                                                                                            <  1:1137114/91‑1 (MQ=255)
                                                     aCTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGTTACAGCAGTTACGCACTGGCGCGCCGGTTTAACGCGTGAGTCGATAAAGAgg                                                                                                            <  2:1428444/103‑1 (MQ=255)
                                                     aCTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGTTACAGCAGTTACGCACTGGCGCGCCGGTTTAACGCGTGAGTCGATAAAGAgg                                                                                                            >  1:1428444/1‑103 (MQ=255)
                                                      cTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGTTACAGCAGTTACGCACTGGCGCGCCGGTTTAACGCGTGAGTCGATAAAGAGgat                                                                                                          <  2:390657/104‑1 (MQ=255)
                                                      cTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGTTACAGCAGTTACGCACTGGCGCGCCGGTTTAACGCGTGAGTCGATAAAGAGgat                                                                                                          >  1:390657/1‑104 (MQ=255)
                                                        tAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGTTACAGCAGTTACGCACTGGCGCGCCGGTTTAACGCGTGAGTCGATAAAGAGGATGATTTATGAGCAGAACGATTTTTTGTACTTTCCTGCa                                                                     >  1:1026807/1‑139 (MQ=255)
                                                        tAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGTTACAGCAGTTACGCACTGGCGCGCCGGTTTAACGCGTGAGTCGATAAAGAGGATGATTTATGAGCAGAACGATTTTTTGTACTTTCCTGCa                                                                     >  1:651098/1‑139 (MQ=255)
                                                          gcgcAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGTTACAGCAGTTACGCACTGGCGCGCCGGTTTAAcgcg                                                                                                                            >  2:1306638/1‑82 (MQ=255)
                                                          gcgcAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGTTACAGCAGTTACGCACTGGCGCGCCGGTTTAAcgcg                                                                                                                            <  1:1306638/82‑1 (MQ=255)
                                                          gcgcAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGTTACAGCAGTTACGCACTGGCGCGCCGGTTTAACGCGTGAGTCGATAAAGAGGATGATTTATGAGCAGAACGATTTTTTGTa                                                                               >  2:971145/1‑127 (MQ=255)
                                                          gcgcAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGTTACAGCAGTTACGCACTGGCGCGCCGGTTTAACGCGTGAGTCGATAAAGAGGATGATTTATGAGCAGAACGATTTTTTGTa                                                                               <  1:971145/127‑1 (MQ=255)
                                                          gcgcAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGTTACAGCAGTTACGCACTGGCGCGCCGGTTTAACGCGTGAGTCGATAAAGAGGATGATTTATGAGCAGAACGATTTTTTGTACTTTCCTGCAAc                                                                   >  1:720162/1‑139 (MQ=255)
                                                          gcgcAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGTTACAGCAGTTACGCACTGGCGCGCCGGTTTAACGCGTGAGTCGATAAAGAGGATGATTTATGAGCAGAACGATTTTTTGTACTTTCCTGCAAc                                                                   >  1:485287/1‑139 (MQ=255)
                                                           cgcAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGTTACAGCAGTTACGCACTGGCGCGCCGGTTTAAcgcg                                                                                                                            <  1:235843/81‑1 (MQ=255)
                                                           cgcAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGTTACAGCAGTTACGCACTGGCGCGCCGGTTTAAcgcg                                                                                                                            >  2:235843/1‑81 (MQ=255)
                                                           cgcAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGTTACAGCAGTTACGCACTGGCGCGCCGGTTTAACGCGTGAGTCGATAAAGAGGATGATTTATGAGCAGAACGATTTTTTGTACTTTCCTGCAACg                                                                  <  2:1130628/139‑1 (MQ=255)
                                                           cgcAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGTTACAGCAGTTACGCACTGGCGCGCCGGTTTAACGCGTGAGTCGATAAAGAGGATGATTTATGAGCAGAACGATTTTTTGTACTTTCCTGCAACg                                                                  >  1:1115182/1‑139 (MQ=255)
                                                           cgcAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGTTACAGCAGTTACGCACTGGCGCGCCGGTTTAACGCGTGAGTCGATAAAGAGGATGATTTATGAGCAGAACGAGTTTTTGTACTTTCCTGCAACg                                                                  >  2:980812/1‑139 (MQ=255)
                                                              aaCCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGTTACAGCAGTTACGCACTGGCGCGCCGGTTTAACGCGTGAGTCGATAAAGAGGATGATTTATGAGCAGAACGATTTTTTGTACTTTCCTGCAACGTGa                                                               <  2:152791/139‑1 (MQ=255)
                                                                   ccgTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGTTACAGCAGTTACGCACTGGCGCGCCGGTTTAACGCGTGAGTCGATAAAGAGGATGATTTATGAGCAGAACGATTTTTTGTACTTTCCTGCAACGTGAAGCAg                                                          <  1:344740/139‑1 (MQ=255)
                                                                     gTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGTTACAGCAGTTACGCACTGGCGCGCCGGTTTAACGCGTGAGTCGATAAAGAGGATGATTTATGAGCAGAACGATTTTTTGTACTTTCCTGCAACGTGAAGCa                                                           >  1:1270306/1‑136 (MQ=255)
                                                                     gTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGTTACAGCAGTTACGCACTGGCGCGCCGGTTTAACGCGTGAGTCGATAAAGAGGATGATTTATGAGCAGAACGATTTTTTGTACTTTCCTGCAACGTGAAGCa                                                           <  2:1270306/136‑1 (MQ=255)
                                                                       cAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGTTACAGCAGTTACGCACTGGCGCGCCGGTTTAACGCGTGAGTCGATAAAGAGGATGATTTATGAGCAGAACGATTTTTTGTACTTTCCTGCAACGTGAAGCAGAAgg                                                      <  2:1115182/139‑1 (MQ=255)
                                                                                    ggCTCCCGTGGAGCGTTAGTTACAGCAGTTACGCACTGGCGCGCCGGTTTAACGCGTGAGTCGATAAAGAGGATGATTTATGAGCAGAACGATTTTTCGTACTTTCCTGCAACGTGAAGCAGAAGGTCAGGATTTTCac                                         <  2:1026807/139‑2 (MQ=255)
                                                                                       tCCCGTGGAGCGTTTGTTACAGCAGTTACGCACTGGCGCGCCGGTTTAACGCGTGAGTCGATAAAGAGGATGATTTATGAGCAGAACGATTTTTTGTACTTTCCTGCAACGTGAAGCAGAAGGTCAGGATTTTCAGCTg                                      >  1:674648/1‑139 (MQ=255)
                                                                                        cccGTGGAGCGTTTGTTACAGCAGTTACGCACTGGCGCGCCGGTTTAAcgcg                                                                                                                            <  2:347137/52‑1 (MQ=255)
                                                                                        cccGTGGAGCGTTTGTTACAGCAGTTACGCACTGGCGCGCCGGTTTAAcgcg                                                                                                                            >  1:347137/1‑52 (MQ=255)
                                                                                        cccGTGGAGCGTTTGTTACAGCAGTTACGCACTGGCGCGCCGGTTTAACGCGTGAGTCGATAAagag                                                                                                             <  2:909919/67‑1 (MQ=255)
                                                                                        cccGTGGAGCGTTTGTTACAGCAGTTACGCACTGGCGCGCCGGTTTAACGCGTGAGTCGATAAagag                                                                                                             >  1:909919/1‑67 (MQ=255)
                                                                                        cccGTGGAGCGTTTGTTACAGCAGTTACGCACTGGCGCGCCGGTTTAACGCGTGAGTCGATAAAGAGGATGATTTATGAg                                                                                                <  1:1052139/80‑1 (MQ=255)
                                                                                        cccGTGGAGCGTTTGTTACAGCAGTTACGCACTGGCGCGCCGGTTTAACGCGTGAGTCGATAAAGAGGATGATTTATGAg                                                                                                >  2:1052139/1‑80 (MQ=255)
                                                                                        cccGTGGAGCGTTTGTTACAGCAGTTACGCACTGGCGCGCCGGTTTAACGCGTGAGTCGATAAAGAGGATGATTTATGAGCAGAACGATTTTTTGTACTTTCCTGCAACGTGAAGCAGAAGGTCAGGATTTTCAGCTGt                                     >  2:1038437/1‑139 (MQ=255)
                                                                                         ccGTGGAGCGTTTGTTACAGCAGTTACGCACTGGCGCGCCGGTTTAACGCGTGAGTCGATAAAGAGGATGATTTATGAGCAGAACGATTTTTTGTACTTTCCTGCAACGTGAAGCAGAAGGTCAGGATTTTCAGCTGTa                                    <  1:304044/139‑1 (MQ=255)
                                                                                            tGGAGCGTTTGTTACAGCAGTTACGCACTGGCGCGCCGGTTTAACGCGTGAGTCGATAAAGAGgat                                                                                                          <  2:798014/66‑1 (MQ=255)
                                                                                            tGGAGCGTTTGTTACAGCAGTTACGCACTGGCGCGCCGGTTTAACGCGTGAGTCGATAAAGAGgat                                                                                                          >  1:798014/1‑66 (MQ=255)
                                                                                            tGGAGCGTTTGTTACAGCAGTTACGCACTGGCGCGCCGGTTTAACGCGTGAGTCGATAAAGAGGATGATTTATGAGCAGAACGATTTTTTg                                                                                 <  2:741793/91‑1 (MQ=255)
                                                                                            tGGAGCGTTTGTTACAGCAGTTACGCACTGGCGCGCCGGTTTAACGCGTGAGTCGATAAAGAGGATGATTTATGAGCAGAACGATTTTTTg                                                                                 >  1:741793/1‑91 (MQ=255)
                                                                                             ggAGCGTTTGTTACAGCAGTTACGCACTGGCGCGCCGGTTTAACGCGTGAGTCGATAAAGAGGATGATTTATGAGCAGAACGATTTTTTGTACTTTCCTg                                                                       >  1:716224/1‑100 (MQ=255)
                                                                                             ggAGCGTTTGTTACAGCAGTTACGCACTGGCGCGCCGGTTTAACGCGTGAGTCGATAAAGAGGATGATTTATGAGCAGAACGATTTTTTGTACTTTCCTg                                                                       <  2:716224/100‑1 (MQ=255)
                                                                                               aGCGTTTGTTACAGCAGTTACGCACTGGCGCGCCGGTTTAACGCGTGAGTCGATAAAGAGGATGATTTATGAGCAGAACGATTTTTTGTACTTTCCTGCAACGTGAAGCAGAAGGTCAGGATTTTCAGCTGTACCCCgg                              <  2:44523/139‑1 (MQ=255)
                                                                                               aGCGTTTGTTACAGCAGTTACGCACTGGCGCGCCGGTTTAACGCGTGAGTCGATAAAGAGGATGATTTATGAGCAGAACGATTTTTTGTACTTTCCTGCAACGTGAAGCAGAAGGTCAGGATTTTCAGCTGTACCCCgg                              <  2:720162/139‑1 (MQ=255)
                                                                                               aGCGTTTGTTACAGCAGTTACGCACTGGCGCGCCGGTTTAACGCGTGAGTCGATAAAGAGGATGATTTATGAGCAGAACGATTTTTTGTACTTTCCTGCAACGTGAAGCAGAAGGTCAGGATTTTCAGCTGTACCCCgg                              <  1:1608006/139‑1 (MQ=255)
                                                                                                   tttGTTACAGCAGTTACGCACTGGCGCGCCGGTTTAACGCGTGAGTCGATAAAGAGga                                                                                                           >  1:1393982/1‑58 (MQ=255)
                                                                                                   tttGTTACAGCAGTTACGCACTGGCGCGCCGGTTTAACGCGTGAGTCGATAAAGAGga                                                                                                           <  2:1393982/58‑1 (MQ=255)
                                                                                                   tttGTTACAGCAGTTACGCACTGGCGCGCCGGTTTAACGCGTGAGTCGATAAAGAGGATGATTTATGAGCAGAACGATTTTTTGTACTTTCCTGCAACGTGAAGCAGAAGGTCAGGAtt                                              >  2:1369760/1‑119 (MQ=255)
                                                                                                   tttGTTACAGCAGTTACGCACTGGCGCGCCGGTTTAACGCGTGAGTCGATAAAGAGGATGATTTATGAGCAGAACGATTTTTTGTACTTTCCTGCAACGTGAAGCAGAAGGTCAGGAtt                                              <  1:1369760/119‑1 (MQ=255)
                                                                                                   tttGTTACAGCAGTTACGCACTGGCGCGCCGGTTTAACGCGTGAGTCGATAAAGAGGATGATTTATGAGCAGAACGATTTTTTGTACTTTCCTGCAACGTGAAGCAGAAGGTCAGGATTTTCAGCTGTACCCCGGCGAg                          >  2:459374/1‑139 (MQ=255)
                                                                                                          cagcagTTACGCACTGGCGCGCCGGTTTAACGCGTGAGTCGATAAAGAGGATGATTTATGAGCAGAACGATTTTTTCTACTTTCCTGCAACGTGAAGCAGAAGGTCAGGATTTTCAGCTGTACCCCGGCGAGCTGGGaa                   <  1:1089248/139‑1 (MQ=255)
                                                                                                           agcagTTACGCACTGGCGCGCCGGTTTAACGCGTGAGTCGATAAAGAGGATGATTTATGAGCAGAACGATTTTTTGTACTTTCCTGCAACGTGAAGCAGAAGGTCAGGATTTTCAGCTGTACCCCGGCGAGCTGGGaaa                  >  2:1386935/1‑139 (MQ=255)
                                                                                                            gcagTTACGCACTGGCGCGCCGGTTTAACGCGTGAGTCGATAAAGAGGATGATTTATGAGCAGAACGATTTTTTGTACTTTCCTGCAACGTGAAGCAGAAGGTCAGGATTTTCAGCTGTACCCCGGCGAGCTGGGaaaa                 <  2:674648/139‑1 (MQ=255)
                                                                                                             cagTTACGCACTGGCGCGCCGGTTTAACGCGTGAGTCGATAAAGAGGATGATTTATGAGCAGAACGATTTTTTGTACTTTCCTGCAACGTGAAGCAGAAGGTCAGGATTTTCAGCTGTACCCCGGCGAGCTGGGAAAAc                >  1:55582/1‑139 (MQ=255)
                                                                                                               gTTACGCACTGGCGCGCCGGTTTAACGCGTGAGTCGATAAAGAGGATGATTTATGAGCAGAACGATTTTTTGTACTTTCCTGCAACGTGAAGCAGAAGGTCAGGATTTTCAGCTGTACCCCGGCGAGCTGGGAAAACGc              <  2:707636/139‑1 (MQ=255)
                                                                                                                ttACGCACTGGCGCGCCGGTTTAACGCGTGAGTCGATAAagag                                                                                                             <  2:663219/43‑1 (MQ=255)
                                                                                                                ttACGCACTGGCGCGCCGGTTTAACGCGTGAGTCGATAAagag                                                                                                             >  1:663219/1‑43 (MQ=255)
                                                                                                                     cACTGGCGCGCCGGTTTAACGCGTGAGTCGATAAAGAGGATGATTTATGAGCAGAACGATTTTTTGTACTTTCCTACAACGTGAAGCAGAAGGTCAGGATTTTCAGCTGTACCCCGGCGAGCTGGGAAAACGCATCtat        <  1:543929/139‑1 (MQ=255)
                                                                                                                       cTGGCGCGCCGGTTTAACGCGTGAGTCGATAAAGAGGATGATTTATGAGCTGAACGATTTTTTGTACTTTCCTGCAACGTGAAGCAg                                                          <  2:1061692/87‑1 (MQ=255)
                                                                                                                       cTGGCGCGCCGGTTTAACGCGTGAGTCGATAAAGAGGATGATTTATGAGCAGAACGATTTTTTGTACTTTCCTGCAACGTGAAGCAg                                                          >  1:1061692/1‑87 (MQ=255)
                                                                                                                        tGGCGCGCCGGTTTAACGCGTGAGTCGATAAAGAGGATGATTTATGAGCAGAACGATTTTTTGTACTTTCCTGCAACGTGAAGCAGAAGGTCAGGATTTTCAGCTGTACCCCGGCGAGCTGGGAAAACGCATCTATAAc     <  2:485287/139‑1 (MQ=255)
                                                                                                                        tGGCGCGCCGGTTTAACGCGTGAGTCGATAAAGAGGATGATTTATGAGCAGAACGATTTTTTGTACTTTCCTGCAACGTGAAGCAGAAGGTCAGGATTTTCAGCTGTACCCCGGCGAGCTGGGAAAACGCATCTATAAc     >  2:905098/1‑139 (MQ=255)
                                                                                                                          gcgcgcCGGTTTAACGCGTGAGTCGATAAAGAGgat                                                                                                          >  2:875461/1‑36 (MQ=255)
                                                                                                                          gcgcgcCGGTTTAACGCGTGAGTCGATAAAGAGgat                                                                                                          <  1:875461/36‑1 (MQ=255)
                                                                                                                           cgcgcCGGTTTAACGCGTGAGTCGATAAAGAGGATGATTTATGAGCAGAACGATTTTTTGTACTTTCCTGCAACGTGAAGCAGAAGGTCAGGATTTTCAGCTGTACCCCGGCGAGCTGGGAAAACGCATCTATAACgag  <  2:372202/139‑1 (MQ=255)
                                                                                                                               ccGGTTTAACGCGTGAGTCGATAAAGAGGATGATTTATGAGCAGAACg                                                                                         >  1:203674/1‑48 (MQ=255)
                                                                                                                               ccGGTTTAACGCGTGAGTCGATAAAGAGGATGATTTATGAGCAGAACg                                                                                         <  2:203674/48‑1 (MQ=255)
                                                                                                                                       |                                                                                                                              
TCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGTTACAGCAGTTACGCACTGGCGCGCCGGTTTAGCGCGTGAGTCGATAAAGAGGATGATTTATGAGCAGAACGATTTTTTGTACTTTCCTGCAACGTGAAGCAGAAGGTCAGGATTTTCAGCTGTACCCCGGCGAGCTGGGAAAACGCATCTATAACGAG  >  NC_000913/3103930‑3104191

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