Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,474,826 +A intergenic (+1/+44) yjgL → / ← argI hypothetical protein/ornithine carbamoyltransferase 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,474,8231.A100.0% 687.1 / NA 182*617K (TAA→AAA) yjgLhypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base A (121/61);  total (121/61)

Only 100 of 200 total aligned reads displayed.

TTAATGAACAGCGTTTAAATAAACTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAATAATGTGATGCCGGGATGGTTTCCGCATCTTATAGCGATCAGATTATTTACTGAGCGTCGCGACCATCACCGCTTTGATAGTAT  >  NC_000913/4474731‑4474910
                                                                                             |                                                                                       
  aaTGAACAGCGTTTAAATAAACTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCataa                                                                                       >  2:1647669/1‑93 (MQ=255)
    tGAACAGCGTTTAAATAAACTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCataaag                                                                                     <  2:6317546/93‑2 (MQ=255)
    tGAACAGCGTTTAAATAAACTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCataaaa                                                                                     >  1:5959578/1‑92 (MQ=255)
    tGAACAGCGTTTAAATAAACTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCataaaa                                                                                     >  1:572458/1‑92 (MQ=255)
    tGAACAGCGTTTAAATAAACTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCataaaa                                                                                     >  1:3605224/1‑92 (MQ=255)
     gAACAGCGTTTAAATAAACTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCataaaat                                                                                    >  2:5283142/1‑91 (MQ=255)
       aCAGCGTTTAAATAAACTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCataaaataa                                                                                  >  1:634406/1‑89 (MQ=255)
          gCGTTTAAATAAACTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAAtgt                                                                                <  2:894572/93‑1 (MQ=255)
          gCGTTTAAATAAACTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAAtgt                                                                                >  1:6730608/1‑93 (MQ=255)
          gCGTTTAAATAAACTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAAtgt                                                                                >  2:5672542/1‑93 (MQ=255)
          gCGTTTAAATAAACTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAAtgt                                                                                >  1:4894377/1‑93 (MQ=255)
           cGTTTAAATAAACTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAAtgtg                                                                               <  2:3892359/93‑1 (MQ=255)
             tttAAATAAACTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGAt                                                                             >  2:6696224/1‑93 (MQ=255)
             tttAAATAAACTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGAt                                                                             >  1:4490496/1‑93 (MQ=255)
               taaataaaCTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGc                                                                           <  2:6055287/93‑1 (MQ=255)
                 aataaaCTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCg                                                                         >  2:4901343/1‑93 (MQ=255)
                  ataaaCTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACGTGGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGcctg                                                                        >  2:779387/1‑91 (MQ=255)
                   taaaCTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCggg                                                                       >  2:4351115/1‑93 (MQ=255)
                    aaaCTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGa                                                                      >  1:944716/1‑93 (MQ=255)
                     aaCTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGAt                                                                     <  1:5225482/93‑1 (MQ=255)
                     aaCTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGAt                                                                     >  2:6299943/1‑93 (MQ=255)
                      aCTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATg                                                                    >  2:6200331/1‑93 (MQ=255)
                       cTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATgg                                                                   >  2:950087/1‑93 (MQ=255)
                        tAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGt                                                                  >  2:1103438/1‑93 (MQ=255)
                         aaaTCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGtt                                                                 <  1:187596/93‑1 (MQ=255)
                         aaaTCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGtt                                                                 <  1:2198921/93‑1 (MQ=255)
                         aaaTCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGtt                                                                 <  1:2811442/93‑1 (MQ=255)
                          aaTCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGttt                                                                >  2:301310/1‑93 (MQ=255)
                           aTCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGGCCTATCATGGCATAAAATAATGTGATGCCGGGCTGGTTTg                                                               >  1:3025965‑M1/1‑92 (MQ=255)
                           aTCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGCTTg                                                               >  1:3508279‑M1/1‑92 (MQ=255)
                            tCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgt                                                              >  2:5821026‑M1/1‑91 (MQ=255)
                             ccACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgta                                                             >  1:3356028‑M1/1‑90 (MQ=255)
                              cACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtat                                                            <  1:4374132‑M1/93‑5 (MQ=255)
                                ccGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtattt                                                          >  1:703696‑M1/1‑87 (MQ=255)
                                  ggAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcc                                                        >  1:5825474‑M1/1‑85 (MQ=255)
                                  ggAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcc                                                        >  1:6672538‑M1/1‑85 (MQ=255)
                                  ggAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcc                                                        <  1:5872213‑M1/93‑9 (MQ=255)
                                  ggAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcc                                                        <  1:4292741‑M1/93‑9 (MQ=255)
                                   gAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttccc                                                       >  2:1707161‑M1/1‑84 (MQ=255)
                                   gAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttccc                                                       >  1:1403375‑M1/1‑84 (MQ=255)
                                     aaaTTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccgg                                                     >  1:3924108‑M1/1‑82 (MQ=255)
                                      aaTTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggc                                                    <  1:386652‑M1/93‑13 (MQ=255)
                                      aaTTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggc                                                    >  2:201147‑M1/1‑81 (MQ=255)
                                      aaTTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggc                                                    >  2:3461942‑M1/1‑81 (MQ=255)
                                           aCGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctt                                               >  1:1300502‑M1/1‑76 (MQ=255)
                                           aCGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctt                                               >  1:4008771‑M1/1‑76 (MQ=255)
                                           aCGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctt                                               >  1:5020136‑M1/1‑76 (MQ=255)
                                             gTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatcttta                                             <  2:4008811‑M1/93‑20 (MQ=255)
                                               ataGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttata                                           <  2:3366368‑M1/93‑22 (MQ=255)
                                               ataGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttata                                           <  2:2483006‑M1/93‑22 (MQ=255)
                                                taGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatag                                          >  1:1698029‑M1/1‑71 (MQ=255)
                                                taGCAATAGAACACTTTGGGTGGAAGAATAAACCTAACACTGCATAAAATAATGTGATGTCGGGATGGTTTgtatttcctggcatctttatag                                          >  1:6204952‑M1/1‑71 (MQ=255)
                                                     aTAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgata                                     >  2:3115256‑M1/1‑66 (MQ=255)
                                                     aTAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgata                                     >  2:2244512‑M1/1‑66 (MQ=255)
                                                      tAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatag                                    >  2:4369887‑M1/1‑65 (MQ=255)
                                                       aGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagc                                   >  2:162904‑M1/1‑64 (MQ=255)
                                                         aaCACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaa                                 >  1:800570‑M1/1‑62 (MQ=255)
                                                            acTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaatta                              >  2:1537133‑M1/1‑59 (MQ=255)
                                                             cTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattat                             >  2:6583799‑M1/1‑58 (MQ=255)
                                                             cTTTGGGGGGAAGAATAGACCTATCACTGCATAAAATAATGGGATGCCGGGGTGGTTTgtatttcccggcatctttatagcgatagcaattat                             >  1:1765140‑M1/1‑58 (MQ=255)
                                                              tttGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatt                            <  2:6641683‑M1/93‑37 (MQ=255)
                                                               ttGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattattt                           >  1:5639913‑M1/1‑56 (MQ=255)
                                                                 gggTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttac                         >  2:2049327‑M1/1‑54 (MQ=255)
                                                                 gggTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttac                         <  1:4431571‑M1/93‑40 (MQ=255)
                                                                 gggTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttac                         >  2:3743369‑M1/1‑54 (MQ=255)
                                                                  ggTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttact                        <  1:4158243‑M1/93‑41 (MQ=255)
                                                                     ggAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgag                     >  1:1964122‑M1/1‑50 (MQ=255)
                                                                     ggAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgag                     <  2:4181421‑M1/93‑44 (MQ=255)
                                                                      gaagaaTAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagt                    <  1:3914657‑M1/93‑45 (MQ=255)
                                                                      gaagaaTAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagt                    <  2:4245162‑M1/93‑45 (MQ=255)
                                                                      gaagaaTAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagt                    >  1:4780679‑M1/1‑49 (MQ=255)
                                                                       aagaaTAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagtg                   <  2:5019215‑M1/93‑46 (MQ=255)
                                                                        agaaTAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagtgt                  >  2:1136167‑M1/1‑47 (MQ=255)
                                                                        agaaTAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagtgt                  >  2:6333427‑M1/1‑47 (MQ=255)
                                                                        agaaTAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagtgt                  >  2:3550093‑M1/1‑47 (MQ=255)
                                                                        agaaTAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagtgt                  <  1:4207600‑M1/93‑47 (MQ=255)
                                                                         gaaTAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagtgtc                 >  1:886312‑M1/1‑46 (MQ=255)
                                                                         gaaTAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagtgtc                 <  2:4650649‑M1/93‑48 (MQ=255)
                                                                         gaaTAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagtgtc                 >  1:3637968‑M1/1‑46 (MQ=255)
                                                                           aTAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGGTTgtatttcccggcatctttatagcgatagcaattatttactgagtgtcgc               >  1:206186‑M1/1‑44 (MQ=255)
                                                                            tAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagtgtcgcg              >  1:4992920‑M1/1‑43 (MQ=255)
                                                                            tAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagtgtcgcg              >  2:4519854‑M1/1‑43 (MQ=255)
                                                                            tAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagtgtcgcg              <  2:5944664‑M1/93‑51 (MQ=255)
                                                                             aGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagtgtcgcga             >  1:5046450‑M1/1‑42 (MQ=255)
                                                                              gACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagtgtcgcgac            >  2:2803713‑M1/1‑41 (MQ=255)
                                                                              gACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagtgtcgcgac            >  2:2140579‑M1/1‑41 (MQ=255)
                                                                               aCCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagtgtcgcgacc           >  1:4570582‑M1/1‑40 (MQ=255)
                                                                                ccTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagtgtcgcgacca          <  2:3304710‑M1/93‑55 (MQ=255)
                                                                                ccTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagtgtcgcgacca          <  1:4400145‑M1/93‑55 (MQ=255)
                                                                                 cTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagtgtcgcgaccat         >  2:4812252‑M1/1‑38 (MQ=255)
                                                                                  tATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagtgtcgcgaccatc        >  1:1510334‑M1/1‑37 (MQ=255)
                                                                                  tATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagtgtcgcgaccatc        >  1:1940803‑M1/1‑37 (MQ=255)
                                                                                   aTCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagtgtcgcgaccatca       >  1:3179406‑M1/1‑36 (MQ=255)
                                                                                    tCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagtgtcgcgaccatcac      <  1:6490696‑M1/93‑59 (MQ=255)
                                                                                    tCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagtgtcgcgaccatcac      >  2:5361806‑M1/1‑35 (MQ=255)
                                                                                     cACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagtgtcgcgaccatcacc     <  1:1210443‑M1/93‑60 (MQ=255)
                                                                                      aCTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagtgtcgcgaccatcaccg    <  2:6516108‑M1/93‑61 (MQ=255)
                                                                                       cTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagtgtcgcgaccatcaccgc   >  1:5598659‑M1/1‑32 (MQ=255)
                                                                                        tGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagtgtcgcgaccatcaccgct  >  1:2273343‑M1/1‑31 (MQ=255)
                                                                                        tGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagtgtcgcgaccatcaccgct  <  1:2643779‑M1/93‑63 (MQ=255)
                                                                                             |                                                                                       
TTAATGAACAGCGTTTAAATAAACTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAATAATGTGATGCCGGGATGGTTTCCGCATCTTATAGCGATCAGATTATTTACTGAGCGTCGCGACCATCACCGCTTTGATAGTAT  >  NC_000913/4474731‑4474910

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: